• Journal of Ship and Ocean Technology
• /
• v.11 no.4
• /
• pp.29-54
• /
• 2007

Systematic sensitivity analysis of smoothed particle hydrodynamics method (SPH), a gridless Lagrangian particle method, was carried out in this study. Unlike traditional grid-based numerical schemes, systematic sensitivity study for computational parameters is very limited for SPH. In this study, the effect of computational parameters in SPH simulation is explored through two-dimensional dam-breaking and sloshing problem. The parameters to be considered are the speed of sound, the type of kernel function, the frequency of density re-initialization, particle number, smoothing length and pressure extraction method. Through a series of numerical test, detailed information was obtained about how SPH solution can be more stabilized and improved by adjusting computational parameters.

• Coupled systems mechanics
• /
• v.10 no.6
• /
• pp.491-507
• /
• 2021

The aim of this paper is to conduct a hydrodynamic analysis of fluid flow over different weir types using the analytical solution, the physical model taken from another article, and numericalsimulations through the Smoothed particle hydrodynamic method (SPH) using the compiled DualSPHysics source code. The paper covers the field of real fluid dynamics that includes a description of different proposed types of weirs in various flow regimes and the optimal solution for the most efficiency structure shape. A detailed presentation of the method, the structure and it's characteristics are included. Apart from the single stepped weir, two other weir types are proposed: a Divided type and a Downstream slopped type. All of them are modeled using the SPH method.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.6 no.4
• /
• pp.89-99
• /
• 2003

Smoothed particle hydrodynamics, SPH, is a gridless Lagrangian technique which is a useful alternative numerical analysis method to simulate high velocity deformation problems as well as astrophysical and cosmological problems. The SPH method brings about some difficulties such as tensile Instability and stress oscillation. A new SPH method, so called normalized algorithm, was introduced to overcome these difficulties. In this paper we aimed to estimate this method and have developed an one-dimensional normalized SPH program. The high velocity impact model of an aluminum bar has been analysed by using the developed program and a commercial hydrocode, LS-DYNA. The obtained numerical results showed good agreement with the results of the same model in reference. The program also showed more stable results than those of LS-DYNA in stress oscillation. We hopefully expect that the developed one-dimensional normalized SPH program can be used to solve hydrodynamic problems especially for explosive detonation analysis.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.85-101
• /
• 2020

The optimal parameters for the fluid-structure interaction analysis using the Smoothed Particle Hydrodynamics (SPH) for fluids and finite elements for structures, respectively, are explored, and the effectiveness of the simulations with those parameters is validated by solving several open surface fluid problems. For the optimization of the Equation of State (EOS) and the simulation parameters such as the time step, initial particle spacing, and smoothing length factor, a dam-break problem and deflection of an elastic plate is selected, and the least squares analysis is performed on the simulation results. With the optimal values of the pivotal parameters, the accuracy of the simulation is validated by calculating the exerted force on a moving solid column in the open surface fluid. Overall, the SPH-FEM coupled simulation is very effective to calculate the fluid-structure interaction. However, the relevant parameters should be carefully selected to obtain accurate results.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• Journal of the Korea Computer Graphics Society
• /
• v.14 no.4
• /
• pp.27-33
• /
• 2008

Physics-based graphic techniques are used when simulating and rendering natural phenomena such as smoke, water and flame with computational physics. We propose novel methods which render simulated particle data fast onto 3D using tetrahedron splat. We calculate the position and the normal vector of splat by SPH(smoothed particle hydrodynamics) method then we reconstruct splat into quadrangular pyramid to reduce seam. We implement this technique for SPH fluid simulation, and animate natural flow of water successfully.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.2
• /
• pp.122-125
• /
• 2011

In this work, the movement pattern of a floating breakwater is numerically analyzed using Smoothed Particle Hydrodynamic (SPH) method as a Lagrangian scheme. At the seaside, the regular incident waves with varying height and period were considered as the dynamic free surface boundary conditions. The smooth and impermeable beach slope was defined as the bottom boundary condition. The effects of various boundary conditions such as incident wave characteristics, beach slope, and water depth on the movement of the floating body were studied. The numerical results are in good agreement with the available experimental data in the literature The results of the movement of the floating body were used to determine the transmitted wave height at the corresponding boundary conditions.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.922-925
• /
• 2011

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method widely used for the modeling fluid flows. Simulations of explosions require, besides the hydrodynamic equations, a realistic equation of state, an energy source term, and a set of chemical kinetic equations to follow the composition changes of the gas during the explosion. The performance of the hydrodynamic equations is investigated in the framework of the Sedov-Taylor blast-wave. The implementation of chemical kinetic equations and equation of state is studied with 1D detonation of TNT slab. Our results are compared to those from analytical and experimental studies.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3B
• /
• pp.261-270
• /
• 2008

Even though there is a great deal of progress in a numerical method of high caliber like SPH, it is very rarely deployed in a water resources community. Despite the great stride in computing environment, depth averaged approach like a nonlinear shallow equation is still efficient tool for flood routing in large watershed, but it can give some misleading information like the inundation height of flood. In this rationale, we numerically simulate the flow into the dry channel, dry channel with an obstacle triggered by the collapse of a two dimensional water column using SPH (Smoothed Particle Hydrodynamics) in order to boost the application of numerical method of high caliber like SPH in a water resources community. As a most severe test of the robustness of SPH, we also carry out the simulation of the flow through a clearance into the wet channel driven by the rapid removal of a water gate. As a hydrodynamic model, we used the Navier-Stokes equation, a numerical integration of which was carried out using SPH. To verify the validity of newly proposed numerical model, we compare the numerically simulated flow with the others in the literature mainly from VOF and MAC, and hydraulic experiments by Martin and Moyce (1952), Koshizuka et al. (1995) and Janosi et al. (2004). It was shown that agreements between the numerical results in this study and hydraulic experiments are remarkable.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.2
• /
• pp.65.2-65.2
• /
• 2014

We present a new hydrodynamic simulation code based on the Voronoi tessellation for estimating the density precisely. The code employs both of Lagrangian and Eulerian description by adopting the movable mesh scheme, which is superior to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for all surfaces of each Voronoi cell so as to update the hydrodynamic states as well as to move current meshes. Besides, the IEM (incremental expanding method) is devised to compute the Voronoi tessellation to desired degree of speed, thereby the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.