Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Accretion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Ag isotopic and chalcophile element evolution of the terrestrial and martian mantles during accretion: new constraints from Bi, Pd, and Ag metal-silicate partitioning.

The Earth’s timing of accretion and acquisition of moderately volatile compounds is uncertain. Hafnium-W and Mn-Cr isotopic data can bracket the timing of early planetary differentiation and core formation. The Ag-Pd system has also been utilized but its application has been limited by a lack of high pressure and temperature metal-silicate partitioning for Pd and Ag. Be-cause Ag (and Bi) are volatile chalcophile siderophile elements, understanding their early distribution can constrain the origin of volatile elements in differentiated bodies and planets. Unfortunately, neither Ag or Bi have been studied across the wide range of pressure and temperature conditions that are relevant to accretion and core-mantle differentiation. Here, new high-pressure and temperature multi-anvil metal-silicate equilibrium experiments for Bi and Ag have been carried out at conditions relevant to planetary accretion and metal silicate differentiation that allow a more refined and complete understanding of element partitioning during core formation. The new metal-silicate partitioning data utilized to predict the distributions of Bi, Pd, and Ag at conditions of accretion for Earth and Mars and show that the Pd/Ag ratio is significantly fractionated during accretion, allowing for the production of detectable 107Ag anomalies produced while 107Pd (half life = 6.5 M.y.) was still extant. Application of the new partitioning results to Earth shows that D(Bi) and D(Ag) (D = metal/silicate concentration ratio) are lowered due to the effect of pressure and Si alloyed in the metallic liquid, resulting in higher predicted mantle Bi and Ag abundances than in the bulk silicate Earth (BSE), as well as high and variable Pd/Ag. The unradiogenic Ag isotopic composition of the BSE could have been generated by early accretion of volatile-poor (high Pd/Ag) precursors, followed by later accretion of volatile–rich (low Pd/Ag) material, in agreement with earlier studies of Pd-Ag and Mn-Cr (Schönbächler et al., 2010). However, these main accretion phases would have to be followed by segregation of a sulfide liquid (at least 1.5% of magma ocean) at high pressures (>30 GPa), to explain the PUM Bi, Pd, and Ag, as well as Au, Pt, Cu and Ni concentrations as proposed previously. If the early accreted bulk Earth was volatile depleted with high Pd/Ag ratios, portions of the mantle may contain ancient domains that developed positive 107Ag isotopic anomalies (as also argued by noble gases, Nd, W, and Os isotopes). In comparison, Bi, Pd, and Ag concentrations in the martian mantle could have been set by simple metal-silicate equilibrium. Mars accreted and differentiated relatively rapidly, while also developing a deep magma ocean with a high Pd/Ag ratio that could have evolved positive 107Ag anomalies, in contrast to Earth. Measurements on shergottites may reveal these predicted Ag isotopic anomalies.

accretion

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.

Ice Accretions and Full-Scale Iced Aerodynamic Performance Data for a Two-Dimensional NACA 23012 Airfoil

