Efficient Simulation of Fluids

Fluid simulation is based on Navier-Stokes equations. Efficient simulation codes may rely on the smooth particle hydrodynamic toolbox (SPH), a method that uses kernel density estimation. Many variants of SPH have been proposed to optimize the simulation, like implicit incompressible SPH (IISPH) or predictive-corrective incompressible SPH (PC-ISPH). This chapter recalls the formulation of SPH and focuses on its effective parallel implementation using the Nvidia common unified device architecture (CUDA), while message passing interface (MPI) is another option. The key to effective implementation is a dedicated accelerating structure, and therefore some well-chosen parallel design patterns are detailed. Using a rough model of the ocean, this type of simulation can be used directly to simulate a tsunami resulting from an underwater earthquake

Fluid Simulation by the Smoothed Particle Hydrodynamics Method: A Survey.

This paper presents a survey of Smoothed Particle Hydrodynamics (SPH) and its use in computational fluid dynamics. As a truly mesh-free particle method based upon the Lagrangian formulation, SPH has been applied to a variety of different areas in science, computer graphics and engineering. It has been established as a popular technique for fluid based simulations, and has been extended to successfully simulate various phenomena such as multi-phase flows, rigid and elastic solids, and fluid features such as air bubbles and foam. Various aspects of the method will be discussed: Similarities, advantages and disadvantages in comparison to Eulerian methods; Fundamentals of the SPH method; The use of SPH in fluid simulation; The current trends in SPH. The paper ends with some concluding remarks about the use of SPH in fluid simulations, including some of the more apparent problems, and a discussion on prospects for future work

Blue Noise Sampling using an SPH-based Method

We propose a novel algorithm for blue noise sampling inspired by the Smoothed Particle Hydrodynamics (SPH) method. SPH is a well-known method in fluid simulation -- it computes particle distributions to minimize the internal pressure variance. We found that this results in sample points (i.e., particles) with a high quality blue-noise spectrum. Inspired by this, we tailor the SPH method for blue noise sampling. Our method achieves fast sampling in general dimensions for both surfaces and volumes. By varying a single parameter our method can generate a variety of blue noise samples with different distribution properties, ranging from Lloyd's relaxation to Capacity Constrained Voronoi Tessellations ({CCVT}). Our method is fast and supports adaptive sampling and multi-class sampling. We have also performed experimental studies of the SPH kernel and its influence on the distribution properties of samples. We demonstrate with examples that our method can generate a variety of controllable blue noise sample patterns, suitable for applications such as image stippling and re-meshing

Experimental and numerical evaluation of the force due to the impact of a dam-break wave on a structure

Flood events caused by the collapse of dams or river levees can have damaging consequences on buildings and infrastructure located in prone areas. Accordingly, a careful prediction of the hydrodynamic load acting on structures is important for flood hazard assessment and potential damage evaluation. However, this represents a challenging task and requires the use of suitable mathematical models. This paper investigates the capability of three different models, i.e. a 2D depth-averaged model, a 3D Eulerian two-phase model, and a 3D Smoothed Particle Hydrodynamics (SPH) model, to estimate the impact load exerted by a dam-break wave on an obstacle. To this purpose, idealised dam-break experiments were carried out by generating a flip-through impact against a rigid squat structure, and measurements of the impact force were obtained directly by using a load cell. The dynamics of the impact event was analyzed and related to the measured load time history. A repeatability analysis was performed due to the great variability typically shown by impact phenomena, and a confidence range was estimated. The comparison between numerical results and experimental data shows the capability of 3D models to reproduce the key features of the flip-through impact. The 2D modelling based on the shallow water approach is not entirely suitable to accurately reproduce the load hydrograph and predict the load peak values; this difficulty increases with the strength of the wave impact. Nevertheless, the error in the peak load estimation is in the order of 10% only, thus the 2D approach may be considered appropriate for practical applications. Moreover, when the shallow water approximation is expected to work well, 2D results are comparable with the experimental data, as well as with the numerical predictions of far more sophisticated and computationally demanding 3D solvers. All the numerical models overestimate the falling limb of the load hydrograph after the impact. The SPH model ensures good evaluation of the long-time load impulse. The 2D shallow water solver and the 3D Eulerian model are less accurate in predicting the load impulse but provide similar results. A sensitivity analysis with respect to the model parameters allows to assess model uncertainty