Flood modelling with hydraTE: 2+1-dimensional smoothed-particle hydrodynamics

We present HydraTE, our own implementation of the smoothed-particle hydrodynamics technique for shallow water that uses the adaptive size of the smoothing kernel as a proxy for the local water depth. We derive the equa- tions of motion for this approach from the Lagrangian before demonstrating that we can model the depth of water in a trough, implement vertical walls, recover the correct acceleration and terminal velocity for water flowing down a slope and obtain a stable hydraulic jump with the correct jump condition. We demonstrate that HydraTE performs well on two of the UK Environ- ment Agency flood modelling benchmark tests. Benchmark EA3 involves flow down an incline into a double dip depression and studies the amount of water that reaches the second dip. Our results are in agreement with those of the other codes that have attempted this test. Benchmark EA6 is a dam break into a horizontal channel containing a building. HydraTE again pro- duces results that are in good agreement with the other methods and the experimetal validation data except where the vertical velocity structure of the flow is expected to be multi-valued, such as the hydralic jump where the precise location is not recovered even though the pre- and post- jump water heights are. We conclude that HydraTE is suitable for a wide range of flood modelling problems as it preforms at least as well as the best available commercial alternatives for the problems we have tested

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

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

Extended plasticity in commercial-purity zinc

90% rolling-reduction of annealed commercial-purity zinc sheet (grain size 100 - 150 μ) results in the fragmentation of the large grains into, finally, stable micro-grains, 1 - 211 in diameter. The stability of the micro-grains is due to the presence of soluble and insoluble impurities which prevent recrystallization. This micro-grain material is strain-rate sensitive, and elongations of 200% have been obtained at room temperature. Although this as-rolled, 90% reduction zinc sheet is not super-plastic according to the current definition, its behaviour has led to the coining of the phase 'extended plasticity'. Evidence of grain-boundary sliding is found on examination of the surface by scanning electron microscopy, while the examination of thin foils and activation energy measurements support the dynamic softening (recovery) theory; thus, both these mechanisms must be operating, to a greater or less extent, to confer on this material the observe mechanical behaviour. It is finally concluded that it is dangerous to draw conclusions regarding the mechanism of plastic deformation from surface observations alone

Optimal traps in graphene

We transform the two-dimensional Dirac-Weyl equation, which governs the charge carriers in graphene, into a non-linear first-order differential equation for scattering phase shift, using the so-called variable phase method. This allows us to utilize the Levinson Theorem to find zero-energy bound states created electrostatically in realistic structures. These confined states are formed at critical potential strengths, which leads to us posit the use of `optimal traps' to combat the chiral tunneling found in graphene, which could be explored experimentally with an artificial network of point charges held above the graphene layer. We also discuss scattering on these states and find the zero angular momentum states create a dominant peak in scattering cross-section as energy tends towards the Dirac point energy, suggesting a dominant contribution to resistivity.Comment: 11 pages, 5 figure

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

The orientation of galaxy dark matter haloes around cosmic voids

Using the Millennium N-body Simulation we explore how the shape and angular momentum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We find that the major and intermediate axes of the haloes tend to lie parallel to the surface of the voids, whereas the minor axis points preferentially in the radial direction. We have quantified the strength of these alignments at different radial distances from the void centres. The effect of these orientations is still detected at distances as large as 2.2 Rvoid from the void centre. Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres we detect, at the 99.9 per cent confidence level, a signal that the angular momentum of those haloes tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, this signal is only detected in shells at the void surface (1 < R < 1.07 Rvoid) and disappears at larger distances. This signal, together with the similar alignment observed using real spiral galaxies, strongly supports the prediction of the Tidal Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, their angular momentum before the moment of turnaround

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$^{-1}$. 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$^{-1}$. The two cores are surrounded by hotter gas ($\mathrm{kT}=9$--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