• Journal of Ocean Engineering and Technology
• /
• v.26 no.5
• /
• pp.55-62
• /
• 2012

In this paper, the rigid lid boundary condition is applied to simulate the influence of floating structures such as ships or pontoons, and the pressure term in both the momentum equations and continuity equation are modified. The pressure of a floating structure under the free surface is dependent on the draft of the structure, generally called a ship. If the free surface is covered by a floating structure, the free surface cannot move freely. The water level should be fixed, using a rigid lid boundary condition. This boundary condition is implemented by reducing the storage area of the grid cell with a factor between zero and one. The numerical model developed by Hong (2009) is verified through a comparison with experimental results, and the influence of the reduction factor is investigated using the verified numerical model.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1037-1048
• /
• 2023

Fishways are constructed to ensure the fish migration because river-crossing structures such as dams and weirs cut off the stream longitudinal connectivity and influence on aquatic ecosystems. However, the passage efficiency of fishes varies depending on flow characteristics in the fishway and fish species. In this study, three-dimensional numerical simulations are carried out using a RANS model and the volume of fluid method for resolving free surface fluctuations to calculate the turbulent flow in the pool and weir type fishway. The Flow velocity and turbulent kinetic energy in the pool of fishway are analyzed according to variation of the upstream water level and the length of pool. The present numerical simulations reasonably well reproduce the stream flow and plunging flow characteristics in the pool. The simulation results show that the stream flow changes to the plunging flow as the length of the pool increases. When the upstream level increases, the stream flow becomes more evident. Key parameters related to the fish migration within the fishway such as the flow velocity and the turbulent kinetic energy are examined to assess the ascending possibility of fishes.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.4
• /
• pp.91-100
• /
• 2008

Numerical implementation with a Cartesian cut-cell method is conducted in this study. A Cartesian cut-cell method is an easy and efficient mesh generation methodology for complex geometries. In this method, a background Cartesian grid is employed for most of computational domain and a cut-cell grid is applied for the peculiar grids where the flow characteristics are changed such as solid boundary to enhance the accuracy, applicability and efficiency. Accurate representation of complex geometries can be obtained by using the cut-cell method. The cut-cell grids are constructed with irregular meshes which have various shape and size. Therefore, the finite volume method is applied to numerical discretization on a irregular domain. The HLLC approximate Riemann solver, a Godunov-type finite volume method, is employed to discretize the advection terms in the governing equations. The weighted average flux method applied on the Cartesian cut cell grid for stabilization of the numerical results. To validate the numerical model using the Cartesian cut-cell grids, the model is applied to the rectangular tank problem of which the exact solutions exist. As a comparison of numerical results with the analytical solutions, the numerical scheme well represents flow characteristics such as free surface elevation and velocities in x-and y-directions in a rectangular tank with the Cartesian and cut-cell grids.

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

