• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.92-101
• /
• 2012

Three-dimensional finite-difference simulation in a small-scale half-sphere basin with planar free-surface is performed for an arbitrary shear-dislocation point source. A new scheme to deal with free-surface boundary condition is presented. Then basin parameters are examined to understand main characteristics on ground-motion response in the basin. To analyze the frequency content of ground motion in the basin, spectral amplitudes are compared with each other for four sites inside and outside the basin. Also particle motions for those sites are examined to find which kind of wave plays a dominant role in ground-motion response. The results show that seismic energy is concentrated on a marginal area of the basin far from the source. This focusing effect is mainly due to constructive interference of the direct Swave with basin-edge induced surface waves. Also, ground-motion amplification over the deepest part of the basin is relatively lower than that above shallow basin edge. In the small-scale basin with relatively simple bedrock interface, therefore, the ground-motion amplification may be more related to the source azimuth or direction of the incident waves into the basin rather than depth of it.

• Journal of Korea Water Resources Association
• /
• v.50 no.1
• /
• pp.55-64
• /
• 2017

This study presents a numerical model for simulating turbidity currents using the ULTIMATE scheme. For this, the layer-averaged model is used. The model is applied to laboratory experiments, where the flume is composed of sloping and flat parts, and the characteristics of propagating turbidity currents are investigated. Due to the universal limiter of the ULTIMATE scheme, the frontal part of the turbidity currents at a sharp gradient without numerical oscillations is computed. Simulated turbidity currents propagate super-critically to the end of the flume, and internal hydraulic jumps occur at the break-in-slope after being affected by the downstream boundary. It is found that the hydraulic jumps are computed without numerical oscillations if Courant number is less than 1. In addition, factors that affect propagation velocity of turbidity currents is studied. The particle size less than $9{\mu}m$ does not affect propagation velocity but the buoyancy flux affects clearly. Finally, it is found that the numerical model computes the bed elevation change due to turbidity currents properly. Specifically, a discontinuity in the bed elevation, arisen from the hydraulic jumps and resulting difference in sediment entrainment, is observed.

• Geomechanics and Engineering
• /
• v.15 no.5
• /
• pp.1113-1124
• /
• 2018

Geological dynamic hazards during deep coal mining are caused by the failure of a composite system consisting of the rock and coal layers, whereas the joint in coal affects the stability of the composite system. In this paper, the compression test simulations for the rock-coal combined body with single joint in coal were conducted using $PFC^{2D}$ software and especially the effects of joint length and joint angle on strength and failure characteristics in a rock-coal combined body were analyzed. The joint length and joint angle exhibit a deterioration effect on the strength and affect the failure modes. The deterioration effect of joint length of L on the strength can be neglected with a tiny variation at ${\alpha}$ of $0^{\circ}$ or $90^{\circ}$ between the loading direction and joint direction. While, the deterioration effect of L on strength are relatively large at ${\alpha}$ between $30^{\circ}$ and $60^{\circ}$. And the peak stress and peak strain decrease with the increase of L. Additionally, the deterioration effect of ${\alpha}$ on the strength becomes larger with the increase of L. With the increase of ${\alpha}$, the peak stress and peak strain first decrease and then increase, presenting "V-shaped" curves. And the peak stress and peak strain at ${\alpha}$ of $45^{\circ}$ are the smallest. Moreover, the failure mainly occurs within the coal and no apparent failure is observed for rock. At ${\alpha}$ between $30^{\circ}$ and $60^{\circ}$, the secondary shear cracks generated in or close to the joint tips, cause the structural instability failure of the combined body. Therefore, their failure models present as a shear failure along partial joint plane direction and partially cutting across the coal body or a shear failure along the joint plane direction. However, at ${\alpha}$ of $60^{\circ}$ and L of 10 mm, the "V-shaped" shear cracks cutting across the coal body cause its final failure. While crack nucleations at ${\alpha}$ of $0^{\circ}$ or $90^{\circ}$ are randomly distributed in the coal, the failure mode shows a V-shaped shear failure cutting across the coal body.

• Journal of the Korean Vacuum Society
• /
• v.10 no.3
• /
• pp.380-386
• /
• 2001