This report documents the data collected during the large wind tunnel campaigns conducted as part of the SUNSET project (StUdies oN Scaling EffecTs due to ice) also known as the Ice-Accretion Aerodynamics Simulation study: a joint effort by NASA, the Office National d'Etudes et Recherches Aérospatiales (ONERA), and the University of Illinois. These data form a benchmark database of full-scale ice accretions and corresponding ice-contaminated aerodynamic performance data for a two-dimensional (2D) NACA 23012 airfoil. The wider research effort also included an analysis of ice-contaminated aerodynamics that categorized ice accretions by aerodynamic effects and an investigation of subscale, low- Reynolds-number ice-contaminated aerodynamics for the NACA 23012 airfoil. The low-Reynolds-number investigation included an analysis of the geometric fidelity needed to reliably assess aerodynamic effects of airfoil icing using artificial ice shapes. Included herein are records of the ice accreted during campaigns in NASA Glenn Research Center's Icing Research Tunnel (IRT). Two different 2D NACA 23012 airfoil models were used during these campaigns; an 18-in. (45.7-cm) chord (subscale) model and a 72-in. (182.9-cm) chord (full-scale) model. The aircraft icing conditions used during these campaigns were selected from the Federal Aviation Administration's (FAA's) Code of Federal Regulations (CFR) Part 25 Appendix C icing envelopes. The records include the test conditions, photographs of the ice accreted, tracings of the ice, and ice depth measurements. Model coordinates and pressure tap locations are also presented. Also included herein are the data recorded during a wind tunnel campaign conducted in the F1 Subsonic Pressurized Wind Tunnel of ONERA. The F1 tunnel is a pressured, high- Reynolds-number facility that could accommodate the full-scale (72-in. (182.9-cm) chord) 2D NACA 23012 model. Molds were made of the ice accreted during selected test runs of the full-scale model in the IRT. From these molds, castings were made that closely replicated the features of the accreted ice. The castings were then mounted on the full-scale model in the F1 tunnel, and aerodynamic performance measurements were made using model surface pressure taps, the facility force balance system, and a large wake rake designed specifically for these tests. Tests were run over a range of Reynolds and Mach numbers. For each run, the model was rotated over a range of angles-of-attack that included airfoil stall. The benchmark data collected during these campaigns were, and continue to be, used for various purposes. The full-scale data form a unique, ice-accretion and associated aerodynamic performance dataset that can be used as a reference when addressing concerns regarding the use of subscale ice-accretion data to assess full-scale icing effects. Further, the data may be used in the development or enhancement of both ice-accretion prediction codes and computational fluid dynamic codes when applied to study the effects of icing. Finally, as was done in the wider study, the data may be used to help determine the level of geometric fidelity needed for artificial ice used to assess aerodynamic degradation due to aircraft icing. The structured, multifaceted approach used in this research effort provides a unique perspective on the aerodynamic effects of aircraft icing. The data presented in this report are available in electronic form upon formal approval by proper NASA and ONERA authorities.

Addy, Harold E., Jr.

Migration of Accreting Planets in Radiative Discs From Dynamical Torques

We present the results of hydrodynamical simulations of the orbital evolution of planets undergoing runaway gas accretion in radiative discs. We consider accreting disc models with constant mass flux through the disc, and where radiative cooling balances the effect of viscous heating and stellar irradiation. We assume that 20–30 M(sub ⨁) giant planet cores are formed in the region where viscous heating dominates and migrate outward under the action of a strong entropy-related corotation torque. In the case where gas accretion is neglected and for an α viscous stress parameter α = 2 × 10(exp -3), we find evidence for strong dynamical torques in accreting discs with accretion rates M ̇ ≳ 7 × 10(exsp -8)M(sub ⨀) yr(exp -1). Their main effect is to increase outward migration rates by a factor of ∼2 typically. In the presence of gas accretion, however, runaway outward migration is observed with the planet passing through the zero-torque radius and the transition between the viscous heating and stellar heating dominated regimes. The ability for an accreting planet to enter a fast migration regime is found to depend strongly on the planet growth rate but can occur for values of the mass flux through the disc of M ̇ ≳ 5 × 10(exp -8)M(sub ⨀) yr(exp -1).We find that an episode of runaway outward migration can cause an accreting planet formed in the 5–10 au region to temporarily orbit at star–planet separations as large as∼60–70 au. However, increase in the amplitude of the Lindblad torque associated with planet growth plus change in the streamline topology near the planet systematically cause the direction of migration to be reversed. Subsequent evolution corresponds to the planet migrating inward rapidly until it becomes massive enough to open a gap in the disc and migrate in the type II regime. Our results indicate that a planet can reach large orbital distances under the combined effect of dynamical torques and gas accretion, but an alternative mechanism is required to explain the presence of massive planets on wide orbits.

Accretion

Accretion flows in binary X-ray systems

