1,996 research outputs found
An improved SPH scheme for cosmological simulations
We present an implementation of smoothed particle hydrodynamics (SPH) with
improved accuracy for simulations of galaxies and the large-scale structure. In
particular, we combine, implement, modify and test a vast majority of SPH
improvement techniques in the latest instalment of the GADGET code. We use the
Wendland kernel functions, a particle wake-up time-step limiting mechanism and
a time-dependent scheme for artificial viscosity, which includes a high-order
gradient computation and shear flow limiter. Additionally, we include a novel
prescription for time-dependent artificial conduction, which corrects for
gravitationally induced pressure gradients and largely improves the SPH
performance in capturing the development of gas-dynamical instabilities. We
extensively test our new implementation in a wide range of hydrodynamical
standard tests including weak and strong shocks as well as shear flows,
turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas
clouds. We jointly employ all modifications; however, when necessary we study
the performance of individual code modules. We approximate hydrodynamical
states more accurately and with significantly less noise than standard SPH.
Furthermore, the new implementation promotes the mixing of entropy between
different fluid phases, also within cosmological simulations. Finally, we study
the performance of the hydrodynamical solver in the context of radiative galaxy
formation and non-radiative galaxy cluster formation. We find galactic disks to
be colder, thinner and more extended and our results on galaxy clusters show
entropy cores instead of steadily declining entropy profiles. In summary, we
demonstrate that our improved SPH implementation overcomes most of the
undesirable limitations of standard SPH, thus becoming the core of an efficient
code for large cosmological simulations.Comment: 21 figures, 2 tables, accepted to MNRA
Using the generalized interpolation material point method for fluid-solid interactions induced by surface tension
This thesis is devoted to the development of new, Generalized Interpolation Material Point Method (GIMP)-based algorithms for handling surface tension and contact (wetting) in fluid-solid interaction (FSI) problems at small scales. In these problems, surface tension becomes so dominant that its influence on both fluids and solids must be considered. Since analytical solutions for most engineering problems are usually unavailable, numerical methods are needed to describe and predict complicated time-dependent states in the solid and fluid involved due to surface tension effects. Traditional computational methods for handling fluid-solid interactions may not be effective due to their weakness in solving large-deformation problems and the complicated coupling of two different types of computational frameworks: one for solid, and the other for fluid. On the contrary, GIMP, a mesh-free algorithm for solid mechanics problems, is numerically effective in handling problems involving large deformations and fracture. Here we extend the capability of GIMP to handle fluid dynamics problems with surface tension, and to develop a new contact algorithm to deal with the wetting boundary conditions that include the modeling of contact angle and slip near the triple points where the three phases -- fluid, solid, and vapor -- meet. The error of the new GIMP algorithm for FSI problems at small scales, as verified by various benchmark problems, generally falls within the 5% range. In this thesis, we have successfully extended the capability of GIMP for handling FSI problems under surface tension in a one-solver numerical framework, a unique and innovative approach.Chapter 1. Introduction -- Chapter 2. Using the generalized interpolation material point method for fluid dynamics at low reynolds numbers -- Chapter 3. On the modeling of surface tension and its applications by the generalized interpolation material point method -- Chapter 4. Using the generalized interpolation material point method for fluid-solid interactions induced by surface tension -- Chapter 5. Conclusions
Hydrodynamic simulations with the Godunov SPH
We present results based on an implementation of the Godunov Smoothed
Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the
GADGET-3 hydrodynamic code. We first review the derivation of the GSPH
discretization of the equations of moment and energy conservation, starting
from the convolution of these equations with the interpolating kernel. The two
most important aspects of the numerical implementation of these equations are
(a) the appearance of fluid velocity and pressure obtained from the solution of
the Riemann problem between each pair of particles, and (b the absence of an
artificial viscosity term. We carry out three different controlled
hydrodynamical three-dimensional tests, namely the Sod shock tube, the
development of Kelvin-Helmholtz instabilities in a shear flow test, and the
"blob" test describing the evolution of a cold cloud moving against a hot wind.
The results of our tests confirm and extend in a number of aspects those
recently obtained by Cha (2010): (i) GSPH provides a much improved description
of contact discontinuities, with respect to SPH, thus avoiding the appearance
of spurious pressure forces; (ii) GSPH is able to follow the development of
gas-dynamical instabilities, such as the Kevin--Helmholtz and the
Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl
structures in the shear-flow test and the dissolution of the cold cloud in the
"blob" test.
We also discuss in detail the effect on the performances of GSPH of changing
different aspects of its implementation. The results of our tests demonstrate
that GSPH is in fact a highly promising hydrodynamic scheme, also to be coupled
to an N-body solver, for astrophysical and cosmological applications.
[abridged]Comment: 19 pages, 13 figures, MNRAS accepted, high resolution version can be
obtained at
http://adlibitum.oats.inaf.it/borgani/html/papers/gsph_hydrosim.pd
Recent advances in the simulation of particle-laden flows
A substantial number of algorithms exists for the simulation of moving
particles suspended in fluids. However, finding the best method to address a
particular physical problem is often highly non-trivial and depends on the
properties of the particles and the involved fluid(s) together. In this report
we provide a short overview on a number of existing simulation methods and
provide two state of the art examples in more detail. In both cases, the
particles are described using a Discrete Element Method (DEM). The DEM solver
is usually coupled to a fluid-solver, which can be classified as grid-based or
mesh-free (one example for each is given). Fluid solvers feature different
resolutions relative to the particle size and separation. First, a
multicomponent lattice Boltzmann algorithm (mesh-based and with rather fine
resolution) is presented to study the behavior of particle stabilized fluid
interfaces and second, a Smoothed Particle Hydrodynamics implementation
(mesh-free, meso-scale resolution, similar to the particle size) is introduced
to highlight a new player in the field, which is expected to be particularly
suited for flows including free surfaces.Comment: 16 pages, 4 figure
- …