• Ocean Systems Engineering
• /
• v.1 no.3
• /
• pp.207-222
• /
• 2011

A stabilized incompressible smoothed particles hydrodynamics (ISPH) method is utilized to simulate free falling rigid body into water domain. Both of rigid body and fluid domain are modeled by SPH formulation. The proposed source term in the pressure Poisson equation contains two terms; divergence of velocity and density invariance. The density invariance term is multiplied by a relaxed parameter for stabilization. In addition, large eddy simulation with Smagorinsky model has been introduced to include the eddy viscosity effect. The improved method is applied to simulate both of free falling vessels with different materials and water entry-exit of horizontal circular cylinder. The applicability and efficiency of improved method is tested by the comparisons with reference experimental results.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.390-403
• /
• 2017

Smoothed Particle Hydrodynamics (SPH) method has a good adaptability for the simulation of free surface flow problems. There are two forms of SPH. One is weak compressible SPH and the other one is incompressible SPH (ISPH). Compared with the former one, ISPH method performs better in many cases. ISPH based on Rankine source solution can perform better than traditional ISPH, as it can use larger stepping length by avoiding the second order derivative in pressure Poisson equation. However, ISPH_R method needs to solve the sparse linear matrix for pressure Poisson equation, which is one of the most expensive parts during one time stepping calculation. Iterative methods are normally used for solving Poisson equation with large particle numbers. However, there are many iterative methods available and the question for using which one is still open. In this paper, three iterative methods, CGS, Bi-CGstab and GMRES are compared, which are suitable and typical for large unsymmetrical sparse matrix solutions. According to the numerical tests on different cases, still water test, dam breaking, violent tank sloshing, solitary wave slamming, the GMRES method is more efficient than CGS and Bi-CGstab for ISPH method.

• Ocean Systems Engineering
• /
• v.12 no.1
• /
• pp.63-86
• /
• 2022

A refined projection-based purely Lagrangian meshfree method is presented towards reliable numerical analysis of fluid flow interactions with saturated/unsaturated porous media of uniform/spatially-varying porosities. The governing equations are reformulated on the basis of two-phase mixture theory with incorporation of volume fraction. These principal equations of mixture are discretized in the context of Incompressible SPH (Smoothed Particle Hydrodynamics) method. Associated with the consideration of governing equations of mixture, a new term arises in the source term of PPE (Poisson Pressure Equation), resulting in modified source term. The linear and nonlinear force terms are included in momentum equation to represent the resistance from porous media. Volume increase of fluid particles are taken into consideration on account of the presence of porous media, and hence multi-resolution ISPH framework is also incorporated. The stability and accuracy of the proposed method are thoroughly examined by reproducing several numerical examples including the interactions between fluid flow and saturated/unsaturated porous media of uniform/spatially-varying porosities. The method shows continuous pressure field, smooth variations of particle volumes and regular distributions of particles at the interface between fluid and porous media.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.93-100
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved successively in the moving least square sense. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as a FVM.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.27-27
• /
• 2021

세굴은 유체와 유사의 상호작용으로 발생하는 중요한 자연 현상 중 하나로, 구조 및 지반 붕괴, 홍수, 생태계 파괴 등의 문제를 야기할 수 있다. 이러한 세굴 현상을 예측하기 위해 많은 수치적 연구가 진행되어왔지만, 대부분의 연구가 기존 격자기반방법인 유한체적법 (FVM)과 개별요소법 (DEM)이 연성된 모델을 이용하였고, 이는 격자 의존도로 인한 정확도와 효율성의 문제점을 보였다. 해결책으로 입자기반 유체해석 방법인 약압축성 SPH (WCSPH)와 개별요소법의 결합모델을 이용한 모의가 연구되어 왔지만, 단순 밀도차를 활용한 유체해석방법이 압력의 불안정성을 야기하여 유사의 운동에도 영향을 주는 결과를 보였다. 또한, 개별요소법의 특성상 모의 입자의 크기를 실제 실험 입자의 크기와 동일하게 설정하면서 입자수가 지나치게 증가해 계산의 효율성이 현저히 낮아지게 되었고, 이로 인해 실제 자연 지형에 적용하는데 어려움을 보여주었다. 본 연구에서는 향상된 세굴 수치모의해석을 위해 반복법을 통해 안정적인 유체 압력을 계산하는 비압축성 SPH (ISPH)와 개별요소법을 연성한 ISPH-DEM 모델을 사용하였다. 또한, 계산속도 향상을 위해 하나의 입자가 다수의 작은 입자의 움직임을 대표하는 Coarse-grained 방법을 적용하여 기존 모델을 개선하였다. 개선된 모델을 NFLOW ISPH PURPLE 소프트웨어를 이용하여 세굴 현상을 수치 모의하였고 실험 결과와 검증을 진행한 결과, 세굴의 깊이, 너비, 형상 등을 비교하였을 때 약 10% 이내의 오차를 보였고, Coarse-grained 방법을 통한 입자 수 감소로 최소 13배 증가된 해석 속도를 보였다. 이를 통해 본 연구에서 제시된 모델이 실제 자연 지형에서의 적용가능성을 확인할 수 있었다.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.16 no.3
• /
• pp.94-99
• /
• 2017

