12,085 research outputs found
Hydra: An Adaptive--Mesh Implementation of PPPM--SPH
We present an implementation of Smoothed Particle Hydrodynamics (SPH) in an
adaptive-mesh PPPM algorithm. The code evolves a mixture of purely
gravitational particles and gas particles. The code retains the desirable
properties of previous PPPM--SPH implementations; speed under light clustering,
naturally periodic boundary conditions and accurate pairwise forces. Under
heavy clustering the cycle time of the new code is only 2--3 times slower than
for a uniform particle distribution, overcoming the principal disadvantage of
previous implementations\dash a dramatic loss of efficiency as clustering
develops. A 1000 step simulation with 65,536 particles (half dark, half gas)
runs in one day on a Sun Sparc10 workstation. The choice of time integration
scheme is investigated in detail. A simple single-step Predictor--Corrector
type integrator is most efficient. A method for generating an initial
distribution of particles by allowing a a uniform temperature gas of SPH
particles to relax within a periodic box is presented. The average SPH density
that results varies by \%. We present a modified form of the
Layzer--Irvine equation which includes the thermal contribution of the gas
together with radiative cooling. Tests of sound waves, shocks, spherical infall
and collapse are presented. Appropriate timestep constraints sufficient to
ensure both energy and entropy conservation are discussed. A cluster
simulation, repeating Thomas andComment: 29 pp, uuencoded Postscrip
Smoothed Particle Hydrodynamics in cosmology: a comparative study of implementations
We analyse the performance of twelve different implementations of Smoothed
Particle Hydrodynamics (SPH) using seven tests designed to isolate key
hydrodynamic elements of cosmological simulations which are known to cause the
SPH algorithm problems. In order, we consider a shock tube, spherical adiabatic
collapse, cooling flow model, drag, a cosmological simulation, rotating
cloud-collapse and disc stability. In the implementations special attention is
given to the way in which force symmetry is enforced in the equations of
motion. We study in detail how the hydrodynamics are affected by different
implementations of the artificial viscosity including those with a
shear-correction modification. We present an improved first-order
smoothing-length update algorithm that is designed to remove instabilities that
are present in the Hernquist and Katz (1989) algorithm.
For all tests we find that the artificial viscosity is the most important
factor distinguishing the results from the various implementations. The second
most important factor is the way force symmetry is achieved in the equation of
motion. Most results favour a kernel symmetrization approach. The exact method
by which SPH pressure forces are included has comparatively little effect on
the results. Combining the equation of motion presented in Thomas and Couchman
(1992) with a modification of the Monaghan and Gingold (1983) artificial
viscosity leads to an SPH scheme that is both fast and reliable.Comment: 30 pages, 26 figures and 9 tables included. Submitted to MNRAS.
Postscript version available at
ftp://phobos.astro.uwo.ca/pub/etittley/papers/sphtest.ps.g
Abell 754: A Non-Head On Collision of Subclusters
We have analyzed spatially resolved spectra of A754 obtained with ASCA. The
combination of spectral and imaging capabilities of ASCA has set unprecedented
constraints on the hydrodynamical effects of a cluster merger. We find
significant gas temperature variations over the cluster face, indicating shock
heating of the atmosphere during the merger. The hottest region, >12 keV, is
located in the region of the Northwest galaxy clump though the entire region
along the cluster axis appears to be hotter than the mean cluster temperature
(~9 keV). The cool, ~5 keV gas originally found with the HEAO1-A2 experiment,
resides in the exterior of the cluster atmosphere and in plume of gas we
identify with a stripped cool atmosphere of the infalling subcluster. We have
also attempted to reconstruct an iron abundance map of this merging system.
Though poorly constrained, no significant deviations of abundance from the mean
value are apparent in the individual regions.
A754 is the only cluster so far which shows the significant temperature
pattern expected in a subcluster merger, in both the ROSAT (Henry & Briel 1995)
and ASCA data, providing the first possibility to compare it with theoretical
predictions. The observed temperature and surface brightness maps suggest that
the two colliding subunits have missed each other by about 1 Mpc, and are now
moving perpendicular to the cluster axis in the image plane (as, e.g., in the
simulations by Evrard etal 1996).Comment: Latex, 10 pages, 3 figures incl. color plate, uses aaspp4.sty,
flushrt.sty and pstricks.sty. Submitted to ApJ Letter
Chandra and XMM-Newton Observations of the Double Cluster Abell 1758
Abell 1758 was classified as a single rich cluster of galaxies by Abell, but
a ROSAT observation showed that this system consists of two distinct clusters
(A1758N and A1758S) separated by approximately 8\arcmin (a projected
separation of 2 Mpc in the rest frame of the clusters). Only a few galaxy
redshifts have been published for these two clusters, but the redshift of the
Fe lines in the Chandra and XMM-Newton spectra shows that the recessional
velocities of A1758N and A1758S are within 2,100 km s. Thus, these two
clusters most likely form a gravitationally bound system, but our imaging and
spectroscopic analyses of the X-ray data do not reveal any sign of interaction
between the two clusters. The Chandra and XMM-Newton observations show that
A1758N and A1758S are both undergoing major mergers.
A1758N is in the late stages of a large impact parameter merger between two 7
keV clusters. The two remnant cores have a projected separation of 800 kpc.
Based on the measured pressure jumps preceding the two cores, they are receding
from one another at less than 1,600 km s. The two cores are surrounded
by hotter gas (--12 keV) that was probably shock heated during
the early stages of the merger. The gas entropy in the two remnant cores is
comparable with the central entropy observed in dynamically relaxed clusters,
indicating that the merger-induced shocks stalled as they tried to penetrate
the high pressure cores of the two merging systems.Each core also has a wake of
low entropy gas indicating that this gas was ram pressure stripped without
being strongly shocked (abridged). (A copy of the paper with higher resolution
images is available at http://asc.harvard.edu/~lpd/a1758.ps).Comment: paper plus 13 figure
A Comparison of Semi-Analytic and Smoothed Particle Hydrodynamics Galaxy Formation
We compare the statistical properties of galaxies found in two different
models of hierarchical galaxy formation: the semi-analytic model of Cole et al.
and the smoothed particle hydrodynamics (SPH) simulations of Pearce et al.
Using a `stripped-down' version of the semi-analytic model which mimics the
resolution of the SPH simulations and excludes physical processes not included
in them, we find that the two models produce an ensemble of galaxies with
remarkably similar properties, although there are some differences in the gas
cooling rates and in the number of galaxies that populate halos of different
mass. The full semi-analytic model, which has effectively no resolution limit
and includes a treatment of star formation and supernovae feedback, produces
somewhat different (but readily understandable) results. Agreement is
particularly good for the present-day global fractions of hot gas, cold dense
(i.e. galactic) gas and uncollapsed gas, for which the SPH and stripped-down
semi-analytic calculations differ by at most 25%. In the most massive halos,
the stripped-down semi-analytic model predicts, on the whole, up to 50% less
gas in galaxies than is seen in the SPH simulations. The two techniques
apportion this cold gas somewhat differently amongst galaxies in a given halo.
This difference can be tracked down to the greater cooling rate in massive
halos in the SPH simulation compared to the semi-analytic model. (abridged)Comment: 19 pages, 13 figures, to appear in MNRAS. Significantly extended to
explore galaxy progenitor distributions and behaviour of models at high
redshift
- âŠ