k-ε 모형을 이용한 불연속 유입 밀도류의 수치모의 = Numerical Simulations of Discontinuous Density Currents using k-ε Model

국문 초록 (Abstract)

본 연구에서는 2-방정식 난류모형인 ${\kappa}-{\varepsilon}$ 모형을 이용하여 2차원적으로 발달하는 하층 밀도류를 모의하기 위한 수치모형을 제시하였다. 타원형 편미분 방정식으로 흐름에 대한 지배방정식을 구성하였으며, 수치기법으로는 유한체적법을 사용하였다. 연구를 통해 개발된 수치모형을 기존의 실험결과와 비교하여 적용성을 검토하였다. 불연속 유입조건 밀도류의 시간별 진행에서 흐름의 불안정에 의한 와(渦)의 발생을 확인하였으며, 진행중인 불연속 유입 밀도류의 진행단계 변화와 선단부 속도의 관계를 살펴보았다. 또한, 무거운 유체 위로 이동하던 주변수체가 벽면을 접하면서 발생된 단파가 밀도류의 진행에 미치는 영향을 수직구조를 통해 살펴보았다. 한편, 개발된 모형을 이용한 수치모의에서 선단부의 진행 속도가 실험보다 느리게 모의되었다. 이는 3차원 현상인 난류의 불안정성을 모의하는 2차원 모형의 한계인 것으로 판단된다.

다국어 초록 (Multilingual Abstract)

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.

참고문헌 (Reference)

1 최성욱, "부유사 밀도류 점화현상의 수치해석" 대한토목학회 21 (21): 193-200, 2001

2 Choi, S.-U., "k-ε turbulence modeling of density currents developing two dimensionally on a slope" ASCE 128 (128): 55-63, 2002

3 Eghbalzadeh, A., "Two- and three-dimensional unsteady RANS modeling of discontinuous gravity currents in rectangular domains" ASCE 2009

4 Ellison, T.H., "Turbulent entrainment in stratified flows" 6 : 423-448, 1959

5 Holmes, P., "Turbulence, coherent structures, dynamical systems and symmetry" Cambridge University Press 1996

6 Huppert, H.E., "The slumping of gravity currents" 99 (99): 785-799, 1980

7 Launder, B.E., "The numerical computation of turbulent flows" 3 : 269-289, 1974

8 Patterson, M.D., "Numerical modelling of two-dimensional and axisymmetric gravity currents" 47 : 1221-1227, 2005

9 Patankar, S.V., "Numerical heat transfer and fluid flow"

10 Cantero, M.I., "Multidimensional CFD simulation of a discontinuous density current" 2003

11 Choi, S.-U., "Modeling of one-dimensional turbidity currents with a dissipative Galerkin finite element method" IAHR 30 (30): 623-648, 1995

12 Ghisalberti, M., "Mixing layers and coherent structures in vegetated aquatic flows" 107 : 1-11, 2002

13 Rottman, J.W., "Gravity currents produced by instantaneous releases of a heavy fluid in a rectangular channel" 135 : 95-110, 1983

14 Simpson, J.E., "Gravity Currents" Cambridge University Press. 1997

15 Chikita, K., "Field study on turbidity currents initiated from spring runoffs" 25 (25): 257-271, 1989

16 Hartel, C., "Analysis and direct numerical simulation of the flow at a gravity-current head. Part 2. The lobe-and-cleft instability" 418 : 213-229, 2000

17 Hartel, C., "Analysis and direct numerical simulation of the flow at a gravity-current head. Part 1. Flow topology and front speed for slip and no-slip boundaries" 418 : 189-212, 2000

18 Versteeg, H.K., "An Introduction to Computational Fluid Dynamics - The Finite Volume Method" Prentice Hall 2007