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

하구에서의 흐름은 하천의 담수와 바다에서부터 유입되는 염수, 조석, 파랑 등으로 인해 복잡한 흐름구조와 혼합 양상을 보인다. 특히 만 내에 하천이 있을 경우 만의 해류특성은 하구에서의 혼합과 이송에 지배적인 영향을 미치며, 하천에서부터 방류되는 입자들은 만의 해류특성 따라 만에서의 체류시간과 이송이 결정된다. 잔차류 특성에 의한 순 물질 플럭스의 방향과 조석비대칭에 따른 하구에서의 퇴적 형태들이 결정되며, 이로 인해 하구에서의 퇴적물 퇴적 및 물질의 집적 위치, 하구 인근과 만에서의 환경변화에 영향을 줄 수 있다. 따라서 만 내에서의 혼합과 입자의 이송확산, 하천 담수의 영향역 등과 같은 만과 하천의 흐름 특성을 이해하는 것은 연안 및 하구의 환경 및 관리에 중요하다. 본 연구에서는 영일만과 형산강을 대상으로 계절변화에 따른 영일만 내 흐름과 형산강 하구에서의 퇴적양상에 대해 수치모의를 통해 수행하였다. 수치모델로는 천수방정식으로 준 3차원 유동해석을 하는 Delft-3D Flow와 파랑모형인 SWAN 모델을 결합하여 형산강 하구와 영일만의 유동을 해석하였다. 상류개방경계는 형산강하구 9 km, 하류개방경계는 영일만 외해 50 km로 설정하였고, 경계조건은 대상지역의 관측소 자료와 전지구 모형자료를 결합하여 구성하였다. 또한, 라그랑쥬 입자추적모델을 통해 형산강 상류에서 유입한 입자들의 영일만 내 체류 시간과 집적 위치를 평가하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.231-231
• /
• 2022

하구에서의 흐름은 하천의 담수와 바다에서부터 유입되는 염수, 조석, 파랑 등으로 인해 복잡한 흐름구조와 혼합 양상을 보인다. 특히 만 내에 하천이 있을 경우 만의 해류특성은 하구에서의 혼합과 이송에 지배적인 영향을 미치며, 하천에서부터 방류되는 입자들은 만의 해류특성 따라 만에서의 체류시간과 이송이 결정된다. 잔차류 특성에 의한 순 물질 플럭스의 방향과 조석비대칭에 따른 하구에서의 퇴적 형태들이 결정되며, 이로 인해 하구에서의 퇴적물 퇴적 및 물질의 집적 위치, 하구 인근과 만에서의 환경변화에 영향을 줄 수 있다. 따라서 만 내에서의 혼합과 입자의 이송확산, 하천 담수의 영향역 등과 같은 만과 하천의 흐름 특성을 이해하는 것은 연안 및 하구의 환경 및 관리에 중요하다. 본 연구에서는 영일만과 형산강을 대상으로 계절변화에 따른 영일만 내 흐름과 형산강 하구에서의 퇴적양상에 대해 수치모의를 통해 수행하였다. 수치모델로는 천수방정식으로 준 3차원 유동해석을 하는 Delft-3D Flow와 파랑모형인 SWAN 모델을 결합하여 형산강하구와 영일만의 유동을 해석하였다. 상류개방경계는 형산강하구 9 km, 하류개방경계는 영일만 외해 50 km로 설정하였고, 경계조건은 대상지역의 관측소 자료와 전지구 모형자료를 결합하여 구성하였다. 또한, 라그랑쥬 입자추적모델을 통해 형산강 상류에서 유입한 입자들의 영일만 내 체류시간과 집적 위치를 평가하였다.

• Journal of Korea Water Resources Association
• /
• v.34 no.1
• /
• pp.19-30
• /
• 2001

A laterally averaged 2-dimensional hydrodynamic model is newly developed. The coordinate system is first transformed to minimize the effects of irregularity of bottom and surface. The advection terms of the governing equations are then discretized by an upwind scheme. By employing an explicit scheme for longitudinal direction and an implicit scheme for vertical direction, the model is free from restriction of temporal step size caused by a relatively small grid ratio. To demonstrate the applicability of the model, calculated time histories of free surface displacements and distributions of velocity and salinity are compared with the field measurements of the Keum River Estuary before construction of the estuary dam. A reasonable agreement is observed between them.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.15-20
• /
• 2014

Accurate simulation of free-surface wave flows around a ship is very important for better hull-form design. In this paper, a computational fluid dynamics (CFD) code, termed SNUFOAM, which is based on the open source libraries, OpenFOAM, was developed to predict the wave patterns around a ship. Additional anti-diffusion source term for minimizing a numerical diffusion, which was caused by convection differencing scheme, was considered in the volume-fraction transport equation. The influence of the anti-diffusion source term was tested by applying it to free-surface wave flow around the Wigley model ship. In results, the band width of the volume fraction contours between 0.1 to 0.9 at the hull surface was narrowed by considering the anti-diffusion term.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.475-483
• /
• 2008

2D finite element model for analysis of transport of accidentally released pollutants in the flow was developed by SUPG method, and the mass balance of this model was checked though two example problems: line source and point source problem in the straight channel and unidirectional 2D flow field, respectively. All the test cases were simulated with both SUPG and conventional Galerkin method to compare the accuraccy of the numerical mass balance. Test results show that the model with SUPG can adequately conserve the released mass though simulation than the model using Galerkin method, so the developed model verified to be appropriate to solve this accidental mass release problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.6C
• /
• pp.267-274
• /
• 2012

