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

This paper describes the development of the stream-tube based dispersion model for modeling contaminant transport in open channels. The operator-splitting technique is employed to separate the 2D contaminant transport equation into the pure advection and pure dispersion equations. Then the total variation diminishing (TVD) schemes are combined with the second-order Lax-Wendroff and third-order QUICKEST explicit finite difference schemes respectively to solve the pure advection equation in order to prevent the occurrence of numerical oscillations. Due to various limiters owning different features, the numerical tests for 1D pure advection and 2D dispersion are conducted to evaluate the performance of different TVD schemes firstly, then the TVD schemes are applied to experimental data for simulating the 2D mixing in a straight trapezoidal channel to test the model capability. Both the numerical tests and model application show that the TVD schemes are very competent for solving the advection-dominated transport problems.

• Water for future
• /
• v.23 no.1
• /
• pp.119-127
• /
• 1990

The concentration of cohesive suspended sediment is determined by the circulation of water and the material dispersion. The equations of the two-dimensional, depth-integrated dispersive transport are the Reynolds equation, continuity equation, and advection-dispersion equation based on the Fick's law. A finite difference method has been applied to two models of circulation and dispersion transport. The circulation model is solved by the explicit scheme and the dispersion transport model is solved by multi-operational scheme. It is investigated wheter advective terms are included when the equation of circulation is applied to the model. For advection-dispersion equation, it was also investigated about variations of suspended sediment concentration with respect to the critical shear stresses.

• Water for future
• /
• v.27 no.2
• /
• pp.155-166
• /
• 1994

Various Eulerian-Lagrangian numerical models for the one-dimensional longitudinal dispersion equation are studied comparatively. In the model studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing adveciton and the other dispersion. The advection equation has been solved using the method of characteristics following fluid particles along the characteristic line and the results are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpolation polynomials are superior to Lagrange interpolation polynomials in reducing dissipation and dispersion errors in the simulation.

• Journal of Environmental Science International
• /
• v.11 no.12
• /
• pp.1321-1326
• /
• 2002

A numerical model is investigated to predict a behavior of the gaseous volatile organic compounds and a subsurface contamination caused by them in the unsaturated zone. Two dimensional advective-dispersion equation caused by a density difference and two dimensional diffusion equation are computed by a finite difference method in the numerical model. A laboratory experiment is also carried out to compare the results of the numerical model. The dimensions of the experimental plume are 1.2m in length, 0.5m in height, and 0.05m in thickness. In comparing the result of 2 methods used in the numerical model with the one of the experiment respectively, the one of the advective-dispersion equation shows better than the one the diffusion equation.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.2
• /
• pp.495-506
• /
• 2013

This study developed the Time-split Mixing Model (TMM) which can represent the pollutant mixing process on a three-dimensional open channel through constructing the conceptual model based on Taylor's assumption (1954) that the shear flow dispersion is the result of combination of shear advection and diffusion by turbulence. The developed model splits the 2-D mixing process into longitudinal mixing and transverse mixing, and it represents the 2-D advection-dispersion by the repetitive calculation of concentration separation by the vertical non-uniformity of flow velocity and then vertical mixing by turbulent diffusion sequentially. The simulation results indicated that the proposed model explains the effect of concentration overlapping by boundary walls, and the simulated concentration was in good agreement with the analytical solution of the 2-D advection-dispersion equation in Taylor period (Chatwin, 1970). The proposed model could explain the correlation between hydraulic factors and the dispersion coefficient to provide the physical insight about the dispersion behavior. The longitudinal dispersion coefficient calculated by the TMM varied with the mixing time unlike the constant value suggested by Elder (1959), whereas the transverse dispersion coefficient was similar with the coefficient evaluated by experiments of Sayre and Chang (1968), Fischer et al. (1979).

• Water for future
• /
• v.28 no.4
• /
• pp.125-136
• /
• 1995