암거는 일반적으로 용수나 배수용의 수로가 도로, 철도, 제방 등의 아래에 매설 된 수로를 지칭한다. 이러한 암거는 산업발전으로 사회기반시설의 신설 및 확충, 재정비 등으로 많이 활용되고 있다. 최근 들어 기후로 인한 재해가 급증하면서 이러한 시설물에 대한 안정성 및 관리에 대한 관심이 높아지고 있는 것이 현실이다. 특히 소하천은 집수면적 및 유로연장이 짧고 하상경사가 급하기 때문에 홍수에 취약하다. 즉, 빨라진 유속으로 인해 구조물 주변의 세굴에 의한 유실, 토사유출로 인한 하상퇴적, 부유물로 인한 차단으로 인해 통수에 지장을 받아 피해가 발생하게 된다. 이러한 암거시설물 피해는 2차 피해로 이어질 수 있으며 사회기반시설 파괴로 도시기능이 마비되고 인근 주변지역에 침수로 인한 재산 및 인명피해까지 발생시킬 수 있는 피해 잠재능력을 보유하고 있다. 그러나 피해에 대한 예방대책은 유지관리를 통해 지속적으로 관리하는 것이 대부분의 지침 등에 소개된 내용들이다. 본 연구에서는 암거를 대상으로 암거의 폐색으로 인해 암거주변에서 변화되는 흐름특성을 축소모형을 통해 검토하고자 하였다. 암거 축소모형실험은 1.5m 폭을 갖는 직선수로에서 수행하였으며, 암거모형은 도로암거표준도(2008)를 참조하여 $3m{\times}3m$ 수로암거를 대상으로 1/10 축소모형을 제작하였다. 암거유입부 퇴적으로 인한 암거의 차단률(차단면적/암거단면적)은 차단이 발생하지 않는 0% 조건에서부터 10%, 20%, 30%, 40%, 50% 조건에 대해 실험을 수행하였다. 실험결과 차단에 따른 암거 상류단의 수위는 차단이 없는 암거의 경우에 비해 차단율이 높아질수록 암거유입부 수위는 20.4% ~ 82.7% 상승하는 것으로 나타났다. 암거의 차단률이 40% 이상일 경우 높아진 수위로 인해 암거통로의 윗상면부까지 다다르고 있으며 50%일 경우 암거를 통과하는 흐름이 자유수면흐름이 아닌 오리피스 흐름이 발생하는 것으로 나타났다. 암거유입부 차단으로 인한 암거주변의 최대유속은 암거 직하류부에서 주로 발생하여 암거 유출부에서의 최대유속은 차단율이 증가할수록 선형적으로 증가하는 것으로 나타났으며 암거 유출부에서의 유속은 차단전의 조건(0%) 대비 4.2% ~ 35.5% 까지 상승하는 것으로 나타났다. 이와 같은 결과를 토대로 고려하였을 때 대부분 산지부에서 설치되는 암거의 경우 유속이 불가피하게 증가하게 됨으로 유속에 따른 유속조절방안(차단 및 우회시설) 및 세굴대책을 세워야 할 것으로 판단된다.

• Journal of Wetlands Research
• /
• v.23 no.3
• /
• pp.242-251
• /
• 2021

Surface water pollution is a serious environmental problem in developing countries, like India, due to the unregulated discharge of untreated wastewater. To overcome this, the constructed wetlands (CWs) have been proven to be an efficient technology for wastewater treatment. In this study, different existing and experimental facilities were reviewed to be able to determine the current status of constructed wetlands in India. Based on the collected data from published literature, industrial wastewater contained the highest average chemical oxygen demand (COD), biochemical oxygen demand (BOD). In terms of total nitrogen (TN), Total phosphorous (TP), the lowest concentration was found on domestic wastewater. Vertical flow constructed wetlands (VFCW) and Horizontal flow constructed wetland (HFCW) were more effective in removing TSS, BOD, TP in domestic and industrial wastewater, whereas hybrid constructed wetlands (HCW) showed the highest removal for COD. The use of constructed wetlands as advanced wastewater treatment facilities in India yielded better water quality. The treatment of wastewater using constructed wetlands also enabled further reuse of wastewater for irrigation and other agricultural purposes. Overall, this study can be beneficial in evaluating and promoting the use of constructed wetlands in India.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.457-457
• /
• 2021

