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

본 연구에서는 고정하상 직하류부에서 발생하는 국부세굴현상을 분석하기 위해서 이동상 수리실험을 수행하였다. 순환식 유량공급이 가능한 개수로 실험장치에 아크릴로 제작된 사각형 모형을 설치하고 하류에 균일한 입경의 모래를 설치하여 다양한 수리조건과 시간변화에 따른 하상고의 변화를 측정하였다. 그리고 측정된 하상고 자료로부터 시간별 최대세굴심의 크기를 산정하여 분석하였다. 그 결과 세굴심의 상류부 경사는 일정한 값으로 유지되었고, 최대세굴심과 세굴공 길이, 최대세굴심 발생거리는 시간변화에 따라 증가폭이 점차적으로 감소하였다. 또한 세굴심의 시간변화식의 기울기와 경계지점에서의 수심평균 상대난류강도값과의 관계를 분석하였다. 세굴공 내에서는 하상근처에서 후류가 발생하고, 수심방향으로 유속편차가 극심하게 증가하는 분리전단층이 발생하였다. 본 연구에서 분석된 최대세굴심의 시간적 변화양상과 세굴심 내부에서 발생하는 흐름 및 난류강도 분포특성을 보다 면밀히 분석하여 세굴심의 발달에 영향을 미치는 인자를 이용하여 향후 적절한 세굴심의 예측식을 개발할 수 있고, 최대세굴심을 저감할 수 있는 방안을 제안할 수 있다.

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

충남서부권은 기상이변으로 강수량이 적고 가뭄이 심화되고 있는 반면, 신규수원 개발 적지 부족 등 용수공급에 어려움을 겪고 있다. 2015년부터 2016년까지 가뭄시 비상용수공급을 위해 비상도수로를 건설하였다. 이 비상도수로는 보령댐의 가뭄단계 상황기준으로 운영중에 있다. 가뭄단계 상황이 경계단계에 진입하면 도수로를 가동시작하여 관심단계 회복시에 중단하게 된다. 동 연구는 현재 운영중인 보령댐도수로의 가동Rule을 수량(이수안전도), 수질(BOD), 전력비용을 고려하여 시나리오별로(상시도수, 총 도수량대비 25%, 50%, 75%, 100%) 도수할 때 어느 도수 Rule이 수량-수질-에너지 넥서스 측면에서 가장 유리한지를 밝혀 내고자 하였다. 연구에 사용된 기초자료는 보령댐도수로 준공이후 2016년~2019년 도수실적 자료를 이용하였다. 수량(이수안전도) 검토는 MODSIM을 활용하였고, 수질개선효과는 실측 BOD를 기준으로 도수조건별로 오염부하량으로 검토하였다. 전력비용은 금강원수 도수년도별로 산정하여 분석하였다. 추가 연구결과를 통해 가뭄단계 상황기준, 월별 물관리 여건을 고려하여 기준향후 비상도수로 운영에 있어서 수량절량, 수질개선, 전력비용 절감을 기대할 수 있을 것이다.

• Journal of Korea Water Resources Association
• /
• v.47 no.3
• /
• pp.297-306
• /
• 2014

In this study, EFDC (Environmental Fluid Dynamics Code) model was used to simulate the salinity distribution for sea water inflow and rainfall runoff. The flowrate was given to the boundary conditions, which can be calculated by areal-specific flowrate method from the measured flowrate of the representative outfall. The boundary condition of the water elevation can be obtained from the hourly tidal elevation. The flowrate from the outfall can be calculated using the condition of the 245 mm raifall. The simulation results showed that at Sites 1~2 and the Mangsan island (Site 4) the salinity becomes 0 ppt after the rainfall. However, the salinity is 30 ppt when there is no rainfall. Time series of the salinity changes were compared with the measured data from January 1 to December 31, 2010 at the four sites (Site 2~5) of Yongwon channel. Lower salinities are shown at the inner sites of Yongwon channel (Site 1~4) and the sites of Songjeong river (Site 7~8). The intensive investigation near the Mangsan island showed that the changes of salinity were 21.9~28.8 ppt after the rainfall of 17 mm and those of the salinity were 2.33~8.05 ppt after the cumulative rainfall of 160.5 mm. This means that the sea water circulation is blocked in Yongwon channel, and the salinity becomes lower rapidly after the heavy rain.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.23 no.1
• /
• pp.32-48
• /
• 2018

