• Journal of Korean Society on Water Environment
• /
• v.32 no.3
• /
• pp.243-252
• /
• 2016

In cases of water pollution accidents, accurate prediction for arrival time and concentration of contaminants in a river is essential to take proper measures and minimize their impact on downstream water intake facilities. It is critical to fully understand the behavior characteristics of contaminants on river surface, especially in case of oil spill accidents. Therefore, in this study, the effects of main parameters of advection and diffusion of contaminants were analyzed and validated by comparing the results of Lagrangian particle tracking (LPT) simulation of Environmental Fluid Dynamic Code (EFDC) model with those of Global Position System (GPS)-equipped drifter experiment. Prevention scenario modeling was accomplished by taking cases of movable weir operation into account. The simulated water level and flow velocity fluctuations agreed well with observations. There was no significant difference in the speed of surface particle movement between 5 and 10 layer modeling. Therefore, 5 layer modeling could be chosen to reduce computational time. It was found that full three dimensional modeling simulated wind effects on surface particle movements more sensitively than depth-averaged two dimensional modeling. The diffusion range of particles was linearly proportional to horizontal diffusivity by sensitivity analysis. Horizontal diffusivity estimated from the results of GPS-equipped drifter experiment was 0.096 m²/sec, which was considered to be valid for applying the LPT module in this area. Finally, the scenario analysis results showed that particle movements could be stagnant when discharge from the upstream weir was reduced, implying the possibility of securing time for mitigation actions such as oil boom installation and wiping oil contaminants. The outcomes of this study can help improve the prediction accuracy of particle tracking simulation to establish the most suitable mitigation plan considering the combination of movable weir operation.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.76-85
• /
• 1993

In order to predict the free surface signature of turbulent ship wakes two things are essential; a basic understanding of the mechanism of turbulent vortical flow/free surface interactions and a mathematical model to accurately predict the signature. The goal of the study described here is both to supplement experimental work to obtain basic understanding, as well as to condense this understanding in a model(or models) that captures the essential phenomena and thus allows predictions. To do so we followed two main paths guided by experimental observations. One is full simulations of the flow using the clavier-Stokes equations. The other is a vortex modeling, where the vortical structures of the flows are approximated by idealized structures, an the interaction assumed to be essentially inviscid. These approaches complement each other. Full simulations are only applicable to small scale phenomena, where the system is simple, and the Reynolds number is low. The vortex modeling, on the other hand, cannot represent essentially viscous aspects of the problem such as the effect of contamination gradient. Obviously, the modeling is what may eventually lead to a prediction method; the full simulations-too limited to mimic all but the simplest circumstances-are to aid and support the construction of realistic models. We address two-dimensional aspects of the vortex/free surface interaction first. Secondly we obtain some basic understanding of the interaction process through an experiment and then talk about several three-dimensional problems hoping to develop a successful prediction model.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the

• Tunnel and Underground Space
• /
• v.22 no.2
• /
• pp.106-119
• /
• 2012

Roughness of rock joint has generally been characterized based upon geometrical aspects of a two-dimensional surface profile. The appropriate description of joint roughness, however, should consider the features of roughness mobilization at contact areas under normal and shear loads. In this study, direct shear tests were conducted on the replicas of tensile fractured gneiss joints and the influence of the shear direction on the shear behavior and effective roughness was examined. In this procedure, a joint surface was represented as a group of triangular planes, and the steepness of each plane was characterized using the concepts of the active and inactive micro-slope angles. The contact areas at peak strength which were estimated by a numerical method showed that the locations of the contact areas were mainly dependent on the distribution of the micro-slope angle and the shear behavior of joint was dominated by only the fractions with active micro-slope angles. Therefore, a three-dimensional coefficient for the quantification of rock joint roughness is proposed based on the distribution of active micro-slope angle: active roughness coefficient, $C_r$. Comparison of the active roughness coefficient and the peak shear strength obtained from the experiment suggests that the active roughness coefficient is the effective parameter to quantify the surface roughness and estimate the shear behavior of rock joint.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.1
• /
• pp.537-544
• /
• 2017

An experimental study was carried out to investigate the characteristics of the propagation distance of a flood wave considering the levee failure speed in a flat inundation area. The Ritter solution for one dimensional flow was considered to formulate the experimental results and a representative form with coefficients of k and m, which consider the three dimensional flow characteristics, was applied. The experiments showed that the propagation velocity of the wave front in the inundation area was influenced by the levee breach speed as well as the initial water level, which is a significant variable representing the flood wave behavior. In addition, coefficients k and m are not constants, but variables that vary with levee breach speed. An empirical formula was also suggested using the experimental results in the form of the relationships between k and m. In this study, a large-scale experiment for flood inundation was carried out to examine the behavior of flooding in the inundated area and the relationships between the levee breach speed and wave-front propagation velocity were suggested based on the experimental results. These research results are expected to be used as the baseline data to draw a flow inundation map, establish an emergency action plan, and verify the two-dimensional numerical model.

• Tunnel and Underground Space
• /
• v.29 no.4
• /
• pp.262-279
• /
• 2019

