• 한국군사과학기술학회지
• /
• 제14권2호
• /
• pp.313-320
• /
• 2011

SPH(Smoothed Particle Hydrodynamics) is a gridless Lagrangian technique that is useful as an alternative numerical analysis method used to analyze high deformation problems as well as astrophysical and cosmological problems. In SPH, all points within the support of the kernel are taken as neighbours. The accuracy of the SHP is highly influenced by the method for choosing neighbours from all particle points considered. Typically a linked-list method or tree search method has been used as an effective tool because of its conceptual simplicity, but these methods have some liability in anisotropy situations. In this study, convex hull algorithm is presented as an improved method to eliminate this artifact. A convex hull is the smallest convex set that contains a certain set of points or a polygon. The selected candidate neighbours set are mapped into the new space by an inverse square mapping, and extract a convex hull. The neighbours are selected from the shell of the convex hull. These algorithms are proved by Fortran programs. The programs are expected to use as a searching algorithm in the future SPH program.

• 대한기계학회논문집B
• /
• 제25권12호
• /
• pp.1802-1812
• /
• 2001

The purpose of this study is to analyze numerically the movement of particles included in turbulent fluid flow characteristics of metallic surfaces. To describe fluid flew, the incompressible Navier-Stokes equation discretized by the finite volume method were solved on the non-orthogonal coordinates with non-staggered variable arrangement, and the k-$\xi$ turbulence model was adapted. After fluid flow was calculated, particle movement was predicted from the Lagrangian approaches. Non-essential complexities were avoided by assuming that the particles had spherical shapes and the Stoke's drag formula only consisted of external farces acting upon them. In order to validate the numerical calculations, the results were compared with the experimental data reported in literature and agreed well with them. The drag force coefficient equation showed better agreement with the experimental data in the prediction of particle movement than the correction factor equation. Impact velocity and impact angle increased as inlet turbulence intensity decreased, relative jet height was lower. or the Reynolds number was larger.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.916-921
• /
• 2003

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method(FVM). This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications.

• Korea-Australia Rheology Journal
• /
• 제21권3호
• /
• pp.193-200
• /
• 2009

We present a direct simulation technique for two-dimensional mold-filling simulations of fluids filled with a large number of circular disk-like rigid particles. It is a direct simulation in that the hydrodynamic interaction between particles and fluid is fully considered. We employ a pseudo-concentration method for the evolution of the flow front and the DLM (distributed Lagrangian multipliers)-like fictitious domain method for the implicit treatment of the hydrodynamic interaction. Both methods allow the use of a fixed regular discretization during the entire computation. The discontinuous Galerkin method has been used to solve the concentration evolution equation and the rigid-ring description has been introduced for freely suspended particles. A buffer zone, the gate region of a finite area subject to the uniform velocity profile, has been introduced to put discrete particles into the computational domain avoiding any artificial discontinuity. From example problems of 450 particles, we investigated the particle motion and effects of particles on the flow for both Newtonian and shear-thinning fluid media. We report the prolonged particle movement toward the wall in case of a shear-thinning fluid, which has been interpreted with the shear rate distribution.

• 한국정밀공학회지
• /
• 제29권6호
• /
• pp.632-639
• /
• 2012

A fully coupled thermo-mechanical model is adopted to study the temperature distribution and the material deformation in friction stir welding(FSW) process. Rotational speed is most important parameters in this research. Three dimension results under different process parameters were presented. Result indicate that the maximum temperature is lower than the melting point of the welding material. The higher temperature gradient occurs in the leading side of the workpiece. The maximum temperature can be increased with increasing the tool angular velocity, rpm in the current numerical modeling. In this research ABAQUS Ver.6.7 is to analyze a fully coupled thermo-mechanical model. ALE(Arbitrary Lagrangian-Eulerian) finite element formulation is used for the large deformation in FSW process and using the Mass scaling for the analysis time efficiency.

• 대한조선학회논문집
• /
• 제30권2호
• /
• pp.86-97
• /
• 1993

