1,065 research outputs found

    Enhanced SPH modeling of free-surface flows with large deformations

    Get PDF
    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

    Get PDF
    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

    Full text link
    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

    Get PDF
    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

    Full text link
    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

    Get PDF
    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
    • …
    corecore