• Journal of the Korean Geotechnical Society
• /
• v.33 no.12
• /
• pp.93-105
• /
• 2017

This study presents direct shear tests, model tests, and numerical simulations to assess the effect of reduction of soil strength because of saturation during formation of ground cave-in caused by damaged sewer pipe lines. The direct shear test results show that the saturation affects the cohesion of soil significantly although it does not influence the friction angle of soil. To experimentally reproduce ground cave-in, the model tests were performed. As ground cave-ins were accompanied with extreme deformation, conventional finite element method has difficulty in simulating them. The present study relies on generalized interpolation material point method, which is one of meshless methods. Although there are differences between the model test and numerical simulation caused by boundary conditions, incomplete saturation, and exclusion of groundwater flow, similar ground deformation characteristics are observed both in the model test and numerical simulation.

• International Journal of Vascular Biomedical Engineering
• /
• v.2 no.2
• /
• pp.27-32
• /
• 2004

The objective of the present study is to visualize the steady and pulsatile flow fields in a branching model by using a high-resolution PIV system. A bifurcated flow system was built for the experiments in the steady and pulsatile flows. Harvard pulsatile pump was used to generate the pulsatile velocity waveforms. Conifer powder as the tracing particles was added to water to visualize the flow fields. CCD cameras($1K{\times}1K$(high resolution camera) and $640{\times}480$(low resolution camera)) captured two consecutive particle images at once for the image processing of several cross sections on the flow system. The range validation method and the area interpolation method were used to obtain the final velocity vectors with high accuracy. The results of the image processing clearly showed the recirculation zones and the formation of the paired secondary flows from the distal to the apex of the branch flow in the bifurcated model. The results also indicated that the particle velocities at the inner wall moved faster than the velocities at the outer wall due to the inertial force effects and the helical motions generated in the branch flows as the flow proceeded toward the outer wall. Even though the PIV images from the high resolution camera were closer to the simulation results than the images from the low resolution camera at some locations, both results of the PIV experiments from the two cameras generally agreed quite well with the results from the computer simulations. Therefore, instead of using the expensive stereoscopic PIV or 3D PIV system, the three-dimensional flow fields in a bifurcated model could be easily and exactly investigated by this study.

• Journal of the Korea Computer Graphics Society
• /
• v.25 no.5
• /
• pp.21-30
• /
• 2019

We propose a physical simulation method based on the shallow water equation(SWE) to represent water surface effectively. In this paper, the water which can be represented has a much larger width compared to the depth does not have a large vertical direction flow. In order to calculate the water flow efficiently, we start with the shallow water equation as the governing equation, which is a simplified version of the Navier-Stokes equation. In order to numerically calculate the advection term of the SWE, we introduce a new conservtive USCIP(CUSCIP) method which improves the Constrained Interpolation Profile (CIP) method to preserve the physical quantity while increasing the numerical accuracy. The proposed method is based on Kim et. al.'s Unsplit Semi-lagrangian CIP[9], and calculates advection term with additional constraints on term that consider integral values. The experimental results show that the CUSCIP method is robust to the loss of physical quantity due to numerical dissipation, which improves wave detail and persistence.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.1
• /
• pp.70-81
• /
• 1999

To predict the viscous boundary layers and wakes around a ship, the CFD techniques are commonly used. For the efficient application of CFD tools to practical hull farms, a 3-D field grid generating system is developed. The present grid generating system utilizes the solution of Poisson equation. Sorenson's method developed for 2-D is extended into 3-D to provide the forcing functions controling grid interval and orthogonality on hull surface, etc. The weighting function scheme is used for the discretization of the Poisson equation and the linear equations are solved by using MSIP salver. The trans-finite interpolation is also employed to assure the smooth transition into boundary surface grids. To rove the applicability, the field grid systems around a container ship and a VLCC with bow and stem bulb are illustrated, and the procedures for generating 3-D field grid system are explained.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.243-253
• /
• 1999