Recent developments in the study of accretion flows in binary X-ray systems are reviewed, concentrating on issues concerning the gas flows rather than the radiation mechanisms. A model for Her X-1 is presented which illustrates the main structures expected to exist in the case where mass transfer is by Roche-lobe overflow and the X-ray source is a rotating magnetized neutron star. A model for Cyg X-1 is given which shows the accretion-flow structures to be expected if the compact object is a black hole, including an accretion wake and a very small accretion disk around the black hole. Physical parameters of the accretion flows in these two instances are examined along with some scale lengths that may be derived from these parameters by dimensional analysis. Detailed attention is given to neutron-star magnetospheres (for slow, intermediate, and fast rotators), accretion-disk structure, accretion from stellar winds, and the relationship between X-rays and the gas dynamics of accretion flows. It is suggested that accretion flows in binary X-ray systems be observed directly by means of high-resolution UV and X-ray spectroscopy.

Mccray, R.

Can accretion onto isolated neutron stars produce gamma-ray bursts?

It is shown here that, if slow accretion onto high-velocity neutron stars from the interstellar medium is to be the origin of gamma-ray bursts, the accretion physics is very different from what applies for local low-velocity neutron stars. For halo neutron stars with high magnetic fields and velocities, electromagnetic dipole radiation pressure prevents accretion unless the period is longer than tens of seconds; the centrifugal barrier will then prevent accretion until the period reaches several thousand seconds. For periods as long as this, accretion may proceed through Kelvin-Helmholtz instability at the magnetopause boundary. At interstellar densities and neutron-star magnetic fields of about 10 exp 12 G, the accretion rate by this process can be much larger than the Bondi-Hoyle accretion rate, but is still well below what is needed for slow-accretion burst models. It is concluded that slow accretion onto high-velocity neutron stars in the halo cannot be the origin of gamma-ray bursts.

Harding, Alice K.

The hard X-ray emission spectra from accretion columns in intermediate polars

We consider the hard (greater than 2 keV) X-ray emission from accretion columns in an intermediate polar system, GK Per, using a simple settling solution. The rate of photon emission per logarithmic energy interval can be fitted with a power law, E(exp -gamma), with gamma approximately 2.0, in agreement with observations. This index is only weakly dependent on the mass accretion rate, dot-M, for dot-M in the range of a few times 10(exp 16-18) g/s. The peak energy of the photon spectra (after photoelectric absorption) is expected to be E(sub p) approximately (5 keV) gamma(exp -1/3) (N(sub H)/10(exp 23)/sq cm)(exp 1/3) where N(sub H) is the hydrogen column density along the line of sight. The observed spectra of GK Per and possibly of V1223 Sgr suggest N(sub H) approximately 10(exp 23)/sq cm. This large N(sub H) may be due to partially ionized preshock column material. Alternatively, we also consider absorption by the cool outer parts of an accretion disk. In this case the photoelectric absorption depth in the disk is a sensitive function of inclination. For GK Per the required inclination is approximately 83 deg. For mass accretion rates larger than a critical rate of approximately 10(exp 18) g/s, X-ray emission from the column accretion is significantly affected by radiation drag. Although the mass accretion rate increases dramatically during outbursts, the observed hard (greater than 2 keV) X-ray luminosity will not rise proportionately. The slope and peak energy of the outburst spectra are only weakly affected. We conclude that the observed X-ray spectra can be explained by this simple analytic solution and that the production of hard X-rays from the accretion shock at the magnetic poles in the intermediate polars is in general agreement with the observations. However, since the X-ray emission and absorption depend on the mass accretion rate in a complicated manner, observed hard X-ray luminosities (greater than 2 keV) are not a good indicator of the mass accretion rate in the X-ray emitting regions.

Yi, Insu

Why Do T Tauri Disks Accrete?

