• Proceedings of the Korea Water Resources Association Conference
• /
• 1998.05a
• /
• pp.395-403
• /
• 1998

오염원의 이송확산에 관한 많은 연구들이 수행되어 왔으나 특히 비포화 영역에서 오염원 이송확산을 측정하는 것은 매우 어려운 것으로 알려지고 있다. 비포화 토양에서의 오염원 이송확산은 매질의 함수량 변화에 영향을 받기 때문에 오염원 거동특성을 이해하려면 비포화 흐름 분석을 선행한 후 오염원의 이송확산 특성을 분석하여야 한다. 본 연구에서는 비포화 영역에서의 오염원 이송특성을 분석하기 위하여 TDR(Time Domain Reflectometry)을 이용하여 비포화 흐름 및 오염원 이송을 측정하였다. 이를 위하여 본 연구에서는 TDR을 이용하여 오염원 이송을 측정하는 방법을 개발하였으며, 이 방법을 이용하여 1차원의 토양기둥시료에서 비포화 흐름 및 오염원 이송확산에 관한 실험을 수행하고 수치모형을 적용함으로써 비포화 영역에서 오염원의 이송확산에 관한 거동특성을 규명하였다. 본 연구에서는 두 종류의 국내 토양시료(SUS, KUS)를 사용하였는데, 토양의 물리적 특성을 예비실험을 통하여 규명한 후 토양기둥시료를 이용한 본실험을 수행하였다. 비포화 천이흐름하의 오염원 이송확산 실험에서는 급격한 습윤전선의 전진에 따른 종형의 함수량변화를 관측할 수 있었고, 이때 오염원의 농도는 함수량의 천이구간의 중심점으로부터 전방영역의 농도분포가 습윤전선에서의 함수량 분포와 유사한 종형을 이루고 있음을 관측할 수 있었다. 비포화 정상흐름하의 오염원 이송확산 실험에서는 오염원이 이송하며 농도 천이구간이 확장되어지는 전형적인 형태를 보였다. 또한 예비실험에서 측정한 매개변수를 입력자료로하여 수행한 수치결과와 실험결과를 비교하였는데 비포화 흐름특성은 실험결과와 수치결과가 정량적으로 일치하는 경향을 보였으나, 오염원 이송확산 특성은 정량적으로 수치결과가 실험결과보다 더 많이 확산되는 경향을 보였다. 따라서 수치모형을 현장에 적용할 경우 확산지수 결정에 주의하여야 할 것으로 판단된다. 즉, 수치모형에 적용할 확산지수는 BTC 실험을 통하여 측정한 확산지수, 수치확산, 흡착계수, 적용영역의 크기 등을 고려하여 결정하여야 한다. 특히 본 논문에서는 TDR을 이용하여 최초로 천이상태의 함수량과 오염원 농도를 측정하였는데 이를 위하여 전기전도도와 함수량관계를 추정하는 식을 제안하였으며, 전기전도도와 토양수 농도, 전기전도도와 함수량의 관계를 이용한 천이상태의 오염원 농도 측정방법을 개발하였다. 특히 제안식에서는 한계함수량의 개념을 도입하여 전기전도도와 함수량관계를 추정하므로 추정식의 실험값 반영 정도를 증가시켰다. 본 연구에서 제안된 식을 이용하여 추정된 전기전도도와 함수량관계는 다른 제안식에 비하여 개선된 결과를 보여 주었고, 본 연구에서 개발한 오염원 농도 측정법을 이용하여 측정한 결과 함수량이 0.15이하에서는 측정오차가 크지만 함수량이 0.15이상일 경우 매우 좋은 결과를 보였는데 질량평형을 검토한 결과 약 5-10%의 오차율을 보였다. 따라서 본 논문에서 개발된 천이상태의 오염원 농도측정법은 용존 오염물질의 이송에 관한 정확한 실험을 제공할 것으로 판단된다.

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

In this study, a 1-D laboratory experiment was conducted to investigate the characteristics of transient unsaturated solute transport by using two kinds of soils of which properties were known by test. Especially the TDR method which is proposed in this study was used to measure water content and solute concentration. As results, in the transient flow, the wetting front moves down rapidly, and the distribution of solute concentration near the wetting front showed the similar type of the water content distribution(semi-bell type). A numerical model HYDRUS was used to compare with the experimental results. Numerical results for the water movement are similar to experimental result. However, numerical results of the distribution of solute concentration are more scattered than experimental results. It means that measured dispersivity, numerical dispersion, adsorption coefficient, and soil sample size etc. should be considered in order to determine the dispersivity used in the numerical model. The present measuring method was proved to be superior to other formula and to be an available method to apply to solute transport test. The measuring error of the developed method is estimated smaller than 10% while water content is larger than 0.15.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.1
• /
• pp.37-44
• /
• 1998

In this study a new hybrid method is developed for solving flow-dominated transport problems accurately and effectively. The method takes the forward-tracking particle method for advection. However, differently from the random-walk Lagrangian approach it solves the diffusion process on the fixed Eulerian grids. Therefore, neither any interpolating algorithm nor a large enough number of particles is required. The method was successfully examined for both cases of instantaneous and continuous sources released at a point. Comparison with a surrounding 5-point Hermite polynomial method (Eulerian-Lagrangian method) and the random-walk pure Lagrangian method shows that the present method is superior in result accuracy and time-saving ability.

