• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.149-149
• /
• 2020

하천에서 물 흐름이 보와 댐과 같은 수공구조물을 지날 때 일반적으로 흐름상태에 다양하고 급진적인 변화가 발생한다. 특히 흐름이 구조물을 지나면서 사류(supercritical flow)로 변하고 다시 상류(subcritical flow)로 복원되면서 일어나는 도수(hydraulic jump) 현상은 수위의 급변화, 흐름 에너지 소산, 변동성이 강한 압력 분포 등이 특징이다. 이러한 흐름 특성들은 보나 여수로와 같은 수공구조물 자체의 성능뿐만 아니라 이들 수공구조물의 하류에서 발생하는 국부세굴로 인해 구조물의 안정성에 부정적인 영향을 줄 수 있다. 따라서 수공구조물을 설계할 때는 이들 구조물을 통과하는 흐름의 비정상 난류 흐름 특성을 정확하게 해석하여 반영하여야 한다. 이 연구에서는 k-omega SST 난류 모형과 자유수면의 급격한 변동을 해석하기 위한 하이브리드-VOF(hybrid volume of fluid)기법을 이용하여 도수현상을 수치적으로 재현하고자 한다. 기존 CFD(computational fluid Dynamics) 모델링에서는 자유수면 변동의 영향을 고려하기 위해 VOF 기법을 많이 사용하였다. 하지면 전통적인 VOF 기법은 다상흐름(multiphase flow)을 오직 부피분율(volume fraction)의 함수로만 고려하며 모의함으로써 강한 수면변동뿐만 아니라 공기연행(air entrainment)를 동반하는 난류 흐름을 모의하는데는 한계가 있다. 이 연구에서 이용하는 Eulerian 기법인 하이브리드 VOF 기법은 물과 공기의 각 상에 대하여 흐름 특성들을 개별적으로 계산하기 때문에 공기연행을 포함한 급변류 흐름에서 전통적인 VOF 기법보다 적용성이 우수하다. 이와 같은 난류모형과 자유수면 포착기법을 이용하여 3차원 비정상 난류 흐름 수치모형을 구축하여 수공구조물 주변에서 발생하는 강한 공기연행과 난류 특성를 보이는 급변류를 수치적으로 재현한다. 이 연구는 계산된 수치해석 결과를 기존 수리실험 결과와 비교하여 수치모형의 적용성을 평가하고 도수 현상에서 발생하는 독특한 흐름 특성을 제시한다.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.13 no.1
• /
• pp.2206-2217
• /
• 1971

Spillway and discharge channel of reservoirs require the Control of Large volume of water under high pressure. The energies at the downstream end of spillway or discharge channel are tremendous. Therefore, Some means of expending the energy of the high-velocity flow is required to prevent scour of the riverbed, minimize erosion, and prevent undermining structures or dam it self. This may be accomplished by Constructing an energy dissipator at the downstream end of spillway or discharge channel disigned to dissipated the excessive energy and establish safe flow Condition in the outlet channel. There are many types of energy dissipators, stilling basins are the most familar energy dissipator. In the stilling basin, most energies are dissipated by hydraulic jump. stilling basins have some length to cover hydraulic jump length. So stilling basins require much concrete works and high construction cost. Flip bucket type energy dissipators require less construction cost. If the streambed is composed of firm rock and it is certain that the scour will not progress upstream to the extent that the safety of the structure might be endangered, flip backet type energy dissipators are the most recommendable one. Following items are tested and studied with bucket radius, $R=7h_2$,(medium of $4h_2{\geqq}R{\geqq}10h_2$). 1. Allowable upstream channel slop of bucket. 2. Adequate bucket lip angle for good performance of flip bucket. Also followings are reviwed. 1. Scour by jet flow. 2. Negative pressure distribution and air movement below nappe flow. From the test and study, following results were obtained. 1. Upstream channel slope of bucket (S=H/L) should be 0.25<H/L<0.75 for good performance of flip bucket. 2. Adequated lip angle $30^{\circ}{\sim}40^{\circ}$ are more reliable than $20^{\circ}{\sim}30^{\circ}$ for the safety of structures.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.4
• /
• pp.329-337
• /
• 2006