Observations of T Tauri stars and young brown dwarfs suggest that the accretion rates of their disks scales roughly with the square of the central stellar mass. No dependence of accretion rate on stellar mass is predicted by the simplest version of the Gammie layered disk model, in which nonthermal ionization of upper disk layers allows accretion to occur via the magnetorotational instability. We show that a minor modification of Gaminie's model to include heating by irradiation from the central star yields a modest dependence of accretion on the mass of the central star. A purely viscous disk model could provide a strong dependence of accretion rate on stellar mass if the initial disk radius (before much viscous evolution has occurred) has a strong dependence on stellar mass. However, it is far from clear that at least the most massive pre-main-sequence disks can be totally magnetically activated by X-rays or cosmic rays. We suggest that a combination of effects are responsible for the observed dependence, with the lowest mass stars having the lowest mass disks, which can be thoroughly magnetically active, while the higher mass stars have higher mass disks that have layered accret,ion and relatively inactive or "dead" central zones at some radii. In such dead zones, we suggest that gravitational instabilities may play a role in allowing accretion to proceed. In this connection, we emphasize the uncertainty in disk masses derived from dust emission and argue that T Tauri disk masses have been systematically underestimated by conventional analyses. Furtlier study of accretion rates, especially in the lowest mass stars, would help to clarify the mechanisms of accretion in T Tauri stars.

Hartmann, Lee

A Model to Assess the Risk of Ice Accretion Due to Ice Crystal Ingestion in a Turbofan Engine and its Effects on Performance

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that were attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was one or more of the following anomalies: degraded engine performance, engine roll back, compressor surge and stall, and flameout of the combustor. The main focus of this research is the development of a computational tool that can estimate whether there is a risk of ice accretion by tracking key parameters through the compression system blade rows at all engine operating points within the flight trajectory. The tool has an engine system thermodynamic cycle code, coupled with a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Assumptions are made to predict the complex physics involved in engine icing. Specifically, the code does not directly estimate ice accretion and does not have models for particle breakup or erosion. Two key parameters have been suggested as conditions that must be met at the same location for ice accretion to occur: the local wet-bulb temperature to be near freezing or below and the local melt ratio must be above 10%. These parameters were deduced from analyzing laboratory icing test data and are the criteria used to predict the possibility of ice accretion within an engine including the specific blade row where it could occur. Once the possibility of accretion is determined from these parameters, the degree of blockage due to ice accretion on the local stator vane can be estimated from an empirical model of ice growth rate and time spent at that operating point in the flight trajectory. The computational tool can be used to assess specific turbine engines to their susceptibility to ice accretion in an ice crystal environment.

Jorgenson, Philip C. E.

The structure and appearance of winds from supercritical accretion disks. I - Numerical models

Equations for the structure and appearance of supercritical accretion disks and the radiation-driven winds which emanate from them are derived and solved by a steady-state hydrodynamic computer code with a relaxation technique used in stellar structure problems. The present model takes into account the mass of the accreting star, the total accretion rate, a generalization of the disk alpha parameter which accounts for heating by processes in addition to viscosity, and the ratio of the total luminosity to the Eddington luminosity. Solutions indicate that for accretion onto a hard-surfaced star, steady, optically thick winds result for even slightly supercritical accretion, and the object will appear as a supergiant star with a high mass loss rate and a nonblackbody spectrum. Winds from black hole accretion disks are expected to depend on the form of the accretion interior to the critical radius, possibly consisting of no ejection at all, a wind similar to that of a hard-surfaced star, or a column of material ejected from a hole in the accretion disk.

Meier, D. L.

CNO abundances resulting from diffusion in accreting nova progenitors

The configurations of potential nova progenitors at the onset of a thermonuclear runaway (TNR) are obtained by evolving a series of white dwarf (WD) models through the quasi-static accretion phase preceding the TNR. Special attention is paid to the effect of diffusion on composition profiles. The initial parameters of the accretion phase are the WD mass, the WD luminosity, and the accretion rate. They lead to a wide range of combinations of accreted mass and heavy element content of the envelope Z(env). It is found that the WD mass is the most important parameter in determining the accreted mass. A lower limit exists for the mass of a WD capable of producing a nova-like outburst which decreases with decreasing WD mass. There is a strong correlation between the accretion rate and Z(env), higher Z(env) values resulting from lower accretion rates. Changes in WD luminosity have no dramatic effect on the outcome of the accretion phase.

