• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.123-136
• /
• 1996

A three-dimensional Galerkin-FEM model which can handle the temporal and spatial variation of density is presented. The hydrostatic approximation is used and density effects are included by means of conservation equation of heat and the equation of state. The finite difference grids are used in the horizontal plane and a set of linear-shape functions is used for the vertical expansion. The similarity transform is introduced to solve resultant matrix equations. The proposed model was first applied to the density-driven circulation in an idealized basin in the presence of the heat exchange between the air and the sea. The advection terms in the momentum equation were ignored, while the convection terms were retained in the heat equation. Coefficients of the vertical eddy viscosity and diffusivity were fixed to be constant. Calculation in a non-rotating idealized basin shows that the difference in heat capacity with depth gives rise to the horizontal gradient of temperature. Consequently, there is a steady new in the upper layer in the direction of increasing depth with compensatory counter flow .in the lower layer. With Coriolis force, geostrophic flow was predominant due to the balance between the pressure gradient and the Coriolis force. As a test in region of irregular topography, the model is applied to the Yellow Sea. Although the resultant flow was very complex, the character of the flow Showed to be geostrophic on the whole.

• Journal of the Korean Society of Groundwater Environment
• /
• v.3 no.2
• /
• pp.51-59
• /
• 1996

In case waste disposal site is to be constructed close to the underground facilities such as LPG storage cavern which is completely maintained by groundwater pressure, it is generally requested that the possibility on leachate contamination of cavern area be reviewed and the countermeasure, if it is estimated cavern area is severely affected by leachate, be taken into consideration. Prediction was performed and leachate control plan was made using by analytical and the numerical analysis on the leachate migration which is likely to happen at the area between the proposed waste disposal site and the underground LPG storage cavern located at the U petrochemical complex. Analytical solutions were obtained by the conservative mass advection-diffusion equation and the effect of advection and dispersion factor on the leachate migration was reviewed through peclet number calculation and the functional relationship between the factors and leachate transport velocity was established, which leads to enable us to predict the leachate transport velocity without difficulties when different parameters (factors) are used for analytical solution. Numerical solutions were obtained by FEM using AQUA2D which is for the simulation of groundwater flow and contaminant transport. 3-D discrete fracture models were simulated and fracture flow analysis was performed and feasibility study on the water-curtain system was conducted through the fracture connectivity analysis in rock mass. As results of those analyses, it was interpreted that the leachate would trespass on the LPG storage cavern area in 30 years from the proposed wate disposal site and the vertical water-curtain system was effective mathod for the prevention of leachate's migration further into the cavern area.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.6B
• /
• pp.551-560
• /
• 2009

This paper describes a newly developed pollutant transport model named ARPTM which was designed to simulate the transport and characteristics of pollutant materials after an accidental spill in upstream of river system up to a given position in the downstream. In particular, the ARPTM incorporated ADCP data to compute longitudinal dispersion coefficient and advection velocity which are necessary to apply one-dimensional advection-dispersion equation. ARPTM was built on top of the geographic information system platforms to take advantage of the technology's capabilities to track geo-referenced processes and visualize the simulated results in conjunction with associated geographic layers such as digital maps. The ARPTM computes travel distance, time, and concentration of the pollutant cloud in the given flow path from the river network, after quickly finding path between the spill of the pollutant material and any concerned points in the downstream. ARPTM is closely connected with a recently developed GIS-based Arc River database that stores inputs and outputs of ARPTM. ARPTM thereby assembles measurements, modeling, and cyberinfrastructure components to create a useful cyber-tool for determining and visualizing the dynamics of the clouds of pollutants while dispersing in space and time. ARPTM is expected to be potentially used for building warning system for the transport of pollutant materials in a large basin.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.6B
• /
• pp.507-514
• /
• 2011

