• Interaction and multiscale mechanics
• /
• v.6 no.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.

• Journal of the Korea Institute of Building Construction
• /
• v.19 no.1
• /
• pp.39-46
• /
• 2019

A three-dimensional flow simulation model for high flow concrete was developed using Incompressible Smoothed Particle Hydrodynamics (ISPH), which can solved Navier-Stokes equation with the assumption of a fluid to be incompressible. For the simulation, a computer program code for ISPH was implemented with MATALB programming code. A piecewise cubic spline function was used for the kernel function of ISPH. Projetion method was used to calculate the velocity and pressure of particles as a function of time. Fixed ghost particle was used for wall boundary condition. Free surface boundaries were determined by using virtual density of particles. In order to validate the model and the code, the simulation results of slump flow test, $T_{500}$ test and L-box test were compared with experimental ones. The simulation results were well matched with the experimental results. The simulation described successfully the characteristics of the flow phenomenon according to the change of the viscosity and yield stress of high flow concrete.

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

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