Kovetz, A.

An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation

An investigation of the ice accretion pattern and performance characteristics of a multi-element airfoil was undertaken in the NASA Lewis 6- by 9-Foot Icing Research Tunnel. Several configurations of main airfoil, slat, and flaps were employed to examine the effects of ice accretion and provide further experimental information for code validation purposes. The text matrix consisted of glaze, rime, and mixed icing conditions. Airflow and icing cloud conditions were set to correspond to those typical of the operating environment anticipated tor a commercial transport vehicle. Results obtained included ice profile tracings, photographs of the ice accretions, and force balance measurements obtained both during the accretion process and in a post-accretion evaluation over a range of angles of attack. The tracings and photographs indicated significant accretions on the slat leading edge, in gaps between slat or flaps and the main wing, on the flap leading-edge surfaces, and on flap lower surfaces. Force measurments indicate the possibility of severe performance degradation, especially near C sub Lmax, for both light and heavy ice accretion and performance analysis codes presently in use. The LEWICE code was used to evaluate the ice accretion shape developed during one of the rime ice tests. The actual ice shape was then evaluated, using a Navier-Strokes code, for changes in performance characteristics. These predicted results were compared to the measured results and indicate very good agreement.

Potapczuk, Mark G.

An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation

An investigation of the ice accretion pattern and performance characteristics of a multi-element airfoil was undertaken in the NASA Lewis 6- by 9-Foot Icing Research Tunnel. Several configurations of main airfoil, slat, and flaps were employed to examine the effects of ice accretion and provide further experimental information for code validation purposes. The text matrix consisted of glaze, rime, and mixed icing conditions. Airflow and icing cloud conditions were set to correspond to those typical of the operating environment anticipated for a commercial transport vehicle. Results obtained included ice profile tracings, photographs of the ice accretions, and force balance measurements obtained both during the accretion process and in a past-accretion evaluation over a range of angles of attack. The tracings and photographs indicated significant accretions on the slat leading edge, in gaps between slat or flaps and the main wing, on the flap leading-edge surfaces, and on flap lower surfaces. Force measurements indicate the possibility of severe performance degradation, especially near C sub Lmax, for both light and heavy ice accretion and performance analysis codes presently in use. The LEWICE code was used to evaluate the ice accretion shape developed during one of the rime ice tests. The actual ice shape was then evaluated, using a Navier-Strokes code, for changes in performance characteristics. These predicted results were compared to the measured results and indicate very good agreement.

Potapczuk, Mark G.

Accretion in wind-driven X-ray sources

The hydrodynamics of axisymmetric accretion flow in stellar wind-fed X-ray sources is investigated, including momentum deposition by radiation, as well as radiative heating and cooling by Compton and bremsstrahlung processes. The results of two-dimensional numerical simulations for mass accretion rates ranging from 0.02 to 0.4 of the Eddington value reveal a variety of behaviors. At low rates radiative effects are unimportant, and the accretion flow is nearly adiabatic. In this case the flow is steady. For intermediate accretion rates, radiative heating and cooling effects become important. The flow remains steady, but the accretion rate is significantly reduced from the Hoyle-Lyttleton estimate. At the highest mass accretion rates, radiative momentum transfer and energy exchange are important and lead to nonsteady behavior. For the parameters relevant to massive X-ray binary systems accreting at high rates, the results reveal that the accretion flows can be complex and time dependent.

Taam, Ronald E.

Plasma physics of accreting neutron stars

