• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.1-10
• /
• 2020

The application of multiphase flows is increasingly being applied to analyze phenomena such as single phase flows where the fluid boundary changes continuously over time or the problem of mixing a liquid phase and a gas phase. In particular, multiphase flow models that take into account incompressible Newtonian fluids for liquid and gas are often applied to solve the problems of the free water surface such as wave fields. In general, multi-phase flow models require time-based the surface tracking of each fluid's phase boundary, which determines the accuracy of the final calculation of the model. This study evaluates the advection performance of representative VOF-type boundary tracking techniques applied to various CFD numerical codes. The effectiveness of the FCT method to control the numerical flux to minimize the numerical diffusion in the conventional VOF-type boundary tracking method and advection calculation was mainly evaluated. In addition, the possibility of tracking performance of free surface using CIP method (Yabe and Aoki, 1991) was also investigated. Numerical results show that the FCT-VOF method introducing an anti-diffusive flux to precent excessive diffusion is superior to other methods under the confined conditions in this study. The results from this study are expected to be used as an important basic data in selecting free surface tracking techniques applied to various numerical codes.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.2
• /
• pp.104-110
• /
• 1997

Increase of marine transpotation in coastal area frequently yields oil spill accidents due to collision or grounding of oil tankers, which affects great deal of damages on ocean environments. Exact prediction of oil pollution area in time domain, which is called oil map, is very important for effective and efficient oil recovery and minimization of environmental damage. The prediction is carried out by considering the two distinct processes which are initial diffusion on the still water surface and advection due to tide, wind wave induced surface currents. In the present paper, only the initial diffusion is dealt with. Somewhat new simulation model and its numerical scheme are proposed to predict it. Simple diffusion experiment is also carried out to check the validity of the present method. Furthermore, some example simulations are performed for virtual oil spill accident. Quite realistic oil map including oil thickness distributions can be obtained by the present model.

• Journal of Radiation Protection and Research
• /
• v.12 no.2
• /
• pp.28-36
• /
• 1987

Gaussian plume model is used to assess environmental dose for abnormal radioactive release in nuclear facility, but there has a problem to use it for complex terrain. In this report, MATTEW and WIND04 Codes which had been verified were used to calculate wind field in the complex terrain. Under the base of these codes principle, wind fields were obtained from the calculation of the finite difference approximation for advection-diffusion equations which satisfy the mass-conservative law. Particle concentrations and external doses were calculated by using PIC model which approximate the particle to radioactive cloud, and atmospheric diffusion of the particles from the random walk method. The results show that the adjusted wind fields and the distributions of the exposure dose vary with the topography of the complex terrain.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.14 no.5
• /
• pp.226-232
• /
• 2004

One predicts the crystal growth rate of ZnSe with a low vapor pressure system in a horizontal configuration based on one dimensional advection-diffusion and two-dimensional diffusion-convection model. The present results show that for the ratios of partial pressures, s = 0.2 and 2.9, the growth rate increases with the temperature differences between the source and crystal. As the ratio of partial pressure approaches the stoichiometric value, s = 2 from s = 1.5 (zinc-deficient case: s < 2) and 2.9 (zinc-rich case: s > 2), the rate increases sharply. For the ranges from 1.5 to 1.999 (zinc-deficient case: s < 2) and from s = 9 to 2.9 (zinc-rich case: s > 2), the rate are slightly varied. From the viewpoint of the order of magnitude, the one-dimensional model for low vapor pressure system falls within the 2D predictions, which indicates the flow fields would be advective-diffusive. For the effects of gravitational accelerations on the rate, the gravitational constants are varied from 1 g to $10^{-6}$ g for $\Delta$T = 50 K and s = 1.5, the rates remain nearly constant, i.e., 211 mg/hr, which indicates Stefan flow is dominant over convection.

• Korean Journal of Computational Design and Engineering
• /
• v.18 no.6
• /
• pp.399-406
• /
• 2013

This paper describes a new simulation technique for advection-diffusion phenomena over the sea surface using the lattice Boltzmann method (LBM), capable of predicting oil dispersion from tankers. The LBM is used to solve the pollutant transport problem within the framework of the ocean environment. The sea space is represented by the lattices, where each lattice has the information on oil transportation. Since dispersed oils (i.e., oil droplets) at sea are transported by convection due to waves, buoyancy, and turbulent diffusion, the conservation of mass and many physical oil transport rules were used in the prediction model. Since the LBM is modeled using the uniform lattices and simple rules, it can be easily accelerated by the parallel mechanism, for example, GPU-accelerated method. The proposed model using the LBM is used to simulate a simple pollution event with the oil pollutants of 10,000 kL. The simulation results indicate that the LBM method accelerated with the GPU is 6 times faster than that without the GPU.

• Nuclear Engineering and Technology
• /
• v.34 no.3
• /
• pp.218-231
• /
• 2002

Using a three-dimensional numerical code, B3R developed for nuclide transport of an arbitrary length of decay chain in the buffer between the canister and adjacent rock in a high- level radioactive waste repository by adopting a finite difference method utilizing the control- volume scheme, some illustrative calculations have been done. A linear sorption isotherm, nuclide transport due to diffusion in the buffer and the rock matrix, and advection and dispersion along thin rigid parallel fractures existing in a saturated porous rock matrix as well as diffusion through the fracture wall into the matrix is assumed. In such kind of repository, buffer and rock matrix are known to be important physico-chemical harriers in nuclide retardation. To show effects of buffer and rock matrix on nuclide transport in HLW repository and also to demonstrate usefulness of B3R, several cases of breakthrough curves as well as three- dimensional plots of concentration isopleths associated with these two barriers are introduced for a typical case of decay chain of $^{234}$ Ulongrightarrow$^{230}$ Thlongrightarrow$^{226}$ Ra, which is the most important chain as far as the human environment is concerned.

• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.

• Journal of Korean Port Research
• /
• v.9 no.1
• /
• pp.45-56
• /
• 1995

Investigation of the physical behavior of dredged material disposal in coastal water includes estimations of water column concentration in the receiving water, exposure time, the initial deposition pattern as well as thickness of material at the dumping fields near the estuary area. Calculation based on vertical setting and horizontal advection of single particles ignore the effects of bulk properties of the disposed material, vertical and horizontal diffusion, and material dilution due to the entrainment of ambient water during descent. This paper focuses on the spatial and temporal changes in the dumping fields for the water column and bottom at a hypothetically confined coastal water, where the ambient time-invariant velocity and density profiles are applied, within the initial time period following the instantaneous release of the dredged material. This model accounts the behavior of material after release divided into three phases: convective descent, dynamic collapse and long-term passive dispersion