• Proceedings of the Korea Water Resources Association Conference
• /
• 1987.07a
• /
• pp.113-122
• /
• 1987

A numerical simulation technique of three-dimensional finite difference model is developed to study the groundwater flow system in Dcjima, an upland area which faces Kasumigaura Lake. For general perspectives of the groundwater flow system, a steady state three-dimentional model is simulated. For the sedimentary mud formations which are found in the representative formation, three situations of hydraulic conductivity are considered, representing an isotropic condition and situations where the horizontal permeability is equal to 10 times and 100times of the vertical one. The finite difference grid used in the simulation has 60x50x30=90,000 nodes. A converged solution with a tolerance of 0.001 meter of hydraulic head is set. Having determined the flow net by using a steady state three-dimensional model. the results for the three cases of hydraulic conductivity are compared with the results of tracer methods (Bae and Kayane 1987) With the aid of four representative vertical cross-sections, groundwater flow systems in the study area are assumed. Water balances for the three cases indicate very good agreement between total recharge and discharge in each case Analyses of groundwater flow system based on the tritium concentrations and water quality measurements (Bae and Kayane 1987) are confirmed by the numerical simulation and the results obtained by these two methods appeared to be in close agreement.

• Journal of Korean Society on Water Environment
• /
• v.29 no.3
• /
• pp.383-392
• /
• 2013

In this research, mixing behavior of the floating pollutant such as oil spill accidents was analyzed by studying the advection-diffusion of GPS floaters at water surface. The LPT (Lagrangian Particle Tracking) model of EFDC (Environmental Fluid Dynamics Computer Code) was used to simulate the motion of the GPS floater tracer. In the field experiment, 35 GPS floaters were injected at the Samun Bridge of Nakdong River. GPS floaters traveled to downstream about 700 m for 90 minutes. The field data by the GPS floater experiments were compared with the simulation in order to calibrate the parameter of LPT model. The turbulent diffusion coefficient of LPT model was determined as $K_H/hu^*$ = 0.17 from the scatter diagram. The arrival time of peak concentration and transverse diffusion from the simulation results were similar with the experiments from the concentration curves. Numerical experiments for anticipation of damage from floating pollutant were conducted in the same reach of the Nakdong River and the results show that the pollutant cloud transported to the left bank where the Hwawon pumping station is located. For this reason, it is suggested that the proper action should be needed to maintain the safety of the water withdrawal at the Hwawon pumping station.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.51 no.6
• /
• pp.53-62
• /
• 2009

The EFDC (Environmental Fluid Dynamics Code), a numerical model for simulating three-dimensional (3D) flow, transport, and biogeochemical processes in surface water systems including rivers, reservoirs, and estuaries, was applied to assess the effect of sea water and fresh water exchange rates ($Q_e$) on the mixing characteristics of a conservative pollutant (tracer) induced from upstreams and salinity in Saemangeum Lake, Korea. The lake has been closed by a 33 km estuary embankment since last April of 2006, and now seawater enters the lake partially through two sluice gates (Sinsi and Garyuk), which is driving the changes of hydrodynamic and water quality properties of the lake. The EFDC was constructed and calibrated with surveyed bathymetry data and field data including water level, temperature, and salinity in 2008. The model showed good agreement with the field data and adequately replicated the spatial and temporal variations of the variables. The validated model was applied to simulated the tracer and salinity with two different gate operation scenarios: RUN-1 and RUN-2. RUN-1 is the case of real operation condition ($Q_e=25,000,000\;m^3$) of 2008, while RUN-2 assumed full open of Sinsi gate to increase $Q_e$ by $120,000,000\;m^3$. Statistical analysis of the simulation results indicate that mixing characteristics of pollutants from upstream can be significantly affected by the amount of $Q_e$.

• Journal of Korean Society for Atmospheric Environment
• /
• v.19 no.4
• /
• pp.397-414
• /
• 2003

