• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.231-237
• /
• 2009

This study presents a numerical model to simulate density currents developing two dimensionally. The ${\kappa}-{\varepsilon}$ model is used for the turbulence closure. Elliptic flow equations are solved by the finite volume method. In order to investigate the applicability of the numerical model, discontinuous density currents are simulated numerically. The vortices due to the instability at the interface are simulated, showing a good agreement with the experimental visualizations in the literature. It is also investigated that the transition from slumping phase to inertial phase occurs when a bore generated at the end wall overtakes the front. However, the propagation of the density current is retarded compared with the experimental results. Two-dimensional modeling seems to have an effect on underestimating the front velocity of the density current.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.159-163
• /
• 2008

본 연구에서는 2-방정식 난류모형을 이용하여 사면을 따라 발달하는 하층밀도류의 시간에 따른 이동현상에 대해 살펴보았다. 이를 위해 타원형의 편미분 방정식을 지배방정식으로 구성하고, 난류 완결을 위해 k-$\varepsilon$ 난류모형을 이용하였다. 개발된 모형을 이용하여 경사의 사면을 따라 진행하는 연속 유입 밀도류를 수치모의 하였다. 완전 발달된 하층밀도류의 거리에 따른 주흐름방향 유속 분포, 체적 농도를 계산하였고, 이를 기존의 실험결과와 비교하였다. 실험과 수치모의 결과가 잘 일치함을 확인하였다. 또한, 불연속 유입 밀도류의 시간에 따른 진행 상황을 수치모의하여 밀도류와 주변수체의 경계부에서 Kelvin-Helmholtz 불안정에 의한 와(渦)가 형성되는 것을 확인하였으며, 밀도류 선단부의 진행 속도와 주변수체의 유입에 대해 고찰하였다.

• Journal of Korea Water Resources Association
• /
• v.50 no.1
• /
• pp.55-64
• /
• 2017

This study presents a numerical model for simulating turbidity currents using the ULTIMATE scheme. For this, the layer-averaged model is used. The model is applied to laboratory experiments, where the flume is composed of sloping and flat parts, and the characteristics of propagating turbidity currents are investigated. Due to the universal limiter of the ULTIMATE scheme, the frontal part of the turbidity currents at a sharp gradient without numerical oscillations is computed. Simulated turbidity currents propagate super-critically to the end of the flume, and internal hydraulic jumps occur at the break-in-slope after being affected by the downstream boundary. It is found that the hydraulic jumps are computed without numerical oscillations if Courant number is less than 1. In addition, factors that affect propagation velocity of turbidity currents is studied. The particle size less than $9{\mu}m$ does not affect propagation velocity but the buoyancy flux affects clearly. Finally, it is found that the numerical model computes the bed elevation change due to turbidity currents properly. Specifically, a discontinuity in the bed elevation, arisen from the hydraulic jumps and resulting difference in sediment entrainment, is observed.