1,065 research outputs found
Enhanced SPH modeling of free-surface flows with large deformations
The subject of the present thesis is the development of a numerical solver to
study the violent interaction of marine flows with rigid structures.
Among the many numerical models available, the Smoothed Particle
Hydrodynamics (SPH) has been chosen as it proved
appropriate in dealing with violent free-surface flows. Due to its
Lagrangian and meshless character it can naturally handle breaking waves and
fragmentation that generally are not easily treated by standard methods. On
the other hand, some consolidated features of mesh-based methods, such as
the solid boundary treatment, still remain unsolved issues in the SPH
context.
In the present work a great part of the research activity has been devoted
to tackle some of the bottlenecks of the method. Firstly, an enhanced SPH
model, called delta-SPH, has been proposed. In this model, a proper numerical diffusive
term has been added in the continuity equation in order to remove the spurious
numerical noise in the pressure field which typically affects the weakly-compressible SPH
models. Then, particular attention has been paid to the development of suitable
techniques for the enforcement of the boundary conditions. As for the free-surface, a
specific algorithm has been designed to detect free-surface particles and
to define a related level-set function with two main targets: to allow the
imposition of peculiar conditions on the free-surface and to analyse and
visualize more easily the simulation outcome (especially in 3D cases).
Concerning the solid boundary treatment, much effort has been spent to
devise new techniques for handling generic body geometries with an adequate
accuracy in both 2D and 3D problems. Two different techniques have been
described: in the first one the standard ghost fluid method has been
extended in order to treat complex solid geometries. Both free-slip and
no-slip boundary conditions have been implemented, the latter being a quite
complex matter in the SPH context. The proposed boundary treatment proved
to be robust and accurate in evaluating local and global loads, though it
is not easy to extend to generic 3D surfaces.
The second technique has been adopted for these cases.
Such a technique has been developed in the context of Riemann-SPH methods
and in the present work is reformulated in the context of the standard SPH scheme.
The method proved to be robust in treating complex 3D
solid surfaces though less accurate than the former.
Finally, an algorithm to correctly initialize the SPH simulation in the case of generic
geometries has been described. It forces a resettlement of the fluid particles
to achieve a regular and uniform spacing even in complex configurations. This
pre-processing procedure avoids the generation of spurious currents due to
local defects in the particle distribution at the beginning of the simulation.
The delta-SPH model has been validated against several problems
concerning fluid-structure interactions. Firstly, the capability of the
solver in dealing with water impacts has been tested by simulating a
jet impinging on a flat plate and a dam-break flow against a vertical
wall. In this cases, the accuracy in the prediction of local loads and of
the pressure field have been the main focus. Then, the viscous flow around
a cylinder, in both steady and unsteady conditions, has been simulated
comparing the results with reference solutions. Finally, the generation
and propagation of 2D gravity waves has been simulated. Several
regimes of propagation have been tested and the results
compared against a potential flow solver.
The developed numerical solver has been applied to several cases of
free-surface flows striking rigid structures and to the problem of the
generation and evolution of ship generated waves. In the former case, the
robustness of the solver has been challenged by simulating 2D and 3D water impacts
against complex solid surfaces. The numerical outcome have been compared
with analytical solutions, experimental data and other numerical results
and the limits of the model have been discussed.
As for the ship generated waves, the problem has been firstly studied
within the 2D+t approximation, focusing
on the occurrence and features of the breaking bow waves. Then, a
dedicated 3D SPH parallel solver has been developed to tackle the simulation
of the entire ship in constant forward motion. This simulation is quite demanding in
terms of complexities of the boundary geometry and computational resources
required. The wave pattern obtained has been compared against experimental
data and results from other numerical methods, showing in both the cases a fair
and promising agreement
A Smoothed Particle Hydrodynamics Method for the Simulation of Centralized Sloshing Experiments
The Smoothed Particle Hydrodynamics (SPH) method is proposed for studying hydrodynamic processes related to nuclear engineering problems. A problem of possible recriticality due to the sloshing motions of the molten reactor core is studied with SPH method. The accuracy of the numerical solution obtained in this study with the SPH method is significantly higher than that obtained with the SIMMER-III/IV reactor safety analysis code
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
Verification and validation in highly viscous fluid simulation using a fully implicit sph method
Catastrophes involving mass movements has always been a great threat to civilizations. We propse to simplify the behavior of the mass movement material as a highly viscous fluid, possibly non-Newtonian. In this context, this study describes the application of two improvements in highly viscous fluid simulations using the smoothed particle hydrodynamics (SPH) method: an implicit time integration scheme to overcome the problem of impractically small
time-step restriction, and the introduction of air ghost particles to fix problems regarding the free-surface treatment. The application of a fully implicit time integration method implies an adaptation of the wall boundary condition, which is also covered in this study. Furthermore, the proposed wall boundary condition allows for different slip conditions, which is usually difficult to adopt in SPH. To solve a persistent problem on the SPH method of unstable pressure
distributions, we adopted the incompressible SPH [1] as a basis for the implementation of these improvements, which guarantees stable and accurate pressure distribution. We conducted non-Newtonian pipe flow simulations to verify the method and a variety of dam break and wave generated by underwater landslide simulations for validation. Finally, we demonstrate the potential of this method with the highly viscous vertical jet flow over a horizontal plate test, which features a complex viscous coiling behavior
A moving least square immersed boundary method for SPH with thin-walled structures
This paper presents a novel method for smoothed particle hydrodynamics (SPH)
with thin-walled structures. Inspired by the direct forcing immersed boundary
method, this method employs a moving least square method to guarantee the
smoothness of velocity near the structure surface. It simplifies thin-walled
structure simulations by eliminating the need for multiple layers of boundary
particles, and improves computational accuracy and stability in
three-dimensional scenarios. Supportive three-dimensional numerical results are
provided, including the impulsively started plate and the flow past a cylinder.
Results of the impulsively started test demonstrate that the proposed method
obtains smooth velocity and pressure in the, as well as a good match to the
references results of the vortex wake development. In addition, results of the
flow past cylinder test show that the proposed method avoids mutual
interference on both side of the boundary, remains stable for three-dimensional
simulations while accurately calculating the forces acting on structure.Comment: 15 pages,11 figure
Smoothed Particle Hydrodynamics for Navier-Stokes Fluid Flow Application
The aim of this publication is to introduce the particle based computational fluid dynamics (CFD) method smoothed particle hydrodynamics (SPH) and introduce an applicable and valid SPH implementation for practical cases. For this purpose, current research approaches are combined regarding performance and numerical stability. The principles of the method, the mathematical basics and the discretization of the Navier-Stokes equations are clarified. Furthermore, the implementation of method-specific boundary conditions, wall, inlet and outlet, as well as several correction procedures and a surface tension setup into the present code framework are described. The advantages and validity of the method are shown based on different cases. The free surface fluid behavior of a dam break is compared to experimental data of the time dependent water level of selected positions. A Karman vortex street is validated by its Strouhal number for different Reynolds numbers. The frequency of an oscillating drop is analysed and compared to the analytical solution. The SPH is utilized for pipe flows influenced by a backward facing step and shows an expected qualitative flow field
- …