• Journal of the Korean Geotechnical Society
• /
• v.28 no.2
• /
• pp.71-78
• /
• 2012

Rainfall infiltration test under one dimensional condition is conducted to evaluate the effect of rainfall intensity on seepage velocity and infiltration characteristics for initial unsaturated sediment. Experimental results are compared with those numerical simulations with respect to variations of pore water pressure, degree of saturation and discharge velocity with time, and both results give good agreement. High rainfall intensity tends to increase seepage velocity almost linearly. But it shows rapid increase as rainfall intensity approaches saturated hydraulic conductivity of the sediment. In addition, the upper part of wetting front depth is partially saturated, not fully. Therefore, actual wetting front depth is considered to advance faster than theoretical prediction, which leads to slope instability of unsaturated slope due to surface rainfall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.577-589
• /
• 2005

The present study has numerically investigated two-dimensional laminar flow over a steadily rotating circular cylinder with a uniform planar shear, where the free-stream velocity varies linearly across the cylinder. It aims to find the combined effect of rotation and shear on the flow. Numerical simulations using the immersed boundary method are performed for the ranges of $-2.5{\le}\alpha{\le}2.5$ and $0{\le}K{\le}0.2$ at a fixed Reynolds number of Re=100, where a and K are respectively the dimensionless rotational speed and velocity gradient. Results show that the positive shear, with the upper side having the higher free-stream velocity than the lower one, favors the effect of the counter-clockwise rotation $(\alpha<0)$ but countervails that of the clockwise rotation $(\alpha>0)$. Accordingly, the absolute critical rotational speed, below which vortex shedding occurs, decreases with increasing K for $(\alpha>0)$, but increases for $\alpha>0$. The vortex shedding frequency increases with increasing \alpha (including the negative) and the variation becomes steeper with increasing K. The mean lift slightly decreases with increasing K regardless of the rotational direction. However, the mean drag and the amplitudes of the lift- and drag-fluctuations strongly depend on the direction. They all decrease with increasing K for $\alpha>0$, but increase for $\alpha<0$. Flow statistics as well as instantaneous flow folds are presented to identify the characteristics of the flow and then to understand the underlying mechanism.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.68-74
• /
• 2012

Near-wall effect on wakes behind particles is one of the important factors in precise tracking of particles in turbulent flows. However, most aerodynamic force models for particles did not fully consider the wall effect. In the present study, we focused on changes of hydrodynamic forces acting on a particle depending on wall proximity. To this end, we developed an immersed boundary method with multi-direct forcing incorporated to a fully implicit decoupling procedure for incompressible flows. We validate the present immersed boundary method through two-dimensional simulations of flow over a circular cylinder. Comprehensive parametric studies on the effect of the wall proximity on the drag and lift forces acting on an immersed circular cylinder in a channel flow are performed in order to investigate general flow patterns behind the circular cylinder for a wide range of Reynolds number (0.01 ${\leq}$ Re ${\leq}$ 200). As the cylinder is closer to the wall, the drag coefficient decreases while the lift coefficient increases with a local maximum. Maximum drag and lift coefficients for different wall proximities decrease with increment of Reynolds number. Normalized drag and lift coefficients by their maximum values show universal correlations between the coefficients and wall proximity in a low Reynolds number regime (Re ${\leq}$ 1).

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.8
• /
• pp.957-964
• /
• 2011

This investigation presents a finite element method to obtain the transmission properties of bulk elastic waves in piezoelectric band gap structures(phonon crystals) for varying frequencies and modes. To this end, periodic boundary conditions are imposed on a three-dimensional model while both in-plane and out-of-plane modes are included. In particular, the mode decoupling characteristics between in-plane and out-of-plane modes are identified for each electric poling direction and the results are incorporated in the finite element modeling. Through numerical simulations, the proposed modeling method was found to be a useful, effective one for analyzing the wave characteristics of various types of piezoelectric phononic band gap structures.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.83-90
• /
• 2013

OpenMP is becoming more and more useful as a simple parallel processing paradigm on SMP (Shared Memory Multi-Processors) computing environment with the development of multi-core processors. However, very few data is available publically regarding the OpenMP performance in CFD (Computational Fluid Dynamics). In the present study a CFD test suite is prepared for the performance evaluation of OpenMP on various multi-core systems. The test suite is composed of two-dimensional numerical simulations for inviscid/viscous and reacting/non-reacting flows using three different levels of grid systems. One to five test runs were carried out on various systems from dual-core dual threads to 16-core 32-threads systems by changing the number of threads engaged for each test up to 80. The results exhibit some interesting results and the lessons learned from the tests would be quite helpful for the further use of OpenMP for CFD studies using multi-core processor systems.

• Journal of the Korean Society for Railway
• /
• v.11 no.3
• /
• pp.300-310
• /
• 2008