A numerical model has been developed by employing a finite element method to simulate the depth-averaged 2-D dispersion of the heat pollutant, which is an important pollutant material in natural streams. Among the finite element methods, the Streamline Upwind/Petrov Galerkin (SUPG) method was applied. Also both linear and quadratic elements can be applied so that irregular river boundaries can be easily represented. To show the movement of heat pollutants, the reaction term describing heat transfer was represented as an equation in which sink/source term is proportional to the difference between the equilibrium temperature and water surface temperature. The equation was expressed so that the water surface temperature changes according to the temperature transfer coefficient and the equilibrium temperature. For the calibration of the model developed, analytic and numerical results from a case of rectangular channel with full width continuous injection have been compared in a steady state. The comparisons showed that the numerical results were in good agreement with analytical solutions. The application site was selected from the downstream of Paldang dam to Jamsil submerged weir, and overall length of this site is about 22.5 km. The change of water temperature caused by the discharge from the Guri sewage treatment plant has been simulated, and results were similar to the observed data. Overall it is concluded that the developed model can represent the water temperature changes due to heat transport accurately. But the verification using observed data will further enhance the validity of the model.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.621-628
• /
• 1997

Excavation works on stream beds have been done for various reasons including aggregate collection, sediment dredging, bridge constructions, or laying pipes under the ground. These activities may cause significant loadings of SS (suspended solids) resulting in water pollution and other detrimental effects to the surrounding environment. This research investigates application potential of a fickian diffustion model, derived from two dimensional advection-diffusion equation through some simplifying assumptions, as a planning tool for the estimation of SS loadings from excavation works and evaluation fo pollution prevention measures in case that sophisticated numerical simulation models are not applicable due to various practical reasons. Through a case study of the Juncheon stream in the Donghae City on the Kangwondo Province, this study demonstrates applicability of the fickian diffustion model as a practical method for the preliminary estimation of Ss loadings from excavation works and evaluation of performance of fabrics made of synthetic fiber for the reduction of downstream SS concentration with deficient field data.

• Journal of Korea Water Resources Association
• /
• v.33 no.2
• /
• pp.181-193
• /
• 2000

This paper presents the flow and the suspended sediment movement over ripples for oscillatory flows. A new numerical model system is developed, and applied to a laboratory experimental condition of regular waves and a fictitious condition of irregular waves. The flow field is obtained from a programme proposed by Kim et. al.(1994), which is a modified version of SOLA based on SMAC scheme. The sub-model solves the continuity and Reynolds momentum equations in the x-z plane. The wave orbital velocities, shear stresses, and pressure are all reasonably reproduced by the model. The model results on the vertical velocity component show good agreement with the measurements. The suspended sediment transport sub-model is newly set up to solve the advection-diffusion equation of suspended sediment using a split method, and involving a special shear entrainment from the whole ripple surface. The calculated suspended sediment concentrations for regular waves show reasonable agreement with measurements at Deltaflume. The model results for random waves show that the suspended sediment concentration is higher than those for regular waves and that the sediment diffuses higher than for regular waves with the significant wave height and the peak wave period of the irregular waves.

• Journal of Korean Society for Atmospheric Environment
• /
• v.22 no.1
• /
• pp.73-84
• /
• 2006

Hazardous Air Pollutants (HAPs) are characterized by being relatively heavier and denser than that of ambient air due to the various reasons such as higher molecular weight, low temperature and other complicated chemical transformations (Witlox, 1994). In an effort to investigate transport and diffusion from instantaneous emission of heavy gas, Lagrangian Particle Dispersion Model (LPDM) coupled with the RAMS output was employed. Both deposition process and buoyancy term were added on the atmospheric diffusion equations of LPDM, and the locations and concentrations of dense gas particle released from instantaneous single point source (emitting initially for 10 minutes only) were analyzed. The result overall shows that adding deposition process and buoyancy terms on the diffusion equation of LPDM has very small but detectable effect on the vertical and horizontal distribution of Lagrangian particles that especially transported for a fairly long traveling time. Also the slumping of dense gas can be found to be ignored horizontally compared to the advection by the horizontal wind suggesting that it was essential to couple the Lagrangian particle dispersion model coupled with the RAMS model in order to explain the dispersion of HAPs more accurately. However, during the initial time of instantaneous emission, buoyancy term play an important role on the vertical locations of dense particles for near surface atmosphere and around source area, indicating the importance of densities of HAPs in the beginning stage or short duration for the risk assessment of HAPs or management of heavy vapors during the explosive accidents.