In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the environment can be exposed to radioactivity, it is very important to clarify gas migration in terms of long-term integrity of the engineered barrier system. In particular, it is necessary to identify the hydro-mechanical mechanism for the dilation flow, which is a very important gas flow phenomenon only in medium containing large amounts of clay materials such as bentonite buffer, and to develop and validate new numerical approach for the quantitative evaluation of the gas migration phenomenon. Therefore, in this study, we developed a two-phase flow model considering the mechanical damage model in order to simulate the gas migration in the engineered barrier system, and validated with 1D gas flow modelling through saturated bentonite under constant volume boundary conditions. As a result of numerical analysis, the rapid increase in pore water pressure, stress, and gas outflow could be simulated when the dilation flow was occurred.

• Journal of Welding and Joining
• /
• v.17 no.1
• /
• pp.88-96
• /
• 1999

The effects of root gap on welding residual stress and deformation are dealt with the multi-pass weldment with three kinds(0, 6, 30mm) of root gap by F.E.M common code, and then compared with experiment data. In this analysis, an 100% ramp heat input model was used to avoid numerical convergence problem due to an instantaneous increase in temperature near the fusion zone, and the effect of a moving arc in a two dimensional plane was also included. During the analysis, a small time increment was applied in a period with instantaneous temperature fluctuation while a large time increment was used in the rest period. The residual stress is distributed as symmetric types and maximum value is also equivalent when the weldment with 0mm and 6mm root gap is welded. In the case of 30mm root gap welding, the distribution of the residual stress extends over a wide range as asymmetric types due to the built-up weld, and most of the residual stress is biased in the side of a built-up weld part. In case of 0mm gap welding and 6mm gap welding, a little angular distortion occurs, but the level of deformation is small. When the weldment with 30mm root gap is welded, the angular deformation of the asymmetric types, however, occurs larger than the other specimens. The experimental and the analytic results show good coincidence and indicate that the welding residual stress and deformation distribution of 30 mm root gap specimen may be asymmetric and the amplitude is larger than those of root gap specimen under standard.

• Ecology and Resilient Infrastructure
• /
• v.8 no.4
• /
• pp.165-178
• /
• 2021

We studied driftwood behaviors including generation and deposition in a tsunami using a numerical simulation. We used an integrated two-dimensional numerical model, which included a driftwood dynamics model. The study area was Sendai, Japan. Observation data collected by Inagaki et al. (2012) were used to verify the simulation results by comparing them with driftwood deposition patterns. A simplified model was developed to consider the threshold of driftwood generation by the drag force of water flows. To consider the volume of driftwood generated, we estimated the total wood number in the study area using Google Earth. Therefore, we simulated more than 13,000 pieces of driftwood that were generated and transported inland from approximately 300,000 trees that were growing in the forest. The final distribution of the driftwood was similar to the observation data. The reproducibility of the generation and deposition patterns of driftwood showed good agreement in terms of longitudinal deposition pattern. In the future, a sensitivity analysis on driftwood parameters, such as the size of the wood, boundary conditions, and grid size, will be implemented to predict the travel patterns of driftwood. Such modeling will be a useful methodology for disaster prediction based on water flow and driftwood.

• Tunnel and Underground Space
• /
• v.17 no.2 s.67
• /
• pp.109-118
• /
• 2007

The excavation damaged zone (EDZ) is an area around an excavation where in situ rock mass properties, stress condition. displacement. groundwater flow conditions have been altered due to the excavation. Various studies have been carried out on EDZ, but most studies have been focused on the mechanical bahavior of EDZ by in situ experiment. Even though the EDZ could potentially form a high permeable pathway of groundwater flow, only a few studies were performed on the analysis of groundwater flow in EDZ. In this study, the' hydraulic EDZ' was defined as the rock Lone adjacent to the excavation where the hydraulic aperture has been changed due to the excavation. And hydraulic EDZ (hydraulic aperture changed zone) estimated by two-dimensional DEM program was considered in three-dimensional DFN model. From this approach the groundwater flow characteristics corresponding to hydraulic aperture change were examined. Together. a parametric study was performed to examine the boundary conditions that frequently used in DFN analysis such as constant head or constant flux condition. According to the numerical analysis, hydraulic aperture change induced by the hydraulic-mechanical interaction becomes one of the most important factors Influencing the hydraulic behavior of jointed rock masses. And also from this study, we suggest the proper boundary condition in three-dimensional DFN model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.609-620
• /
• 1992

A numerical solution of temperature and thermally induced stress in a wafer during rapid thermal processing (R.T.P) is obtained, and an analysis of onset and propagation of slip is performed and compared with experiment. In order to calculate temperature distribution of a wafer in R.T.P system, heat conduction equation that incorporated with radiative and convective heat transfer model is solved, and the solution of the equation is calculated numerically using alternating direction implicit (A.D.I) method. In dealing with radiative heat transfer, a partially transparent body that absorbs the radiation energy is assumed and this transparent body undergoes multiple internal reflections and absorptions. Two dimensional (assuming plane stress) thermoelastic constitutive equation is used to calculate thermal stress induced in a wafer and finite element method is employed to solve the equation numerically. The stress resolved in the slip directions on the slip planes of silicon is compared with the yield stress of silicon in order to predict the slip. The result of the analysis shows that the wafer temperature at which slip occurs is affected by the heating rate of the R.T.P system. It is observed that once slip occurs in the wafer, the slip grows.