Through this study, a conceptual design for the next generation high-speed EMU has been derived to meet the crash worthiness requirements of the Korean rollingstock safety regulation. The crashworthiness regulations require some performance requirements for two heavy collision accident scenarios; a train-to-train collision at the relative speed of 36 km/h, and a collision against a standard deformable obstacle of 15 ton at 110km/h. The complete train set will be composed of 2TC-6M with 13 ton axle load, which is different from KTX with the power car of 17 ton axle load. Using theoretical and numerical analyses, a crashworthy conceptual design was derived in terms of mean crush forces and energy absorptions for principal crushable structures and devices. The derived conceptual design was evaluated and improved using one dimensional dynamic simulations for the bar-spring-damper-mass model. It is shown from the simulation results that the suggested conceptual design can easily satisfy domestic crashworthiness requirements.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.35 no.4
• /
• pp.356-365
• /
• 2002

One of the most important factors that cause eutrophication is nutrient materials containing nitrogen and phosphorus which stem from excreation of terrestial sources and release from sediment. Therefore, to improve water quality, the reduction of these nutrients loads should be indispensible. At this study, the three-dimensional numerical hydrodynamic and ecosystem model, which was developed by Institute for Resources and Environment of Japan, were applied to analyze the processes affecting the eutrophication. The residual currents, which were obtained by integrating the simulated tidal currents over 1 tidal cycle, showed the presence of a typical counterclockwise eddies between Gyewha and Garyuk island. Density driven currents were generated westward at surface and eastward at the bottom in Saemankeum area where the fresh waters are flowing into, The ecosystem model was calibrated with the data surveyed in the field of the study area in annual average. The simulated results were fairly good coincided with the observed values within relative error of $30\%$. The simulations of DIN and DIP concentrations were performed using ecosystem model under the conditions of $40\~100\%$ pollution load reductions from pollution sources. In study area, concentration of DIN and DIP were reduced to $59\%$ and $28\%$ in case of the $80\%$ reduction of the input loads from fresh water respectively. But pollution loads from sediment had hardly affected DIN and DIP concentration, The $95\%$ input load abatement is necessary to meet the DIN and DIP concentration of second grade of ocean water quality criteria.

• Journal of Welding and Joining
• /
• v.34 no.6
• /
• pp.42-46
• /
• 2016

Hydrogen assisted cracking (HAC) is one of the most complicated problem in welding. Huge amount of studies have been done for decades. Based on them, various standards have been established to avoid HAC. But it is still a chronic problem in industrial field. It is well known that the main causes of the hydrogen crack are residual stress, crack susceptible micro structures and a certain critical level of hydrogen concentration. Even though the exact generating mechanism is unclear till today, it has been reported that the hydrogen level in the weld metal should be managed less than a certain amount to prevent it. Matsuda studied that the residual hydrogen level in the weld metal can be varied even if the initial hydrogen content is same. It is also insisted in this report that the residual hydrogen concentration is in stronger correlation with hydrogen crack than the initial hydrogen content. But, in practical point of view, the residual hydrogen is still hard to consider because measuring hydrogen level is time and cost consuming process. In this regard, numerical analysis is the only solution for considering the residual hydrogen content. Meanwhile, Takahashi showed the possibility of predicting the residual hydrogen by a rigorous FE analysis. But, few commercial software suitable for solving the weld metal hydrogen has been reported yet. In this study, two dimensional thermal - hydrogen coupled analysis was developed by using the commercial FE software MARC. Since the governing equation of the hydrogen diffusion is similar to the heat transfer, it is shown that the heat transfer FE analysis in association with hydrogen diffusion property can be used for hydrogen diffusion analysis. A series of simulation was performed to verify the accuracy of the model. For BOP (Bead-On-Plate) and the multi-pass butt welding simulations, remaining hydrogen contents in the weld metal is well matched with measurements which are referred from Kim and Masamitsu.

• Proceedings of KOSOMES biannual meeting
• /
• 2007.05a
• /
• pp.169-174
• /
• 2007

The purpose in having a control rod an a buoy system is to control the motion of it. The system may be composed entirely af a single circular cylinder finder and a lang mooring anchor cable. A control rod has one function to perform in meeting its purpose, and that is to develop a control force in consequence of its orientation and movement relative to the water. The forces and moments generated as a result of the effects of mutual interference, then determine the stability characteristics of the body. In this paper, the study of control-rod-attached buoy's 2-dimensional section was accomplished. model tests and numerical simulations had been carried out with different diameters of control rods. and varying the Reynolds number $Re=5,000{\sim}25,000$ based an the cylinder diameter(D=50mm) to. predict the performance af the body and the 2 frame particle tracking method had been used to obtain the velocity distribution in the flaw field 50mm circular cylinder had been used during the whale experiments and measured results had been compared with each other.