인공습지는 자연이 가진 정화기작을 인위적으로 증가시키기 위하여 조성한 자연기반해법에 해당한다. 인공습지는 습지 내 식물, 미생물, 토양 등의 상호기작에 의하여 오염물질이 제거된다. 인공습지의 오염물질저감효율은 시설의 규모와 유량, 유입물질의 부하량 수리학적 부하량, 체류시간 등의 영향을 받게 된다. 일반적으로 인공습지 적정 규모는 유역 및 기상인자의 특성과 조성목적에 고려하여 산정된다. 본 연구는 전국 35개 지역에 설치된 54개 인공습지를 선정하여 모니터링을 수행하였으며, 2011년부터 2018년에 설치된 시설이다. 54개 시설 중 도심지역에 13개, 농업지역 25개, 공업지역 3개, 상업지역 3개, 축산 10개가 설치되어있다. 습지형태는 Cell형 자유수면형 인공습지(Free Water Surface, Cell-FWS), 유로형(Flow) 자유수면형 인공습지(Cell-FWS), Cell과 Flow형이 결합된 Hybrid-FWS, 수직흐름형 인공습지(vertical flow constructed wetland)와 수평지하흐름형 인공습지(vertical flow constructed wetland)가 결합된 HYBIRD 형 습지로 구분된다. 연구결과, 일반적으로 SA/CA 비율이 클수록 오염물질의 저감효율은 증가하는 것으로 나타났다. 오염 물질별 인공습지 규모를 비교할 경우 저감효율 60%에서 인공습지의 규모는 유기물>영양염류>입자상물질 순으로 나타났다. 목표 제거효율 60%에서 SA/CA 비는 BOD에서 약 3.2%, COD에서 2.5%, SS에서 1.9%, TN 2.5%, TP 2.3%로 나타났다. 입자상물질인 SS는 유기물 및 영양염류에 비하여 유역면적 대비 시설면적이 가장 적게 나타났으며, 유기물질 제거에 큰 시설규모가 필요한 것으로 나타났다. 따라서 인공습지 설계시 유역 토지이용 및 강우특성을 고려하여 적정한 수질과 유량모니터링이 필요하며, 이를 토대로 목표 오염물질 선정이 중요한 것으로 나타났다. 또한, 농업지역의 최적화된 인공습지 위치는 임야가 20~30%, 밭이 20% 이하, 논이 10~50%를 포함하는 곳이 적정한 것으로 평가되었다. 도시지역 인공습지는 도시면적이 증가할수록 효율이 크게 변하지 않기에 가용위치가 적정한 위치로 평가된다. 인공습지의 효율은 유역의 세부 토지이용에 크게 의존하는 것으로 평가되었다. 따라서 인공습지 설계시 농업지역에서는 임야, 밭 및 논의 적정면적을 고려하여 인공습지 위치가 결정되어야 하는 것으로 나타났다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.643-651
• /
• 2008

In the present paper, turbulent open-channel flows over longitudinal bedforms are numerically simulated. The Reynolds- averaged Navier-Stokes equations in curvilinear coordinates are solved with the non-linear $k-{\varepsilon}$ model by Speziale( 1987). First, the developed model is applied to rectangular open channel flows for purposes of model validation and parameter sensitivity studies. It is found that the parameters $C_D$ and $C_E$ are important to the intensity of secondary currents and the level of turbulent anisotropy, respectively. It is found that the non-linear $k-{\varepsilon}$ model can hardly reproduce the turbulence anisotropy near the free surface. However, the overall pattern of the secondary currents by the present model is seen to coincide with measured data. Then, numerical simulations of turbulent flows over longitudinal bedforms are performed, and the simulated results are compared with the experimental data in the literature. The simulated secondary currents clearly show upflows and downflows over the ridges and troughs, respectively. The numerical results of secondary currents, streamwise mean velocity, and turbulence structures compare favorably with the measured data. However, it is observed that the secondary currents towards the troughs were significantly weak compared with the measured data.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.743-752
• /
• 2019

The flow in the meandering channel is characterized by the spiral motion of secondary currents that typically cause the erosion along the outer bank. Hydraulic structures, such as spur dike and groyne, are commonly installed on the channel bottom near the outer bank to mitigate the strength of secondary currents. This study is to investigate the effects of submerged vanes installed in a $90^{\circ}$ meandering channel on the development of secondary currents through three-dimensional numerical modeling using the hybrid RANS/LES method for turbulence and the volume of fluid method, based on OpenFOAM open source toolbox, for capturing the free surface at the Froude number of 0.43. We employ the second-order-accurate finite volume methods in the space and time for the numerical modeling and compare numerical results with experimental measurements for evaluating the numerical predictions. Numerical results show that the present simulations well reproduce the experimental measurements, in terms of the time-averaged streamwise velocity and secondary velocity vector fields in the bend with submerged vanes. The computed flow fields reveal that the streamwise velocity near the bed along the outer bank at the end section of bend dramatically decrease by one third of mean velocity after the installation of vanes, which support that submerged vanes mitigate the strength of primary secondary flow and are helpful for the channel stability along the outer bank. The flow between the top of vanes and the free surface accelerates and the maximum velocity of free surface flow near the flow impingement along the outer bank increases about 20% due to the installation of submerged vanes. Numerical solutions show the formations of the horseshoe vortices at the front of vanes and the lee wakes behind the vanes, which are responsible for strong local scour around vanes. Additional study on the shapes and arrangement of vanes is required for mitigate the local scour.