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

• Journal of Korea Game Society
• /
• v.11 no.6
• /
• pp.193-199
• /
• 2011

Producing high density crowd simulation is time-consuming task as increasing the number of individuals in the crowds. In this paper, we propose a new control technique that can create realistic high density crowd simulation by using Smoothed Particle Hydrodynamics (SPH) method from fluid simulation field. Equations in SPH method are modified for evacuation, distance maintenance, and group maintenance forces for individual behaviors in the crowds. Experimental results showed that the proposed system could enable natural high density crowd simulation efficiently.

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

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.752-762
• /
• 2021

The Smoothed Particle Hydrodynamics is one of the most widely used mesh-free numerical method for thermo-fluid dynamics. Due to its Lagrangian nature and simplicity, it is recently gaining popularity in simulating complex physics with large deformations. In this study, the 3D single/two-phase numerical simulations are performed on the Liquid Metal Reactor (LMR) centralized sloshing benchmark experiment using the SPH parallelized using a GPU. In order to capture multi-phase flows with a large density ratio more effectively, the original SPH density and continuity equations are re-formulated in terms of the normalized-density. Based upon this approach, maximum sloshing height and arrival time in various experimental cases are calculated by using both single-phase and multi-phase SPH framework and the results are compared with the benchmark results. Overall, the results of SPH simulations show excellent agreement with all the benchmark experiments both in qualitative and quantitative manners. According to the sensitivity study of the particle-size, the prediction accuracy is gradually increasing with decreasing the particle-size leading to a higher resolution. In addition, it is found that the multi-phase SPH model considering both liquid and air provides a better prediction on the experimental results and the reality.

• Journal of Mechanical Science and Technology
• /
• v.16 no.3
• /
• pp.395-402
• /
• 2002

This paper shows how to formulate the transient analysis of 2-dimensional Hagen-Poiseuille flow using smoothed particle hydrodynamics (SPH). Treatments of viscosity, particle approximation and boundary conditions are explained. Numerical tests are calculated to examine effects caused by the number of particles, the number of particles per smoothing length, artificial viscosity and time increments for 2-dimensional Hagen-Poiseuille flow. Artificial viscosity for reducing the numerical instability directly affects the velocity of the flow, though effects of the other parameters do not produce as much effect as artificial viscosity. Numerical solutions using SPH show close agreement with the exact ones for the model flow, but SPH parameter must be chosen carefully Numerical solutions indicate that SPH is also an effective method for the analysis of 2-dimensional Hagen-Poiseuille flow.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.756-759
• /
• 2011

일반적으로 컴퓨터를 이용한 수치 해석에는 격자 수치 해석 방법인 유한요소법 또는 유한차분법이 주로 사용되어 왔다. 그러나 이러한 방법들은 해석하고자 하는 영역을 요소나 격자 등으로 분할해야 하기 때문에 복잡한 현상들을 다루는 데 어려움을 갖게 된다. 이를 극복하기 위해 개발된 방법이 무요소법(Meshfree Method)이며 본 논문에서는 다양한 무요소법들 중 SPH(Smoothed Particle Hydrodynamics)가 고려되어진다. SPH는 라그랑지안 수치 근사 기법을 사용하는 입자법(Particle Method)으로 SPH를 정확하게 실행하기 위해서는 적절한 경계 처리법이 요구된다. 그러나 기존의 경계 처리법은 유체 입자의 침투현상 및 커널(Kernel) 끊김 현상이 발생하기 때문에 적합하지 않다. 따라서 지금까지 SPH의 경계 처리법을 향상시키기 위해 다양한 접근법들이 제안되었으며 본 논문에서는 이러한 접근법들 중 정반사(Specular Reflection), 재회복(Bounce-back), 재도입(Reintroduce) 방법 및 경계 반발력(Repulsive Force)과 가상 입자(Ghost Particle)의 적용이 분석되고 현상 접목을 통해 적절한 경계 처리법이 제안되어진다.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Structural Engineering and Mechanics
• /
• v.78 no.4
• /
• pp.497-518
• /
• 2021

In this study, a new method for treating the wall boundary in smoothed particle hydrodynamics (SPH) is proposed to simulate free surface flows effectively. Unlike conventional methods of wall boundary treatment through boundary particles, in the proposed method, the wall boundary condition is directly imposed by adding boundary truncation terms to the mass and momentum conservation equations. Thus, boundary particles are not used in boundary modeling. Doing so, the wall boundary condition is accurately imposed, boundary modeling is simplified, and computation is made efficient without losing stability in SPH. Performance of the proposed method is demonstrated through several numerical examples: dam break, dam break with a wedge, sloshing, inclined bed, cross-lever rotation, pulsating tank and sloshing with a flexible baffle. These results are compared with available experimental results, analytical solutions, and results obtained using the boundary particle method.

• Nuclear Engineering and Technology
• /
• v.53 no.10
• /
• pp.3264-3274
• /
• 2021

In a severe accident of light water reactor (LWR), molten core material (corium) can be released into the wet cavity, and a fuel-coolant interaction (FCI) can occur. The molten jet with high speed is broken and fragmented into small debris, which may cause a steam explosion or a molten core concrete interaction (MCCI). Since the premixing stage where the jet breakup occurs has a large impact on the severe accident progression, the understanding and evaluation of the jet breakup phenomenon are highly important. Therefore, in this study, the jet breakup simulations were performed using the Smoothed Particle Hydrodynamics (SPH) method which is a particle-based Lagrangian numerical method. For the multi-fluid system, the normalized density approach and improved surface tension model (CSF) were applied to the in-house SPH code (single GPU-based SOPHIA code) to improve the calculation accuracy at the interface of fluids. The jet breakup simulations were conducted in two cases: (1) jet breakup without structures, and (2) jet breakup with structures (control rod guide tubes). The penetration depth of the jet and jet breakup length were compared with those of the reference experiments, and these SPH simulation results are qualitatively and quantitatively consistent with the experiments.