We modify an existing magnetohydrodynamics algorithm to make it more compatible with a dimensionally-split (DS) framework. It is based on the standard reconstruct-solve-average strategy (using a Riemann solver), and relies on constrained transport to ensure that the magnetic field remains divergence-free (div B = 0). The DS approach, combined with the use of a single, cell-centred grid (for both the fluid quantities and the magnetic field), means that the algorithm can be easily added to existing DS hydrodynamics codes. This makes it particularly useful for mature astrophysical codes, which often model more complicated physical effects on top of an underlying DS hydrodynamics engine, and therefore cannot be restructured easily. Several test problems have been included to demonstrate the accuracy of the algorithm, and illustrative source code has been made freely available online.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

(shortened) We perform 3D hydrodynamic simulations of gas flowing around a planetary core of mass \mplan=10\me embedded in a near Keplerian background flow, using a modified shearing box approximation. We employ a nested grid hydrodynamic code with as many as six nested grids, providing spatial resolution on the finest grid comparable to the present day diameters of Neptune and Uranus. We find that a strongly dynamically active flow develops such that no static envelope can form. The activity is not sensitive to plausible variations in the rotation curve of the underlying disk. It is sensitive to the thermodynamic treatment of the gas, as modeled by prescribed equations of state (either `locally isothermal' or `locally isentropic') and the temperature of the background disk material. The activity is also sensitive to the shape and depth of the core's gravitational potential, through its mass and gravitational softening coefficient. The varying flow pattern gives rise to large, irregular eruptions of matter from the region around the core which return matter to the background flow: mass in the envelope at one time may not be found in the envelope at any later time. The angular momentum of material in the envelope, relative to the core, varies both in magnitude and in sign on time scales of days to months near the core and on time scales a few years at distances comparable to the Hill radius. We show that material entering the dynamically active environment may suffer intense heating and cooling events the durations of which are as short as a few hours to a few days. Peak temperatures in these events range from $T \sim 1000$ K to as high as $T \sim 3-4000$ K, with densities $\rho\sim 10^{-9}-10^{-8}$ g/cm$^3$. These time scales, densities and temperatures span a range consistent with those required for chondrule formation in the nebular shock model.

It has been known since the 1950's that the observed gas content of Galactic globular clusters (GCs) is 2-3 orders of magnitude less than the mass lost by stars between Galactic disk crossings. In this work we address the question: What happens to this stellar gas? Using an Eulerian nested grid code, we present 3D simulations to determine how stellar wind material evolves within the GC environment. We expand upon work done in the 70's and move a single-mass King-model GC through the Galactic halo medium, stripping a 10^5 Msun GC of its intra-cluster medium but predicting a detectable medium for a 10^6 Msun cluster. We find from new multi-mass King model simulations, the first to incorporate empirical mass-loss formulae, that the single-mass King model underestimates the retention of intra-cluster gas in the cluster. Lastly, we present a simple discretised multi-mass GC model, which yields lower levels of intra-cluster medium compared to the continuous single- and multi-mass King models. Our results show that there is still an issue with the predicted intra-cluster gas content of massive GCs. We conclude that by modelling GC systems more accurately, in particular the stellar structure and description of mass loss, we will be able to work towards resolving this issue and begin to fill in some of the gaps in our understanding of the evolution of globular clusters.

Context: Mergers of neutron stars (NS) and black holes (BH) are among the strongest sources of gravitational waves and are potential central engines for short gamma-ray bursts. Aims: We aim to compare the general relativistic (GR) results by other groups with Newtonian calculations of models with equivalent parameters. We vary the mass ratios between NS and BH and the compactness of the NS. The mass of the NS is 1.4 M_sol. We compare the dynamics in the parameter-space regions where the NS is expected to reach the innermost stable circular orbit (ISCO) before being tidally disrupted (mass shedding, MS) and vice versa. Methods: The hydrodynamics is evolved by a Newtonian PPM scheme with four levels of nested grids. We use a polytropic EoS (Gamma=2), as was done in the GR simulations. However, instead of full GR we use a Newtonian potential supplemented by a Paczynski-Wiita-Artemova potential for the BH, both disregarding and including rotation of the BH. Results: If the NS is compact (C=0.18) it is accreted by the BH more quickly, and only a small amount of mass remains outside the BH. If the mass ratio is small (Q=2 or 3) or the NS is less compact (C=0.16 or less) the NS is tidally torn apart before being accreted. Although most of the mass is absorbed by the BH, some 0.1 M_sol remain in a tidal arm. For small mass ratios the tidal arm can wrap around the BH to form a thick disk. When including the effects of BH spin-up or spin-down by the accreted matter, more mass remains in the surroundings (0.2-0.3 M_sol). Conclusions: Although details and quantitative results differ, the general trends of our Newtonian calculations are similar to the GR calculations. A clear delimiting line that separates ISCO from the MS cases is not found. Inclusion of BH rotation as well as sufficient numerical resolution are extremely important.

We investigate the efficiency and time-dependence of thermally and cosmic ray driven galactic winds for the metal enrichment of the intra-cluster medium (ICM) using a new analytical approximation for the mass outflow. The spatial distribution of the metals are studied using radial metallicity profiles and 2D metallicity maps of the model clusters as they would be observed by X-ray telescopes like XMM-Newton. Analytical approximations for the mass loss by galactic winds driven by thermal and cosmic ray pressure are derived from the Bernoulli equation and implemented in combined N-body/hydrodynamic cosmological simulations with a semi-analytical galaxy formation model. Observable quantities like the mean metallicity, metallicity profiles, and 2D metal maps of the model clusters are derived from the simulations. We find that galactic winds alone cannot account for the observed metallicity of the ICM. At redshift $z=0$ the model clusters have metallicities originating from galactic winds which are almost a factor of 10 lower than the observed values. For massive, relaxed clusters we find, as in previous studies, a central drop in the metallicity due to a suppression of the galactic winds by the pressure of the ambient ICM. Combining ram-pressure stripping and galactic winds we find radial metallicity profiles of the model clusters which agree qualitatively with observed profiles. Only in the inner parts of massive clusters the observed profiles are steeper than in the simulations. Also the combination of galactic winds and ram-pressure stripping yields too low values for the ICM metallicities. The slope of the redshift evolution of the mean metallicity in the simulations agrees reasonably well with recent observations.

We consider hot accretion disk outflows from black hole - neutron star mergers in the context of the nucleosynthesis they produce. We begin with a three dimensional numerical model of a black hole - neutron star merger and calculate the neutrino and antineutrino fluxes emitted from the resulting accretion disk. We then follow the element synthesis in material outflowing the disk along parameterized trajectories. We find that at least a weak r-process is produced, and in some cases a main r-process as well. The neutron-rich conditions required for this production of r-process nuclei stem directly from the interactions of the neutrinos emitted by the disk with the free neutrons and protons in the outflow.

We investigate the efficiency of galactic mass loss, triggered by ram-pressure stripping and galactic winds of cluster galaxies, on the chemical enrichment of the intra-cluster medium (ICM). We combine N-body and hydrodynamic simulations with a semi-numerical galaxy formation model. By including simultaneously different enrichment processes, namely ram-pressure stripping and galactic winds, in galaxy-cluster simulations, we are able to reproduce the observed metal distribution in the ICM. We find that the mass loss by galactic winds in the redshift regime z>2 is ~10% to 20% of the total galactic wind mass loss, whereas the mass loss by ram-pressure stripping in the same epoch is up to 5% of the total ram-pressure stripping mass loss over the whole simulation time. In the cluster formation epochs z<2 ram-pressure stripping becomes more dominant than galactic winds. We discuss the non-correlation between the evolution of the mean metallicity of galaxy clusters and the galactic mass losses. For comparison with observations we present two dimensional maps of the ICM quantities and radial metallicity profiles. The shape of the observed profiles is well reproduced by the simulations in the case of merging systems. In the case of cool-core clusters the slope of the observed profiles are reproduced by the simulation at radii below ~300 kpc, whereas at larger radii the observed profiles are shallower. We confirm the inhomogeneous metal distribution in the ICM found in observations. To study the robustness of our results, we investigate two different descriptions for the enrichment process interaction.

We investigate the long-term evolution and observability of remnants originating from the merger of compact binary systems and discuss the differences to supernova remnants. Compact binary mergers expel much smaller amounts of mass at much higher velocities, as compared to supernovae, which will affect the dynamical evolution of their remnants. The ejecta of mergers consist of very neutron rich nuclei. Some of these neutron rich nuclei will produce observational signatures in form of gamma ray lines during their decay. The composition of the ejecta might even give interesting constraints about the internal structure of the neutron star. We further discuss the possibility that merger remnants appear as recently discovered 'dark accelerators' which are extended TeV sources which lack emission in other bands.

We assess the importance of AGN outflows with respect to the metal enrichment of the intracluster medium (ICM) in galaxy clusters. We use combined N-body and hydrodynamic simulations, along with a semi-numerical galaxy formation and evolution model. Using assumptions based on observations, we attribute outflows of metal-rich gas initiated by AGN activity to a certain fraction of our model galaxies. The gas is added to the model ICM, where the evolution of the metallicity distribution is calculated by the hydrodynamic simulations. For the parameters describing the AGN content of clusters and their outflow properties, we use the observationally most favorable values. We find that AGNs have the potential to contribute significantly to the metal content of the ICM or even explain the complete abundance, which is typically ~0.5 Z_sun in core regions. Furthermore, the metals end up being inhomogeneously distributed, in accordance with observations.

(abridged) Three-dimensional hydrodynamic simulations with an Eulerian PPM code are presented for the time-dependent evolution of accretion tori around nonrotating and rotating stellar-mass black holes (BHs), using a pseudo-Newtonian gravitational potential to approximate the effects of general relativity. The initial configurations are assumed to be remnants of binary neutron star (NS) or NS+BH mergers and consist of a 4 solar mass BH with varied spin, girded by a torus with a mass between 0.01 and 0.2 solar masses. The evolution of tori without and with physical shear viscosity is simulated, using a realistic equation of state and following the energy loss and lepton number change due to neutrino emission by a neutrino-trapping scheme. The time-dependent efficiency of converting rest-mass energy to neutrinos is found to reach 10 percent, the efficiency of converting neutrino energy into a pair-photon fireball by neutrino annihilation can reach several percent. The rate of the latter process declines with time much less steeply than the total neutrino luminosity, because the ongoing protonization of the torus ensures a rather stable product of neutrino and antineutrino luminosities. The neutrino emission increases steeply with higher viscosity, larger torus mass, and larger BH spin in corotation with the torus. For rotation rates as expected for post-merger BHs (a > 0.5) and reasonable values of the alpha viscosity (alpha ~ 0.1), the considered tori release sufficient energy in neutrinos to account for the energetics of the well-localized short gamma-ray bursts recently detected by Hete and Swift, if collimation is invoked as predicted by hydrodynamic jet simulations.

We investigate the impact of galactic mass loss triggered by ram-pressure stripping of cluster galaxies on the evolution of the intra-cluster medium (ICM). We use combined N-body and hydrodynamic simulations together with a phenomenological galaxy formation model and a prescription of the effect of ram-pressure stripping on the galaxies. We analyze the effect of galaxy -- ICM interaction for different model clusters with different masses and different merger histories. Our simulations show that ram-pressure stripping can account for ~ 10% of the overall observed level of enrichment in the ICM within a radius of 1.3 Mpc. The efficiency of metal ejection of cluster galaxies depends at the first few Gyr of the simulation mainly on the cluster mass and is significantly increased during major merger events. Additionally we show that ram-pressure stripping is most efficient in the center of the galaxy cluster and the level of enrichment drops quite fast at larger radii. We present emission weighted metallicity maps of the ICM which can be compared with X-ray observations. The resulting distribution of metals in the ICM shows a complex pattern with stripes and plumes of metal rich material. The metallicity maps can be used to trace the present and past interactions between the ICM and cluster galaxies.

We present an investigation of the metal enrichment of the intra-cluster medium (ICM) by galactic winds and merger-driven starbursts. We use combined N-body/hydrodynamic simulations with a semi-numerical galaxy formation model. The mass loss by galactic winds is obtained by calculating transonic solutions of steady state outflows, driven by thermal, cosmic ray and MHD wave pressure. The inhomogeneities in the metal distribution caused by these processes are an ideal tool to reveal the dynamical state of a galaxy cluster. We present surface brightness, X-ray emission weighted temperature and metal maps of our model clusters as they would be observed by X-ray telescopes like XMM-Newton. We show that X-ray weighted metal maps distinguish between pre- or post-merger galaxy clusters by comparing the metallicity distribution with the galaxy-density distribution: pre-mergers have a metallicity gap between the subclusters, post-mergers a high metallicity between subclusters. We apply our approach to two observed galaxy clusters, Abell 3528 and Abell 3921, to show whether they are pre- or post-merging systems. The survival time of the inhomogeneities in the metallicity distribution found in our simulations is up to several Gyr. We show that galactic winds and merger-driven starbursts enrich the ICM very efficiently after z=1 in the central (~ 3 Mpc radius) region of a galaxy cluster.

We investigate the long-term evolution and observability of remnants originating from the merger of compact binary systems and discuss the differences to supernova remnants. Compact binary mergers expel much smaller amounts of mass at much higher velocities, as compared to supernovae and therefore the free expansion phase of the remnant will be short (~ 1 - 10 yr). In general the remnants will be observable for a considerable time (~ 10^6 - 10^7 yr). Events releasing large amounts of kinetic energy may be responsible for a subsample of observed giant HI holes of unknown origin as compact binaries merge far away from star forming regions. If the ejecta consist primarily of actinides, on long timescales the expelled material will contain mainly the few quasi-stable nuclei in the actinides range. Consequently the abundance of each isotope in the ejecta might be of the order of a few percent. During their decay some actinides will produce observational signatures in form of gamma ray lines. We particularly investigate the gamma ray emission of Am 243, Cm 247, Cm 248 and Bi 208 and estimate their observability in nearby remnants. Detections of the gamma ray lines with INTEGRAL will be possible only in very advantageous cases but these remnants are promising targets for future instruments using focusing optics for soft gamma rays. Due to the low mass expelled in mergers and due to the lack of free electrons in the ejecta, the merger remnants might be significantly fainter in bremsstrahlung and synchrotron radiation than comparable supernova remnants. Hence merger remnants might represent a candidate for very recently discovered 'dark accelerators' which are hard gamma ray sources with no apparent emission in other bands.

We propose that the nebular shocks currently favored as a model to form chondrules and other annealed silicates in the solar nebula originate in the dynamical activity present in the envelope of forming Jovian planets. In contrast to the classic `core accretion model', our 3D hydrodynamic simulations show that this envelope is not a 1D hydrostatic structure but is instead vigorously active and contains densities and temperatures that appear similar in magnitude and spatial extent to those thought to be responsible for the production of chondrules.

We present numerical simulations of galaxy clusters which include interaction processes between the galaxies and the intra-cluster gas. The considered interaction processes are galactic winds and ram-pressure stripping, which both transfer metal-enriched interstellar medium into the intra-cluster gas and hence increase its metallicity. We investigate the efficiency and time evolution of the interaction processes by simulated metallicity maps, which are directly comparable to those obtained from X-ray observations. We find that ram-pressure stripping is more efficient than quiet (i.e. non-starburst driven) galactic winds in the redshift interval between 1 and 0. The expelled metals are not mixed immediately with the intra-cluster gas, but inhomogeneities are visible in the metallicity maps. Even stripes of higher metallicity that a single galaxy has left behind can be seen. The spatial distribution of the metals transported by ram-pressure stripping and by galactic winds are very different for massive clusters: the former process yields a centrally concentrated metal distribution while the latter results in an extended metal distribution.

The instability of Bondi-Hoyle-Lyttleton accretion, observed in numerical simulations, is analyzed through known physical mechanisms and possible numerical artefacts. The mechanisms of the longitudinal and transverse instabilities, established within the accretion line model, are clarified. They cannot account for the instability of BHL accretion at moderate Mach number when the pressure forces within the shock cone are taken into account. The advective-acoustic instability is considered in the context of BHL accretion when the shock is detached from the accretor. This mechanism naturally explains the stability of the flow when the shock is weak, and the instability when the accretor is small. In particular, it is a robust proof of the instability of 3D accretion when gamma=5/3 if the accretor is small enough, even for moderate shock strength (M sim 3). The numerical artefacts that may be present in existing numerical simulations are reviewed, with particular attention paid to the advection of entropy/vorticity perturbations and the artificial acoustic feedback from the accretor boundary condition. Several numerical tests are proposed to test these mechanisms.

We present numerical simulations of the dynamical and chemical evolution of galaxy clusters. X-ray spectra show that the intra-cluster medium contains a significant amount of metals. As heavy elements are produced in the stars of galaxies material from the galaxies must have been expelled to enrich the ambient medium. We have performed hydrodynamic simulations investigating various processes. In this presentation we show the feedback from gas which is stripped from galaxies by ram-pressure stripping. The efficiency, resulting spatial distribution of the metals and the time dependency of this enrichment process on galaxy cluster scale is shown.

We present three-dimensional hydrodynamic simulations of the evolution of selfgravitating, thick accretion discs around hyperaccreting stellar-mass black holes. The black hole-torus systems are considered to be remnants of compact object mergers, in which case the disc is not fed by an external mass reservoir and the accretion is non-stationary. Our models take into account viscous dissipation, described by an alpha-law, a detailed equation of state for the disc gas, and an approximate treatment of general relativistic effects on the disc structure by using a pseudo-Newtonian potential for the black hole including its possible rotation and spin-up during accretion. Magnetic fields are ignored. The neutrino emission of the hot disc is treated by a neutrino-trapping scheme, and the neutrino-antineutrino annihilation near the disc is evaluated in a post-processing step. Our simulations show that the neutrino emission and energy deposition by neutrino-antineutrino annihilation increase sensitively with the disc mass, with the black hole spin in case of a disc in corotation, and in particular with the alpha-viscosity. We find that for sufficiently large alpha-viscosity neutrino-antineutrino annihilation can be a viable energy source for gamma-ray bursts.

(Abridged) In this paper we present a compilation of results from our most advanced neutron star merger simulations, including a description of the employed numerical procedures and a more complete overview over a large number of computed models. The three-dimensional hydrodynamic simulations were done with a code based on the Piecewise Parabolic Method with up to five levels of nested Cartesian grids. The simulations are basically Newtonian, but gravitational-wave emission and the corresponding back-reaction are taken into account. The use of a physical nuclear equation of state allows us to follow the thermodynamic history of the stellar medium and to compute the energy and lepton number loss due to the emission of neutrinos. The computed models differ concerning the neutron star masses and mass ratios, the neutron star spins, the numerical resolution expressed by the cell size of the finest grid and the number of grid levels, and the calculation of the temperature from the solution of the entropy equation instead of the energy equation. Our simulations show that the details of the gravitational-wave emission are still sensitive to the numerical resolution, even in our highest-quality calculations. The amount of mass which can be ejected from neutron star mergers depends strongly on the angular momentum of the system. Our results do not support the initial conditions of temperature and proton-to-nucleon ratio assumed in recent work for producing a solar r-process pattern for nuclei around and above the A approx 130 peak. The improved models confirm our previous conclusion that gamma-ray bursts are not powered by neutrino emission during the dynamical phase of the merging of two neutron stars.

Merging neutron stars and neutron star-black holes binaries are powerful sources of gravitational waves. They have also been suggested as possible sources of cosmic gamma-ray bursts and are discussed as sites for the formation of r-process elements. Whereas the first aspect is undoubtable, the latter two are very uncertain. The current status of our knowledge and of the numerical modeling is briefly reviewed and the results of simulations are critically discussed concerning their significance and implications for potentially observable signals.

Hydrodynamic simulations of the merger of stellar mass black hole - neutron star binaries (BH/NS) are compared with mergers of binary neutron stars (NS/NS). The simulations are Newtonian, but take into account the emission and backreaction of gravitational waves. The use of a physical nuclear equation of state allows us to include the effects of neutrino emission. For low neutron star to black hole mass ratios the neutron star transfers mass to the black hole during a few cycles of orbital decay and subsequent widening before finally being disrupted, whereas for ratios near unity the neutron star is already distroyed during its first approach. A gas mass between about 0.3 and about 0.7 solar masses is left in an accretion torus around the black hole and radiates neutrinos at a luminosity of several 10^{53} erg/s during an estimated accretion time scale of about 0.1 s. The emitted neutrinos and antineutrinos annihilate into electron-positron pairs with efficiencies of 1-3% percent and rates of up to 2*10^{52} erg/s, thus depositing an energy of up to 10^{51} erg above the poles of the black hole in a region which contains less than 10^{-5} solar masses of baryonic matter. This could allow for relativistic expansion with Lorentz factors around 100 and is sufficient to explain apparent burst luminosities of up to several 10^{53} erg/s for burst durations of approximately 0.1-1 s, if the gamma emission is collimated in two moderately focussed jets in a fraction of about 1/100-1/10 of the sky.

The adiabatic shock produced by a compact object moving supersonically relative to a gas with uniform entropy and no vorticity is a source of entropy gradients and vorticity. We investigate these analytically. The non-axisymmetric Rayleigh-Taylor and axisymmetric Kelvin-Helmholtz linear instabilities are potential sources of destabilization of the subsonic accretion flow after the shock. A local Lagrangian approach is used in order to evaluate the efficiency of these linear instabilities. However, the conditions required for such a WKB type approximation are fulfilled only marginally: a quantitative estimate of their local growth rate integrated along a flow line shows that their growth time is at best comparable to the time needed for advection onto the accretor, even at high Mach number and for a small accretor size. Despite this apparently low efficiency, several features of these mechanisms qualitatively match those observed in numerical simulations: in a gas with uniform entropy, the instability occurs only for supersonic accretors. It is nonaxisymmetric, and begins close to the accretor in the equatorial region perpendicular to the symmetry axis. The mechanism is more efficient for a small, highly supersonic accretor, and also if the shock is detached. We also show by a 3-D numerical simulation an example of unstable accretion of a subsonic flow with non-uniform entropy at infinity. This instability is qualitatively similar to the one observed in 3-D simulations of the Bondi-Hoyle-Lyttleton flow, although it involves neither a bow shock nor an accretion line.

The hydrodynamics of a variant of classical Bondi-Hoyle-Lyttleton accretion is investigated: a totally absorbing sphere moves at various Mach numbers (3 and 10) relative to a medium, which is taken to be an ideal gas having a density gradient (of 3%, 20% or 100% over one accretion radius) perpendicular to the relative motion. Similarly to the 3D models published previously, both with velocity gradients and without, the models with a density gradient presented here exhibit non-stationary flow patterns, although the Mach cone remains fairly stable. The accretion rates of mass, linear and angular momenta do not fluctuate as strongly as published previously for 2D models. No obvious trend of the dependency of mass accretion rate fluctuations on the density gradient can be discerned. The average specific angular momentum accreted is roughly between zero and 70% of the total angular momentum available in the accretion cylinder in the cases where the average is prograde. Due to the large fluctuations during accretion, the average angular momentum of some models is retrograde by up to 25%. Small gradients hardly influence the average accretion rates as compared to accretion from a homogeneous medium, while very large ones succeed to dominate and form an accretion flow in which the sense of rotation is not inverted.

Hydrodynamic simulations were performed of the dynamical phase of the merging of binary neutron stars (NS-NS) and of neutron star black hole binaries (NS-BH), using the physical nuclear equation of state of Lattimer and Swesty (1991) and taking into account the emission of gravitational waves and neutrinos.

By means of three-dimensional hydrodynamic simulations with a Eulerian PPM code we investigate the formation and the properties of the accretion torus around the stellar mass black hole which originates from the merging of two neutron stars. The simulations are performed with four nested cartesian grids which allow for both a good resolution near the central black hole and a large computational volume. They include the use of a physical equation of state as well as the neutrino emission from the hot matter of the torus. The gravity of the black hole is described with a Newtonian and alternatively with a Paczynski-Wiita potential. In a post-processing step, we evaluate our models for the energy deposition by nu-nubar annihilation around the accretion torus. Our models show that nu-nubar annihilation can yield the energy to account for weak, short gamma-ray bursts, if moderate beaming is involved. In fact, the barrier of the dense baryonic gas of the torus suggests that the low-density pair-photon-plasma is beamed as axial jets into a fraction 2 delta Omega/ (4 pi) between 1/100 and 1/10 of the sky, corresponding to opening half-angles of roughly ten to several tens of degrees. Thus gamma-burst energies of 10^{50}--10^{51} erg seem within the reach of our models (if the source is interpreted as radiating isotropically), corresponding to luminosities around 10^{51} erg/s for typical burst durations of 0.1--1 s. Gravitational capture of radiation by the black hole, redshift and ray bending do not reduce the jet energy significantly. Effects associated with the Kerr character of the rapidly rotating black hole, however, could increase the gamma-burst energy considerably, and effects due to magnetic fields might even be required to get the energies for long complex gamma-ray bursts.

Three-dimensional hydrodynamical simulations are presented for the direct head-on or off-center collision of two neutron stars, employing a basically Newtonian PPM code but including the emission of gravitational waves and their back-reaction on the hydrodynamical flow. A physical nuclear equation of state is used that allows us to follow the thermodynamical evolution of the stellar matter and to compute the emission of neutrinos. Predicted gravitational wave signals, luminosities and waveforms, are presented. The models are evaluated for their implications for gamma-ray burst scenarios. We find an extremely luminous outburst of neutrinos with a peak luminosity of more than 4E54 erg/s for several milliseconds. This leads to an efficiency of about 1% for the annihilation of neutrinos with antineutrinos, corresponding to an average energy deposition rate of more than 1E52 erg/s and a total energy of about 1E50 erg deposited in electron-positron pairs around the collision site within 10ms. Although these numbers seem very favorable for gamma-ray burst scenarios, the pollution of the $e^\pm$ pair-plasma cloud with nearly 0.1$M_{\odot}$ of dynamically ejected baryons is 5 orders of magnitude too large. Therefore the formation of a relativistically expanding fireball that leads to a gamma-ray burst powered by neutrino emission from colliding neutron stars is definitely ruled out.

We report on the first three-dimensional hydrodynamic simulations of secular and dynamical non-axisymmetric instabilities in collapsing, rapidly rotating stellar cores which extend well beyond core bounce. The resulting gravitational radiation has been calculated using the quadrupole approximation. We find that secular instabilities do not occur during the simulated time interval of several 10 ms. Models which become dynamically unstable during core collapse show a strong nonlinear growth of non-axisymmetric instabilities. Both random and coherent large scale initial perturbations eventually give rise to a dominant bar-like deformation ($\exp(\pm im\phi$) with $m=2$). In spite of the pronounced tri-axial deformation of certain parts of the core no considerable enhancement of the gravitational radiation is found. This is due to the fact that rapidly rotating cores re-expand after core bounce on a dynamical time scale before non-axisymmetric instabilities enter the nonlinear regime. Hence, when the core becomes tri-axial, it is no longer very compact.

The dynamics, time evolution of the mass distribution, and gravitational wave signature of coalescing neutron stars described by polytropes are compared with three simulations published previously: (a) ``Run 2'' of Zhuge et al. (1994), (b) ``Model III'' of Shibata et al. (1992), and (c) ``Model A64'' of Ruffert et al. (1996). We aim at studying the differences due to the use of different numerical methods, different implementations of the gravitational wave backreaction, and different equations of state. Comparison (a) confronts the results of our grid-based PPM scheme with those from an SPH code. We found that due to the lower numerical viscosity of the PPM code, the post-merging oscillations and pulsations could be followed for a longer time and lead to larger secondary and tertiary maxima of the gravitational wave luminosity. In case (b) two grid based codes with the same backreaction formalism but differing hydrodynamic integrators and different numerical resolution are compared. Satisfactory agreement of the amplitude of the gravitational wave luminosity is established, although due to the different initial conditions a small time delay develops in the onset of the dynamical instability. In (c) we find that using a polytropic equation of state instead of the high-density equation of state of Lattimer & Swesty (1991) does not change the overall dynamical evolution of the merger and yields agreement of the gravitational wave signature to within 20% accuracy. However, differences of the structure and evolution of the outer layers of the neutron stars are present, which has important implications for questions like mass loss and disk formation during the merging of binary neutron stars.

Three-dimensional hydrodynamical, Newtonian calculations of the coalescence of equal-mass binary neutron stars are performed, including a physical high-density equation of state and a treatment of the neutrino emission of the heated matter. The total neutrino luminosity climbs to a maximum value of 1--$1.5\cdot 10^{53}$~erg/s of which 90--95\% originate from the toroidal gas cloud surrounding the very dense core formed after the merging. When the neutrino luminosities are highest, $\nu\bar\nu$-annihilation deposits about 0.2--0.3\% of the emitted neutrino energy in the immediate neighborhood of the merger, and the maximum integral energy deposition rate is 3--$4\cdot 10^{50}$~erg/s. Since the $3\,M_{\odot}$ core of the merged object will most likely collapse into a black hole within milliseconds, the energy that can be pumped into a pair-photon fireball is insufficient by a factor of about 1000 to explain $\gamma$-ray bursts at cosmological distances with an energy of the order of $10^{51}/(4\pi)$~erg/steradian. Analytical estimates show that the additional energy provided by the annihilation of $\nu\bar\nu$ pairs emitted from a possible accretion torus of $\sim 0.1\,M_{\odot}$ around the central black hole is still more than a factor of 10 too small, unless focussing of the fireball into a jet-like expansion plays an important role. About $10^{-4}$--$10^{-3}$~$M_\odot$ of material lost during the neutron star merging and swept out from the system in a neutrino-driven wind might be a site for nucleosythesis. Aspects of a possible r-processing in these ejecta are discussed.

We investigate the hydrodynamics of a variant of classical Bondi-Hoyle-Lyttleton accretion: a totally absorbing sphere moves at various Mach numbers (3 and 10) relative to a medium, which is taken to be an ideal gas having a velocity gradient (of 3% or 20% over one accretion radius) perpendicular to the relative motion. We examine the influence of the Mach number of the flow and the strength of the gradient upon the physical behaviour of the flow and the accretion rates of the angular momentum in particular. The hydrodynamics is modeled by the ``Piecewise Parabolic Method'' (PPM). The resolution in the vicinity of the accretor is increased by multiply nesting several grids around the sphere. Similarly to the 3D models without gradients published previously, models exhibit non-stationary flow patterns, although the Mach cone remains fairly stable. The accretion rates of mass, linear and angular momenta do not fluctuate as strongly as published previously for 2D models, but similarly to the 2D models, transient disks form around the accretor that alternate their direction of rotation with time. The average specific angular momentum accreted is roughly between 7% and 70% of the total angular momentum available in the accretion cylinder and is always smaller than the value of a vortex with Kepler velocity around the surface of the accretor. The fluctuations of the mass accretion rate in the models with small gradients (2%) are similar to the values of the models without gradients, while the models with large gradients (20%) exhibit larger fluctuations. The mass accretion rate is maximal when the specific angular momentum is zero, while the specific entropy tends to be smaller when the disks are prograde.

We prove that the sonic surface of axisymmetric meridional stationary flows is always attached to the accretor, however small, if the adiabatic index of the gas is gamma=5/3. Using local expansions near a point-like accretor, we extend Bondi's classification of spherically symmetric flows to axisymmetric flows, introducing the possibility of angular sectors reached by no flow lines, and singular directions of infinite mass flux, in addition to the angular regions of subsonic and supersonic accretion. For gamma<5/3, we show the impossibility of subsonic accretion onto a point-like accretor when the entropy of the flow is not uniform. The special case gamma=5/3 is treated separately. We analyse the influence of the adiabatic index and Mach number of the flow at infinity on the mass accretion rate of shocked spherical flows. We propose an interpolation formula for the mass accretion rate of axisymmetric flows as a function of the Mach number and the adiabatic index, in the range 9/7<gamma<5/3.

We perform three-dimensional hydrodynamical simulations of the coalescence of binary neutron stars. We include the emission and backreaction of gravitational waves into the Newtonian ``Piecewise Parabolic Method''. The use of the physical equation of state of Lattimer & Swesty (1991) allows us to take into account the production of neutrinos. We evaluate our models for the efficiency of neutrino-antineutrino annihilation in the surroundings of the coalescing neutron stars. The corresponding energy deposition prior to and during merging turns out to be 2--3 orders of magnitude too small to power a typical gamma-ray burst with an energy output of about 10^{51}/(4 pi) erg/sterad at cosmological distances. Analytical estimates of the subsequent evolution of the disk which possibly surrounds the central black hole show that even under the most favorable conditions the energy in a fireball of electron-positron pairs and photons that is created by neutrino-antineutrino annihilation falls short by at least an order of magnitude.

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. A totally absorbing sphere of different sizes (1, 0.1 and 0.02 accretion radii) moves at different Mach numbers (0.6, 1.4, 3.0 and 10) relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal, nearly isothermal, gas ($\gamma=1.01$). The hydrodynamics is modeled by the ``Piecewise Parabolic Method'' (PPM). The resolution in the vicinity of the accretor is increased by multiply nesting several $32^3$-zone grids around the sphere, each finer grid being a factor of two smaller in zone size than the next coarser grid. grids. For small Mach numbers (0.6 and~1.4) the flow patterns tend towards a steady state, while in the case of supersonic flow (Mach~3 and~10) and small enough accretors (radius of~0.1 and~0.02 accretion radii), an unstable Mach cone develops, destroying axisymmetry. The shock cones in the supersonic models never clear the surface of the accretors (they are tail shocks, not bow shocks) and the opening angle is smaller (compared to models with larger $\gamma$) especially for the highly supersonic models. The densities in the shock cone is larger for models with smaller $\gamma$. The fluctuations of the accretion rates and flow structures are weaker than in the corresponding models with larger $\gamma$. The hydrodynamic drag of all models with accretor sizes of 0.1~$R_{\rm A}$ or smaller acts in an accelerating direction, while the gravitational drag is always decelerating and larger than the hydrodynamic drag (thus the net force is decelerating).

We investigate the dynamics and evolution of coalescing neutron stars. Although the code (Piecewise Parabolic Method) is purely Newtonian, we do include the emission of gravitational waves and their backreaction on the hydrodynamic flow. The properties of neutron star matter are described by the physical equation of state of Lattimer \& Swesty (1991). Energy loss by all types of neutrinos and changes of the electron fraction due to the emission of electron neutrinos and antineutrinos are taken into account by an elaborate ``neutrino leakage scheme''. We simulate the coalescence of two identical, cool neutron stars with a baryonic mass of $\approx\!1.6\,M_\odot$ and a radius of $\approx\!15$~km and with an initial center-to-center distance of 42~km. The initial distributions of density and electron concentration are given from a model of a cold neutron star in hydrostatic equilibrium (central temperature about $8\,{\rm MeV}$). We investigate three cases which differ by the initial velocity distribution in the neutron stars, representing different cases of the neutron star spins relative to the direction of the orbital angular momentum vector. Within about 1~ms the neutron stars merge into a rapidly spinning ($P_{\rm spin}\approx 1$~ms), high-density body ($\rho\approx 10^{14}$~g/cm$^3$) with a surrounding thick disk of material with densities $\rho\approx 10^{10}-10^{12}$~g/cm$^3$ and orbital velocities of~0.3--0.5~c. In this work we evaluate the models in detail with respect to the gravitational wave emission using the quadrupole approximation. In a forthcoming paper we will concentrate on the neutrino emission and implications for gamma-ray bursters. A maximum luminosity in excess of $10^{55}$~erg/s is reached for about 1~ms.

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. A totally absorbing sphere of different sizes (1, 0.1 and 0.02 accretion radii) exerts gravity on and moves at different Mach numbers (0.6, 1.4, 3.0 and 10) relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas ($\gamma=4/3$). We examine the influence of Mach number of the flow and size of the accretor upon the physical behaviour of the flow and the accretion rates. The hydrodynamics is modeled by the ``Piecewise Parabolic Method'' (PPM). The resolution in the vicinity of the accretor is increased by multiply nesting several $32^3$-zone grids around the sphere, each finer grid being a factor of two smaller in zone dimension than the next coarser grid. 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. For small Mach numbers (0.6 and~1.4) the flow patterns tend towards a steady state, while in the case of supersonic flow (Mach~3 and~10) and small enough accretors, (radius of~0.1 and~0.02 accretion radii) an unstable Mach cone develops, destroying axisymmetry. Our 3D models do not show the highly dynamic flip-flop flow so prominent in 2D calculations performed by other authors. In the gamma=4/3 models, the shock front remains closer to the accretor and the mass accretion rates are higher than in the gamma=5/3 models, whereas the rms of the specific angular momentum accreted does not change.

We investigate the dynamics and evolution of coalescing neutron stars. The three-dimensional Newtonian equations of hydrodynamics are integrated by the `Piecewise Parabolic Method' However, we do include the effects of the emission of gravitational waves on the hydrodynamics. The properties of neutron star matter are described by the equation of state of Lattimer & Swesty. In addition to the fundamental hydrodynamic quantities, density, momentum, and energy, we follow the time evolution of the electron density in the stellar gas. Energy losses and changes of the electron abundance due to the emission of neutrinos are taken into account by an elaborate ``neutrino leakage scheme'', which employs a careful calculation of the lepton number and energy source terms of all neutrino types. The grid is Cartesian and equidistant with a resolution of 64**3 or 128**3, which allows us to calculate the self-gravity via fast Fourier transforms.