• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.13-20
• /
• 2024

In this paper, a parallel analysis algorithm for Smoothed Particle Hydrodynamics (SPH), one of the numerical methods for fluidic materials, is introduced. SPH, which is a meshless method, can represent the behavior of a continuum using a particle-based approach, but it demands substantial computational resources. Therefore, parallel analysis algorithms are essential for SPH simulations. The domain decomposition algorithm, which divides the computational domain into partitions to be independently analyzed, is the most representative method among parallel analysis algorithms. In Discrete Element Method (DEM) and Molecular Dynamics (MD), the Cartesian coordinate-based domain decomposition method is popularly used because it offers advantages in quickly and conveniently accessing particle positions. However, in SPH, it is important to share particle information among partitioned domains because SPH particles are defined based on information from nearby particles within the smoothing length. Additionally, maintaining CPU load balance is crucial. In this study, a highly parallel efficient algorithm is proposed to dynamically minimize the size of orthogonal domain partitions to prevent excess CPU utilization. The efficiency of the proposed method was validated through numerical analysis models. The parallel efficiency of the proposed method is evaluated for up to 30 CPUs for fluidic models, achieving 90% parallel efficiency for up to 28 physical cores.

• Korean Chemical Engineering Research
• /
• v.60 no.4
• /
• pp.535-543
• /
• 2022

A CPFD (Computational particle fluid dynamics) model of solar fluidized bed receiver of silicon carbide (SiC: average d_p=123 ㎛) particles was established, and the model was verified by comparing the simulation and experimental results to analyze the effect of particle behavior on the performance of the receiver. The relationship between the heat-absorbing performance and the particles behavior in the receiver was analyzed by simulating their behavior near bed surface, which is difficult to access experimentally. The CPFD simulation results showed good agreement with the experimental values on the solids holdup and its standard deviation under experimental condition in bed and freeboard regions. The local solid holdups near the bed surface, where particles primarily absorb solar heat energy and transfer it to the inside of the bed, showed a non-uniform distribution with a relatively low value at the center related with the bubble behavior in the bed. The local solid holdup increased the axial and radial non-uniformity in the freeboard region with the gas velocity, which explains well that the increase in the RSD (Relative standard deviation) of pressure drop across the freeboard region is responsible for the loss of solar energy reflected by the entrained particles in the particle receiver. The simulation results of local gas and particle velocities with gas velocity confirmed that the local particle behavior in the fluidized bed are closely related to the bubble behavior characterized by the properties of the Geldart B particles. The temperature difference of the fluidizing gas passing through the receiver per irradiance (∆T/I_DNI) was highly correlated with the RSD of the pressure drop across the bed surface and the freeboard regions. The CPFD simulation results can be used to improve the performance of the particle receiver through local particle behavior analysis.

• Journal of Biomedical Engineering Research
• /
• v.20 no.5
• /
• pp.601-608
• /
• 1999

동맥의 일부분의 팽창하는 동맥류는 높은 사망률을 야기하는 혈관계 질환이다. 동맥류의 발생 및 파열에는 동맥류 내부의 혈류의 유동에 의한 혈관벽 전단 응력 및 압력이 주용한 원인 중 하나로 의심되고 있다. 복부대동맥류 내부의 혈류 유동 특성을 밝히기 위해서 동맥류의 최대 확장부가 복부동맥의 1.5배, 2배인 유리 모델을 제작하였다. 정상류 상태에서 다양한 레이놀즈수에 대해서 속도 및 난동도를 입자영상속도계를 이용하여 측정하였다. 경계층 박리로 인한 재순환 부분이 끝나는 재부착점은 동맥류 최대 확장부 후부에서 발생하였으며, 이 위치는 레이놀즈수의 변화에 따라 바뀌었다. 축방향 속도의 난동은 최대 확장부 후부에서 크게 나타났으며, 이 위치에서 난동에 의한 부가적 응력이 크며 혈관벽 구조변화가 발생하리라 예측된다. 동맥류 내부의 압력분포는 수치해석에 의해 계산되었다. 동맥류 내부 압력은 크기가 증가함에 따라 커졌으며 압력은 동맥류 최대 확장부 후부에서 발생하는 재부착점에서 최대값을 나타내었다. 동맥류 최대확장부 후부는 압력이 최대값을 가지며, 전단력의 변화 및 난동이 큰 지역이므로 동맥류의 파열이 발생하기 쉬운 지역으로 예측된다.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.4
• /
• pp.8-14
• /
• 2007

A particle method recognized as one of gridless methods has been developed to investigate incompressible viscous flaw. The method is more feasible and effective than conventional grid methods for solving the flaw field with complicated boundary shapes or multiple bodies. The method is consists of particle interaction models representing pressure gradient, diffusion, incompressibility and the boundary conditions. In the present study, the models in case of various simulation condition were checked with the analytic solution, and applied to the two-dimensional Poiseuille flow in order to validate the developed method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.12
• /
• pp.1079-1089
• /
• 2013

