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

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

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

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.3
• /
• pp.329-347
• /
• 2018

The Smoothed Particle Hydrodynamics (SPH) method has proved to have great potentials in dealing with the wave-structure interactions. Compared with the Weakly Compressible SPH (WCSPH) method, the ISPH approach solves the pressure by using the pressure Poisson equation rather than the equation of state. This could provide a more stable and accurate pressure field that is important in the study of wave-structure interactions. This paper improves the solid boundary treatment of ISPH by using a high accuracy Simplified Finite Difference Interpolation (SFDI) scheme for the 2D wave-structure coupling problems, especially for free-moving structure. The proposed method is referred as the ISPH_BS. The model improvement is demonstrated by the documented benchmark tests and laboratory experiment covering various wave-structure interaction applications.

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

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

• 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배 증가된 해석 속도를 보였다. 이를 통해 본 연구에서 제시된 모델이 실제 자연 지형에서의 적용가능성을 확인할 수 있었다.

• Proceedings of the Plant Resources Society of Korea Conference
• /
• 2018.10a
• /
• pp.114-114
• /
• 2018

Among the various volatile organic compounds (VOCs) emitted by the plant, floral scent plays a key role in attracting pollinators for reproduction and mediates ecological interactions. Floral scent is an important trait and industry drives the competition for flowers with novel scents. Chrysanthemum is one of the well-known ornamental plants and is a popular cut flower across the world. Floral scent and the genes responsible for the floral scent emission are poorly studied in chrysanthemum. In the present study, floral scent and the expression pattern of floral scent genes were analyzed in two chrysanthemum cultivars 'Golden Egg' and 'Gaya Glory'. Initially, intensity of the floral scent in five developing stages of flower including 'budding (B), bud developing (BD), initial blooming (IB), almost open (AO) and open flower (OF)' was analyzed using electronic nose (E-nose) with six metal oxide sensors. Based on the distance analysis, different stages of flower showed different relative intensity of scent according to the sensory evaluation. Although the scent pattern differed by stage, scent intensity was strongest in the OF stage in the completely opened flower in both the cultivars. Further, expression pattern of six genes in the floral scent pathway including FDS, IDI, ISPH, TPS2, TPS5 and TPS6 was observed in all the five stages of the flower in both the cultivars. The expression pattern of all the six genes differed by stage and the terpene synthase genes TPS2, TPS5 and TPS6 showed good expression levels in the $5^{th}$ flower stage compared to other stages. This study provides a preliminary data for understanding the regulation of floral scent in chrysanthemum.

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

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