• Journal of Korea Water Resources Association
• /
• v.56 no.10
• /
• pp.603-617
• /
• 2023

This study aims to provide a better understanding of the turbulent flow characteristics in swash zone. A double dam-break method is employed to generate the swash zone flow. Comparing with the conventional single dam-break method, a delay between two gate opening can be controlled to reproduce various interactions between uprush and backwash. For numerical simulations, overInterDyMFoam based on OpenFOAM is adopted. Using overInterDyMFoam, interface between two immiscible fluids having different densities (i.e., air and water phases) can be tracked in a moving mesh with multiple layers. Two-dimensional Reynolds-Averaged Navier-Stokes equations are solved with a standard 𝜅-𝜖 turbulence model for momentum and continuity. Numerical model results are validated with laboratory experiment data for the time series of water depth and streamwise velocity. Turbulent kinetic energy distribution is further investigated to identify the turbulence evolution for each flow regime (i.e., uprush, backwash, and swash-swash interaction).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.302-302
• /
• 2015

본 연구에서는 PMP 조건하에서 농업용 목적의 시례저수지의 가상붕괴에 따른 홍수파 해석을 통해서 하류부에 대한 수리학적 분석을 실시하였다. 그리고, 사력댐으로 구성된 댐 제체의 붕괴로 인한 유사 이송을 모의하고 이로 인한 하류부에서의 영향을 분석하였다. 이를 위하여 경상남도 김해에 위치한 시례저수지 유역을 대상으로 수리 수문자료 및 하천현황을 상세히 조사하고 댐의 가상 붕괴로 인한 하류부 홍수파 해석 및 댐 본체의 유실에 따른 유사의 이송으로 인해 발생 가능한 하류부에서의 영향을 분석하였다. 시례저수지의 댐 형식은 흙댐으로써, 구성방식은 침하량이 적은 중심코아형이며 댐 높이는 34.3 m, 유역면적은 3.13 km이다. 본 연구에서는 하천에서의 흐름과 유사이송 계산이 가능한 1차원 수치 모형인 CCHE1D를 이용하였다. 적용 모형은 확산파 및 동역학적 방정식을 모두 고려한다. 또한, 비평형 이송 모형을 사용하여 유사의 발생과 퇴적을 산정하고 비균질 상태의 유사 이송을 계산 하는데 있어서 하천 단면형의 변화와 하상물질의 구성 상태, 제방침식, 수로 확폭의 과정과 함께 산정된다. CCHE1D는 하상물질의 공극률과 non-equilibrium adaptation length, mixing layer thickness 등의 여러 변수들을 현재까지 개발된 식들을 제공하며, 토사이동해석을 위한 공식은 SEDTRA 모듈, Wu et al(2000) 공식, 수정 Ackers and White(1973) 공식, 수정 Engelund and Hansen(1967)이 사용된다. 비균질 유사이송과 bed deformation, bed-material sorting은 완전연계과정(coupling)로, 유사와 흐름방정식의 계산은 반연계과정(semi coupling)을 적용하여 계산하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.67-67
• /
• 2016