The dual reciprocity boundary element method (DRBEM) is used to solve the Graetz problem of laminar flow inside circular duct. In this method the domain integral tenn of boundary integral equation resulting from source term of governing equation is transformed into equivalent boundary-only integrals by using the radial basis interpolation function, and therefore complicate domain discretization procedure Is completely removed. Velocity profile is obtained by solving the momentum equation first and then, using this velocities as Input data, energy equation Is solved to get the temperature profile by advancing from duct entrance through the axial direction marching scheme. DRBEM solution is tested for the uniform temperature and heat flux boundary condition cases. Local Nusselt number, mixed mean temperature and temperature profile inside duct at each dimensionless axial location are obtained and compared with exact solutions for the accuracy test Solutions arc in good agreement at the entry region as well as fully developed region of circular duct, and their accuracy are verified from error analysis.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.1-8
• /
• 2007

A code is developed using the hybrid Cartesian/immersed boundary method and it is applied to simulate flows around a three-dimensional deforming body. A new criterion is suggested to distribute the immersed boundary nodes based on edges crossing a body boundary. Velocities are reconstructed at the immersed boundary nodes using the interpolation along a local normal line to the boundary. Reconstruction of the pressure at the immersed boundary node is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other experimental and numerical results for the velocity profiles around a circular cylinder under the forced in-line oscillation and the pressure coefficient distribution on a sphere. The code is applied to simulate the flow fields around a plate whose tail is periodically flapping under a translation. The effects of the velocity and acceleration due to the deformation on the periodic shedding of pairs of tip vortices are investigated.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.1
• /
• pp.1-14
• /
• 2009

When the stagnation temperature increases, the specific heat does not remain constant and start to vary with this temperature. The gas is perfect, it's state equation remains always valid, except, it was called by gas calorically imperfect or gas at high temperatures. The purpose of this work is to develop a mathematical model for a normal shock wave normal at high temperature when the stagnation temperature is taken into account, less than the dissociation of the molecules as a generalisation model of perfect for constant heat specific. A study on the error given by the perfect gas model compared to our model is presented in order to find a limit of application of the perfect gas model. The application is for air.

• Transactions of Materials Processing
• /
• v.10 no.3
• /
• pp.206-213
• /
• 2001

Dynamic deformation of metallic materials mostly accompanies substantial amounts of deformation heat. Since the flow stress of deformation is sensitive to temperature, implication of heat due to plastic work is essential to the evaluation of constitutive relations. In this study, a series of compression tests were conducted for SAF 2507 super duplex stainless steel at various temperatures and strain rates. The accumulation of plastic work was calculated through numerical integration and converted into the elevation of temperature. Subsequent logarithmic interpolation deduced isothermal flow surfaces, which were primary input data of finite element analysis. Simple closed die forging process was analyzed and optimized with commercial FEM code applying both raw and calibrated material database. The effect of accounting deformation heat was more noticeable in high-speed forming process.

• Journal of Power System Engineering
• /
• v.12 no.4
• /
• pp.18-25
• /
• 2008

This paper investigated the new calibration algorithm of a straight-type five-hole pressure probe necessary for calculating three-dimensional flow velocity components. The new data reduction method Includes a look-up, a geometry transformation such as the translation and reflection of nodes, and a binary search algorithm. This new calibration map was applied up to the application angle, ${\pm}55^{\circ}$ of a probe. As a result, this data reduction method showed a perfect performance without any kind of interpolation errors In calculating yaw and pitch angle from the calibration map.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.44-50
• /
• 2015

The computational efficiency of flow simulations with Multi-resolution analysis (MRA) was enhanced via the boundary treatment of the computational domain. In MRA, an adaptive dataset to a solution is constructed through data decomposition with interpolating polynomial and thresholding. During the decomposition process, the basis points of interpolation should exceed the boundary of the computational domain. In order to resolve this problem, the weight coefficients of interpolating polynomial were adjusted near the boundaries. By this boundary treatment, the computational efficiency of MRA was enhanced while the numerical accuracy of a solution was unchanged. This modified MRA was applied to two-dimensional steady Euler equations and the enhancement of computational efficiency and the maintenance of numerical accuracy were assessed.