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

여름철 강우시 저수지로 유입하는 고탁도의 하천수는 저수지 주변 수체에 비해 낮은 수온과 높은 부유물질 농도(SS)로 인해 상대적으로 높은 밀도를 가지며, 저수지 내에서 동일 밀도층을 따라 분포하며 밀도류를 형성하게 된다. 탁수는 대부분 장기간 저수지에 체류하며 수질과 수생태계에 큰 영향을 주게 된다. 따라서 저수지로 유입한 탁수의 밀도류 거동특성을 파악하는 것은 저수지 운영, 수질관리 및 수생태계 보호를 위한 관리대책 수립에 있어서 중요한 요소이다. 본 연구의 목적은 횡방향 평균 2차원 수치모델을 이용하여 초기수위별 다양한 유량규모에서 소양호로 유입하는 밀도류의 거동특성(유입경계지점 수심($h_0$), 침강점 수심($d_p$)과 거리($X_p$), 분리점 수심($d_s$), 중층류 관입두께($h_i$), 댐 축 도달시간($t_a$), 희석율(${\beta}$))을 분석하여 저수지 수질관리를 위한 의사결정지원 기초정보를 제공하는데 있다. 밀도류의 거동분석을 위해 사용된 유량조건은 그동안 소양호에서 발생한 홍수크기를 바탕으로 7개의 등급으로 나누었으며, 각각의 유량조건별 수문곡선은 2007년 수문사상 중 7월 30일부터 8월 30일까지의 수문사상을 바탕으로 수정가우시안 공식을 사용하여 산정하였다. 탁수 거동 특성의 모의는 서로 다른 초기 수위 및 유량조건하에서의 탁수거동을 나타냈고, SS농도 25 mg/L 이상을 기준으로 하여 탁수층의 $d_p$, $X_p$, $d_s$, $h_i$, $t_a$, ${\beta}$를 산정하였다. 일반적으로 유량규모가 커질수록 $d_p$와 $d_s$는 증가하였고, $X_p$는 댐 축으로부터 가까워짐이 확인되었으나, 동일 유량규모에 대해 초기 수위가 증가함에 따라 $X_p$는 대체로 댐으로부터 멀어졌다. $h_i$는 유량규모가 증가함에 따라 완만하게 증가하는 경향을 보였고, $t_a$는 초기수위가 EL. 165 m일 때 유량이 2,000 CMS 미만인 경우 댐 앞까지 도달하지 못하는 것으로 나타났으며, 나머지 수위조건에서는 유량이 3,000 CMS 미만인 경우 댐 앞까지 도달하지 못하는 것으로 나타났다. 유량 규모에 따라 $X_p$가 0 km인 지점과 19 km인 지점에서의 ${\beta}$ 값을 산정한 결과 일반적으로 유량규모가 커질수록 유입수의 희석효과는 작아지며 초기수위가 커질수록 증가하는 경향을 나타냈다. 연구 결과는 탁수 발생 초기 저수지 운영 실무자들이 유량규모 및 초기 수위조건에 따른 밀도류의 거동특성을 신속히 예측하는 목적으로 사용될 수 있다.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.6 no.1
• /
• pp.40-48
• /
• 2001

The vertical distributions of temperature, salinity, dissolved oxygen, nutrients, chlorophyll, and primary production were investigated within the top 200m water depth in the south Pacific Ocean in February,2000. The study area ($24^{\circ}-41^{\circ}S,\;81^{\circ}-168^{\circ}W$) can be hydrologically divided into two regions. Upwelling was actively occurring in the eastern region of the $110^{\circ}S$ line, meanwhile it was not active in the western region. Accordingly, chemical properties in the surface waters were different between the two regions; nitrate+nitrite and phosphate concentrations were much higher in the eastern region than in the western region due to the active upwelling, but silicate concentration was higher in the western region. Among the nutrients, the major element influencing primary production was also different between the two regions; silicon would be a major element influencing primary production in the eastern region, but nitrogen may act as a major element for primary production in the western region. Primary production showed similar values in the two regions in spite of the large differences of nutrient concentrations in the surface waters, but the total chlorophyll integrated within the 200 m water depth was almost twice as much as in the western region than that in the eastern legion.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.6 no.1 s.20
• /
• pp.59-68
• /
• 2006

Flow separation of recirculation zone by increasing of flow and change of its direction at confluence results in backwater due to conveyance reduction. The hydraulic characteristics of flow separation are analysed by experimental results of flow ratios of tributary and main streams and approaching angles. The boundary of flow separation by dimensionless length and width is defined by the streamline of zero and this definition agrees well to the existing investigation. Because flow separation doesn't appear in small flow ratio and approaching angle of $30^{\circ}$, the equation of flow separation with flow ratio and approaching angle is provided. In flow separation consideration and comparing with previous results, the existing equations of dimensionless length and width ratios by function of approaching angle, flow ratio, and downstream Froude number are modified and also contraction coefficient and shape factor are analysed. Dimensionless length and width ratios are proportional to the flow ratio and approaching angle. In analysis of water surface profiles, the backwater effects are proportional to the flow ratio and approaching angle and the magnitude at outside wall is greater than that of inside wall of main stream. The length, $X_l$ from the beginning of confluence to downstream of uniform flow, where the depth is equal to uniform depth, is characterized by width of stream, flow ratio, approaching angle, and contraction coefficient. The ratios between maximum water depth by backwater and minimum depth at separation are analysed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2011