• Atmosphere
• /
• v.23 no.1
• /
• pp.1-11
• /
• 2013

The angular momentum transport of an idealized axisymmetric vortex in the developing stage was investigated using the Weather Research and Forecast (WRF) model. The balanced axisymmetric vortex was constructed based on an empirical function for tangential wind, and the temperature, geopotential, and surface pressure were obtained from the balanced equation. The numerical simulation was carried out for 6 days on the f-plane with the Sea Surface Temperature (SST) set as constant. The weak vortex at initial time was intensified with time, and reached the strength of tropical cyclone in a couple of days. The Absolute Angular Momentum (AAM) was transported along with the secondary circulation of the vortex. Total AAM integrated over a cylinder of radius of 2000 km decreased with simulation time, but total kinetic energy increased rapidly. From the budget analysis, it was found that the surface friction is mainly responsible for the decrease of total AAM. Also, contribution of the surface friction to the AAM loss was about 90% while that of horizontal advection was as small as 8%. The trajectory of neutral numerical tracers following the secondary circulation was presented for the Lagrangian viewpoint of the transports of absolute angular momentum. From the analysis using the trajectory of tracers it was found that the air parcel was under the influence of the surface friction continuously until it leaves the boundary layer near the core. Then the air parcel with reduced amount of angular momentum compared to its original amount was transported from boundary layer to upper level of the vortex and contributed to form the anti-cyclone. These results suggest that the tropical cyclone loses angular momentum as it develops, which is due to the dissipation of angular momentum by the surface friction.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.1
• /
• pp.1-11
• /
• 2010

The definition of the speed of sound is reexamined since it is crucial in the numerical analysis of compressible real gas flows. The thermodynamic speed of sound (TSS), $a_{th}$, and the characteristic speed of sound (CSS), $a_{ch}$, are derived using generalized equation of state (EOS). It is found that the real gas EOS, for which pressure is not linearly dependent on density and temperature, results in slightly different TSS and CSS. in this formalism, Roe's approximate Riemann solver was derived again with corrections for real gases. The results show a little difference when the speeds of sound are applied to the Roe's scheme and Advection Upstream Splitting Method (AUSM) scheme, but a numerical instability is observed for a special case using AUSM scheme. It is considered reasonable to use of CSS for the mathematical consistency of the numerical schemes. The approach is applicable to multi-dimensional problems consistently.

• Journal of Korean Society of Water and Wastewater
• /
• v.26 no.5
• /
• pp.609-617
• /
• 2012

A dynamic water quality model is presented in order to simulate water quality under slowly varying flow conditions over time. To improve numerical accuracy, the proposed model uses a lumped system approach instead of extended period simulation, unlike the other available models. This approach can achieve computational efficiency by assuming liquid and pipe walls to be rigid, unlike the method of characteristics, which has been successfully implemented in rapidly varying flows. The discrete volume method is applied to resolve the advection and reaction terms of the transport equation for water quality constituents in pipes. Numerical applications are implemented to the pipe network examples under steady and unsteady conditions as well as hydraulic and water quality simulations. The numerical results are compared with EPANET2, which is a widely used simulation model for a water distribution system. The model results are in good agreement with EPANET2 for steady-state simulation. However, the hydraulic simulation results under unsteady flows differ from those of EPANET2, which causes a deviation in water quality prediction. The proposed model is expected to be a component of an integrated operation model for a water distribution system if it is combined with a computational model for rapidly varying flows to estimate leakage, pipe roughness, and intensive water quality.