5,016 research outputs found
A turbulence model for smoothed particle hydrodynamics
The aim of this paper is to devise a turbulence model for the particle method
Smoothed Particle Hydrodynamics (SPH) which makes few assumptions, conserves
linear and angular momentum, satisfies a discrete version of Kelvin's
circulation theorem, and is computationally efficient. These aims are achieved.
Furthermore, the results from the model are in good agreement with the
experimental and computational results of Clercx and Heijst for two dimensional
turbulence inside a box with no-slip walls. The model is based on a Lagrangian
similar to that used for the Lagrangian averaged Navier Stokes (LANS)
turbulence model, but with a different smoothed velocity. The smoothed velocity
preserves the shape of the spectrum of the unsmoothed velocity, but reduces the
magnitude for short length scales by an amount which depends on a parameter
. We call this the SPH- model. The effectiveness of the
model is indicated by the fact that the second order velocity correlation
function calculated using the smoothed velocity and a coarse resolution, is in
good agreement with a calculation using a resolution which is finer by a factor
2, and therefore requires 8 times as much work to integrate to the same time.Comment: 34 pages, 11 figure
Smoothed Particle Magnetohydrodynamics II. Variational principles and variable smoothing length terms
In this paper we show how a Lagrangian variational principle can be used to
derive the SPMHD (smoothed particle magnetohydrodynamics) equations for ideal
MHD. We also consider the effect of a variable smoothing length in the SPH
kernels after which we demonstrate by numerical tests that the consistent
treatment of terms relating to the gradient of the smoothing length in the
SPMHD equations significantly improves the accuracy of the algorithm. Our
results complement those obtained in a companion paper (Price and Monaghan
2003a, paper I) for non ideal MHD where artificial dissipative terms were
included to handle shocks.Comment: 14 pages, 4 figures, accepted to MNRA
Smoothed Particle Magnetohydrodynamics (some shocking results...)
There have been some issues in the past in attempts to simulate magnetic
fields using the Smoothed Particle Hydrodynamics (SPH) method. SPH is well
suited to star formation problems because of its Lagrangian nature. We present
new, stable and conservative methods for magnetohydrodynamics (MHD) in SPH and
present numerical tests on both waves and shocks in one dimension to show that
it gives robust and accurate results.Comment: Kluwer latex, 6 pages, 3 figures; Proceedings of the International
Workshop "Magnetic Fields and Star Formation: Theory vs Observations",
Madrid, 21-25 April 2003. Revised version accepted to proceedings (exact
solutions added, other minor changes
SPH simulations of turbulence in fixed and rotating boxes in two dimensions with no-slip boundaries
In this paper we study decaying turbulence in fixed and rotating boxes in two
dimen- sions using the particle method SPH. The boundaries are specified by
boundary force particles, and the turbulence is initiated by a set of gaussian
vortices. In the case of fixed boxes we recover the results of Clercx and his
colleagues obtained using both a high accuracy spectral method and experiments.
Our results for fixed boxes are also in close agreement with those of Monaghan1
and Robinson and Monaghan2 obtained using SPH. A feature of decaying turbulence
in no-slip, square, fixed boundaries is that the angular momentum of the fluid
varies with time because of the reaction on the fluid of the viscous stresses
on the boundary. We find that when the box is allowed to rotate freely, so that
the total angular momentum of box and fluid is constant, the change in the
angular momentum of the fluid is a factor ~ 500 smaller than is the case for
the fixed box, and the final vorticity distribution is different. We also
simulate the behaviour of the turbulence when the box is forced to rotate with
small and large Rossby number, and the turbulence is initiated by gaussian
vortices as before. If the rotation of the box is maintained after the
turbulence is initiated we find that in the rotating frame the decay of kinetic
energy, enstrophy and the vortex structure is insensitive to the angular
velocity of the box. On the other hand, If the box is allowed to rotate freely
after the turbulence is initiated, the evolved vortex structure is completely
different
Lateralised sleep spindles relate to false memory generation
The file attached to this record is the author's final peer reviewed version. The Publisher's final version can be found by following the DOI link.Sleep is known to enhance false memories: After presenting participants with lists of semantically related words, sleeping before recalling these words results in a greater acceptance of unseen “lure” words related in theme to previously seen words. Furthermore, the right hemisphere (RH) seems to be more prone to false memories than the left hemisphere (LH). In the current study, we investigated the sleep architecture associated with these false memory and lateralisation effects in a nap study. Participants viewed lists of related words, then stayed awake or slept for approximately 90 min, and were then tested for recognition of previously seen-old, unseen-new, or unseen-lure words presented either to the LH or RH. Sleep increased acceptance of unseen-lure words as previously seen compared to the wake group, particularly for RH presentations of word lists. RH lateralised stage 2 sleep spindle density relative to the LH correlated with this increase in false memories, suggesting that RH sleep spindles enhanced false memories in the RH
Direct Numerical Simulation of decaying two-dimensional turbulence in a no-slip square box using Smoothed Particle Hydrodynamics
This paper explores the application of SPH to a Direct Numerical Simulation
(DNS) of decaying turbulence in a two-dimensional no-slip wall-bounded domain.
In this bounded domain, the inverse energy cascade, and a net torque exerted by
the boundary, result in a spontaneous spin up of the fluid, leading to a
typical end state of a large monopole vortex that fills the domain. The SPH
simulations were compared against published results using a high accuracy
pseudo-spectral code. Ensemble averages of the kinetic energy, enstrophy and
average vortex wavenumber compared well against the pseudo-spectral results, as
did the evolution of the total angular momentum of the fluid. However, while
the pseudo-spectral results emphasised the importance of the no-slip boundaries
as generators of long lived coherent vortices in the flow, no such generation
was seen in the SPH results. Vorticity filaments produced at the boundary were
always dissipated by the flow shortly after separating from the boundary layer.
The kinetic energy spectrum of the SPH results was calculated using a SPH
Fourier transform that operates directly on the disordered particles. The
ensemble kinetic energy spectrum showed the expected k-3 scaling over most of
the inertial range. However, the spectrum flattened at smaller length scales
(initially less than 7.5 particle spacings and growing in size over time),
indicating an excess of small-scale kinetic energy
Evaluations of topological Tutte polynomials
We find new properties of the topological transition polynomial of embedded
graphs, . We use these properties to explain the striking similarities
between certain evaluations of Bollob\'as and Riordan's ribbon graph
polynomial, , and the topological Penrose polynomial, . The general
framework provided by also leads to several other combinatorial
interpretations these polynomials. In particular, we express , ,
and the Tutte polynomial, , as sums of chromatic polynomials of graphs
derived from ; show that these polynomials count -valuations of medial
graphs; show that counts edge 3-colourings; and reformulate the Four
Colour Theorem in terms of . We conclude with a reduction formula for the
transition polynomial of the tensor product of two embedded graphs, showing
that it leads to additional relations among these polynomials and to further
combinatorial interpretations of and .Comment: V2: major revision, several new results, and improved expositio
Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
In this paper we discuss the treatment of discontinuities in Smoothed
Particle Hydrodynamics (SPH) simulations. In particular we discuss the
difference between integral and differential representations of the fluid
equations in an SPH context and how this relates to the formulation of dissip
ative terms for the capture of shocks and other discontinuities.
This has important implications for many problems, in particular related to
recently highlighted problems in treating Kelvin-Helmholtz instabilities across
entropy gradients in SPH. The specific problems pointed out by Agertz et al.
(2007) are shown to be related in particular to the (lack of) treatment of
contact discontinuities in standard SPH formulations which can be cured by the
simple application of an artificial thermal conductivity term. We propose a new
formulation of artificial thermal conductivity in SPH which minimises
dissipation away from discontinuities and can therefore be applied quite
generally in SPH calculations.Comment: 31 pages, 10 figures, submitted to J. Comp. Phys. Movies + hires
version available at http://www.astro.ex.ac.uk/people/dprice/pubs/kh/ . v3:
modified as per referee's comments - comparison with Ritchie & Thomas
formulation added, quite a few typos fixed. No major change in metho
- …