Finite element grid system using h-refinement on the Yellow Sea was constructed based on previous study (Suh, 1999b) from 14 K to 210 K and special attention was concentrated on refining the coastal zone. In grid generation, depth change between adjacent points and non-dimensional tidal wave length ratio were considered. As a result approximately a quarter of the total nodes are located nearby 5 m of shallow area. Accurate bathymetry data using 30's and ETOPO1 with open boundary conditions of 8 major tidal constituents extracted automatically from FES2004 have been applied. In tidal simulation a 3-dimensional nonlinear harmonic model was setup and tidal amplification due to changes in vertical turbulent and bottom friction were simulated. In this study not only 8 major tidal constituents but also nonlinear shallow tides $M_4,$, $MS_4$ and long period $M_f,$, $M_{sf}$ were reproduced. It is found that implication of spatial variation of friction coefficient plays a very important role in reproduction of astronomical and shallow tides which are computed by iterative computation of nonlinear terms. Also it should be considered differently with respect to tidal periods. To understand the distribution of tidal asymmetry, amplitude ratio of $M_4/M_2$ and phase differences $2g(M_2)-g(M_4)$ were calculated. Tidal distortion ratio marks up to 0.2 on the west coast showing shallow coastal characteristics and somewhat wide range of ebb-dominances in front of Mokpo area are reproduced.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.20-26
• /
• 2019

Where hydraulic structures are to be installed over the entire width of a river or stream, usually a bed protection structure is to be installed. However, a local scour occurs in which the river bed downstream of the river protection system is eroded due to the influence of the upstream flow characteristics. This local scour is dominant in the flow and turbulence characteristics at the boundary of the flow direction and in the material of the bed materials, and may gradually become dangerous over time. Therefore, in this study, we compared the turbulence patterns in the local scour hole at the downstream of the river bed protection with the results of the analysis of the mobile bed experiment, and compared with the application of OpenFoam, a three dimensional numerical analysis model. The distribution of depth-averaged relative turbulence intensities along the flow direction was analyzed. In addition to this result, the stabilization of scour hole was compared with the bed shear stress and Shields parameter, and the results were compared by changing the initial turbulent flow conditions. From the results, it was confirmed that the maximum depth of generation of the three-stage was dominantly developed by the magnitude of depth-averaged relative turbulence intensity rather than the mean flow velocity. This result also suggests that design, construction or gate control are needed to control the depth-averaged relative turbulence intensities in order to reduce or prevent the local scour faults that may occur in the downstream part of the bed protection.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2010.04a
• /
• pp.389-391
• /
• 2010

The Free surface influenced the wake behind a circular cylinder and its effects were investigated experimentally in a circulating water channe with the variation of water depth Instantaneous velocity fields were measured in this paper. the measured results has been compared with each other to investigate the flow characteristics of the circular cylinder's 2-dimensional section at $Re=1.0\times10^3$ using 2-frame grey level cross correlation PIV method. The flow around the circular cylinder with free surface affected the wake structure. especially, in case of d=1.0D, the boundary layer was measured in the whole area. The separation point and boundary layer of the circular cylinder could be controlled by the water depth.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.3 no.2
• /
• pp.53-58
• /
• 1998

Ekman layer equation with rough-bottom boundary condition has been solved to determine the effect of roughness on the magnitude of Ekman veering. The bottom boundary condition and the flow field were expanded in a power of roughness (h) which is always smaller than the Ekman layer thickness (${\delta}_E$). By changing the magnitude of roughness parameter (h/${\delta}_E$), the magnitude of the veering, which rotates counterclockwise from the interior geostrophic flow, has been computed. At a fixed depth within the Ekman layer, the magnitude of veering increases as the roughness parameter increases. However, the cross-isobar flux turns out to decrease with increasing roughness. To verify the analytic solution, laboratory experiments were carried out. Rough-bottom cylinderical container filled with homogeneous fluid was sit on a rotating table. The flow pattern during the period of steady spin-up shows that the degree of veering coincides well with the analytic results for various roughness parameters.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.2
• /
• pp.103-111
• /
• 2014

The wave power, extracted from a circular array of small power buoys, is investigated under the potential theory. It is assumed that the buoy's radius, the draft, and the separation distance are much smaller than the water depth, the wave length, and the radius of a circular deployment area. The boundary value problem involving the macro-scale boundary condition on the mean surface covered by buoys is solved using the eigenfunction expansion method. The capture width, which is defined as the ratio of the extracted power to the wave power per unit length of the incident wave crest, is assessed for various combinations of packing ratio, radius of a circular array, and PTO damping coefficient. It is found that the circular array deployment is more effective in the viewpoint of efficiency than the single large buoy of the same total displaced volume.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.205-214
• /
• 2011

Finite element grid refinements with different intensities having 14 K, 52 K and 211 K on the Yellow Sea (YS) have been constructed to make precise tidal simulations. In the meanwhile 57 K grid was made to the extended North Western Pacific (NWP) sea. Numerical simulation were done based on 32 parallel processors by using pADCIRC v 49.21 model. In the YS tidal simulation on YS-G52K and YS-G211K grid structure, KorBathy30s and ETOPO1 bathymetry data are used and 4 major tidal constituents are prescribed from FES2004. Computed results are in good agreement within 0.138 meter in RMS error for amplification and 14.80 degree of phase compared to observed tidal records. Similar error bounds are acquired in the extended NWP tidal simulation on NWP-G57K grid with 8 tidal constituent prescription on the open boundary.