The buoy of the wave energy converter moves by direct contact with the fluid. In order to design a buoy by using the numerical method, it is necessary to analyze not only the contact with the fluid but also the exact behavior of the fluid. In this paper, differences between weakly compressible smoothed particle hydrodynamics (WCSPH) and incompressible smoothed particle hydrodynamics (ISPH) are compared and analyzed for two-dimensional dam breaking simulation. ABAQUS, which is a commercial analysis program, is used for WCSPH analysis. A laboratory code is developed for ISPH analysis. The surface shape, the velocity, and the pressure pattern of the fluid are compared. The results of the laboratory code show the similar tendencies with those of ABAQUS, and there is a little difference in the pressure result.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.239-244
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.6
• /
• pp.560-568
• /
• 2011

In SPH(Smoothed Particle Hydrodynamics) method, the fluid has been assumed that it is weakly compressible to solve the basic equations composed of Navier-Stokes equations and continuity equation. That leads to some drawbacks such as non-physical pressure fluctuations and a restriction as like small time steps in computation. In this study, to improve these problems we assume that the fluid is incompressible and the velocity-pressure coupling problem is solved by a projection method(that is, by ISPH method). The two-dimensional computation results of dam breaking and gravitational wave generation are respectively compared with the results of finite volume method and analytical method to confirm the accuracy of the present numerical computation technique. And, the agreements are comparatively acceptable. Subsequently, the green water simulations of a two-dimensional fixed barge are carried out to inspect the possibility of practical application to ship hydrodynamics, those correspond to one of the violent free surface motions with impact loads. The agreement between the experimental data and the present computational results is also comparatively good.

• Journal of the Korea Computer Graphics Society
• /
• v.26 no.5
• /
• pp.1-13
• /
• 2020

In this paper, we propose a new method to stable simulation the directional ice shape by coupling of freezing solver and viscous water flow. The proposed ice modeling framework considers viscous fluid flow in the direction of ice growth, which is important in freezing simulation. The water simulation solution uses the method of applying a new viscous technique to the IISPH(Implicit incompressible SPH) simulation, and the ice direction and the glaze effect use the proposed anisotropic freezing solution. The condition in which water particles change state to ice particles is calculated as a function of humidity and new energy with water flow. Humidity approximates a virtual water film on the surface of the object, and fluid flow is incorporated into our anisotropic freezing solution to guide the growth direction of ice. As a result, the results of the glaze and directional freezing simulations are shown stably according to the flow direction of viscous water.

• Journal of the Korea Computer Graphics Society
• /
• v.18 no.3
• /
• pp.17-27
• /
• 2012

In this paper, we describe a case study of the film 'Sector 7' which was produced by technologies applied fluid simulation. For the CG scenes in the movie which include highly detailed fluid motions, we used smoothed particle hydrodynamics(SPH) technique to express subtle movements of seawater from a crashed huge tank, and used hybrid simulation method of particles and levelsets to describe bursting water from a submarine's broken canopy. We also used detonation shock dynamics(DSD) technique for detailed flame simulations to produce a burning monster, the film"s main character. At this point, the divergence-free vortex particle method was applied to conserve the incompressible property of fluids. In addition, we used an upsampling method to achieve more efficient video production. Consequently, we could produce the high-quality visual effects by using the domestic technologies.