Plasma concepts and phenomena that are needed to understand X- and gamma-ray sources are discussed. The capture of material from the wind or from the atmosphere or envelope of a binary companion star is described and the resulting types of accretion flows discussed. The reasons for the formation of a magnetosphere around the neutron star are explained. The qualitative features of the magnetospheres of accreting neutron stars are then described and compared with the qualitative features of the geomagnetosphere. The conditions for stable flow and for angular and linear momentum conservation are explained in the context of accretion by magnetic neutron stars and applied to obtain rough estimates of the scale of the magnetosphere. Accretion from Keplerian disks is then considered in some detail. The radial structure of geometrically thin disk flows, the interaction of disk flows with the neutron star magnetosphere, and models of steady accretion from Keplerian disks are described. Accretion torques and the resulting changes in the spin frequencies of rotating neutron stars are considered. The predicted behavior is then compared with observations of accretion-powered pulsars. Magnetospheric processes that may accelerate particles to very high energies, producing GeV and, perhaps, TeV gamma-rays are discussed. Finally, the mechanisms that decelerate and eventually stop accreting plasma at the surfaces of strongly magnetic neutron stars are described.

Ghosh, Pranab

The evidence for clumpy accretion in the Herbig Ae star HR 5999

Analysis of IUE high- and low-dispersion spectra of the young Herbig Ae star HR 5999 (HD 144668) covering 1978-1992 has revealed dramatic changes in the Mg II h and k (2795.5, 2802.7 A) emission profiles, changes in the column density and distribution in radial velocity of accreting gas, and flux in the Ly(alpha), O I and C IV emission lines, which are correlated with the UV excess luminosity. We also observe variability in the spectral type inferred from the UV spectral energy distribution, ranging from A5 IV-III in high state to A7 III in the low state. The trend of earlier inferred spectral type with decreasing wavelength and with increasing UV continuum flux has previously been noted as a signature of accretion disks in lower mass pre-main sequence stars (PMS) and in systems undergoing FU Orionis-type outbursts. Our data represent the first detection of similar phenomena in an intermediate mass (M equal to or greater than 2 solar mass) PMS star. Recent IUE spectra show gas accreting toward the star with velocities as high as +300 km/s, much as is seen toward beta Pic, and suggest that we also view this system through the debris disk. The absence of UV lines with the rotational broadening expected given the optical data (A7 IV, upsilon sin i = 180 plus or minus 20 km/s) for this system also suggests that most of the UV light originates in the disk, even in the low continuum state. The dramatic variability in the column density of accreting gas, consistent with clumpy accretion, such as has been observed toward beta Pic, is a hallmark of accretion onto young stars, and is not restricted to the clearing phase, since detectable amounts of accretion are present for stars, and is not restricted to the clearing phase, since detectable amounts of accretion are present for stars with 0.5 less than t(sub age) less than 2.8 Myr. The implications for models of beta Pic and similar systems are briefly discussed.

Perez, M. R.

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian

Observations of Accreting Pulsars

We summarize 5 years of continuous monitoring of accretion-powered pulsars with the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory. Our 20-70 keV observations have determined or refined the orbital parameters of 13 binaries, discovered five new transient accreting pulsars, measured the pulsed flux history during outbursts of 12 transients (GRO J1744-28, 4U 0115+634, GRO J1750-27, GS 0834-430, 2S 1417-624, GRO J1948+32, EXO 2030+375, GRO J1008-57, A0535+26, GRO J2058+42, 4U 1145-619, and A1118-616), and also measured the accretion torque history during outbursts of six of those transients whose orbital param- eters were also known. We have also continuously measured the pulsed flux and spin frequency for eiaht persistently accreting pulsars (Her X-1, Cen X-3, Vela X-1, OAO 1657-415, GX 301-2, 4U 1626-67, 4U 1538-52, and GX 1+4). Because of their continuity and uniformity over a long baseline, BATSE observations have provided new insights into the long-term behavior of accreting magnetic neutron stars. We have found that all accreting pulsars show stochastic variations in their spin frequencies and luminosities, including those displaying secular spin-up or spin-down on long timescales, which blurs the con- ventional distinction between disk-fed and wind-fed binaries. Pulsed flux and accretion torque are strongly correlated in outbursts of transient accreting pulsars but are uncorrelated, or even anti- correlated, in persistent sources. We describe daily folded pulse profiles, frequency, and flux measurements that are available through the Compton Observatory Science Support Center at NASA/Goddard Space Flight Center.

Bildsten, Lars