The performance of one-particle stochastic Lagrangian models for passive tracer dispersion are evaluated against measurements in horizontally-homogeneous neutrally-stratified atmospheric surface layer. State-of-the-technology models as well as classical Langevin models, all in class of well mixed models are numerically implemented for inter-model comparison study. Model results (far-downstream asymptotic behavior and vertical profiles of the time averaged concentrations, concentration fluxes, and concentration fluctuations) are compared with the reported measurements. The results are: 1) the far-downstream asymptotic trends of all models except Reynolds model agree well with Garger and Zhukov's measurements. 2) profiles of the average concentrations and vertical concentration fluxes by all models except Reynolds model show good agreement with Raupach and Legg's experimental data. Reynolds model produces horizontal concentration flux profiles most close to measurements, yet all other models fail severely. 3) With temporally correlated emissions, one-particle models seems to simulate fairly the concentration fluctuations induced by plume meandering, when the statistical random noises are removed from the calculated concentration fluctuations. Analytical expression for the statistical random noise of one-particle model is presented. This study finds no indication that recent models of most delicate theoretical background are superior to the simple Langevin model in accuracy and numerical performance at well.

• Journal of the Korean GEO-environmental Society
• /
• v.9 no.2
• /
• pp.39-46
• /
• 2008

Due to the growth in industrialization, potential hazards in subsurface environments are becoming increasingly significant. The extraction of the contaminant from the soil and movement of the water are restricted due to the low permeability and adsorption characteristics of the reclaimed soils. There are a number of approaches to in-situ remediation that are used in contaminated sites for removing contaminants. These include soil flushing, dual phase extraction, and soil vapor extraction. Among these techniques, soil flushing was the focus of the investigation in this paper. Incorporated technique with PVDs has been used for dewatering from fine-grained soils for the purpose of ground improvement by means of prefabricated vertical drain systems. The laboratory model tests were performed by using the flushing tracer solutions for silty soils and recorded the tracer concentration changes with the elapsed time and flow rates. The modeling was intended to predict the effectiveness and time dependence of the remediation process. Modeling has been performed on the extraction, considering tracer concentration and laboratory model test characteristics. The computer model used herein are SEEP/W and CTRAN/W, this 2-D finite element program allows for modeling to determine hydraulic head and pore water pressure distribution, efficiency of remediation for the subsurface environment. It is concluded that the coefficient of permeability of contaminated soil is related with vertical velocity and extracted flow rate. The vertical velocity and extracted flow rate have an effect on dispersivity and finally are played an important role in-situ soil remediation.

• Journal of Environmental Impact Assessment
• /
• v.7 no.1
• /
• pp.35-48
• /
• 1998

Gaussian type models have limitations on predicting a detailed description of the near flow and pollution leads over complex terrains under neutral atmospheric conditions. Also, most models used recently have lack of ability to include atmospheric reactions. The model based on the numerical solution of the time-averaged Navier-Stokes equations and conservation equations needs to be developed to improve the limitations mentioned above. When the model was applied to the Shihwa area where the tracer experiment had been carried out, the simulation results have a great difference from the experimental results. There are two reasons that make the difference between the results by the model and the experiment. First, the Shihwa area is not a complex terrain. Second, meteorological data is insufficient. Therefore, the model should be applied to predict the dispersion of air pollutants over complex terrain rather than flat terrain in order that the model could be verified because the model was developed for the prediction of the dispersion over a complex terrain.

• Journal of Advanced Marine Engineering and Technology
• /
• v.35 no.5
• /
• pp.608-615
• /
• 2011

