• Journal of Wetlands Research
• /
• v.15 no.1
• /
• pp.59-69
• /
• 2013

The numerical simulations were performed to analyze the advection-diffusion processes of dredging-induced turbidity plume using RMA2 and RMA4 models in Bunam reservoir, Seosan, Chungnam. Field survey was also performed to measure the turbidity using the multi water quality monitoring system (YSI6600EDS). In the field survey, the vertical and horizontal distributions of the turbidity were measured during the dredging operation in Bunam reservoir. RMA2 model was used to simulate the velocity distributions in both the whole domain and the 2nd part of Bunam reservoir. RMA4 model was also used to simulate the concentration distribution in only the 2nd part of Bunam reservoir, where the dredging work were conducted. The comparison of the simulation results with the field data for the advection-diffusion of the turbidity plume using the concentration ratio concepts shows that the numerical model can be applied to analyze the environmental impact of dredging works.

• Water for future
• /
• v.28 no.5
• /
• pp.151-162
• /
• 1995

A two-dimensional stream tube dispersion model is developed to simulate accurately transverse dispersion processes of pollutants in natural channels. Two distinct features of the stream tube dispersion model derived herein are that it employs the transverse cumulative discharge as an independent variable replacing the transverse distance and that it is developed in a natural coordinate system which follows the general direction of the channel flow. In the model studied, Eulerian-Lagrangian method is used to solve the stream tube dispersion equation. The stream tube dispersion equation is decoupled into two components by the operator-splitting approach; one is governing advection and the other is governing dispersion. The advection equation has been solved using the method of characteristics and the results are interpolated onto Eulerian grid on which the dispersion equation is solved by centered difference method. In solving the advection equation, cubic spline interpolating polynomials is used. In the present study, the results of the application of this model to a natural channel are compared with a steady-state flow measurements. Simulation results are in good accordance with measured data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.749-753
• /
• 2010

본 논문에서는 하천에 유입되는 오염물질 중 대부분을 차지하는 비보존성 오염물질의 확산거동을 분석하기 위해 2차원 수심 평균된 이송분산방정식에 유한요소법을 적용하였다. 수치모형 구성을 위해 SUPG(Streamline-upwind Petrov-Galerkin)법을 이용한 가중잔차법을 사용하였다. 모의대상 수질인자는 BOD와 DO이며, BOD 농도 결과가 DO 농도 계산에서의 입력 자료로 이용되도록 상호 연계를 형성하였다. 모형의 검증을 위하여 직사각형 수로에 선원으로 연속주입하여 얻은 수치해와 해석해를 비교하였다. 비교결과 수치해와 해석해의 결과가 서로 일치하는 것을 볼 수 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.991-995
• /
• 2005

본 연구에서는 1차원 이송-분산 과정을 연구하고 전단류 흐름 및 분산거동에 있어 Taylor 이론의 핵심이라 할 수 있는 '종방향 이송과 횡방향 확산의 균형'을 기본 개념으로 하여, 이송과 확산을 분리하여 이 두 과정이 순차적으로 발생한다는 가정에 의거한 순차혼합모형을 제시하였다. 본 모형에서는 가상의 하천을 여러 개의 행과 종방향 거리를 길이가 일정한 구획으로 나누어 연속적인 분산과정을 이산적인 형태로 나타낼 수 있게 하고, 횡방향 유속분포에 따라 각 행에 각기 다른 유속을 할당한다. 오염물질은 하폭방향 선오염원으로 원점에 순간주입되며, 주어진 혼합시간 $t_m$ 동안 각 행의 오염물질들이 각자에 할당된 유속을 따라 진행하고 진행이 끝난 후 횡방향 확산이 순간적으로 이루어진다. 횡방향 확산은 횡방향으로 완전하게 일어남을 가정하여, 횡방향 확산이 끝나면 각 열에서의 농도 평균값이 할당된다. 이러한 혼합시간 $t_m$ 동안의 순차적인 이송-확산 과정이 반복되면서 오염물질의 분산이 일어나며, 농도 분포 그래프를 그릴 수 있게 된다. 순차혼합모형을 가상의 직선하천에 적용하여 종분산계수를 유도하였는데, 본 연구에서 유도된 종분산계순식은 Fischer.가 제안한 식과 유사한 형태로 나타남을 알 수 있었다. 본 모형에서 계산된 농도분포 곡선을 해석해와 비교한 결과,두 곡선이 적절히 일치함을 확인할 수 있었으며 해석해와의 비교를 통해 종분산계수 K가 혼합시간 $t_m$과 선형관계임을 확인할 수 있었다. 수심대하폭비에 따라 각기 다른 유속분포에 적용하여 종분산계수 K가 유속편차강도의 제곱에 비례관계에 있음이 밝힐 수 있었다. 수압은 $4.69kg/cm^2$으로 나타났다. 밸브 개폐도가 $100\%$일 때가 밸브를 $60\%$와 $80\%$ 개폐시켰을 때보다 $0.3kg/cm^2,\;0.29kg/cm^2$ 낮게 나타나 밸브를 전체 개방 했을 때 관로내의 수압이 상수설계기준에 적합한 수압을 유지함을 알 수 있다. 상수관로 설계 기준에서는 관로내 수압을 $1.5\~4.0kg/cm^2$으로 나타내고 있는데 $6kg/cm^2$보다 과수압을 나타내는 경우가 $100\%$로 밸브를 개방하였을 때보다 $60\%,\;80\%$ 개방하였을 때가 더 빈번히 발생하고 있으므로 대상지역의 밸브 개폐는 $100\%$ 개방하는 것이 선계기준에 적합한 것으로 나타났다. 밸브 개폐에 따른 수압 변화를 모의한 결과 밸브 개폐도를 적절히 유지하여 필요수량의 확보 및 누수방지대책에 활용할 수 있을 것으로 판단된다.8R(mm)(r^2=0.84)$로 지수적으로 증가하는 경향을 나타내었다. 유거수량은 토성별로 양토를 1.0으로 기준할 때 사양토가 0.86으로 가장 작았고, 식양토 1.09, 식토 1.15로 평가되어 침투수에 비해 토성별 차이가 크게 나타났다. 이는 토성이 세립질일 수록 유거수의 저항이 작기 때문으로 생각된다. 경사에 따라서는 경사도가 증가할수록 증가하였으며 $10\% 경사일 때를 기준으로 $Ro(mm)=Ro_{10}{\times}0.797{\times}e^{-0.021s(\%)}$로