최근 기후온난화로 인한 이상기후와 지진의 발생가능성으로 인해 직렬 혹은 병렬로 위치한 2개 이상의 댐(저수지)들의 연속 붕괴 가능성도 점점 커지고 있어 연속 댐 붕괴에 대한 비상대처계획도 수립에 대한 관심도 증대되고 있다. 국내에서는 댐(저수지) 붕괴로 인한 극한홍수해석이나 비상대처계획 수립시 국내에서는 DAMBRK, FLDWAV, HEC-RAS와 같은 1차원 수리동역학 모형이 주로 사용되어지고 있다. 하지만 1차원 모형은 흐름을 하나의 방향으로만 한정하여 해석하고, 각각의 적용 횡단면에서 동일 수위를 가지다는 가정으로 홍수범람해석을 수행하므로 정확성뿐만 아니라 실제 적용성에서도 한계를 가질 수밖에 없다. 댐(저수지) 붕괴로 인한 홍수범람해석에서 2차원 이상의 고차원 모형은 앞서 언급한 1차원 모형에 적용된 비현실적인 가정을 포함하지 않으므로 더욱 정교하고 신뢰성 있는 결과를 얻을 수 있다. 하지만, 홍수범람해석을 위한 상용 프로그램들은 단일 댐(저수지) 붕괴의 적용에는 큰 어려움이 없으나, 지형단면을 초기에 한번만 고려하는 문제로 인하여 연속 댐(저수지) 붕괴의 고려에는 한계를 가진다. 따라서 본 연구에서는 댐(저수지)의 연속 붕괴를 해석할 수 있는 2차원 수리동역학 모형을 개발하고자 한다. 각 댐(저수지)의 붕괴 전과 후의 지형 단면을 여러 번 반영할 수 있는 모형을 개발함으로써 시간차를 두고 붕괴되는 댐(저수지)의 연속 붕괴를 해석할 수 있도록 하였다. 개발모형은 2002년 태풍 루사로 인해 실제로 붕괴된 저수지에 적용되었으며, 이 저수지들은 붕괴시 시간차이를 두고 붕괴가 이루어져 개발모형의 적용 유역으로 선택하였다. 저수지의 연속붕괴 모델링을 위해 지형자료로는 저수지 단면, 댐 제체 및 여수로 붕괴 단면, 하천 단면 그리고 홍수터 지형 반영을 위한 수치지도, 경계조건으로는 저수지로의 유입유량과 하류단 조위조건이 고려되었다. 그리고 태풍 루사 당시의 기록적인 강우를 반영하기 위해 연구유역 인근에 관측된 강우를 모형에서 하천 및 홍수터에서 고려할 수 있도록 하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.379-379
• /
• 2015

최근 몇 년간 기록적인 집중호우와 태풍으로 인해 홍수피해가 점점 증가하고 있다. 국내 대부분의 하천은 고수부지 또는 제방으로 정비되어있다. 국지적으로 발생되는 집중호우는 수위를 급격하게 증가시켜 하천의 고수부지 또는 제방범람을 야기한다. 집중호우로 인한 흐름이 고수부지 또는 제방을 지나면 흐름양상은 매우 복잡하게 변화한다. 그러나 급격히 증가하는 흐름에 의해 하천의 제방, 옹벽 또는 계단지형으로 형성된 고수부지 등을 지나는 흐름특성을 분석하는 것은 매우 어렵고 이러한 불연속 단면을 고려하여 예측하기는 불가능하다. 이와 같이 급변하는 흐름이 계단형 구조물을 지날 때 잠김/드러남 특성을 조사하기위해 수리 실험을 실시하였다. 수리실험은 길이 6 m, 폭 0.4 m로 제작된 직선수로에서 수행되었다. 급격히 증가하는 댐 붕괴 흐름을 재현하기위해 수로중간에 게이트를 설치하고 상류에 물을 채운 후 빠른 속도로 게이트를 개방하였다. 하류에는 수로를 따라 계단형 구조물을 설치하고 하류경계는 물이 빠져나갈 수 있게 하였다. 하류 초기수위는 첫 번째 계단보다 약간 낮게 설정하였다. 계단형 구조물과 그 주변에 파고계를 설치하여 시간에 따른 수위변화를 측정하였다. 다양한 수리조건을 적용하여 수리실험을 수행하였고 측정된 데이터를 이용하여 급변한 비정상흐름이 계단형 구조물을 지날 때 잠김/드러남 특성을 분석하였다. 차후 이 실험 결과는 하천의 제방 및 계단지형을 고려할 수 있는 고정확도 수치모형을 개발하는 검증자료로 유용하게 사용될 수 있을 것이다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.4B
• /
• pp.333-341
• /
• 2011

This study analyzed the propagation of dam-break waves in an area directly downstream of a dam by using 3D numerical modeling with RANS as the governing equation. In this area, the flow of the waves has three dimensional characteristics due to the instantaneous dam break. In particular, the dam-break flows are characterized by a highly unsteady and discontinuous flow, a mixture of the sharp flood waves and their reflected waves, a mixture of subcritical and supercritical flow, and propagation in a dry and movable bed. 2D numerical modeling, in which the governing equation is the shallow water equation, was regarded as restricted in terms of dealing with the sharp fluctuation of the water level at the dam-breaking point and water level vibration at the reservoir. However, in this 30 analysis of flood wave propagation due to partial dam breaking and dam-break in channels with $90^{\circ}$ bend, those phenomena were properly simulated. In addition, the flood wave and bed profiles in a movable bed with a flat/upward/downward bed step, which represents channel aggradation or degradation, was also successfully simulated.