The objective of this study is to investigate the flow characteristics of polluted air behavior in rectangular tunnels using a PIV system and a commercial CFD program. The PIV experiments are simulated by using the olive oil as the tracer particles in scaled rectangular tunnels. Each model has one of four different outlet vents, each dimensionless L/H ratio of which is 0, 0.375, 0.75 and 1.125, respectively as the locations of each outlet are away from the vertical centerline through the inlet. A commercial CFD program, ANSYS CFX, was used to examine the velocity fields and the pressure distributions in numerical simulations. The kinematic viscosity of the air flow of $1.51{\times}10^{-5}m^2/s$ and the flow velocity of 0.3 m/s at the inlet are given under the same conditions in order to analyze the polluted air flow characteristics experimentally and computationally. This study is considered to examine the effect of the outlet locations in the naturally ventilated tunnel models.

• Nuclear Engineering and Technology
• /
• v.39 no.4
• /
• pp.327-336
• /
• 2007

The influence of density differences on the mixing of the primary loop inventory and the Emergency Core Cooling (ECC) water in the downcomer of a Pressurised Water Reactor (PWR) was analyzed at the ROssendorf COolant Mixing (ROCOM) test facility. ROCOM is a 1:5 scaled model of a German PWR, and has been designed for coolant mixing studies. It is equipped with advanced instrumentation, which delivers high-resolution information for temperature or boron concentration fields. This paper presents a ROCOM experiment in which water with higher density was injected into a cold leg of the reactor model. Wire-mesh sensors measuring the tracer concentration were installed in the cold leg and upper and lower part of the downcomer. The experiment was run with 5% of the design flow rate in one loop and 10% density difference between the ECC and loop water especially for the validation of the Computational Fluid Dynamics (CFD) software ANSYS CFX. A mesh with two million control volumes was used for the calculations. The effects of turbulence on the mean flow were modelled with a Reynolds stress turbulence model. The results of the experiment and of the numerical calculations show that mixing is dominated by buoyancy effects: At higher mass flow rates (close to nominal conditions) the injected slug propagates in the circumferential direction around the core barrel. Buoyancy effects reduce this circumferential propagation. Therefore, density effects play an important role during natural convection with ECC injection in PWRs. ANSYS CFX was able to predict the observed flow patterns and mixing phenomena quite well.

• Journal of Korea Water Resources Association
• /
• v.37 no.12
• /
• pp.979-991
• /
• 2004

In this study, RMA2 and RMA4, the 2-D depth-averaged models, were employed to simulate the two-dimensional mixing characteristics of the pollutants in the natural streams. The velocity and depth were first calculated using RMA2, 2-D hydrodynamic model, and then the resulting flow field was inputted to RMA4, 2-D water quality model, to compute the concentration field. RMA models were verified using the velocity and concentration data measured in S-curved meandering channel. The results showed that the RMA2 model simulated well the phenomenon that the maximum velocity line is located at the Inner bank of meandering channel, and the RMA4 model was well adapted to reproduce the general mixing behavior and the separation of tracer clouds. Comparing model simulations with measured data in the field experiments, RMA2 model simulated well general flow field and tendency that the maximum velocity line skewed toward the outer bank which were found in field experiments. The simulations of RMA4 model showed that the center of the tracer cloud tends to follow the path in which the maximum velocity occurs. In this study, the dispersion coefficients are fine-tuned based on the measured coefficients calculated using field concentration data, and the results show reasonable agreement with predictive equations.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2540-2543
• /
• 2008

Liquid crystal displays (LCD's) are continuously coated with some chemicals in the clean room of a factory. Spreading of these chemicals is causing serious problems both in controlling clean room quality as well as to the workers inside the factory. It is required to alleviate or properly control the offensive odor which is mainly composed of propylene glycol mono ethyl acetate, novolak resin and photo active compound. The control strategy employed is to bleed the offensive odor gas out the clean room. A full scale 3D CFD model was created with anisotropic porous media, chemical species transport with no volumetric reaction, and thermal diffusion with propane gas (tracer gas) to simulate the odor spreading. A segregated implicit solver with standard k-$\varepsilon$ model is employed. The detailed CFD analysis made it possible to develop an effective method of ventilating the coater room and optimizing their capacities.