• 한국방재학회:학술대회논문집
• /
• 2009.02b
• /
• pp.76-76
• /
• 2009

• 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
• /
• 2015.05a
• /
• pp.379-379
• /
• 2015

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

• Journal of Korea Water Resources Association
• /
• v.41 no.7
• /
• pp.669-679
• /
• 2008

The present study examines similarity of behavior between an overtopping wave generated by a plunging wave and a dam-break flow through hydraulic model tests. The dam-break flow has been employed to estimate the overtopping effect on the basis of the dam-break flow's behavior similar to the overtopping. In this study, the overtopping velocity was measured by a modified image technique using bubble and bubble texture images called bubble image velocitmetry. From the measurements, the vertical profiles of horizontal overtopping velocity at cross-sections along the deck were presented and discussed. Maximum velocity and depth-averaged velocity at each cross-section were compared with an analytical solution solving the dam-break flow, Ritter's solution. The initial water depth of importance for the solution was determined from the tested wave condition and the overtopping measurements. The comparison shows that the solution with the initial water depth estimated using the front velocity of the overtopping wave is in good agreements with the measurements.

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

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

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.550-554
• /
• 2009

본 연구에서는 댐 붕괴시 발생하는 홍수파가 건물군으로 이루어진 가상의 도시지역을 관통하면서 이동될 경우 시간에 따라 변화되는 특성을 수치적으로 분석하였다. 분석에 적용된 수치모형은 2차원 유한체적법과 불연속 흐름을 모의하기 위해 시간과 공간상으로 1차 정확도를 가진 HLLC 기법을 적용한 모형이다. 본 연구에서 적용된 가상의 도시지역은 Soares-Frazao와 Zech(2008)에 의해 수행된 수리모형실험에 근거하여 세 가지 case에 대해 구성하였다. 첫 번째 case(Case 1)는 25개의 건물로 구성된 정방형의 도시지역이 흐름에 직각으로 배열되어 있는 경우이며, 두 번째와 세 번째 case는 각각 흐름에 대해 $22.5^{\circ}$(Case 2)와 $45^{\circ}$(Case 3)로 경사진 경우이다. 세 경우에 대해 시간에 따른 홍수파의 변화양상을 수치적으로 분석하였으며, 첫 번째와 두번째 case에 대해 모의된 수위변화결과를 수리모형실험결과와 비교하여 모형의 적용성을 검증하였다. 검증결과 수치모형을 통해 계산된 수위변화는 수리모형실험결과와 비교적 일치하는 결과를 얻었다. 또한 본 연구에서는 도시지역의 조도계수(Manning계수: 0.01,015,$0.025sm^{-1/3}$)와 비구조적 격자시스템을 구성하는 격자수(T1: 19759개 요소, T2: 193859)에 따른 홍수파의 변화양상을 분석하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1383-1393
• /
• 2014

Recently, with rapid improvement in computer hardware and theoretical development in the field of computational fluid dynamics, high-order accurate schemes also have been applied in the realm of computational hydraulics. In this study, numerical solutions of 1D shallow water equations are presented with TVD Runge-Kutta discontinuous Galerkin (RKDG) finite element method. The transcritical flows such as dam-break flows due to instant dam failure and transcritical flow with bottom elevation change were studied. As a formulation of approximate Riemann solver, the local Lax-Friedrichs (LLF), Roe, HLL flux schemes were employed and MUSCL slope limiter was used to eliminate unnecessary numerical oscillations. The developed model was applied to 1D dam break and transcritical flow. The results were compared to the exact solutions and experimental data.