A small hydro-power plant operated by cooling water discharged from the power plant is under construction. In this study, the flow characteristics of the effluent channel and the outfall coastal zone in which the facilities are constructed have been measured and analysed. The flow pattern is highly dependent on the effluent discharge and clearly classified as these typical areas; the upstream and downstream areas of the weir, and the outfall coastal zone. The discharge and the width of the channel in the upstream area of the weir are increased step by step, so the water level fluctuation is small. The flow overtopping the weir is rapidly changing and has highly vertical fluctuation patterns after hydraulic jump just below the weir. The flow pattern in the outfall zone is directed toward the seaward direction and the velocity is dominated by the tidal level fluctuation. The mean tidal range in this area is about 10% greater than that of the Tongyeong tidal gauging station and the wave effects are negligible because of the sheltering effects of this area.

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.11a
• /
• pp.177-184
• /
• 2003

Numerical simulation on marine wind and sea surface elevation was carried out using both three-dimensional hydrostatic and non-hydrostatic models and a simple oceanic model from 0900 LST, August 13 to 0900 LST, August 15, 1995. As daytime easterly meso-scale sea-breeze from the eastern sea penetrates Kangnung city in the center part as basin and goes up along the slope of Mt. Taegullyang in the west, it confronts synoptic-scale westerly wind blowing over the top of the mountain at the mid of the eastern slope and then the resultant wind produces an upper level westerly return flow toward the East Sea. In a narrow band of weak surface wind within 10km of the coastal sea, wind stress is generally small, less than l${\times}$10E-2 Pa and it reaches 2 ${\times}$ 10E-2 Pa to the 35 km. Positive wind stress curl of 15 $\times$ 10E-5Pa $m^{-1}$ still exists in the same band and corresponds to the ascent of 70 em from the sea level. This is due to the generation of northerly wind driven current with a speed of 11 m $S^{-1}$ along the coast under the influence of south-easterly wind and makes an intrusion of warm waters from the southern sea into the northern coast, such as the East Korea Warm Current. On the other hand, even if nighttime downslope windstorm of 14m/s associated with both mountain wind and land-breeze produces the development of internal gravity waves with a hydraulic jump motion of air near the coastal inland surface, the surface wind in the coastal sea is relatively moderate south-westerly wind, resulting in moderate wind stress. Negative wind stress curl in the coast causes the subsidence of the sea surface of 15 em along the coast and south-westerly coastal surface wind drives alongshore south-easterly wind driven current, opposite to the daytime one. Then, it causes the intrusion of cold waters like the North Korea Cold Current in the northern coastal sea into the narrow band of the southern coastal sea. However, the band of positive wind stress curl at the distance of 30km away from the coast toward further offshore area can also cause the uprising of sea waters and the intrusion of warm waters from the southern sea toward the northern sea (northerly wind driven current), resulting in a counter-clockwise wind driven current. These clockwise and counter-clockwise currents much induce the formation of low clouds containing fog and drizzle in the coastal region.

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

