17,445 research outputs found
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
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
Integrals of Motion for Critical Dense Polymers and Symplectic Fermions
We consider critical dense polymers . 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
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 \%. 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. 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
- …