• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.211-235
• /
• 2013

An incompressible smoothed particle hydrodynamics (ISPH) method based on the incremental pressure projection method is developed in this study. The Rayleigh-B$\acute{e}$nard convection in a square enclosure is used as a validation case and the results obtained by the proposed ISPH model are compared to the benchmark solutions. The comparison shows that the established ISPH method has a good performance in terms of accuracy. Subsequently, the proposed ISPH method is employed to simulate natural convection from a heated cylinder in a square enclosure. It shows that the predictions obtained by the ISPH method are in good agreements with the results obtained by previous studies using alternative numerical methods. A rotating and heated cylinder is also considered to study the effect of the rotation on the heat transfer process in the enclosure space. The numerical results show that for a square enclosure at, the addition of kinetic energy in the form of rotation does not enhance the heat transfer process. The method is also applied to simulate forced convection from a circular cylinder in an unbounded uniform flow. In terms of results, it turns out that the proposed ISPH model is capable to simulate heat transfer problems with the complex and moving boundaries.

• 한국건축시공학회지
• /
• 제19권1호
• /
• pp.39-46
• /
• 2019

본 연구에서는 비압축성 Navier-Stokes 방정식을 적용한 ISPH 기법을 이용하여 3차원 유동 수치해석 모델을 개발하였다. 수치해석을 위해 MATLAB을 사용하여 ISPH 프로그램을 구현하였다. ISPH의 커널 함수로 piecewise cubic spline 함수를 사용하였다. 벽 경계조건으로 고정 가상 입자를 사용하였으며, 가상 밀도를 적용하여 자유 표면 경계 부근의 입자들을 결정하였다. 수치해석 모델과 코드의 정도를 확인하기 위해 $T_{500}$ 시험, 슬럼프 플로우 시험, L-box 시험의 수치해석 결과와 실험 결과를 비교하였다. 수치해석 결과 고유동 콘크리트의 점성계수 및 항복응력 변화에 따른 유동 현상의 특성을 잘 묘사하였으며, 기존의 실험값과 비교적 잘 일치함을 확인할 수 있었다.

• 대한조선학회논문집
• /
• 제48권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.

• 한국기계가공학회지
• /
• 제16권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.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제9권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
• /
• 제12권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.