The mechanisms of $CH_4$ uptake into and release from the ocean are not well understood due mainly to complexity of the biogeochemical cycle and to lack of regional-scale and/or process-scale observations in the marine boundary layers. Without complete understanding of oceanic mechanisms to control the carbon balance and cycles on a various spatial and temporal scales, however, it is difficult to predict future perturbation of oceanic carbon levels and its influence on the global and regional climates. High frequency, high precision continuous measurements for carbon isotopic compositions from dissolved $CH_4$ in the surface ocean and marine atmosphere can provide additional information about the flux pathways and production/consumption processes occurring in the boundary of two large reservoirs. This paper introduces recent advances on optical instruments for real time $CH_4$ isotope analysis to diagnose potential applications for in situ, continuous measurements of carbon isotopic composition of dissolved $CH_4$. Commercially available, three laser absorption spectrometers - quantum cascade laser spectroscopy (QCLAS), off-axis integrated cavity output spectrometer (OA-ICOS), and cavity ring-down spectrometer (CRDS) are discussed in comparison with the conventional isotope ratio mass spectrometry (IRMS). Details of functioning and performance of a CRDS isotope instrument for atmospheric ${\delta}^{13}C-CH_4$ are also given, showing its capability to detect localized methane emission sources.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.6
• /
• pp.587-601
• /
• 2020

The ocean circulation was simulated in the East Sea and Ulleungdo-Dokdo region using ROMS (Regional Ocean Modeling System) model. By adopting the East Sea 3 km model and the HYCOM 9 km data, Ulleungdo 1 km model and Ulleungdo-Dokdo 300 m model were constructed with one-way grid nesting method. During the model development, a correction method was proposed for the distortion of the open boundary data which may be caused by the bathymetry data difference between the mother and child models and the interpolation/extrapolation method. Using this model, a super-high resolution ocean circulation with a horizontal resolution of 300 m near the Ulleungdo and Dokdo region was simulated for year 2018. In spite of applying the same conditions except for the initial and boundary data, the numerical models result indicated significantly different characteristics in the study area. Therefore, these results were compared and verified by using the surface current data estimated by satellites altimeter data and temperature data from NIFS (National Institute of Fisheries Science). They suggest that in general, the improvement of the one-way grid nesting with the HYCOM data on RMSE, Mean Bias, Pattern correlation and Vector correlation is greater in 300 m model than in the 1 km model. However, the nesting results of using East Sea 3 km model showed that simulations of the 1 km model were better than 300 m model. The models better resolved distinct ridge/trough structures of isotherms in the vertical sections of water temperature when using the higher horizontal resolution. Furthermore, Karman vortex street was simulated in Ulleungdo-Dokdo 300 m model due to the terrain effect of th islands that was not shown in the Ulleungdo 1 km model.

• Journal of the Korea Concrete Institute
• /
• v.18 no.1 s.91
• /
• pp.109-116
• /
• 2006

An incremental approach to predict the time dependent flexural behavior of composite girder is presented in the framework of incremental finite element method. Age dependent nature of creep, shrinkage, and maturing of elastic modulus of concrete is prescribed in the incremental tangent description of constitutive relation derived based on the first order Taylor series expansion applying to the total from of stress-strain relation. The loop phenomenon in which age dependent nature of concrete causes stress redistribution and it causes creep in turn is taken into account in the formulation through the incremental representation of constitutive relation. The developed algorithm predicts the time dependent deflections of 4.8m long two span double composite box girder subjected to shrinkage, maturing of elastic modulus, and creep initially induced by self weight. Comparison shows a good agreement between the predicted and measured results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2B
• /
• pp.179-189
• /
• 2010