The flow characteristics in a solid particle incinerator are investigated numerically for high burning rate of wastes. The studied incinerator employs both a swirl flow used in the furnace of powerplants and a design concept applied to a rocket combustor. As the first step, the non-reactive flow field is analyzed in the incinerator with primary and secondary injectors through which solid fuel and air are injected. The deflection angle of a primary injector, inclination angle of a secondary injector, and gap between the two types of injectors are selected as design parameters. The swirl number is adopted for evaluating the degree of swirl flow and estimated over wide ranges of three parameters. The swirl number increases with deflection angle, but it is affected little by inclination angle. Recirculation zones are formed near the injectors, and their size affects the swirl number. The swirl number decreases with the zonal size of recirculation. From the numerical results, the design points can be found with strong swirl flow.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.21 no.10
• /
• pp.144-149
• /
• 2020

Bubble reduction devices are intended to solve problems related to the quantitative supply of oil. Therefore, in this study, numerical analysis was conducted to verify the flow characteristics of bubble particles during the operation of a bubble reduction device. As a result of the basic analysis, the area where the rise and fall of bubbles were most active was found, and numerical calculations were performed focusing on the points. Before the numerical calculations, a non-dimensional derivation was performed to secure homogeneity among the variables. Based on the data obtained from non-dimension derivation, 25 variable conditions for each particle size and fluid velocity were set. Through separate calculations, the equation for bubble rise and fall was derived. By calculating the ratio of drag and buoyancy for each variable, if the drag force acting on the bubble was greater than buoyancy, the bubble falls, and bubbles are not reduced. If the buoyancy is larger than drag, the bubble rises, and the bubble is reduced. Through the analysis, the rise and fall of the bubble were confirmed, and the results were consistent with the separate numerical calculations.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.2
• /
• pp.110-116
• /
• 2005

This paper deals with interaction between two bubbles or particles and flow around them, visualized by a moving object flow image analyzer(MOFIA) consisting of a three-dimensional (3D) moving object image analyzer(MOIA) and two-dimensional particle image velocimetry(PIV). The experiments were carried out for rising bubbles or particles of various densities, sizes, and/or shapes in stagnant water in a vertical pipe. In the MOFIA employed, 3D-MOIA was used to measure particles or bubbles motion and PIV was used to measure fluid flow, The experimental results showed that the interaction was characterized by the shape, size and density of two particles or bubbles.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.05a
• /
• pp.173-176
• /
• 2009

The Arbitrary Lagrangian-Eulerian(ALE, in short) method is the new description of continum motion, which combines the advantages of the classical kinematical descriptions, i.e. Lagrangian and Eulerian description, while minimizing their respective drawbacks. In this paper, the ALE description is adapted to simulate fluid-structure interaction problems. An automatic re-mesh algorithm and a fluid-structure coupling process are included to analyze the interaction and moving motion during the 2-D axisymmetric solid rocket interior FSI phenomena simulation.

• Journal of the Korea Society of Computer and Information
• /
• v.18 no.1
• /
• pp.149-156
• /
• 2013

In this paper, Commercial program, STAR-CCM+, was used for computer simulation. And also Multi Inner Stage(MIS) cyclone which especially designed for the real experiments of particle removal efficiency. Under negative pressure condition of outlet, computer simulation was performed to predict the removal efficiency for $5{\mu}m$ and $10{\mu}m$ particles by using a turbulence model and lagrangian method. The simulation results are 55.7% and 64.1% for $5{\mu}m$ and $10{\mu}m$ particles, respectively. To compare the simulation results with the actual test of MIS cyclone, we generated the $SiO_2$ particles by heat reactor. Although removal efficiency of actual tests is 63~76% at different flow rate, the size of $SiO_2$ particles which confirmed by SEM(scanning electron microscope) and WAPS(wide range aerosol particle spectrometer) is too small(15~30nm) to compare each results. And so the alternative experiments were performed by using commercial alumina particles ($5{\mu}m$, $10{\mu}m$ and $20{\mu}m$). It was shown that the actual removal efficiency, 76~95%, from MIS cyclone is higher than simulated one.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.1
• /
• pp.73-79
• /
• 2010

Propellant tank pressurizing gas injector is used in the pressurization system of liquid propellant rocket to reduce incoming gas velocity and distribute the gas in the tank. Temperature distribution in the propellant tank ullage is varied according to the gas injector shape, and it has influence on the required pressurant gas and thermal phenomena in the tank. In this paper, diffuser type gas injector was studied to make the ullage have stratified temperature distribution. Injected gas flow at the outlet of prototype diffuser was visulized using particle image velocimetry method and it was compared with the results of calculation. Calculation was well agreed with measurement and was used as an inlet condition of propellant tank ullage calculation.