• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.72-72
• /
• 2018

우리나라의 대형 저수지에 담수된 수체는 여름에 온도가 높은 표수층과 온도가 낮은 심수층, 그리고 두 층 사이에서 온도가 크게 변화되는 변온층으로 나뉘게 된다. 여름 홍수기에 상류 하천에서 저수지로 유입되는 부유사는 큰 밀도를 가지기 때문에 저수지의 바닥을 타고 하층 밀도류의 형태로 저수지 하류로 전파된다. 그러나 밀도류가 성층화된 저수지의 변온층에 도달하면 심수층과 변온층 하층의 낮은 온도로 인해 발생하는 수체의 큰 밀도로 하층으로 더 침투하지 못하고 변온층이나 심수층 상층에서 남아 하류에 중층밀도류의 형태로 전파된다. 대량의 탁수를 하류에 방류하면 하류 수질에 문제가 발생되며 저수시키면 댐 수질 및 심수층 생물체의 태양광 차단 문제가 발생되므로 저수지에서 밀도류의 조절은 대형 저수지 운영에 매우 중요한 사항이다. 따라서 본 연구의 목적은 층적분 모형을 사용하여 대형 저수지에서 발생하는 중층밀도류의 전파 및 선택취수시설을 이용한 밀도류의 차단과 방류에 대해 수치모의하는 것이다. 이를 위하여 중층 밀도류를 수치모의하기 위한 1차원 수치모형을 제시하였으며 소양강댐 상류 소양호에 적용하여 적용성을 검토하고 밀도류의 전파 및 방류에 대해 분석하였다.

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

