scholarly journals A MULTI-NODE APPROACH TO SIMULATE THIN COASTAL STRUCTURES IN THE SPH CONTEXT

Author(s):  
Francesco Aristodemo ◽  
Domenico Davide Meringolo ◽  
Paolo Veltri

We propose an improvement in modeling solid boundary conditions for 2D weakly-compressible Smoothed Particle Hydrodynamics (SPH) simulations for cases in which the thickness of the body is small compared to the desired particle size and the fluid surrounds the body from more than one side. Specifically, the fixed ghost particles technique developed by Marrone et al. (2011), based on interpolation nodes located within the fluid domain, is here extended to a multi-node approach. The fluid domain is thus divided into various sub-areas and an interpolation node for the considered solid particle is associated to every sub-area. Consequently, the solid particles present an array of values interpolated at different sub-areas for the same physical quantity. When a fluid particle located in a specific region interacts with a multi-node fixed ghost particle, the last assumes the field values interpolated in the reference area through the associated node. The present modeling allows to adopt a coarser spatial resolution to model the same physical problem, resulting in a reduction of the computational cost. The proposed solid boundary treatment is applied to horizontal decks and perforated wall-caisson breakwaters subjected to regular waves. In this context, an automatic hybrid diffusive formulation is introduced in order to prevent shock waves during water impacts and preserve the hydrostatic pressure. The formulation is obtained by defining a variable parameter detecting the occurrence of relevant density gradients induced by fluid impacts, resulting in an automatic switch between the two formulations.

Author(s):  
Rubens A. Amaro Junior ◽  
Liang-Yee Cheng ◽  
Sergei K. Buruchenko

Lagrangian particle-based methods have opened new perspectives for the investigation of complex problems with large free-surface deformation. Some well-known particle-based methods adopted to solve non-linear hydrodynamics problems are the smoothed parti- cle hydrodynamics (SPH) and the moving particle semi-implicit (MPS). Both methods model the continuum by a system of Lagrangian particles (points), but adopting distinct approaches for the numerical operators, pressure calculation, and boundary conditions. Despite the ability of the particle-based methods in modeling highly nonlinear hydrodynamics, some shortcomings, such as unstable pressure computation and high computational cost remain. In order to assess the performance of these two methods, the weakly-compressible SPH (WCSPH) parallel solver, DualSPHysics, and an in-house incompressible MPS solver are adopted in this work. Two test cases consisting of three-dimensional (3D) dam-break problems are simulated, and wave heights, pressures and forces are compared with the available experimental data. The influence of the artificial viscosity on the accuracy of WCSPH is investigated. Computational times of both solvers are also compared. Finally, the relative benefits of the methods for solving free-surface problems are discussed, therefore providing directions of their applicability.


2018 ◽  
Vol 19 (2) ◽  
pp. 172-181
Author(s):  
Muhammad Zahir Ramli

ABSTRACT: An implementation of the weakly compressible smoothed particle hydrodynamics (WCSPH) method is demonstrated to determine the hydrodynamics coefficients through radiation problem of an oscillating 2D rectangular box. Three possible modes of motion namely swaying, heaving, and rolling are carried out to establish the influence of oscillating motions in predicting the added mass and damping. Both solid boundary and fluid flow are modelled by WCSPH and validated by the potential flow and experimental results. Discrepancies observed at lower frequencies are further investigated using different particle resolutions, different time steps, and extending the domain with longer runtime to demonstrate the performance of WCSPH. Finally, flow separation and vortices are discussed and compared with experimental results. ABSTRAK: Bagi fenomena yang melibatkan radiasi dalam air, segiempat kotak 2D diosilasikan dengan menggunakan simulasi WCSPH untuk memperoleh pekali hidrodinamik. Mod osilasi terbahagi kepada 3 iaitu sway, heave dan roll. Osilasi dengan mengguna pakai kotak akan mempengaruhi pergerakan air dalam menentukan nilai penambahan jisim dan rendaman. Keseluruhan domain air dan sempadan telah dimodelkan dengan menggunakan WCSPH. Semua model tersebut kemudiannya akan dibandingkan melalui keputusan eksperimen dan teori. Jika keputusan melalui kaedah WCSPH ini berbeza, terutama pada frekuensi rendah, penyelidikan lanjut akan dilakukan dengan menggunakan zarah resolusi yang berbeza, langkah masa yang berbeza dan menambah masa domain ujikaji bagi menilai keputusan WCSPH. Akhirnya, kriteria aliran dan kadar pusaran yang terhasil di sekeliling kotak akan dibincang dan dibandingkan bersama keputusan eksperimen.


2019 ◽  
Vol 16 (07) ◽  
pp. 1850114 ◽  
Author(s):  
Jiandong He ◽  
Juanmian Lei