물리적 녹조저감 기술인 커튼형 수류차단막은 유입 하천과 저수지 천이부에서 높은 영양염류와 조류를 포함한 표층 수류의 차단 또는 우회를 통해 본류 수역의 녹조발생을 저감하는 대책이다. 본 연구에서는 2009년 5월 회남대교 약 2 km 상류에 시범 설치된 대청호의 수류차단막 효과를 분석하고자 선행연구에서 보정한 2차원 횡방향 평균 수리 및 수질 모델을 최근의 수문사상인 2009년 6~8월까지 적용하여 검보정하고 수치모의를 실시하였다. 저수지 수위와 실측수위를 비교한 결과, 7월 중순 유입량 증가에 따른 수위 상승을 잘 반영하였고, 결정계수값($R^2$)이 0.997로 나타나 모델은 저수지 물수지 계산에 있어서 높은 신뢰도를 보였다. 댐앞과 회남수역에서 수심별 수온예측 오차는 AME 0.258~$1.584^{\circ}C$, RMSE 0.393~$2.548^{\circ}C$의 범위로 실측값을 잘 반영하는 것으로 나타났다. 회남, 댐, 추동, 문의 수역의 표층에서 $PO_4$-P 및 Chl-a 농도에 대한 모의값과 실측값의 시계열 비교 결과, 모델은 저수지내 각 측정지점에서 실측값의 시계열 변화를 잘 모의하였고, 회남수역에서 7월 중순 홍수 유입에 따라 증가하였다. 특히, $PO_4$-P 농도가 0.06 mg/L까지 증가하는 것으로 나타났다. 이는 홍수기에 높은 영양염류를 포함한 탁수가 수류차단막 하단을 통과하여 수리학적 도수현상(Hydraulic jump)을 일으킨 것이 원인이라 판단된다. 수류차단막 설치에 따라 회남수역의 표층에서 Chl-a 농도의 저감 효과가 두드러졌으나, 댐, 추동 및 문의수역에서의 제어 효과는 미미하였다. 이와 같이 회남수역에서 효과가 큰 이유는 2009년 수문사상의 영향과 저수지 지형특성상 유입수의 영향을 직접받기 때문으로 판단된다. 또한, 회남수역에서 수류차단막 설치에 따른 T-N 및 T-P 평균 저감 효율(수류차단막이 설치되지 않은 경우에 대한 설치 후 농도 저감 비)은 각각 10.8% 및 19.1%이었으며, 평균 저감농도는 모의값을 기준으로 각각 1.637 mg/L에서 1.461 mg/L 및 0.047 mg/L에서 0.038 mg/L로 나타났다. DIN 및 $PO_4$-P 평균 저감 효율은 각각 6.4% 및 24.6%이었고, Chl-a 평균 저감 효율은 25.5%이었으며, 평균 저감농도는 모의값을 기준으로 0.025 mg/L에서 0.018 mg/L로 나타났다. 모의결과를 종합해 볼 때, 대청호에 시범 설치된 조류제어용 수류차단막은 2009년의 수문사상에서 회남수역의 녹조발생 저감에 기여한 것으로 판단된다. 또한 대청호에서 유사한 수문사상을 보인 2006년(62일간)에 비해 2009년(28일간)에 조류주의보 발령 일수가 대폭 줄었다는 사실도 차단막의 효과를 간접적으로 확인해 준다.

• Journal of Korea Water Resources Association
• /
• v.50 no.3
• /
• pp.201-210
• /
• 2017

This study proposes a numerical model which is able to simulate turbidity currents intruding into a reservoir and resulting sediment depositions. The proposed model is applied to laboratory experiments by Toniolo and Schultz (2005), and propagation of turbidity currents, morphological change, and trap of suspended sediment are simulated. It is simulated that the turbidity current after plunging at the foreset of the model delta, propagates along the bottom. The thickness of the turbidity current increases significantly after being blocked by the dam, and this effect is propagated in the upstream direction. In addition, it is simulated that the foreset moves in the downstream direction due to both the bedload and suspended load and the thickness of the bottom set increases due to the suspended load. It is found that the height of the intake affects the thickness of the turbidity current and the location of the internal hydraulic jump. The impact of the height of the intake on the trap efficiency is not clear in the experimental results, however, overall trap efficiency is predicted quite successfully by the model. Also, sensitivity analysis is carried out, and the results indicates that the particle size affects the trap efficiency most.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5B
• /
• pp.447-457
• /
• 2006

This paper presents a numerical model for the simulations of 2D depth-averaged flows. The shallow water equations are solved numerically by the Characteristic Dissipative Galerkin (CDG) finite element method. For validation, the developed model is applied to the hydraulic jump. The computed results are compared with the analytical solution, revealing good agreement. In addition, flow in a contracting channel showing standing waves is simulated. The calculated water surface profile appears to be qualitatively consistent with the observed data. The foregoing results indicate that the model is capable of simulating the abrupt change in flow field. Next, the model is applied to the flow in a $180^{\circ}$ curved channel. The simulated results show that the velocity near the inner bank is faster than that near the outer bank and the water depth near the inner bank is shallower than that near the outer bank. However, the simulated results show that the velocity distribution across the channel is almost uniform in the bend except the reach close to the end of the bend. This is due to the limitation of the governing equations in which the transverse convection of momentum by the secondary flows along a channel bend is not taken into account.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.