• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.575-583
• /
• 2017

Various types of flow patterns around the stepped weir and spillway, such as the skimming flow over such structures and the wave-type flow with a standing undular hydraulic jump and roller downstream of the structures, are developed in open channels. Unsteady three-dimensional numerical simulations are carried out using a hybrid RANS-LES turbulence modeling approach and the volume of fluid method for resolving free surface fluctuations to represent the turbulent flow including the skimming flow and wave-type flow over a stepped weir installed in a rectangular channel. The comparison of numerical results with an existing experimental measurement reveals that the present numerical simulations reasonably well reproduce the turbulent flow passing the stepped weir, in terms of time-averaged velocity profiles at selected locations downstream of the weir, flow topology characterized by the wave-type and skimming flows, the maximum height and length of the standing wave and the length of reattachment of recirculating zone. The numerical result further elucidates the distinct flow behaviors of the wave-type and skimming flow by presenting instantaneous intense variations of free surface and velocity vectors, the distributions of Reynolds shear stress and turbulent kinetic energy and three-dimensional complex features of coherent structures and total pressure distribution.

• Proceedings of the Korean Society of Agricultural Engineers Conference
• /
• 2004.11a
• /
• pp.103-103
• /
• 2004

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

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

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

Various types of flow conditions are developed in the region just downstream of hydraulic structures such as weir and drop structures. One of distinct flow conditions occurred downstream of drop structures is the wave type flow with undular hydraulic jump formation. We present three-dimensional numerical simulations of a wave type flow formed downstream of a stepped weir which were experimentally investigated by Kang et al. (2010). The turbulent flow over the weir structure is modeling using the unsteady Reynolds-averaged Navier-Stokes (URANS) simulation employing the Spalart-Allmaras one equation model and the detached eddy simulation. Numerical modeling and the performance of turbulence modeling approaches are evaluated by comparing with the experimental measurements in terms of the free surface variation, the shapes and sizes of undular wave, roller near at free surface, recirculation zone near the channel bottom downstream of the structures, and streamwise velocity profiles at selected longitudinal locations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.6
• /
• pp.1037-1048
• /
• 2016

In this study, nonlinear wave interaction in the presence of a uniform current is studied using numerical model, named CADMAS-SURF which is based on the Navier-Stokes equations coupled with Volume of Fluid for tracking free surface deformation. The original CADMAS-SURF developed for interaction of wave with structure is modified/extended to simulate nonlinear fluid dynamic motions within wave-current coexisting field. The capability of Numerical Wave-Current Tank (NWCT) in this study is validated by comparing with available existing laboratory experiments for both wave-following and wave-opposing current. The numerical results for interaction between wave and current are shown to be in good agreement with experimental data. Then, this study focused on the dynamic motions of the water velocity, surface elevation and vorticity within combined wave-current field in demonstrating complex nonlinear physical phenomena due to interaction between wave and current. In addition, NWCT is applied to simulate a more complex wave-current-structure field for wave propagating over a submerged breakwater associated with current. Detailed discussion including characteristics of velocity and vorticity fields and the relation between free surface and vorticity are given.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2003.08a
• /
• pp.76-86
• /
• 2003

부방파제는 파랑에너지가 집중되어 있는 자유수면에서 수면하 20%정도에 부체를 띄워 부체에 의한 반사파와 투과파 및 부체동요에 의한 발산파에 의해 파랑을 제어한다. 부방파제의 장점은 무엇보다 대수심해역에 적합하고, 해수흐름을 방해하지 않아 항내오염을 줄일 수 있을 뿐만 아니라, 지진의 피해를 최소한으로 줄일 수 있고, 연약지반상에서도 쉽게 건설될 수 있다는 점이다. (중략)

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.21-32
• /
• 1992

Numerical solutions are presented for solving the free surface flow created by a three-dimensional body moving beneath the free surface with constant velocity at an angle of attack. The solution is obtained using a panel method based on the perturbation potential, which employs Havelock sources and normal dipoles distributed on the body surface and Havelock normal dipoles in the wake downstream of the trailing edge. A pressure Kutta condition with an iterative solution procedure is implemented to satisfy equal pressure condition on the upper and lower surfaces at the trailing edge. Numerical calculation examples in the present paper include an ellipsoid at zero angle of attack, a rectangular planform wing at a small angle of attack in the limit of zero Froude number and then free surface flows and hydrodynamic forces acting on the submerged spheroid and parabolic strut are calculated. Discussions are made about the validity of the present method.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.2
• /
• pp.95-106
• /
• 1999

A numerical model is represented to calculate the wave fields such as the reflected waves, the transmitted waves and the depth-averaged velocities over submerged breakwaters for the normally incident wave trains of nonlinear mono-chromatic wave and solitary wave. The finite amplitude shallow water equations with the effects of bottom friction are solved numerically in time domain using an explicit dissipative Lax-Wendroff finite difference method. The numerical model is verified by comparisons with the other numerical results and the measured data. It is found that the submerged breakwater may be more useful for protecting the energies of monochromatic waves rather than solitary waves. Finally, the armor stability on submerged breakwater is indirectly analyzed using the hydrodynamic characteristics of flow fields.

• Journal of Korea Water Resources Association
• /
• v.48 no.4
• /
• pp.281-290
• /
• 2015

Numerical simulations of free surface flow over a broad-crested rectangular weir are conducted by using the volume of fraction (VOF) method and three different turbulence models, the k-${\varepsilon}$, RNG k-${\omega}$ and k-${\omega}$ SST models. The governing equations are solved by a second-order accurate finite volume method and the grid sensitivity study of solutions is carried out. The numerical results are evaluated by comparing the solutions with experimental and numerical results of Kirkgoz et al. (2008) and some non-dimensionalized experimental results obtained by Moss (1972) and Zachoval et al. (2012). The results show that the present numerical model can reasonably reproduce the experimental results, while three turbulent models yield different numerical predictions of two distinct zones of flow separation, the first zone is in front of the upstream edge of the weir and the second is created immediately behind the upstream edge of the weir where the flow is separated to form the separation bubble. The standard k-${\varepsilon}$ model appears to significantly underestimate the size of both separation zones and the k-${\omega}$ SST model slightly over-estimates the first separation zone in front of the weir. The RNG k-${\varepsilon}$ model predicts both separation zones in overall good agreement with the experimental measurement, while the k-${\omega}$ SST model yields the best numerical prediction of separation bubble at the upstream edge of the weir.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.40 no.2
• /
• pp.145-156
• /
• 2020

The climate change occurring all over the world increases the risk, specially in urban area, Accordingly, rainwater pumping station expansion is required than before. In order to increase the efficiency of the rainwater pumping station, the analysis of flow characteristics in the pump sump is needed for vortex control. Many efforts have been made to illuminate the vortex behavior using PIV, but any reliable results have not been obtained yet, because of the limitations in image capturing and dependency of measured velocity values on the interrogation area and time interval used for velocity calculation. In this study, therefore, experiments were carried out to find out the limitation of PIV and estimate the validation of the velocity values associated with the analysis parameters such as interrogation area, time interval, grid size. For the experimental condition used in this study, the limitation of PIV and the effects of parameters on the velocity estimation are presented.