The conventional approach to evaluate the contaminant transport in soils adopts the macro-scale implementation while the pore configuration and network is a dominant factor to determine the fate of contaminant. However, the observation of fate and transport at pore scale may not be readily approachable because of the computational expenses to solve Navier-Stokes equation. We herein present the 2D Lattice-Boltzmann method that enables to assess the local fluid velocity and density efficiently for the case of single phase and multi-components. The solute fate spatio-temperal space is explicitly determined by the advection of fluid flow. Two different types of idealized pore space provides the path of fluid. Also, solute transport, the velocity field and average concentration of solute are computed in steady state. Results show that the pore geometry such as tortuosity mainly affect the solute fate. It highlights the significance of the pore configuration and shape in granular soils and rock discontinuity in spite of the equivalent porosity.

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

하천에 유입된 오염물질의 2차원 혼합거동은 하천 주흐름에 의한 이송현상과 유속 성분의 수심평균 값에 대한 공간적 편차로부터 야기되는 분산현상으로 설명 할 수 있다. 이는 3차원 이송확산 방정식으로부터 수심 적분된 2차원 이송-분산 방정식으로 수학적 유도가 가능하며, 수심방향으로 적분하는 과정에서 발생되는 농도의 분산항은 Taylor Dispersion 개념에 기초하여 종방향 및 횡방향의 2차원 분산계수로 표현된다. Fischer(1978)는 연직방향 유속분포로부터 2차원 분산계수를 추정하는 해석해를 수학적으로 유도하였으나, 실제 하천에서 정밀한 연직방향 유속분포를 계측하는 것은 많은 비용 및 노동력을 초래한다. 따라서 선행 연구자들은 2차원 혼합모형의 분산계수를 산정하고자 실험적 방법으로써 추적자실험을 수행하였다. 추적자실험은 추적자 물질을 수체에 주입한 후 농도의 변화를 관측함으로써 추적자물질이 하천에서 이송 및 분산되는 과정을 이해하는데 유용하다. 기존의 추적자실험은 고정된 위치에서 농도를 계측하여 시계열적인 농도의 변화를 관측한 후, 오염운 동결가정을 통해 종,횡방향 분산계수의 산정이 가능하지만, 오염물질 농도의 공간적 분포를 얻기에는 한계가 있다. 본 연구에서는 기존의 추적자실험법의 한계를 극복하고자 형광물질을 이용한 추적자실험을 수행함과 동시에 드론에 장착된 디지털카메라를 이용하여 항공영상을 취득 및 분석하여, 하천에 주입된 형광물질의 농도분포를 시공간적으로 추출하는 기법을 개발하고, 이를 바탕으로 오염물질의 2차원 혼합거동을 분석하였다. 본 실험은 한국건설기술연구원의 안동하천실험센터의 A3실험수로에서 수행되었으며, 실험수로는 평균 하폭 5 m, 평균 수심 0.44 m, 유량 $0.96m^3/s$의 실제 소규모 하천과 유사한 축척을 가지고 있다. 추적자물질은 Rhodamine WT 용액이 사용되었으며, 실험수로 내 설치된 15개의 형광광도계(YSI-600OMS)를 이용하여 농도를 측정하였다. 항공영상의 취득을 위해 이용된 드론은 DJI-Phantom 3 Professional 이며, 3840x2160의 해상도로 초당 30 frame의 동영상으로 취득되었다. 영상의 정합 및 좌표화를 위해 RTK-GPS를 이용하여 12개의 지상 기준점의 좌표를 취득한 후, 사영변환을 통해 영상좌표를 지상좌표로 변환하였다. 영상의 픽셀값을 농도장으로 변환하기 위해 각 RGB 밴드의 픽셀값을 통계적으로 분석하여 농도장으로 변환하였으며, 영상으로부터 얻은 농도장은 형광광도계에 의해 실측된 농도와 결정계수 0.9이상의 수준으로 정확도를 나타냈다.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.474-484
• /
• 2018

In this study, a new numerical modeling system was proposed to predict oil spills, which increasingly occur at sea as a result of abnormal weather conditions such as global warming. The hydrodynamic conditions such as the flow velocity needed to calculate oil dispersion were estimated using a three dimensional hydrodynamic model based on the Navier-Stokes equation, which considered all of the physical variations in the vertical direction. This improved the accuracy compared to those estimated by the conventional shallow water equation. The advection-diffusion model for the spilled oil was combined with the hydrodynamic model to predict the movement and fate of the oil. The effects of absorption, weathering, and wind were also considered in the calculation process. The combined model developed in this study was then applied to various test cases to identify the characteristics of oil dispersion over time. It is expected that the developed model will help to establish initial response and disaster prevention plans in the event of a nearshore oil spill.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Water for future
• /
• v.26 no.3
• /
• pp.137-148
• /
• 1993

A hybrid finite difference method for the longitudinal dispersion equation was developed. The method is based on combining the Holly-Preissmann scheme with the fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme. Longitudinal dispersion of an instantaneously-loaded pollutant source was simulated by the model and other characteristics-based numerical methods. Computational results were compared with the exact solution. The present method was free from wiggles regardless of the Courant number, and exactly reproduced the location of the peak concentration. Overall accuracy of the computation increased for smaller value of the weighting factor, $\theta$ of the model. Larger values of $\theta$ overestimated the peak concentration. Smaller Courant number gave better accuracy, in general, but the sensitivity was very low, especially when the value of $\theta$ was small. From comparisons with the hybrid method using the third-degree interpolating polynomial and with split-operator methods, the present method showed the best performance in reproducing the exact solution as the advection becomes more dominant.