• Ecology and Resilient Infrastructure
• /
• v.1 no.1
• /
• pp.1-7
• /
• 2014

In this study, a numerical simulation technique for coastal area where wave and current interactions are observed is proposed. Considering the spatial scale of coastal area and the coastal processes such as wave, current, shoaling, wave breaking, and inundation processes, boussinesq equation model is used. A depth-integrated transport model based on the consistent assumption with the boussinesq equation model is used for the prediction of solute transport. To solve the equations, finite volume method with an approximate riemann solver is used. The proposed model is applied to a coastal area and reasonable computational results are obtained.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.106-106
• /
• 2019

많은 도심의 하천들은 오염물질의 유입에 취약하다. 최근 신소재 공학 등 첨단산업이 발전하게 되면서 유해화학물질의 유입문제는 더욱 대두되고 있으며, 실제로 최근 유해화학물질 유입사고 발생건수가 늘어나고 있다. 특히 국내 취수량의 90%는 지표수에서 취수하고 있어, 하천오염사고는 직접적인 피해로 이어지게 된다. 따라서 이러한 사고에 대응하기 위하여 수환경에 유입된 유해물질의 거동 매커니즘을 반영한 수질해석이 필요하다. 수체 내에 유입된 유해화학물질은 기본적으로 흐름에 따른 이송 확산을 하며 흡 탈착, 휘발, 침전 부유, 생화학 반응과 같은 다양한 반응과 함께 혼합거동을 한다. 특히 소수성물질의 경우 용해된 상태뿐만 아니라, 유사에 흡착된 상태로 수체에 존재하게 된다. 결국 유해화학물질의 거동을 해석하기 위해서는 유체의 흐름 해석뿐만 아니라 수체에 존재하는 유사의 이송 또한 해석해야한다. 본 연구에서는 흐름해석을 위하여 서울대에서 개발한 흐름모형(HDM-2D)을 사용하였으며, 부유사 거동모의를 위해 부유사거동모형(STM-2D)을 개발하였다. 또한 유해화학물질의 거동모의를 위해 서울대에서 개발한 수질모형(CTM-2D)에 생성/소멸항을 추가하였으며 흐름모형과 부유사모형과의 연계를 통해 유해화학물질의 혼합거동 수치모형을 개발하였다. 각 반응항(흡 탈착, 휘발, 침전 부유, 생화학 반응)을 수치모형에 반영 시에는 보통 두 계(물-토양, 물-공기) 사이의 선형 물질교환으로 이해된다. 따라서 물질의 각 반응 별 평형농도와 물질교환속도계수를 추정식을 통해 산정하여 사용하게 된다. 하지만 각 기작이 반영유무에 따라 계산시간 및 필요입력변수가 늘어나게 되므로, 유해화학물질 유입사고와 같은 빠른 대처가 필요한 경우 각 반응 텀의 유의성을 판단하여 모형에 반영여부를 결정을 통해 경제적인 모의를 할 수 있어야 한다. 이에 따라 본 연구에서는 개발된 모형의 각 매개변수들의 민감도를 분석하고, 흐름조건 및 물질의 특성에 따른 반응항의 유의성을 판단하였다. 본 연구에서는 개발된 모형(부유사거동모형, 유해화학물질의 혼합거동모형)은 해석해 및 현장 데이터와 비교검증을 통해 개발을 완료하였으며, 각 반응항의 민감도 분석을 통해 매개변수의 임계값을 결정하였다.