자연하천에서 다양한 목적으로 설치된 수공구조물을 볼 수 있으며 대표적으로 교량의 교각을 들 수 있다. 교각 주위에서 발생하는 국부세굴은 교량붕괴의 가장 큰 원인이므로 이에 대한 연구는 필수적이다. 수리실험 또는 수치모형을 이용하여 교각주변의 국부세굴에 관한 많은 연구들이 수행되었지만 세굴인자의 특성 및 메커니즘에 대한 연구는 여전히 부족하다. 본 연구에서는 LES(large-eddy simulation)에 유사이송 및 하상변동 모형과 결합하여 2개의 원형교각 주변의 흐름 및 세굴을 수치모의 하였다. LES와 유사이송 및 하상변동 모형의 결합은 난류의 영향을 직접 모형에 고려할 수 있기 때문에 교각 주변에서 발생하는 말굽형 와 구조와 같은 복잡한 흐름에 의한 영향이 반영된다. 계산영역은 흐름방향으로 10 m, 횡 방향으로 2.4 m로 하였고, 지름(D) 0.16 m를 가지는 원형교각을 유입부로부터 2.4 m 떨어진 위치에 배치하였다. 이때 두 개 교각사이의 각은 $0^{\circ}{\sim90^{\circ}$이고 원형교각 사이의 거리는 5D로 하였다. 수치모의에 사용된 조건은 이전의 수리실험(Khosronejad et al., 2012)을 참고하여 접근평균유속은 0.25 m/s, 수심은 0.15m를 사용하였다. 수치모의는 원형교각 주변의 최대세굴심이 평형상태에 이를 때 까지 수행하였다. 접근각도 변화에 따른 원형교각 주변의 세굴과정 및 특성을 분석하였으며 최대 세굴심 결과를 Hannah(1978)의 수리실험결과와 비교하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.308-308
• /
• 2017

국내 하천에 설치된 횡단구조물은 2009년을 기준으로 약 20,753개로 알려져 있으며 수위 유지, 하상고 유지 등의 목적을 가지고 설치된다. 그러나 낙차를 가진 빠른 유속의 흐름을 형성시켜, 하류 하상에서 국소 세굴을 발생시킨다. 이를 방지하기 위해 횡단구조물 하류에 굴요성 구조(flexible structure)인 돌망태, 블록공, 사석 등으로 이루어진 바닥보호공(bed protection)이 설치되나 유실, 침하 등의 문제가 빈번히 발생되고 있어, 문제를 해결하기 위한 연구가 필요하다. 본 연구는 이러한 피복 대책에서 일어날 수 있는 파괴 기구인 전단파괴, 흡출 파괴, 경계 파괴, 하상 형태 변화에 따른 하부 침식 중 흡출 파괴(winnowing failure)를 유발하는 흐름을 검토하기 위한 수치모의를 수행하였다. 이때 흡출 파괴는 바닥보호공의 공극으로 미세한 하상 재료가 난류와 침투류의 작용에 의해 침식되어 바닥보호공이 침하되는 것을 말한다. 수치모의는 전산유체동역학(Computational Fluid Dynamics, CFD) 모형인 FLOW-3D 모형을 이용하였으며, 난류 모형으로 LES 모형을 적용하고 조밀한 격자를 부여하여 바닥보호공의 공극에서 발생되는 비교적 작은 척도의 와(vortex)를 해상할 수 있도록 하였다. 수치모의에 적용된 횡단구조물은 보, 물받이공, 바닥보호공으로 구성하였으며 특히, 바닥보호공의 형상은 구체(sphere)로 가정하여 다층으로 배치하였다. 바닥보호공의 공극 또는 구체 사이에서 발생되는 유속, 압력 등의 흐름특성을 분석한 결과, 바닥보호공 두께가 두꺼울수록 흡출 파괴에 대해 안정적인 것으로 나타났다. 이는 바닥보호공 설계를 위한 기초자료로 활용될 수 있으며, 향후 입자영상유속계(Particle Image Velocimeter, PIV)와 같이 공극에서 흐름을 측정할 수 있는 방법과 병행한 연구를 수행할 수 있을 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.204-204
• /
• 2017

WRF 모형은 실제 자연에서 나타나는 대기 현상의 원인을 물리적 동적 방정식들의 항으로 표현한 수치예보모형으로 전세계의 상업적 비상업적인 수치예보모형 안에서 성능이 뛰어나다고 평가되어지고 있다. WRF 모형은 오픈소스 기반의 비상업적 모형으로 사용 및 수정이 자유로운 특징이 있으며, 위성 및 레이더와 같은 고도화된 다양한 기상관측자료를 입력자료로 활용할 수 있는 장점이 있다. WRF-HYDRO 모형은 WRF 모형이 갖는 공간적인 저해상도 문제를 해결할 수 있는 고해상도의 격자를 구축할 수 있으며 유출량과 수문 변량을 추정할 수 있는 추적 모형을 추가하여 수문학적 예측 능력을 향상하고자 개발되었다. 기존 모형과의 차별성으로는 기상인자로 인하여 도출된 지표면의 수문인자들이 시간의 변동에 따라서 다음 시간의 기상인자에 영향을 미치는 피드백 구조로 구성되어 기상과 지표면이 양방향으로 연결되는 특징이 있다. 기존 모형에 비하여 향상된 구조적인 특징은 수문학적 순환과정을 자연스럽게 재현함으로서 신뢰성 있는 결과를 도출할 수 있을 것으로 판단된다. 본 연구에서는 만경강 유역의 실제 유출 사상에 대하여 WRF-HYDRO 모형을 적용하고, 홍수통제소 관할 만경강 유역내 수문 관측소 자료와의 비교를 통해 WRF-HYDRO 모형의 적용성을 검토하였다. 수문 관측소를 통한 검토 결과를 기반으로 WRF-HYDRO 모형에서 제시된 수문-기상 정보를 통하여 만경강 유역의 홍수 사상의 발생 과정에 대한 추적 및 미계측 변량의 추정에 유용하게 사용할 수 있을 것으로 판단된다.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1997.10a
• /
• pp.147-150
• /
• 1997

일반적으로 공학적인 문제를 해결하기 위하여는 수치모형을 이용한 접근방법도 중요하지만 실험, 실측 및 관측분야도 도외시 할 수는 없다. 그 주된 이유는 이들 두 분야가 서로 별개로 독립적으로 존재하는 것이 아니라 서로 보완적이기 때문이며 따라서 이들 두 분야의 균형적인 발전없이는 아무리 최신의 모형을 사용하더라도 그 한계가 있다. 실제로 수치모형의 신뢰성을 확보하는데는 물리적 현상에 기초한 실측값들과의 비교ㆍ검증절차가 필수적 사안이다. (중략)

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.1429-1433
• /
• 2004

본 연구에서는 쇄파대에서 정현파의 쇄파에 대해 수리모형실험과 수치모형실험을 수행하였다. 수치해석 모형에서는 Reynolds 방정식을 지배방정식으로 사용하고 난류해석을 위해 $k-\varepsilon$모델을 적용하였으며, 자유수면변위를 추적하기 위해 VOF기법을 사용하져다. 사면 및 평탄지형상에서 발생하는 쇄파양상을 서로 다르게 설정하기 위해 수심과 입사파의 주기와 파고를 변화시킨다. 발생된 정현파의 파형은 해석해와 잘 일치하였으며, 입사파와 파고계가 설치된 위치에서 측정된 파고비 $H/H_0$는 관측값과 비교해 본 길과 놀은 정확도를 나타내었다.

• Journal of Korea Water Resources Association
• /
• v.31 no.5
• /
• pp.515-528
• /
• 1998

A quasi-two-dimensional model for simulating the flood plain flow is developed. The model consists, in general, of a multiply-connected network which combines the main channel and two-dimensional flood plain cells. The main channel flow is described by the Saint Venant equations for one-dimensional unsteady flow, and the flood plain flow by the cell continuity and river-or weir-type stage-discharge relations between flood plain cells. The implicit algorithm for unsteady flow in looped channel network is extended to incorporate the flood plain flow. To verify the performance of the model, it is applied to three test problems, and sensitivities to various model parameters are analyzed. It turns out that the present model gives more accurate result than that by Cunge (1975) as the shape of cross section becomes more complex and irregular. Not only the inundation of water from the main channel but the return flow from the flood plain is successfully simulated.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Journal of Korea Water Resources Association
• /
• v.45 no.6
• /
• pp.545-555
• /
• 2012

A numerical modeling has become increasingly popular and more important to the study of water waves with a rapid advancement of computer technology. However, different types of problems are induced during simulating wave motion. One of the key problems is re-reflection to a computation domain at the incident boundary. The internal wave generating-absorbing boundary conditions have been commonly used in numerical wave models to prevent re-reflection. For the Navier-Stokes equations model, the internal wave maker using a mass source function of the continuity equation has been used to generate various types of waves. Nonetheless, almost every numerical experiment is performed in two dimensions and only a few tests have been expanded to three dimensions. More recently, a momentum source function of the Boussinesq equations is applied to generate essentially directional waves in the three dimensional Navier-Stokes equations model. In this study, the internal wave maker using a momentum source function is employed to generate targeted linear waves in the three-dimensional LES model.

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

In the present paper, turbulent open-channel flows with alternate vegetated zones are numerically simulated using threedimensional model. The Reynolds-averaged Navier-Stokes Equations are solved with the ${\kappa}-{\varepsilon}$ model. The CFD code developed by Olsen(2004) is used for the present study. For model validation, the partly vegetated channel flows are simulated, and the computed depth-averaged mean velocity and Reynolds stress are compared with measured data in the literature. Comparisons reveal that the present model successfully predicts the mean flow and turbulent structures in vegetated open-channel. However, it is found that the ${\kappa}-{\varepsilon}$ model cannot accurately predict the momentum transfer at the interface between the vegetated zone and the non-vegetated zone. It is because the ${\kappa}-{\varepsilon}$ model is the isotropic turbulence model. Next, the open channel flows with alternate vegetated zones are simulated. The computed mean velocities are compared well with the previously reported measured data. Good agreement between the simulated results and the experimental data was found. Also, the turbulent flows are computed for different densities of vegetation. It is found that the vegetation curves the flow and the meandering flow pattern becomes more obvious with increasing vegetation density. When the vegetation density is 9.97%, the recirculation flows occur at the locations opposite to the vegetation zones. The impacts of vegetation on the flow velocity and the water surface elevation are also investigated.