ELM, a characteristic line based method, was applied to advection-dispersion equation, and the results obtained were compared with those of Eulerian schemes(Stone-Brian and QUICKEST). The calculation methods consisted of Lagrangian interpolation scheme and cubic spline interpolation scheme for the advection calculation, and the Crank-Nicholson scheme for the dispersion calculation. The results of numerical methods were as follows: (1) for Gaussian hill: ELM, using Lagrangian interpolation scheme, gave the most accurate computational result, ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical damping for Peclet number 50. Stone-Brian scheme gave phase shift introduced in the numerical solution for Peclet number 10 and 50. (2) for advanced front: All schemes gave accurate computational results for Peclet number 1 and 4. ELM, Lagrangian interpolation scheme, and Stone,Brian scheme gave dissipation error and ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical oscillation for Peclet number 50.

• Journal of Environmental Science International
• /
• v.8 no.3
• /
• pp.355-361
• /
• 1999

The aim of this study is to develop the water quality simulation model (BAYQUAL) that deal with the physical, chemical and biological aspects of fate/behavior of pollutants in the bay. BAYQUAL is a two dimensional, time-variable finite element water quality model based on the flow simulation model in bay(BAYFLOW). The algorithm is composed of a hydrodynamic module which solves the equations of motion and continuity, a pollutnat dispersion module which solves the dispersion-advection equation. The applicability and feasibility of the model are discussed by applications of the model to the Kwangyang bay of south coastal waters of Korea. Based on the field data, the BAYQUAL model was calibrated and verified. The results were in good agreement with measured value within relative error of 14% for COD, T-N, T-P. Numerical simulations of velocity components and tide amplitude(M2) were agreed closely with the actual data.

• Geomechanics and Engineering
• /
• v.2 no.2
• /
• pp.141-156
• /
• 2010

In the framework of meshfree methods, a new methodology is developed based on radial point interpolation method (RPIM). This methodology is applied to a one-dimensional contaminant transport modelling in the saturated porous media. The one-dimensional form of advection-dispersion equation involving reactive contaminant is considered in the analysis. The Galerkin weak form of the governing equation is formulated using 1D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Numerical examples representing various phenomena, which occur during migration of contaminants, are presented to illustrate the applicability of the proposed method and the results are compared with those obtained from the analytical and finite element solutions. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. In order to test the practical applicability and performance of the RPIM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2005.11a
• /
• pp.321-331
• /
• 2005

Recently. three-dimensional models have been used for aquatic dispersion of radioactive effluents in relation to nuclear power plant siting based on the Notice No. 2003-12 'Guideline for investigating and assessing hydrological and aquatic characteristics of nuclear facility site' of the Ministry of Science and Technology (MOST) in Korea. Several nuclear power plants have been under construction or planed. which are Shin-Korl Unit 1 and 2, Shin-Wolsong Unit 1 and 2, and Shln-Ulchin Unit 1 and 2. For assessing the aquatic dispersion of radionuclides released from the above nuclear power plants, it is necessary to know the coastal currents around sites which are affected by circulation of East Sea. In this study, a three dimensional hydrodynamic model for the circulation of the East Sea of Korea has been developed as the first Phase, which Is based on the RIAMOM. The model uses the primitive equation with hydrostatic approximation, and uses Arakawa-B grid system horizontally and Z-coordinate vertically. Model domain is $126.5^{\circ}E\;to\;142.5^{\circ}E$ of east longitude and $33^{\circ}N\;and\;52^{\circ}N$ of the north latitude. The space of the horizontal grid was $1/12^{\circ}$ to longitude and latitude direction and vortical level was divided to 20. This model uses Generalized Arakawa Scheme. Slant Advection, and Mode-Splitting Method. The input data were from JODC, KNFRDI, and ECMWF. The model ing results are in fairly good agreement with schematic patterns of the surface circulation in the East Sea The local current model and aquatic dispersion model of the coastal region will be developed as the second phase. The oceanic dispersion experiments will be also tarried out by using ARGO Drifter around a nuclear pelter plant site.