• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.1741-1745
• /
• 2006

본 연구에서는 하천에 유입되는 오염물질 중 대부분을 차지하는 비보존성 오염물질의 확산거동을 예측하기 위해 2차원 수심평균된 이송분산방정식에 유한요소모형을 적용하였다. 수치모형 구성을 위해 Galerkin법을 이용한 가중잔차법을 사용하였으며, 복잡한 하천경계를 보다 정확히 재현할 수 있도록 삼각 및 사각요소망의 혼용이 가능하도록 하였다. 모의대상 수질인자는 BOD, DO, 질소화합물, 인화합물, 수온, pH 및 대장균군수이며, 이 가운데 BOD와 DO는 상호 쌍을 이루는 방정식을 풀어야 하는 특수한 형태이므로 별도로 취급하였다. 순간주입 및 연속주입에 의한 비보존성 오염물질의 확산거동을 모의하였으며, 시간에 따른 민감도를 분석하기 위해 보존성 오염물질의 확산거동과 비교하였다. 해석해를 이용해 순간주입된 오염물질의 오염운을 구하는 식을 유도하고 해석해와 본 연구에서 개발한 수치모형에 의한 수치해를 비교하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.867-874
• /
• 1994

The complex nature of low flow transport and transformation of nonconservative pollutants in natural streams has been investigated using a numerical solution of a proposed mathematical model that is based on a pair of mass balance equations describing the advection, dispersion, decay and mass exchange mechanisms in streams and in storage zones. In the present study, a mathematical model (named "Storage-Transformation Model") has been developed to predict adequately the non-Fickian nature of mixing and transformation mechanisms for decaying substances in natural streams under low flow conditions. Comparisons of the computed concentration-time curves with the measured data show that the Storage-Transformation Model yields better agreements in general shape, peak concentration and time to peak than the conventional 1-D dispersion model. The proposed model shows significant improvement over the 1-D dispersion model in predicting natural transport and transformation processes in streams through pools and riffles.

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

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

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.607-616
• /
• 1999

A fractional step finite difference model for the longitudinal dispersion of nonconservative contaminants is developed. It is based on splitting the longitudinal dispersion equation into a set of three equations each to be solved over a one-third time step. The fourth-order Holly-Preissmann scheme, an analytic solution, and the Crank-Nicholson scheme are used to solve the equations for the pure advection, the first-order decay, and the diffusion, respectively. To test the model, it is applied to simulate the longitudinal dispersion of continuous source released into a nonuniform flow field as well as the dispersion of an instantaneous source in a uniform flow field. The results are compared with the exact solution and those computed by an existing model. Compared to the existing model which uses Euler method for the first-order decay equation, the present model yield more accurate results as the decay coefficient increases.

• Journal of Korea Water Resources Association
• /
• v.34 no.4
• /
• pp.289-301
• /
• 2001

A fractional step finite difference model for the longitudinal dispersion of nonconservative pollutants is applied to the Nakdong River to simulate the phenol spill accident which occurred on March, 1971. Prior to the dispersion calculation, the flow conditions are simulated to provide inputs to the dispersion model. An unsteady flow model based on Preissmann's four-point scheme is used for this purpose. Sensitivities of the dispersion calculation to empirical equations for dispersion coefficient and to the first-order decay coefficient are analyzed. The time to peak concentration at a downstream location is significantly different depending on the formula for the dispersion coefficient. Although the decay coefficient does not affect the shape of the temporal concentration distribution, the concentration values depend on the decay coefficient verb significantly. An optimization technique is used to calibrate the dispersion model as well as the flow model. The time to the peak concentration is simulated for major positions of water intake along the Nakdong River.

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

하천으로 유입된 오염물질의 거동을 정확하게 예측하는 것은 하천 시스템의 수질 유지관리에 매우 중요하다. 본 연구에서는 1차 감쇠율(decay rate)을 가진 비보존성 오염물의 비정상 이송해석 방정식의 해를 위해 유한체적기법이 개발되어졌다. 하천 흐름 해석을 위해 자연형 단면에서의 마름-젖음 해석이 가능한 기법이 도입되었다. 이 기법은 2차-정도의 정확성와 Courant 수가 1 까지 안정함을 보장한다. 도입된 기법은 Godnov 형의 유한체적기법을 이용하여 St. Venant 방정식들을 해석하였고 질량 및 운동량 플럭스는 Roe 형의 Riemann Solver 를 사용하여 연산하였다. 오염물의 이송 해석은 추가적인 이송-확산 방정식을 도입을 통해 기존의 St. Venant 방정식과 함께 풀려질 수 있다. 추가된 방정식과 St. Venant 식은 3×3 eigenstructure를 구성하였고 이는 2차원 흐름해석 기법과 유사하게 해석될 수 있었다. 본 연구 모형의 검증을 위해 오염물의 계속적 주입을 가정한 가상 및 실제 하천에 적용되었다. 연구된 기법은 모든 적용에서 합리적 정확도를 가지고 오염물질의 연속적인 특성을 잘 모의하고 있었다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.221-230
• /
• 2009