수치방법은 포텐셜 유동의 가정하에서 Semi-Lagrangian 기법을 사용하여 2차원 쇄기의 비선형운동과 축대칭 물체의 강제 상하동요 운동에 대해서 개발되었다. 2차원에서 Cauchy 이론은 경계를 따라서 복소포텐셜과 그것의 미분치를 계산하기 위해 적용되었고, 3차원에서 Rankinering 쏘오스가 사용되고 대수방정식을 풀기위해서 그린 제2정리를 이용하였다. 해는 완전한 사유표면 조건을 수치적분함으로서 시간전진시킨다. 수치계산 예는 정속도로 입수하는 쇄기형 주상체와 정지 상태로 부터 강제상하동요하는 문제를 택하였다. 쇄기입수 문제는 Chapman [4], Kim[11]의 계산결과와 비교된다. 위에서 적용된 기법을 이용하여 구한 시간영역에서 힘을 Fourier 변환함으로서 부가질량계수, 감쇄계수, 2차조화력등이 얻어지고 Yamashita[5]의 실험치와 비교된다.

• Wind and Structures
• /
• 제16권5호
• /
• pp.411-431
• /
• 2013

A computational study of vortex-induced transverse vibrations of a cylinder with low mass-damping is presented. An Arbitrary Lagrangian-Eulerian (ALE) formulation of the Unsteady Reynolds-Averaged Navier-Stokes equations (URANS), along with the Spalart-Allmaras (SA) one-equation turbulence model, are coupled conservatively with rigid body motion equations of the cylinder mounted on elastic supports in order to study the amplitude and frequency response of a freely vibrating cylinder, its flow-induced motion, Vortex Street, near-wake flow structure, and unsteady loading in a moderate range of Reynolds numbers. The time accurate response of the cylinder from rest to its limit cycle is studied to explore the effects of Reynolds number on the start of large displacements, motion amplitude, and frequency. The computational results are compared with published physical experiments and numerical studies. The maximum amplitudes of displacements computed for various Reynolds numbers are smaller than the experimental values; however, the overall agreement of the results is quite satisfactory, and the upper branch of the limit-cycle displacement amplitude vs. reduced velocity response is captured, a feature that was missed by other studies. Vortex shedding modes, lock-in phenomena, frequency response, and phase angles are also in agreement with experiments.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제4권2호
• /
• pp.83-95
• /
• 2012

In the present study, the structural response of breakwaters installed on container carriers against green water impact loads was numerically investigated on the basis of the fluid-structure interaction analysis. A series of numerical studies is carried out to induce breakwater collapse under such conditions, whereby a widely accepted fluid-structure interaction analysis technique is adopted to realistically consider the phenomenon of green water impact loads. In addition, the structural behaviour of these breakwaters under green water impact loads is investigated simultaneously throughout the transient analysis. A verification study of the numerical results is performed using data from actual collapse incidents of breakwaters on container carriers. On the basis of the results of a series of numerical analyses, the pressure distribution of green water was accurately predicted with respect to wave mass and velocity. It is expected that the proposed analytical methodology and predicted pressure distribution could be used as a practical guideline for the design of breakwaters on container carriers.

• 대한건축학회논문집:구조계
• /
• 제34권12호
• /
• pp.73-82
• /
• 2018

The objective of this study was to evaluate particle collection performance of electrostatic precipitator (ESP) integrated with double skin façade in naturally ventilated residential buildings using numerical method. To evaluate the efficiency, computational fluid dynamics (CFD) simulation based on electric potential and Lagrangian method was applied. The CFD model was validated by comparing the simulated results with the experimental data including thermal characteristic of double skin façade (DSF) and particle removal characteristic of electrostatic precipitator. The validation results showed that the root mean square error (RMSE) between predicted values and measured values of velocity and temperature in intermediate space of DSF was 1.2%. The adequacy of ion space charge density and turbulent model were determined. The RMSE between predicted values and measured values of collection efficiency of ESP was 9.2%. In addition, the case study was performed to present the application of the numerical method based on validation results of ESP integrated with façade.

• 한국물환경학회지
• /
• 제32권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.