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 $\sim\pm1.3$\%. 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

Integrals of Motion for Critical Dense Polymers and Symplectic Fermions

We consider critical dense polymers ${\cal L}(1,2)$. We obtain for this model the eigenvalues of the local integrals of motion of the underlying Conformal Field Theory by means of Thermodynamic Bethe Ansatz. We give a detailed description of the relation between this model and Symplectic Fermions including the indecomposable structure of the transfer matrix. Integrals of motion are defined directly on the lattice in terms of the Temperley Lieb Algebra and their eigenvalues are obtained and expressed as an infinite sum of the eigenvalues of the continuum integrals of motion. An elegant decomposition of the transfer matrix in terms of a finite number of lattice integrals of motion is obtained thus providing a reason for their introduction.Comment: 53 pages, version accepted for publishing on JSTA

The effect of radiative cooling on scaling laws of X-ray groups and clusters

We have performed cosmological simulations in a ΛCDM cosmology with and without radiative cooling in order to study the effect of cooling on the cluster scaling laws. Our simulations consist of 4.1 million particles each of gas and dark matter within a box size of 100 h-1 Mpc, and the run with cooling is the largest of its kind to have been evolved to z = 0. Our cluster catalogs both consist of over 400 objects and are complete in mass down to ~1013 h-1 M☉. We contrast the emission-weighted temperature-mass (Tew-M) and bolometric luminosity-temperature (Lbol-Tew) relations for the simulations at z = 0. We find that radiative cooling increases the temperature of intracluster gas and decreases its total luminosity, in agreement with the results of Pearce et al. Furthermore, the temperature dependence of these effects flattens the slope of the Tew-M relation and steepens the slope of the Lbol-Tew relation. Inclusion of radiative cooling in the simulations is sufficient to reproduce the observed X-ray scaling relations without requiring excessive nongravitational energy injection

Simulation of associative learning with the replaced elements model

Associative learning theories can be categorised according to whether they treat the representation of stimulus compounds in an elemental or configural manner. Since it is clear that a simple elemental approach to stimulus representation is inadequate there have been several attempts to produce more elaborate elemental models. One recent approach, the Replaced Elements Model (Wagner, 2003), reproduces many results that have until recently been uniquely predicted by Pearce’s Configural Theory (Pearce, 1994). Although it is possible to simulate the Replaced Elements Model using “standard” simulation programs the generation of the correct stimulus representation is complex. The current paper describes a method for simulation of the Replaced Elements Model and presents the results of two example simulations that show differential predictions of Replaced Elements and Pearce’s Configural Theor

DISSOLVED OXYGEN CHARACTERISTICS OF THE GAMTOOS ESTUARY, SOUTH AFRICA

The fall and recovery of dissolved oxygen (DO) is documented throughout the Gamtoos Estuary, South Africa during dry conditions and following light and heavy rainfall over a 13-month period from November 1992 to November 1993. Hypoxic conditions generally occurred in the near-bottom waters of the upper estuary. Localized fluctuations in DO levels are related to the natural diurnal fluctuation associated with photosynthesis of aquatic flora. The drop in DO levels following light rainfall is associated with the volume of oxygen-consuming compounds entering the estuary via runoff from adjacent agricultural fields. This hypoxia was short-lived. Following an extreme rainfall event, however, almost immediate hypoxia was recorded throughout the estuary, and DO levels deteriorated for some time thereafter as a result of the substantial input of organic matter into the estuary. The area of hypoxia and recovery was governed by the freshwater input at the tidal head, estuarine hydrometry and hydrodynamics. Tidal processes were identified as a source of replenishment of oxygen when, during high tide, seawater with a higher DO content penetrated the estuary.Afr. J. mar. Sci. 25: 99–10

Refined conformal spectra in the dimer model

Working with Lieb's transfer matrix for the dimer model, we point out that the full set of dimer configurations may be partitioned into disjoint subsets (sectors) closed under the action of the transfer matrix. These sectors are labelled by an integer or half-integer quantum number we call the variation index. In the continuum scaling limit, each sector gives rise to a representation of the Virasoro algebra. We determine the corresponding conformal partition functions and their finitizations, and observe an intriguing link to the Ramond and Neveu-Schwarz sectors of the critical dense polymer model as described by a conformal field theory with central charge c=-2.Comment: 44 page

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

Fusion algebra of critical percolation

We present an explicit conjecture for the chiral fusion algebra of critical percolation considering Virasoro representations with no enlarged or extended symmetry algebra. The representations we take to generate fusion are countably infinite in number. The ensuing fusion rules are quasi-rational in the sense that the fusion of a finite number of these representations decomposes into a finite direct sum of these representations. The fusion rules are commutative, associative and exhibit an sl(2) structure. They involve representations which we call Kac representations of which some are reducible yet indecomposable representations of rank 1. In particular, the identity of the fusion algebra is a reducible yet indecomposable Kac representation of rank 1. We make detailed comparisons of our fusion rules with the recent results of Eberle-Flohr and Read-Saleur. Notably, in agreement with Eberle-Flohr, we find the appearance of indecomposable representations of rank 3. Our fusion rules are supported by extensive numerical studies of an integrable lattice model of critical percolation. Details of our lattice findings and numerical results will be presented elsewhere.Comment: 12 pages, v2: comments and references adde