In this paper, we developed a GPU parallelized Total Lagrangian Formation of Smoothed Particle Hydrodynamics (TLSPH) algorithm for 3D geometrical nonlinear structure analysis. The code was developed using NVDIA CUDA C++. Both the TLSPH and GPU parallelization algorithms are described in detail. Compared to the traditional FEM method for structure analysis, TLSPH method is much easier to be implemented and parallelized. In addition, as a meshless based method, there is no need to mesh the domain for TLSPH method. Also, the computational cost of TLSPH is much lower than the Weakly Compressible Smoothed Particle (WCSPH) method. By introducing GPU acceleration, we have significantly improved the code performance. Two benchmark test cases for 3D geometrical nonlinear structure analysis are carried out. The simulation results are compared with analysis results and the data obtained by Abaqus, which is a popularly-used software for structure analysis based on FEM method. In order to show the efficiency of GPU parallelization, a serial code based on the same TLSPH method is also developed as a reference. Results show GPU parallelization accelerates the code obviously. In summary, the GPU parallelized TLSPH method shows the potential to become an alternative way to deal with 3D geometrical nonlinear structure analysis.


Author(s):  
Steven J. Lind ◽  
Benedict D. Rogers ◽  
Peter K. Stansby

This paper presents a review of the progress of smoothed particle hydrodynamics (SPH) towards high-order converged simulations. As a mesh-free Lagrangian method suitable for complex flows with interfaces and multiple phases, SPH has developed considerably in the past decade. While original applications were in astrophysics, early engineering applications showed the versatility and robustness of the method without emphasis on accuracy and convergence. The early method was of weakly compressible form resulting in noisy pressures due to spurious pressure waves. This was effectively removed in the incompressible (divergence-free) form which followed; since then the weakly compressible form has been advanced, reducing pressure noise. Now numerical convergence studies are standard. While the method is computationally demanding on conventional processors, it is well suited to parallel processing on massively parallel computing and graphics processing units. Applications are diverse and encompass wave–structure interaction, geophysical flows due to landslides, nuclear sludge flows, welding, gearbox flows and many others. In the state of the art, convergence is typically between the first- and second-order theoretical limits. Recent advances are improving convergence to fourth order (and higher) and these will also be outlined. This can be necessary to resolve multi-scale aspects of turbulent flow.


Author(s):  
Shuai Meng ◽  
Qian Wang ◽  
Rui Yang

The phenomenon of impaction between liquid droplets and solid particles is involved in many scientific problems and engineering applications, such as impaction between sprayed droplet and solid particles in limestone injection desulfurization system and the collision between a droplet of the liquid to be granulated and a seed particle in fluidized bed spray granulation process. There are a lot of factors affected this phenomenon: droplet and particle size, momentum of both liquid droplet and solid particles, materials, surface conditions of the solid particles and so on. However the experimental or numerical researches have been done mostly pay attention to Specific application or process, so the impaction phenomenon has not been through studied, for example how different factors affected the impaction process with its effect on different applications. This paper focuses on the basic issue of interaction between droplet and solid particles. Three main factors were considered: ratio of diameter between the droplet and solid particle, relative velocity and the surface tension (including the contact angle between droplet and solid particle). All the study is based on simulation using SPH (smoothed particle hydrodynamics) method, and the surface tension is simulated by particle-particle interaction.


2014 ◽  
Author(s):  
Zhen Chen ◽  
Li Zou ◽  
Zhi Zong

In this paper, the impact pressures of two different base forms are comparatively studied using Smoothed Particle Hydrodynamics (SPH) method. It is summarized from previous works that the improved weakly compressible SPH model shows better performances than incompressible SPH model in numerical simulations of free surface flows accompany with large deformations and strong discontinuities. Such advantages are observed in numerical accuracy, stability and efficiency. The weakly compressible SPH model used in this paper is equipped with some new correction algorithms, among which include the density reinitialization algorithm and a new coupled dynamic Solid Boundary Treatment (SBT) on solid boundaries. The new boundary treatment combines the advantages of both the repulsive boundary treatment and the dynamic boundary treatment, intending to obtain more stable and accurate numerical results. A benchmark test of dam breaking is conducted to prove the reliability of the numerical model used in this paper. Two representative cases, among which one has one cavity and the other one has three cavities, are numerically investigated and compared to support the conclusion that the base form with cavities generally experience lower local and overall impact pressures than the base form of flat plate. It is found that with the application of cavities on the bottom, the peak values of the boundary pressure near central bottom significantly decrease, leading to smaller force load and better structural stability. The mechanisms of such phenomenon might be the pressure absorption effect conducted by the cavities.


Sign in / Sign up

Export Citation Format

Share Document