• 한국전산유체공학회지
• /
• 제16권2호
• /
• pp.1-10
• /
• 2011

A new computational program, which is based on the CIP/CCUP(Constraint Interpolation Profile/CIP Combined Unified Procedure) method, has been developed to numerically analyse sloshing phenomena dealt as multiphase-flow problems. For the convection terms of Navier-Stokes equations, the RCIP(Rational function CIP) method was adopted and the THINC-WLIC(Tangent of Hyperbola for Interface Capturing-Weighted Line Interface Calculation) method was used to capture the air/water interface. To validate the present numerical method, two-dimensional dam-breaking and sloshing problems in a rectangular tank were solved by the developed method in a stationary Cartesian grid system. In the case of sloshing problems, simulations by using a improved MPS(Moving Particle Simulation) method, which is named as PNU-MPS(Pusan National University-MPS), were also carried out. The computational results are compared with those of experiments and most of the comparisons are reasonably good.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집B
• /
• pp.535-540
• /
• 2000

The objective of the present study is to visualize the pulsatile flow fields by using three-dimensional computer simulation and the PIV system. A closed flow loop system was built for the steady and unsteady experiments. The Harvard pulsatile pump was used to generate the pulsatile pressure and velocity waveforms. Conifer powder as the tracing particles was added to water to visualize the flow field. Two consecutive particle images were captured by a CCD camera for the image processing. The cross-correlation method in combination with the moving searching area algorithm was applied for the image processing of the flow visualization. The pulsatile flow fields were visualized effectively by the PIV system in conjunction with the applied algorithm. The range validation and the area interpolation methods were used to obtain the final velocity vectors with high accuracy. The finite volume predictions were used to analyze three-dimensional flow patterns in the bifurcation model. The results of the PIV experiment and the computer simulation are in good agreement and the results show the recirculation zones and formation of the paired secondary flow distal to the apex of the bifurcated model. The results also show that the branch flow is pushed strongly to the inner wall due to the inertial force effect and helical motions are generated as the flow proceeds toward the outer wall.

• 한국컴퓨터그래픽스학회논문지
• /
• 제28권4호
• /
• pp.13-22
• /
• 2022

본 논문에서는 최소 입자 단위의 역동적인 회전 움직임을 나타낼 수 있는 MPM(Material Point Method) 기반 단일 프레임워크를 소개한다. 우리가 표현하고자 하는 입자는 다양한 형상(Shape)을 가질 수 있음과 동시에, 선형(Linear momentum), 회전(Angular momentum) 운동을 함께 묘사할 수 있다. 그 결과 기존 구형 입자의 선형 움직임만을 나타내던 입자 기반 시뮬레이션과는 달리, 시각적으로 단일 입자의 역동적인 모습을 표현할 수 있다. 제안하는 프레임워크는 회전 운동을 큰 변형(Large Deformation)으로부터 분해 및 추출 할 수 있다는 점에서 MPM을 활용하였다. 본 기법은 MPM 적분 과정 중 계산되는 변형 구배 텐서(Deformation Gradient Tensor)를 극 분해(Polar Decomposition)하는 과정을 통해 회전 텐서(Rotation Tensor)를 추출하고, 각 입자의 선형 운동과 함께 이를 적용하여 결과적으로 입자 자체의 회전, 선형 운동을 동시에 표현 하는 것이 가능하다. 본 연구에서는 제안하는 기법의 검증을 위해 바람에 흩날리며 회전하는 입자의 모습 및 움직이는 물체와 정지한 입자간의 상호작용 시뮬레이션을 기존 MPM을 이용한 시뮬레이션과의 비교를 통해 진행하였다.

• 천문학회보
• /
• 제39권2호
• /
• pp.65.2-65.2
• /
• 2014

We present a new hydrodynamic simulation code based on the Voronoi tessellation for estimating the density precisely. The code employs both of Lagrangian and Eulerian description by adopting the movable mesh scheme, which is superior to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for all surfaces of each Voronoi cell so as to update the hydrodynamic states as well as to move current meshes. Besides, the IEM (incremental expanding method) is devised to compute the Voronoi tessellation to desired degree of speed, thereby the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• 상하수도학회지
• /
• 제24권6호
• /
• pp.713-722
• /
• 2010

This paper has assessed the flow patterns and settling efficiency in the sedimentation basin using the particle tracking method of the CFD code and has reached the following conclusions: In the original design where no baffle is installed in the sedimentation basin, a large recirculating area where the flow stagnates is created in the right side of the sedimentation basin, with most of the particles moving to the left side of the sedimentation basin following the flow. This biased flow structure in the sedimentation basin reduces the residence time of particles and thereby undermines settling efficiency. The biased flow toward the left side of the sedimentation basin is alleviated by installing a baffle in the sedimentation basin, promptly reducing the fast flow of over 0.7 m/s in the inlet of the sedimentation basin to the rate below 0.2 m/s. In this paper's simulation conditions, if a one-sided baffle is to be installed in the sedimentation basin, placing it 15 meters away from the basin's inlet leads to the best settling efficiency; it has also been analyzed that installing a two-sided baffle-rather than a one-sided one-is a better option in terms of settling efficiency. The highest settling efficiency of 96.2% is achieved when the underwater length of the two-sided baffle is set at 8 meters.

• 한국융합학회논문지
• /
• 제10권2호
• /
• pp.201-206
• /
• 2019