• Journal of Korea Water Resources Association
• /
• v.47 no.11
• /
• pp.1077-1086
• /
• 2014

The flood wave generated due to dam break is affected by initial depth upstream since it is related with hydraulic characteristics propagating downstream, and flow resistance stress has influence on the celerity, travel distance, and approaching depth of shock wave in implementing numerical simulation. In this study, a shallow water flow model employing SU/PG scheme was developed and verified by analytic solutions; propagation characteristics of dam break according to flow resistance and initial depth were analyzed. When bottom frictional stress was applied, the flow depth was relatively higher while the travel distance of shock wave was shorter. In the case of Coulomb stress, the flow velocity behind the location of dam break became lower compared with other cases, and showed values between no stress and turbulent stress at the reach of shock wave. The value of Froude number obtained by no frictional stress at the discontinuous boundary was the closest to 1.0 regardless of initial depth. The adaption of Coulomb stress gave more appropriate results compared with turbulent stress at low initial depth. However, as the initial depth became increased, the dominance of flow resistance terms was weakened and the opposite result was observed.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.381-387
• /
• 2013

In compiling flood hazard maps for the case of dam-failure, a scenario-based numerical modeling approach is commonly used, involving the modeling of important parameters that capture peak discharge, such as breach formation and progress. In this study, an earth-dam-break model is constructed assuming an identical mechanism and hydraulic process for all dam-break processes. A focus of the analysis is estimation of the hydrograph at the outlet as a function of time. The constructed hydrograph then serves as an upper boundary condition in running the flood routing model downstream, although flood routing is not considered here. Validation was performed using the record of the Tangjishan dam-break in China. The results were satisfactory, with a coefficient of determination of 0.974, Nash-Sutcliffe Coefficient of Efficiency (NSC) of 0.94, and Root Mean Square Error (RMSE) of $610m^3/sec$. The proposed model will contribute to assessments of potential flood hazards caused by dam-break.

• Water for future
• /
• v.41 no.5
• /
• pp.75-80
• /
• 2008

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.222-222
• /
• 2023

오픈 소스 유체역학 소프트웨어인 OpenFOAM은 다양한 유체 흐름에 적용 가능한 프로그램들로 구성되어 있다. 이 중 interFoam은 밀도가 다른 두 유체(i.e., 물, 공기) 간의 경계를 추적하는 기법을 기반으로 한 프로그램으로, 파랑의 거동 모의에 주로 쓰이고 있다. 파생형 프로그램으로는 동격자(dynamic mesh) 및 중첩 격자 기법(overset grid method)을 interFoam에 추가한 overInterDyMFoam이 있다. 두 기법을 사용하면 각각 여러 영역에서 유체흐름과 다중 물체 간의 상호작용을 효율적으로 모의할 수 있다. 본 연구에서는 overInterDyMFoam을 사용하여 두 개 수문의 개방 움직임을 구현하고 생성된 파랑이 포말대(swash zone)에 접근하였을 때의 흐름 특성을 조사하였다. 수치모형실험 결과 수문 개방 속도가 댐 붕괴 파랑 흐름 전파속도에 영향을 미치는 사실을 발견하였다. 또한, 처오름과 처내림의 상호작용에 의한 난류 운동 특성을 조사하기 위해 수문 개방시간 간격을 0초~3초로 설정하였다. 수치모형실험 결과는 수리모형실험의 수면 변동 시계열과 속도 시계열 결과와 비교하여 모형의 정확성이 검증되었다.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.3
• /
• pp.69-78
• /
• 2007

In this study, a time-dependent aspect of an embankment failure is considered to simulate a flood inundation map and calculate overtopping discharge induced by an embankment failure. A numerical model has been developed by solving the two dimensional nonlinear shallow water equations with a finite volume method on unstructured grids. To analyze a Riemann problem, the HLLC approximate Riemann solver and the Weighted Averaged Flux method are employed by using a TVD limiter and the source term treatment is also employed by using the operator splitting method. Firstly, the numerical model is applied to a dam break problem and a sloping seawall. Obtained numerical results show good agreements with experimental data. Secondly, the model is applied to a flow induced by an embankment failure by assuming that the width and elevation of embankment are varied with time-dependent functions. As a result of the comparison with each numerical overtopping discharge, established flood inundation discharges in the previous studies are overestimated than the result of the present numerical model.