• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.474-484
• /
• 2018

In this study, a new numerical modeling system was proposed to predict oil spills, which increasingly occur at sea as a result of abnormal weather conditions such as global warming. The hydrodynamic conditions such as the flow velocity needed to calculate oil dispersion were estimated using a three dimensional hydrodynamic model based on the Navier-Stokes equation, which considered all of the physical variations in the vertical direction. This improved the accuracy compared to those estimated by the conventional shallow water equation. The advection-diffusion model for the spilled oil was combined with the hydrodynamic model to predict the movement and fate of the oil. The effects of absorption, weathering, and wind were also considered in the calculation process. The combined model developed in this study was then applied to various test cases to identify the characteristics of oil dispersion over time. It is expected that the developed model will help to establish initial response and disaster prevention plans in the event of a nearshore oil spill.

• Atmosphere
• /
• v.14 no.4
• /
• pp.3-18
• /
• 2004

Neglecting the vertical transport from the surface, most of the previous studies on the long-range transport of pollutants have only considered the horizontal transport caused by the free atmosphere wind. I used a three dimensional numerical model, MM5 (The fifth generation Penn State Univ./NCAR Mesoscale Model) for the simulation of vertical transport of pollutants and investigated the mechanism of the vertical transport of atmospheric pollutants between planetary boundary layer(PBL) and free atmosphere by fronts. From the three dimensional simulation of MM5, the amount of pollutants transport from PBL to free atmosphere is 48% within 18 hour after the development of front, 55% within 24 hour, and 53% within 30 hour. The ratios of the vertically transported pollutant for different seasons are 62%, 60%, 54%, and 43% for spring, summer, fall, and winter, respectively. The most active areas for the vertical transport are the center of low pressure and the warm sector located east side of cold front, in which the strong upward motion slanted northward occurs. The horizontal advection of pollutants at the upper level is stronger than at the lower level simply because of the stronger wind speed. The simulation results shows the well known plum shape distribution of pollutants. The high concentration area is located in the center and north of the low pressure system, while the second highest concentration area is in the warm sector. It is shown that the most important mechanism for the vertical transport is vertical advection, while the vertical diffusion process plays an important role in the redistribution of pollutants in the PBL.

• Journal of Korea Water Resources Association
• /
• v.45 no.5
• /
• pp.473-480
• /
• 2012

A two-dimensional depth-integrated hydrodynamic and a depth-averaged passive scalar transport models were developed by using a Compact Finite Volume Method (CFVM) which can assure a higher order accuracy. A typical wave current interaction experimental data set was compared with the computed results by the proposed CFVM model, and resonable agreements were observed from the comparisons. One and two dimensional scalar advection tests were conducted, and very close agreements were observed with very little numerical diffusion. Finally, a turbulent mixing simulation was done in an open channel flow, and a reasonable similarity with LES data was observed.

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

• Journal of the Korean earth science society
• /
• v.39 no.4
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

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

• 한국가시화정보학회:학술대회논문집
• /
• 2005.12a
• /
• pp.99-102
• /
• 2005

This study was conducted for obtaining the optimized data to build the mixer or micro fluid device as measuring the three dimensions flow field in micro mixer. To acquire the rapid diffusion on the region of low Reynolds (Re < 100), the staggered herringbone mixer using chaotic advection was selected in this case. At first, by conducting the numerical analytical virtual experiment using CFD-ACE+, three dimensions flow field in the micro mixer was estimated As this flow field was proven using defocusing particle tracing method, the behavior of micro flow with three dimensional aspects could be analyzed. Numerical analysis and flow pattern in the micro mixer by experimental verification made to be able to analyze the chaotic advection. These can be important sources for building more optimized form. Verifying the information of three dimensional flow structure, these information can be used as the data for developing and improving the $\mu$ -TAS.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO method.

• Journal of Korean Society of Water Science and Technology
• /
• v.26 no.6
• /
• pp.3-12
• /
• 2018

The objective of this study was to evaluate the radius of influence of effluent of water treatment system developed for the purpose of improvement of reservoir water quality using fluorescent dye (Rhodamine-WT) tracer experiment and 3-D numerical model. The tracer experiment was carried out in a medium-sized agricultural reservoir with a storage capacity of $227,000m^3$ and an average depth of 1.6 m. A guideline with a total length of 160 m was installed at intervals of 10 m in the horizontal direction from the discharge part, and a Rhodamine measurement sensor (YSI 6130, measurement range $0-200{\mu}g/L$) was used to measure concentration changes in time, distance, and depth. Experimental design was established in advance through Jet theory and the diffusion process was simulated using ELCOM, a three dimensional hydraulic dynamics model. As a result of the study, the direct effect radius of the jet emitted from the applied water treatment system was about 50-70 m, and the radius of physical effect by the advection diffusion was judged to be 100-120 m. The numerical simulations of effluent advection-diffusion of the water treatment system using ELCOM showed very similar results to those of the impact radius analysis using the tracer experiment and jet flow empirical equations. The results provide valuable information on the spatial extent of the water quality improvement devices installed in the reservoir and the facility layout design.

• Journal of Environmental Science International
• /
• v.10 no.5
• /
• pp.343-350
• /
• 2001

To investigate air quality away from the coastal urban source region, we used a hybrid Eulerian - Lagrangian method which can describe the formation, transport, transform and deposition processes in complex terrain, with inclusion of shipping sources that were considered to be important emission in the coastal urban region. The result of the Eulerian advection - diffusion prediction was quite similar to that of the Lagrangian particle diffusion prediction. It showed that pollutants emitted from Sasang and Janglim industrial complexes can affect Hwamyeong and the coastal, respectively. During the daytime the concentration was low due to large deposition flux and terrain effect.