해안 영역에서의 해수의 범람은 재산 및 인명피해를 야기한다. 이러한 피해를 최소화 하고자 방파시설을 설치하지만, 고정식 방파시설의 경우, 변화하는 해상상태에 대응할 수 없으며, 과도한 설정에 의한 설치는 환경적 피해를 일으킨다. 본 연구에서는 환경적 피해를 줄이고 선박의 항행에 지장이 없는 가변식 초기 쇄파 장치에 대한 연구를 실험적 방법과 입자법을 이용한 수치 해석적 방법으로 연구하였다. 초기 쇄파 장치의 유무에 따라 해안절벽으로 넘어오는 파의 범람은 큰 차이가 있었으며, 이는 실험과 수치 해석 모두 확인되었다. 또한 초기 쇄파 장치의 위치에 따른 영향 역시 고려되었으며, 위치 역시 중요한 요소임이 확인되었다.

• Tribology and Lubricants
• /
• 제34권6호
• /
• pp.279-283
• /
• 2018

In this study, we investigate the behavior of abrasive particles and change of the stick-slip pattern according to chemical mechanical polishing (CMP) process parameters when a large number of abrasive particles are fixed on a pad. The CMP process is simulated using the finite element method. In the simulation, the abrasive grains are composed of those used in the actual CMP process. Considering the cohesion of the abrasive grains with the start of the CMP process, abrasive particles with various sizes are fixed onto the pad at different intervals so that stick-slip could occur. In this analysis, we determine that when the abrasive particle size is relatively large, the stick-slip period does not change as the pressure increases while the moving speed is constant. However, if the size of the abrasive grains is relatively small, the amount of deformation of the grains increases due to the elasticity of the pad. Therefore, the stick-slip pattern may not be observed. As the number of abrasive particles increases, the stick-slip period and displacement decrease. This is consistent with the decrease in the von Mises yield stress value on the surface of the wafer as the number of abrasive grains increases. We determine that when the number of the abrasive grains increases, the polishing rate, and characteristics are improved, and scratches are reduced. Moreover, we establish that the period of stick-slip increases and the change of the stick-slip size was not large when the abrasive particle size was relatively small.

• Nuclear Engineering and Technology
• /
• 제54권6호
• /
• pp.2204-2212
• /
• 2022

Recently, the demand for alpha imaging detectors for quantifying the distributions of alpha particles has increased in various fields. This study aims to reconstruct a high-resolution image from an alpha imaging detector by applying a super-spatial resolution method combined with the maximum-likelihood expectation maximization (MLEM) algorithm. To perform the super-spatial resolution method, several images are acquired while slightly moving the detector to predefined positions. Then, a forward model for imaging is established by the system matrix containing the mechanical shifts, subsampling, and measured point-spread function of the imaging system. Using the measured images and system matrix, the MLEM algorithm is implemented, which converges towards a high-resolution image. We evaluated the performance of the proposed method through the Monte Carlo simulations and phantom experiments. The results showed that the super-spatial resolution method was successfully applied to the alpha imaging detector. The spatial resolution of the resultant image was improved by approximately 12% using four images. Overall, the study's outcomes demonstrate the feasibility of the super-spatial resolution method for the alpha imaging detector. Possible applications of the proposed method include high-resolution imaging for alpha particles of in vitro sliced tissue and pre-clinical biologic assessments for targeted alpha therapy.

• Water Engineering Research
• /
• 제3권2호
• /
• pp.75-84
• /
• 2002

An Euler-Lagrange two-phase flow model is presented fur simulation sheet-flow transport under wave and current. The flow is computed by solving the Reynolds Averaged Navier-Stokes equation in conjunction with the k-$\varepsilon$ turbulence model for turbulence closure. The sediment transport is introduced as a motion of granular media under the action of unsteady flow from the Lagragian point of view. In other word, motion of every single particle is numerically traced with Movable Bed Simulator (MBS) code based on the Distinct Element Method (DEM), in which the frequent interparticle collision of the moving particles during the sheet-flow transport is sophisticatedly taken into account. The particle diameter effect on time-dependent developing process of sheet-flow transport is investigated, by using three different diameter sizes of sediment. The influence of an imposed current on oscillatory sheet-flow transport is also investigated. It is concluded that the sediment transport rate increases due to the relaxation process related to the time-lag between flow velocity and sediment motion.

• Nuclear Engineering and Technology
• /
• 제53권3호
• /
• pp.833-841
• /
• 2021

During the postulated severe accident of nuclear reactor, eutectic reaction leads to low-temperature melting of fuel cladding and early failure of core structure. In order to model eutectic melting with the moving particle semi-implicit (MPS) method, the eutectic reaction model is developed to simulate the eutectic reaction phenomenon. The coupling of mass diffusion and phase diagram is applied to calculate the eutectic reaction with the uniform temperature. A heat transfer formula is proposed based on the phase diagram to handle the heat release or absorption during the process of eutectic reaction, and it can combine with mass diffusion and phase diagram to describe the eutectic reaction with temperature variation. The heat transfer formula is verified by the one-dimensional melting simulations and the predicted interface position agrees well with the theoretical solution. In order to verify the eutectic reaction models, the eutectic reaction of uranium and iron in two semi-infinite domains is simulated, and the profile of solid thickness decrease over time follows the parabolic law. The modified MPS method is applied to calculate Transient Reactor Test Facility (TREAT) experiment, the penetration rate in the simulations are agreeable with the experiment results. In addition, a hypothetical case based on the TREAT experiment is also conducted to validate the eutectic reaction with temperature variation, the results present continuity with the simulations of TREAT experiment. Thus the improved method is proved to be capable of simulating the eutectic reaction in the severe accident.