A nuclear microprobe system with adjustable precision object slits and a magnetic quadrupole doublet was designed by the beam optics simulation using a first order matrix formalism, and installed in a $30^{\circ}$ beam line connected with KIGAM 1.7 MV Tandem VDG Accelerator. Demagnification factors for x and y axis are calculated to be 25 and 4.9, respectively, and a minimum beam spot side is expected to be about 5 $\mu\textrm{m}$ for 3 MeV proton beams with a current of about 1 nA. A multi-purpose octagonal target chamber has been built to facilitate MeV ion-beam analytical techniques of PIXE, RBS, ERDA, and ion beam micro-machining. It contains X-ray and particle detectors, a zoom microscope, a Faraday cup, a 4-axis sample manipulator and a high vacuum pumping system. The system performance of the nuclear microprobe is now being tested, and automatic manipulator control and data acquisition system will be installed for routine applications of micro ion-beam analytical techniques.

• The Journal of Engineering Geology
• /
• v.19 no.1
• /
• pp.81-88
• /
• 2009

Two dimensional Monte-Carlo simulations of a non-reactive solute plume in isotropic porous media which are physically and chemically heterogeneous are conducted to determine the variations of moment. Retardation factors of 1, 2 and 5 are given to ascertain how the second moments are changed as adsorption increased. Retarded longitudinal second spatial moment, ${Z_{11}}^{'R}(t',l')$, increased during the transport process and as the dimensionless lengths of line plume source, $l_2'$, increased. ${Z_{11}}^{'R}(t',l')$ decreased as the retardation factors increased, and the simulated moments fit well to the first-order analytical results. Retarded longitudinal plume centroid variance, ${Z_{11}}^{'R}(t',l')$, decreased as the dimensionless lengths of line plume source, $l_2'$, increased and as the retardation factor increased. The result indicates that the uncertainty about the plume center decreased, and the ergodic condition for the second spatial moments is far from reaching. Simulated longitudinal one particle displacement covariance, ${Z_{11}}^{'R}(t')$, well consistent with the first-order analytical results for the three degrees of retardation factors of 1, 2 and 5 respectively. It is, consequently, concluded that the retarded longitudinal second moments could be produced by stochastic simulation, and that the first-order analytical results definitely provides very close values of the longitudinal retarded moments.

• Journal of the Korean Vacuum Society
• /
• v.10 no.4
• /
• pp.431-439
• /
• 2001

In this paper, the a-Diamond films were synthesized using filtered vacuum arc source (FVAS), FVAS was composed of a torus structure with bending angle of 60 degree. The radius of torus was 266 mm, the radius of plasma duct was 80 mm and the total length was 600 mm. The magnet parts were composed of one permanent magnet and five solenoid magnets. The plasma duct was electrically isolated from the ground so that a bias voltage could be applied. The baffles inside plasma duct were installed in order to prevent the recoil effect of macro-particles. Cathode was made of graphite with 80 mm in diameter. The effects of solenoid magnet on plasma extraction were investigated by computer simulation and experiment using Taguchi's methode. The source and extraction magnet affected the arc stabilization. The extraction beam current was maximized with low value of the source magnet current and high value of the filtering magnet current. The beam current density was 3.2 mA/$\textrm{cm}^2$ and average deposition rate was 5 $\AA$/sec when the currents of arc discharge, source, extraction, bending, deflection and outlet magnet were 30 A, 1 A, 3 A, 5 A, 5 A, and 5 A, respectively. The beam current density and the efficiency of beam transportation were increased with the positive bias voltage of the plasma duct.

• Journal of Korea Water Resources Association
• /
• v.52 no.12
• /
• pp.1047-1055
• /
• 2019

Various types of numerical modeling techniques are used to predict the behavior of pollutants under various water environmental conditions in the event of a water pollutant accident. Among them, a hydraulic model that can consider water flow characteristics is the most basic and very important. The process of evaluating whether the hydraulic model accurately predicts the applied river characteristics is very important. In the verification of the modeling result, the measuring data are often used in the river. Currently, ADCP and FlowTrackers are widely used to measure the flow velocity of rivers. However, ADCP is not accurate when the depth is less than 0.6 m and also when the ratio of irreversibility near the surface is more than 40%. Futhermore, FlowTracker has a limitation in measuring at high depth and high velocity due to the direct measurement method in rivers. Simuation results, which are validated by these methods, are not reliable for low depth conditions of low flowrate and high velocity conditions of high flowrate. In this study, Lagrangian GPS floaters which measures physical quantity of water according to particle movement is used without the conventional method measured by Eulerian technique. The verification method of the model results was studied by comparing the simulation results of the hydraulic model with the velocities measured using the GPS floaters. When comparing the traveling distance of the GPS floaters with the traveling distance of the LPT simulations, the average error rate was 13.6% on distances, and the average error rate was 3.2% on velocities except for the stagnant section. Therefore, GPS floaters can be used for a correction and verification method of hydraulic model simulations.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.3
• /
• pp.8-13
• /
• 2018

A fluid-structure interaction analysis should be performed to evaluate the behavior of the internal fuel and its influence in order to confirm the structural soundness of the fuel tank against external impacts. In the past, fluid-structure interaction analyses have been limited to the obtention of numerical simulation results due to the need for considerable computational resources and excessive computation time. However, recently, computer performance has been dramatically improved, enabling complex numerical analyses such as fluid-structure interaction analysis to be conducted. Lagrangian and Euler coupling methods and Lagrangian based analysis methods are mainly used for fluid-structure interaction analysis. Since both of these methods have their advantages and disadvantages, it is necessary to select the more appropriate one when conducting a numerical analysis. In this study, a numerical analysis of a crash impact test for a fuel tank is performed using ALE. The purpose of the numerical analysis is to estimate the possibility of failure of the fuel tank mounted inside the container when it is subjected to a crash impact. As a result of the numerical analysis, the fluid behavior inside the fuel tank is investigated and the stress generated in the fuel tank and the container structure is calculated, thereby enabling the possibility of fuel tank failure and leakage of the internal fluid to be evaluated.

• Journal of Korea Water Resources Association
• /
• v.33 no.5
• /
• pp.581-592
• /
• 2000

Turbulence structure and suspended sediment transport capacity in vegetated open-channel flows are investigated numerically in the present paper. The $\textsc{k}-\;\varepsilon$ model is employed for the turbulence closure. Mean velocity and turbulence characteristics including turbulence intensity, Reynolds stress, and production and dissipation of turbulence kinetic energy are evaluated and compared with measurement data available in the literature. The numerical results show that mean velocity is diminished due to the drag provided by vegetation, which results in the reduction of turbulence intensity and Reynolds stress. For submerged vegetation, the shear at the top of vegetation dominates turbulence production, and the turbulence production within vegetation is characterized by wakes. For emergent condition, it is observed that the turbulence generation is dominated by wakes within vegetation. In general, simulated profiles compares favorably to measured data. Computed values of eddy viscosity are used to solve the conservation equation for suspended sediment, yielding sediment concentration more uniform over the depth compared with the one in the plain channel. The simulation reveals that the suspended load decreases as the vegetation density increases and the suspended load increases as the particle diameter decreases for the same vegetation density.

• The KIPS Transactions:PartA
• /
• v.16A no.6
• /
• pp.419-426
• /
• 2009

To realistically simulate fluid, the Navier-Stokes equations are generally used. Solving these Navier-Stokes equations on the Eulerian framework, the non-linear advection terms invoke heavy computation and thus Semi-Lagrangian methods are used as an approximated way of solving them. In the Semi-Lagrangian methods, the locations of advection sources are traced and the physical values at the traced locations are interpolated. In the case of Stam's method, there are relatively many chances of numerical losses, and thus there have been efforts to correct these numerical errors. In most cases, they have focused on the numerical interpolation processes, even simultaneously using particle-based methods. In this paper, we propose a new approach to reduce the numerical losses, through improving the tracing method during the advection calculations, without any modifications on the Eulerian framework itself. In our method, we trace the grids with the velocities which will let themselves to be moved to the current target position, differently from the previous approaches, where velocities of the current target positions are used. From the intuitive point of view, we adopted the simple physical observation: the physical quantities at a specific position will be moved to the new location due to the current velocity. Our method shows reasonable reduction on the numerical losses during the smoke simulations, finally to achieve real-time processing even with enhanced realities.