In Korea, many water intake plants are easily affected by effluents of sewage treatment plants because sewage treatment plants are usually located upstream or nearby the plants of the same riverine area. Furthermore, the inflow of harmful contaminants owing to pollutant spills or transportation accidents of vehicles using the roads and bridges intersecting the river causes significant impact on the management of water intake plants. Paldang lake, the main water intake plants in Korea, is especially exposed to various water pollution accidents, because the drainage basin area is significantly large compared to the water surface area of the lake. Therefore it is necessary to predict the possible pollutant spill in advance and consider measurements in case of water pollution. In this study, water quality prediction was performed in Paldang Lake in Korea durig the dry season using two-dimensional numerical models. In order to represent the cases of pollutant accidents, the difference of pollutant transport patterns with varying injection points was analyzed. Numerical simulations for hydrodynamics of water flow and water quality predictions were performed using RMA-2 and RAM4 respectively. As a result of simulation, the difference of pollutant transport with the injection points was analyzed. As a countermeasure against the pollutant accident, the augmentation of the flow rate is proposed. In comparison with the present state, the rapid dilution and flushing effects on the pollutant cloud could be expected with increase of flow rate. Thus, increase of flow rate can be used for operation of water intake plants in case of pollutant spill accidents.

• Journal of Korea Water Resources Association
• /
• v.32 no.1
• /
• pp.61-70
• /
• 1999

The two-dimensional Random-Walk model in which fluid and pollutant particles are tracked using statistical concept was developed to simulate dispersion processes in natural streams. The calibration of the model shows that the error decreases as the number of grid increases, and/or the number of particles in each grid increases. The proposed model is tested against the dispersion data collected in the Grand River, Canada. The simulation results show that the 2-D Random-Walk model describes two-dimensional mixing phenomena occurred in the irregular meandering stream very accurately.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1219-1223
• /
• 2007

수온은 수계에서 가장 중요한 물리적 특성 중 하나로서, 생물 군집, 특히 어류와 무척추 동물에 관련된 많은 수질 인자에 영향을 미친다. 하천의 생태학적 모습을 개선하기 위한 하천 복원 사업 수행에 있어서, 서식처 및 산란처 조건으로서의 수온 조사 및 모델링의 필요성이 점차 증가하고 있는 상황이다. 본 연구는 낙동강의 중상류에 위치한 반변천 합류부에 평면 이차원 비정상 수치모형인 KU-RLMS 모형을 적용하여 수온의 변화 특성을 규명할 목적으로 수행하였다. KU-RLMS 모형은 하천 및 저수지의 국부적인 수리, 수질, 유사이동 해석을 위해 개발된 평면 이차원 비정상 수치모형이다. 직사각형 격자를 사용하는 유한차분법의 단점을 보완하기 위해, 수심적분된 2차원 연속방정식, 운동량방정식, 이송확산방정식을 불규칙한 경계를 현실적으로 모사할 수 있는 직교곡선 좌표계로 변환한 방정식을 사용한다. 이 모형은 흐름, 농도, 지형변화를 조합하여 계산할 수 있는 모형으로서 점착성 및 비점착성 유사의 이동, 보존성 및 비보존성 오염물질의 이동, 수온 변화를 모의할 수 있다. 수치모형 적용을 위한 현황분석으로 안동 및 임하 조정지댐의 방류량, 안동 수위관측소의 수위, 법흥교 및 포진교 지점의 수온 자료를 분석하였다. 이송확산모형의 보정을 위해, 안동대교 지점의 수온 횡분포 측정자료를 사용하여 확산계수에 대한 매개변수 추정 및 검증을 수행하였다. 또한, 안동조정지댐과 임하조정지댐의 방류량 및 방류수온을 고려하여 수치모의조건을 결정하였으며, 각 조건에 대한 수온 변화 특성을 분석하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.731-741
• /
• 2008

Previous studies on the numerical simulation at the tidal reach of Han River tend to restrict downstream boundary as Jeon-ryu station due to difficulties in gaining cross section data and tidal elevation values at Yu-do. But, in this study, geometries beyond the confluence of Gok-reung stream and Im-jin River are constructed based on the numerical sea map; tidal elevation at the downstream boundary, Yu-do is estimated by harmonic analysis of In-cheon tide gage station so that hydrodynamic and diffusion behavior have been analyzed. The domain ranging from Shin-gok submerged weir to Yu-do is selected (which is 36.8 km in length). RMA-2 and RAM4 developed by Il Won Seo (2008) are applied to simulate flow and diffusion behavior, respectively. Numerical results of flow characteristic are compared with the measured data at Jeon-ryu station. Simulation is carried out from June 23 to 25 in 2006 on the ground that hydrologic data is satisfactory and tidal difference is huge during that period. The result shows that reverse flow occurs 5 times according to the tidal elevation at Yu-do and the maximum reverse flow is observed up to Jang-hang IC, which is 32.9 km in length. Also analysis is focused on the process of generation and disappearance of reverse flow, the distribution of water surface elevation and velocity along the maximum velocity line, and the transport of nonconservative pollutant. Pollutant injected from Gul-po stream spreads widely across the river; however, the size of BOD cloud entering from Gok-reung stream is relatively small because water depth at the mid and left side becomes deeper and maximum velocity occurs along the right bank so that transverse mixing is completed quickly. Finally, mixing characteristic of horizontal salinity distribution is obtained by estimating the salinity input with analytical solution of 1D advection-dispersion equation.