• Journal of Korea Water Resources Association
• /
• v.45 no.2
• /
• pp.179-189
• /
• 2012

To simulate the transport of pollutant, a numeric model for the advection-diffusion equation with the decay term is developed. This is finite-difference model using the implicit method (with the weight factor ${\alpha}$) and Gauss-Seidel SOR(successive over-relaxation). This model is compared to the analytical solutions (of simpler dimensional or boundary conditions), and in the condition of Peclet number < 5~20, the result shows stable condition, and Crank-Nicolson method (${\alpha}$=0.5) shows the more accurate results than fully-implicit method (${\alpha}$=1). The mass of advection, diffusion and decay is calculated and the error of mass balance is less than 3%. This model can evaluate the 3-D concentrations of the advection-diffusion and decay problems, but this model uses only the finite-difference method with the fixd grid system, so it can be effectively used in the problems with small Peclet numbers like the pollutant transport in groundwater.

• 한국방재학회:학술대회논문집
• /
• 2010.02a
• /
• pp.98.1-98.1
• /
• 2010

본 연구에서는 한강 본류 팔당댐부터 잠실수중보까지의 22.5 km 구간에서 발생할 수 있는 수질 오염사고에 대한 오염물질 혼합거동 모의를 위해 RAMS를 이용하였다. 이를 위해서 2차원 흐름모형인 RAM2를 이용하여 유속장을 모의하고, 2차원 수질해석 모형인 RAM4를 이용하여 사고주입에 의한 오염물질의 시간에 따른 거동을 분석하였다. 최종적으로 잠실수중보에 위치한 주요 취수장에 미치는 영향을 분석하였다. 모의를 진행한 결과 오염물질의 이송 및 확산 거동은 투입 위치에 따라 크게 달라진다는 것을 알 수 있었다. 특히 팔당대교에서 투입된 오염물질이 풍납 취수장에 도착하는데 까지 걸리는 시간이 좌안주입과 우안주입의 경우 60시간 가까이 차이가 나는 것으로 나타났다. 따라서 오염 물질이 투입되었을 시에 일괄적으로 취수를 통제할 것이 아니라 상황에 따라 오염물질 투입 지점에 맞춰서 적절한 대책을 세워야 하는 것이 필요하다. 이 때 본 연구가 기초적인 정보를 제공할 수 있을 것이다.

• 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.

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

해양으로 유출된 5 mm 이하의 크기로 분해된 미세플라스틱이 해양 환경 오염의 주요 원인으로 자리잡았다. 최근에는 청정해역으로 알려진 남극해에서도 발견되고 있어 남극해에 잔류하는 미세플라스틱 오염 수준을 이해하기 위해 노력하고 있다. 하지만, 파랑의 효과를 고려한 남극해의 해수 순환 구조와 미세플라스틱의 고유 특성을 반영한 미세플라스틱의 거동 및 공간적 분포에 대한 복합적 이해는 상대적으로 부족하다. 남극해에서 발견된 미세플라스틱은 과학기지들의 방류수나, 조사선 등과 같은 인위적인 활동으로 인해 집적될 수 있으며, 특히 영구적으로 거주하는 과학기지에서 흘려보내는 방류수에 포함된 미세플라스틱은 과학기지 주변 해수 오염에 직접적인 영향을 줄 것으로 예상된다. 따라서, 본 연구에서는 파랑 효과에 따른 남극 킹조지 섬(King George Island)에 위치한 세종과학기지의 방류수에 포함된 미세플라스틱의 이송에 대해 모의하였다. 세종과학기지가 위치한 킹조지 섬과 넬슨 섬(Nelson Island) 사이의 멕스웰 만(Maxwell Bay)의 해수 흐름을 재현하기 위하여 해수 유동 모델(Delft3D-FLOW)이 사용되었다. 또한, 해수 유동 모델에 파랑 모델(Delft3D-WAVE)을 결합하여 파랑의 효과가 미세플라스틱의 이송에 미치는 영향을 확인하였다. 세종과학기지의 방류수가 흘러나가는 마리안 소만(Marian Cove)의 유속장을 바탕으로 이송, 확산, 입자의 침강 속도를 고려하여서, 세종과학기지에서 밀물 시 방출한 입자를 라그랑지안 입자 추적(Lagrangian Particle Tracking) 방법을 이용해 추적하였다. 해수의 밀도보다 가벼운 플라스틱의 경우 해수 표층의 흐름을 따라 소만 내부로 이송되어 해안선에 도달하고, 해수의 밀도보다 무거운 플라스틱의 경우 소만 내부로 이송되나 입자의 침강 속도로 인해 방출 위치 근처에서 집적된다. 파랑의 효과를 고려하게 되면, 고려하기 전보다 두 종류의 미세플라스틱 모두 소만 내부로 더 멀리 이송되는데, 이는 파랑으로 인한 힘(wave-induced force)이 해수 유동 모델의 운동방정식에 추가되며 파랑 에너지 분산으로 인해 해수 흐름에 변화를 준 것으로 보인다.