• Journal of Korea Water Resources Association
• /
• v.30 no.2
• /
• pp.143-154
• /
• 1997

A numerical model for solving advection-diffusion equation is presented by splitoperator method combining the Holly-Preissmann scheme with a fifth-degree interpolating polynomial for advection operator and the explicit scheme porposed by Hobson et al. for diffusion operator. To examine the developed model, the obtained numerical solutions are compared with both the analytic solution and those from the existing models for the instantaneous source (Gaussian hill) and the continuous source (advanced front) at upstream boundary with constant velocity and diffusivity condition. For the various cases having different Courant and Peclet numbers, it is shown that the present study provides stable solutions even for Courant numbers exceeding one. The result obtained by the present study also agree well with existing analytical solutions for both cases. The proposed explicit scheme somewhat releases the conventional restriction of explicit schemes for determining the time step size and provides satisfactory results for relatively large time step size.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.39 no.6
• /
• pp.751-758
• /
• 2019

This study tested various assumptions that simplified the configuration of the numerical model for unsaturated zone's contaminant transport to simulate the pathway to exposed point. This study investigated the contaminant migration through in the pollutant exposure pathway of vadoze zone for risk assessment of the contaminated site. For the purpose, generic scenarios as well as contaminant-based scenarios were simulated using the numerical code for transport of the contaminant in the pathway. The finite-difference one-dimensional transport with adsorption and biodegradation were considered, and it also assumed that the initial concentration was also depleted over time. The results of the generic-scenario show that as the groundwater infiltration rate decreases, the longer the path from the source to the groundwater level, the lower the concentration at the point of inflow into the groundwater level. In particular, in the case of high biodegradation rate and rapid depletion of pollutant sources, statistically outliers were found in the simulated results and generic scenarios was good at prediction.

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

This study is to develop a method of measuring the soil water concentration by using TDR, which is based on the relationships between the bulk soil electrical conductivity of soil and the reflected wave of TDR. The proposed monitoring method is combined with two important relationships. One is that between the bulk soil electrical conductivity and the solute concentration, which is known to be linear at a constant volumetric soil water content and the other is that between the relative bulk soil electrical conductivity and the water content at a constant concentration. Some formulas have been proposed to solve the second relationship, but a new formula and the critical water content are proposed to improve the accuracy of measurement. This proposed formula estimates the relative bulk soil electrical conductivity for water contents which is divided to two regions, linear and nonlinear, by the critical water content. As the result of the comparison with other formulas, the proposed formula is proved to be superior to other formulas and to be an available method to apply to the unsaturated transient solute transport.

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

본 연구에서는 영산강의 수질이 가장 심각한 광주하수처리장 방류구 직하류부에 대한 오염물질 확산에 시.공간적 분포양상을 파악하기 위하여 대상하천 영산강 본류의 25km구간에 RMA4 모델을 적용하여 오염 물질의 농도변화를 예측하였다. RMA4 모델의 적용을 위하여 RMA2 모델을 수행하여 구한 출력자료를 입력자료 사용하였으며, 수치모델의 적용을 위하여 연구 대상하천의 25km구간에 대하여 유한요소망은 2,304개의 절점과 7,635개의 요소로 구성하였다. 광주하수처리장 방류구 직하류부에 대한 오염물질 확산에 시.공간적 분포양상을 파악하기 위하여 광주하수종말처리장의 방류수에 의한 오염물질유입에 대한 4개의 시나리오를 작성하여 RMA4 모델을 수행하였으며, 광주하수종말처리장의 평균방류량 6.13cms와 BOD 수질기준 110mg/L(2004년 하수종말처리장 유입원수의 평균농도)의 농도를 유입 오염물질로 가정하여 영산강의 저수기의 유출량을 대상으로 오염물질 유입시간에 따라 4개의 시나리오에 대한 오염물질의 확산의 시.공간적 분포양상, 유하길이에 대한 오염물질의 최대 오염원농도 도달시간을 제시하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.1
• /
• pp.131-141
• /
• 1994

Three characteristics-based split-operator methods were applied to a longitudinal pollutant dispersion problem, and the results were compared with those of several Eulerian schemes. The split-operator methods consisted of generalized upwind, two-point fourth-order and sixth-order Holly-Preissmann schemes, respectively, for the advection calculation, and the Crank-Nicholson scheme for the diffusion calculation. Compared with the Eulerian schemes tested, split-operator methods using the Holly-Preissmann schemes gave much more accurate computational results. Eulerian schemes using centered difference approximations for the advection term resulted in numerical oscillations, and those using backward difference resulted in numerical diffusion, both of which were more severe for smaller value of the longitudinal dispersion coefficient.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.191-199
• /
• 1993

Diffusion of a steady horizontal line source in a turbulent shear flow is simulated by numerically solving a steady two-dimensional advective diffusion equation. The computational result is compared with the analytic solution for uniform velocity and diffusivity distributions over the depth. The analytic solution for constant velocity and diffusivity overestimates the degree of vertical mixing. The normalized equation indicates that friction factor is the only physical parameter that governs the vertical diffusion process. Sensitivities of the diffusion process to the friction factor and initial source position are analyzed. The rate of vertical mixing varies approximately as the square root of the friction factor. The optimal source position, which gives the most rapid mixing, lies above the mid-depth and moves toward the water surface as the friction factor increases.

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