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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.1253-1254
• /
• 2005

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09a
• /
• pp.749-750
• /
• 2005

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.213-217
• /
• 2010

댐붕괴파 (dam-break flow)나 지진해일에 의해 발생하는 undular bore와 충격파 (shock) 현상을 동수압 및 분산효과를 고려하여 수치모의를 수행하였다. 완전비선형 Boussinesq-type equations 모형을 이용하여, 동수압 및 분산 효과를 고려하였다. 방정식은 4차 정확도의 유한체적법을 이용하여 해석하였고, 시간적으로도 4차정확도의 기법을 이용하여 고차미분항에 대한 수치분산을 억제하였다. 다양한 경우의 1차원과 2차원 공간에서의 수치모의를 수행하고 검증을 수행하였다. 그 결과, 완전비선형 Boussinesq-type equations 모형은 천수방정식 (shallow water equations) 기반의 모형에서 재현이 불가능한 undular bore 등을 재현 하는 등, 전반적으로 천수방정식 기반의 모형 보다 물리적으로도 타당하고 정량적으로도 실험결과와 잘 일치하는 경향을 보였다. 즉, 댐붕괴파나 지진해일 등에 의한 범람 모의에 있어 동수압과 분산 효과의 중요성이 공학적으로도 매우 중요한 고려사항 임이 나타났다.

• Journal of Korea Water Resources Association
• /
• v.46 no.2
• /
• pp.139-153
• /
• 2013

For analyzing shallow-water flows over the uneven bottom, the HLLL scheme and the divergence form for bed slope source term (DFB) technique, respectively were applied to the flux gradient and the bottom gradient source terms in a finite-volume model for the shallow water equations. And also the model incorporated the volume/free-surface relationship (VFR) to consider the partially submerged cells (PSC). It was identified that a simpler version of the weighted surface-depth gradient method in the MUSCL was equivalent to the original one in the accuracy for 1D steady flows. It was verified that the flux gradient term and the bottom gradient source term were well-balanced exactly by the VFR for the 1D PSC. The VFR for the triangular PSC settled the problem which the governing equations were not well-balanced by the DFB technique for the 2D PSC. There were good agreements in simulations and experiments for 2D dam-break flows over a triangular sill and a round bump. In addition, the partial dam-break flow was successfully simulated for flooding of roughnesses in an irregular bottom as well as a sloping one. Therefore, this model is expected to be applied to the real river with uneven topography.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• Journal of Korea Water Resources Association
• /
• v.36 no.1
• /
• pp.1-11
• /
• 2003

This paper describes a finite volume, two-dimensional model. It adopts a recently developed essentially non-oscillatory(ENO) schemes based on the Lax-Friedrichs solver, which was modified for a finite volume grid, and employs a modified MUSCL(Monotonic Upstream centered Scheme for Conservation Law) for second-order accuracy in space. To demonstrate the applications of the model, it is applied to solve the 1-D and 2-D dam-break problems. The model in conjunction with the modified MUSCL showed a better agreement with analytical solutions than the minmod function in 1-D dam-break problems and is satisfactorily validated with documented published data in 2-D dam-break problems. The model was applied to tidal wane entering channel at one end, and the results showed a good agreement with analytical solutions. In the channel with reflective boundary conditions specified at the extremities, the model was capable of accurately simulating the wave propagation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.2B
• /
• pp.121-129
• /
• 2009

In this study, dam-break flows are simulated numerically by using an efficient and accurate Cartesian cut-cell mesh system. In the system, most of the computational domain is discretized by the Cartesian mesh, while peculiar grids are done by a cutcell mesh system. The governing equations are then solved by the finite volume method. An HLLC approximate Riemann solver and TVD-WAF method are employed to calculation of advection flux of the shallow-water equations. To validate the numerical model, the model is applied to some problems such as a steady flow convergence on an ideal bed, a steady flow over an irregular bathymetry, and a rectangular tank problem. The present model is finally applied to a simulation of dam-break flow on an experimental channel. The predicted water surface elevations are compared with available laboratory measurements. A very reasonable agreement is observed.

• The Journal of the Korea Contents Association
• /
• v.10 no.10
• /
• pp.451-455
• /
• 2010

Hydraulics usually deals with flows with free surface. When the surface curvature is small, the assumption of hydrostatic pressure distribution is enough. However, in the case, when the curvature is big, the non-hydrostatic pressure distribution should be taken into account and the Navier-Stokes equations should be employed instead of the depth-averaged shallow water equations. For the simulation of two immiscible fluids with different characteristics (e.g. water and air, water and oil), the level set method is selected for this purpose. The developed model is applied to classical dam break problem and the computational results are compared with the experimental data. The effectiveness of the developed model is confirmed.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.2
• /
• pp.90-96
• /
• 2021

In this study, the method of determining the state of grid points in the adaptable surface particle method based on grid system developed as a free-surface tracing method was improved. The adaptable surface particle method is a method of determining the state of the grid point according to the shape of the free-surface and obtaining the intersection of the given free-surface and grid line where the state of the grid point changes. It is difficult to determine the state of grid points in the event of rapid flow, such as collision or separation of free-surfaces, and this study suggests a method for determining the state of current grid points using the state of surrounding grid points where the state of grid point are known. A grid layer value was assigned sequentially to a grid away from the free-surface, centering on the boundary cell where the free-surface exists, to identify the connection information that the grid was separated from the free-surface, and to determine the state of the grid point sequentially from a grid away from the free-surface to a grid close to the free-surface. To verify the improved method, a numerical analysis was made on the problem of dam break in which a sudden collision of free-surface occurred and the results were compared, and the results were relatively reasonable.