• 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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.1
• /
• pp.9-15
• /
• 2014

It is not easy to analyze the behavior of a structural body composed of particles such as rocks using the finite element method facilitating typical element meshes because we cannot ignore the interactions among particles. In the study, we investigated the applicability of smooth particle hydrodynamics (SPH) element method for collision analysis of rock-berm by comparison with the conventional Lagrange method. As the result, SPH technique is expected to be capable of realistic simulation under collision analysis of material composed of particles.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.13-20
• /
• 2024

In this paper, a parallel analysis algorithm for Smoothed Particle Hydrodynamics (SPH), one of the numerical methods for fluidic materials, is introduced. SPH, which is a meshless method, can represent the behavior of a continuum using a particle-based approach, but it demands substantial computational resources. Therefore, parallel analysis algorithms are essential for SPH simulations. The domain decomposition algorithm, which divides the computational domain into partitions to be independently analyzed, is the most representative method among parallel analysis algorithms. In Discrete Element Method (DEM) and Molecular Dynamics (MD), the Cartesian coordinate-based domain decomposition method is popularly used because it offers advantages in quickly and conveniently accessing particle positions. However, in SPH, it is important to share particle information among partitioned domains because SPH particles are defined based on information from nearby particles within the smoothing length. Additionally, maintaining CPU load balance is crucial. In this study, a highly parallel efficient algorithm is proposed to dynamically minimize the size of orthogonal domain partitions to prevent excess CPU utilization. The efficiency of the proposed method was validated through numerical analysis models. The parallel efficiency of the proposed method is evaluated for up to 30 CPUs for fluidic models, achieving 90% parallel efficiency for up to 28 physical cores.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.796-808
• /
• 2019

In this paper, a Smoothed Particle Hydrodynamics (SPH) method is extended to simulate the ice failure process and ice-ship interactions. The softening elastoplastic model integrating Drucker-Prager yield criterion is embedded into the SPH method to simulate the failure progress of ice. To verify the accuracy of the proposed SPH method, two benchmarks are presented, which include the elastic vibration of a cantilever beam and three-point bending failure of the ice beam. The good agreement between the obtained numerical results and experimental data indicates that the presented SPH method can give the reliable and accurate results for simulating the ice failure progress. On this basis, the extended SPH method is employed to simulate level ice interacting with sloping structure and three-dimensional ice-ship interaction in level ice, and the numerical data is validated through comparing with experimental results of a 1:20 scaled Araon icebreaker model. It is shown the proposed SPH model can satisfactorily predict the ice breaking process and ice breaking resistance on ships in ice-ship interaction.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.1
• /
• pp.8-22
• /
• 2015

Most jet spray and atomization simulations are done with the Eulerian method which has inherent disadvantage in representing jet breakups and droplets. Full Lagrangian particles method called Smoothed Particle Hydrodynamics(SPH) is used in this work. We develop the SPH code and perform validations that confirm the suitability of our SPH method for simulating liquid jet atomization problem. Then, we conduct the simulation of liquid-liquid swirl coaxial injector for comparison against the experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.4
• /
• pp.398-407
• /
• 2007

A smoothed particle hydrodynamics (SPH) method is applied for simulating two-dimensional free-surface problems. The SPH method based on the Lagrangian formulation provides realistic flow motions with violent surface deformation, fragmentation and reunification. 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 frequency of density re-initialization, the number of particle and smoothing length. Through a series of numerical test. detailed information was obtained about how SPH solution can be more stabilized and improved by adjusting computational parameters. Finally, some numerical simulations for various fluid flow problem were carried out based on the parameters chosen through the sensitivity study.

• 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.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.2
• /
• pp.313-320
• /
• 2011

SPH(Smoothed Particle Hydrodynamics) is a gridless Lagrangian technique that is useful as an alternative numerical analysis method used to analyze high deformation problems as well as astrophysical and cosmological problems. In SPH, all points within the support of the kernel are taken as neighbours. The accuracy of the SHP is highly influenced by the method for choosing neighbours from all particle points considered. Typically a linked-list method or tree search method has been used as an effective tool because of its conceptual simplicity, but these methods have some liability in anisotropy situations. In this study, convex hull algorithm is presented as an improved method to eliminate this artifact. A convex hull is the smallest convex set that contains a certain set of points or a polygon. The selected candidate neighbours set are mapped into the new space by an inverse square mapping, and extract a convex hull. The neighbours are selected from the shell of the convex hull. These algorithms are proved by Fortran programs. The programs are expected to use as a searching algorithm in the future SPH program.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.21-28
• /
• 2024

In this study, a polyhedral domain decomposition method for Smoothed Particle Hydrodynamics (SPH) analysis is introduced. SPH which is one of meshless methods is a numerical analysis method for fluid flow simulation. It can be useful for analyzing fluidic soil or fluid-structure interaction problems. SPH is a particle-based method, where increased particle count generally improves accuracy but diminishes numerical efficiency. To enhance numerical efficiency, parallel processing algorithms are commonly employed with the Cartesian coordinate-based domain decomposition method. However, for parallel analysis of complex geometric shapes or fluidic problems under dynamic boundary conditions, the Cartesian coordinate-based domain decomposition method may not be suitable. The introduced polyhedral domain decomposition technique offers advantages in enhancing parallel efficiency in such problems. It allows partitioning into various forms of 3D polyhedral elements to better fit the problem. Physical properties of SPH particles are calculated using information from neighboring particles within the smoothing length. Methods for sharing particle information physically separable at partitioning and sharing information at cross-points where parallel efficiency might diminish are presented. Through numerical analysis examples, the proposed method's parallel efficiency approached 95% for up to 12 cores. However, as the number of cores is increased, parallel efficiency is decreased due to increased information sharing among cores.

• Interaction and multiscale mechanics
• /
• v.1 no.2
• /
• pp.251-278
• /
• 2008

In this study, we investigate the discontinuous-derivative treatment at the gas-liquid interface in underwater explosion (UNDEX) problems by using the Moving Least Squares-Smoothed Particle Hydrodynamics (MLS-SPH) method, which is known as one of the particle methods suitable for problems where large deformation and inhomogeneity occur in the whole domain. Because the numerical oscillation of pressure arises from derivative discontinuity in the UNDEX analysis using the standard SPH method, the MLS shape function with Discontinuous-derivative Basis Function (DBF) that is able to represent the derivative discontinuity of field function is utilized in the MLS-SPH formulation in order to suppress the nonphysical pressure oscillation. The effectiveness of the MLS-SPH with DBF is demonstrated in comparison with the standard SPH and conventional MLS-SPH though a shock tube problem and benchmark standard problems of UNDEX of a trinitrotoluene (TNT) charge.