This paper presents a numerical investigation of vorticity generation in fully vegetated open-channel flows. The Reynolds stress model is used for the turbulence closure. Open-channel flows with rough bed-smooth sidewalls and smooth bed-rough sidewalls are simulated. The computed vectors show that in channel flows with rough bed and rough sidewalls, the free-surface secondary currents become relatively smaller and larger, respectively, compared with that of plain channel flows. Also, open-channel flows over vegetation are simulated. The computed bottom vortex occupies the entire water depth, while the free-surface vortex is reduced. The contours of turbulent anisotropy and Reynolds stress are presented with different density of vegetation. The budget analysis of vorticity equation is carried out to investigate the generation mechanism of secondary currents. The results of the budget analysis show that in plain open-channel flow, the production by anisotropy is important in the vicinity of the wall and free-surface boundaries, and the production by Reynolds stress is important in the region away from the boundaries. However, this rule is not effective in vegetated channel flows. Also, in plain channel flows, the vorticity is generated mainly in the vicinity of the free-surface and the bottom, while in vegetated channel flows, the regions of the bottom and vegetation height are important to generate the vorticity.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.192-205
• /
• 2009

The unexpected wind over the Mt. Hwawang on 9 February 2009 was deadly when many spectators were watching a traditional event to burn dried grasses and the fire went out of control due to the wind. We analyzed the fatal wind based on wind flow simulations over a digitized complex terrain of the mountain with a localized heating area using a three dimensional computational fluid dynamics model, CFD_NIMR_SNU (Computational Fluid Dynamics_National Institute of Meteorological Research_Seoul National University). Three levels of fire intensity were simulated: no fire, $300^{\circ}C$ and $600^{\circ}C$ of surface temperature at the site on fire. The surface heat accelerated vertical wind speed by as much as $0.7\;m\;s^{-1}$ (for $300^{\circ}C$) and $1.1\;m\;s^{-1}$ (for $600^{\circ}C$) at the center of the fire. Turbulent kinetic energy was increased by the heat itself and by the increased mechanical force, which in turn was generated by the thermal convection. The heating together with the complex terrain and strong boundary wind induced the unexpected high wind conditions with turbulence at the mountain. The CFD_NIMR_SNU model provided valuable analysis data to understand the consequences of the fatal mountain fire. It is suggested that the place of fire was calm at the time of the fire setting due to the elevated terrain of the windward side. The suppression of wind was easily reversed when there was fire, which caused updraft of hot air by the fire and the strong boundary wind. The strong boundary wind in conjunction with the fire event caused the strong turbulence, resulting in many fire casualties. The model can be utilized in turbulence forecasting over a small area due to surface fire in conjunction with a mesoscale weather model to help fire prevention at the field.

• Journal of Environmental Science International
• /
• v.6 no.4
• /
• pp.323-333
• /
• 1997

An objective method for the generation of velocity streamfunction is presented for dealing with discretely sampled oceauc data. The method treats a Poisson equation (forced by vorticity) derived from Helmholtz theorem In which streamfunction is obtained by isolating the non-divergent part of the two-dimensional flow field. With a mixed boundary condition and vorticity field estimated from observed field, the method Is Implemented over the Texas-Louisiana show based on the current meter data of the Texas-Louisiana Shelf Circulation and Transport Processes Study (LATEX) measured at 31 moorings for 32 months (April 1992 - November 1994). The resulting streamfunction pattern is quote consistent with observations. The streamfunction field by this method presents an opportunity to initiauze and to verier computer models for local forecasts of enoronmental flow conditions for ell spill, nutrient and plankton transports as well as opportuuty to understand shelf-wlde low-frequency currents.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.2
• /
• pp.287-298
• /
• 1989

A numerical work was performed to study the flow behaviors of the open channel type flow with its geometric boundary conditions being similar to that of the Multi-Stage-Flash evaporator with and without a baffle. For the analysis, two-dimensional steady turbulent flow was assumed and the widely known k-.epsilon. turbulence model was usded. SIMPLE algorithm and the power difference scheme were used for the numerical approach. Numerical results generally agree with the previous experimental results though there are some uncertainties at far downstream and near the free surface due to the three dimensionality of the flow and surface waves. Without a baffle, the flow has basically the shape of the submerged plane wall jet with its upper boundary at downstream being sharply curved toward the free surface. For the flow with a baffle, recirculation flow patterns are observed at the upper inlet portion and at the backside of the baffle. For the case without a baffle, it was also confirmed that the ratio between the liquid level and the gate opening height is the most important parameter to determine the flow behavior.