Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

We present two new Lagrangian methods for hydrodynamics, in a systematic comparison with moving-mesh, smoothed particle hydrodynamics (SPH), and stationary (non-moving) grid methods. The new methods are designed to simultaneously capture advantages of both SPH and grid-based/adaptive mesh refinement (AMR) schemes. They are based on a kernel discretization of the volume coupled to a high-order matrix gradient estimator and a Riemann solver acting over the volume ‘overlap’. We implement and test a parallel, second-order version of the method with self-gravity and cosmological integration, in the code gizmo:1 this maintains exact mass, energy and momentum conservation; exhibits superior angular momentum conservation compared to all other methods we study; does not require ‘artificial diffusion’ terms; and allows the fluid elements to move with the flow, so resolution is automatically adaptive. We consider a large suite of test problems, and find that on all problems the new methods appear competitive with moving-mesh schemes, with some advantages (particularly in angular momentum conservation), at the cost of enhanced noise. The new methods have many advantages versus SPH: proper convergence, good capturing of fluid-mixing instabilities, dramatically reduced ‘particle noise’ and numerical viscosity, more accurate sub-sonic flow evolution, and sharp shock-capturing. Advantages versus non-moving meshes include: automatic adaptivity, dramatically reduced advection errors and numerical overmixing, velocity-independent errors, accurate coupling to gravity, good angular momentum conservation and elimination of ‘grid alignment’ effects. We can, for example, follow hundreds of orbits of gaseous discs, while AMR and SPH methods break down in a few orbits. However, fixed meshes minimize ‘grid noise’. These differences are important for a range of astrophysical problems.

A new class of accurate, mesh-free hydrodynamic simulation methods

In this survey we consider the development and mathematical analysis of numerical methods for kinetic partial differential equations. Kinetic equations represent a way of describing the time evolution of a system consisting of a large number of particles. Due to the high number of dimensions and their intrinsic physical properties, the construction of numerical methods represents a challenge and requires a careful balance between accuracy and computational complexity. Here we review the basic numerical techniques for dealing with such equations, including the case of semi-Lagrangian methods, discrete-velocity models and spectral methods. In addition we give an overview of the current state of the art of numerical methods for kinetic equations. This covers the derivation of fast algorithms, the notion of asymptotic-preserving methods and the construction of hybrid schemes.

https://hal.archives-ouvertes.fr/hal-00986714/document

Numerical methods for kinetic equations

Adaptive Mesh and Algorithm Refinement Using Direct Simulation Monte Carlo

Unified solver for rarefied and continuum flows with adaptive mesh and algorithm refinement

We consider a self-propelled interacting particle system for the collective behavior of swarms of animals, and extend it with an attraction term called roosting force, as it has been suggested in Ref. 30. This new force models the tendency of birds to overfly a fixed preferred location, e.g. a nest or a food source. We include roosting to the existing individual-based model and consider the associated mean-field and hydrodynamic equations. The resulting equations are investigated analytically looking at different asymptotic limits of the corresponding stochastic model. In addition to existing patterns like single mills, the inclusion of roosting yields new scenarios of collective behavior, which we study numerically on the microscopic as well as on the hydrodynamic level.

/pdf/self-propelled-interacting-particle-systems-with-roosting-5a268zbrr6.pdf

Self-propelled interacting particle systems with roosting force

A Lagrangian particle scheme is applied to the projection method for the incompressible Navier-Stokes equations. The approximation of spatial derivatives is obtained by the weighted least squares method. The pressure Poisson equation is solved by a local iterative procedure with the help of the least squares method. Numerical tests are performed for two dimensional cases. The Couette flow, Poisseulle flow and the driven cavity flow are presented. The numerical solutions are obtained for stationary as well as instationary cases and are compared with the analyatical solutions for channel flows. Finally, the driven cavity flow in a unit square is considered and the stationary solution obtained from this scheme is compared with that from the finite element method!

/pdf/finite-pointset-method-based-on-the-projection-method-for-4fo2a3ba9p.pdf

Finite Pointset Method Based on the Projection Method for Simulations of the Incompressible Navier-Stokes Equations

Modeling of two-phase flows with surface tension by finite pointset method (FPM)

https://kluedo.ub.uni-kl.de/files/625/gruen_189.pdf

Sudarshan Tiwari

Papers

Self-propelled interacting particle systems with roosting force

Finite Pointset Method Based on the Projection Method for Simulations of the Incompressible Navier-Stokes Equations

Modeling of two-phase flows with surface tension by finite pointset method (FPM)

An adaptive domain decomposition procedure for Boltzmann and Euler equations

Finite pointset method for simulation of the liquid–liquid flow field in an extractor