• 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.

• Korea-Australia Rheology Journal
• /
• 제18권4호
• /
• pp.171-181
• /
• 2006

We present a short review for authors' previous work on direct numerical simulations for inertialess hard particle suspensions formulated either with a Newtonian fluid or with viscoelastic polymeric fluids to understand the microstructural evolution and the bulk material behavior. We employ two well-defined bi-periodic domain concepts such that a single cell problem with a small number of particles may represent a large number of repeated structures: one is the sliding bi-periodic frame for simple shear flow and the other is the extensional bi-periodic frame for planar elongational flow. For implicit treatment of hydrodynamic interaction between particle and fluid, we use the finite-element/fictitious-domain method similar to the distributed Lagrangian multiplier (DLM) method together with the rigid ring description. The bi-periodic boundary conditions can be effectively incorportated as constraint equations and implemented by Lagrangian multipliers. The bulk stress can be evaluated by simple boundary integrals of stresslets on the particle boundary in such formulations. Some 2-D example results are presented to show effects of the solid fraction and the particle configuration on the shear and elongational viscosity along with the micro-structural evolution for both particles and fluid. Effects of the fluid elasticity has been also presented.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제1권1호
• /
• pp.20-28
• /
• 2009

In thus paper we validate a numerical model for wave-structure interaction by comparing numerical results with laboratory data. The numerical model is based on the Navier-Stokes (N-S) equations for an incompressible fluid. The N-S equations are solved by a two-step projection finite volume scheme and the free surface displacements are tracked by the volume of fluid (VOF) method The numerical model is used to simulate solitary waves and their interaction with a group of slender vertical piles. Numerical results are compared with the laboratory data and very good agreement is observed for the time history of free surface displacement, fluid particle velocity and wave force. The agreement for dynamic pressure on the cylinder is less satisfactory, which is primarily caused by instrument errors.

• 한국전산유체공학회지
• /
• 제13권1호
• /
• pp.28-34
• /
• 2008

A particle method recognized as one of the gridless methods has been developed to investigate the nonlinear free-surface motions interacting to the structures. The method is more feasible and effective than convectional grid methods for solving the non-linear free-surface motion with complicated boundary shapes. The right-handed side of the governing equations for incompressible fluid, which includes gradient, viscous and external force terms, can be replaced by the particle interaction models. In the present study, the developed method is applied to the dam-broken problem on dried- and wet-floor and its adequacy will be discussed by the comparison with the experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2007년도 추계 학술대회논문집
• /
• pp.258-263
• /
• 2007

A particle method recognized as one of the gridless methods has been developed to investigate the nonlinear free-surface motions interacting to the structures. The method is more feasible and effective than convectional grid methods for solving the non-linear free-surface motion with complicated boundary shapes. The right-handed side of the governing equations for incompressible fluid, which includes gradient, viscous and external force terms, can be replaced by the particle interaction models. In the present study, the developed method is applied to the dam-broken problem on dried- and wet-floor and its adequacy will be discussed by the comparison with the experimental results.

• 항공우주시스템공학회지
• /
• 제12권5호
• /
• pp.24-31
• /
• 2018

항공기의 블레이드 안테나에 대해 더미와 실 안테나의 조류충돌 해석과 시험을 수행하였다. 해석에서 조류는 SPH(Smooth Particle Hydrodynamics) 방법을 이용하여 모델링하였으며, 유체-구조 연성해석 (FSI, Fluid-Structure Interaction) 기법으로 조류와 안테나, 기체 체결부의 거동을 시뮬레이션 하였다. 실제 조류를 사용한 시험을 수행하여 안테나와 동체 사이의 체결부 손상 및 이탈여부를 확인하였으며, 항공기 기체의 구조건전성과 해석 및 시험 결과 사이의 상관성이 있음을 입증하였다.

• 한국컴퓨터정보학회논문지
• /
• 제29권6호
• /
• pp.23-29
• /
• 2024

본 논문에서는 옷감과 입자 기반 유체 해법인 SPH(Smoothed particle hydrodynamics)를 이용한 액체 간의 상호작용으로 표현되는 다공성 흐름(Pore flow), 흡수, 방출 그리고 확산 효과를 GPU 기반으로 빠르게 표현할 수 있는 방법을 제안한다: 1) 옷감-액체의 상호작용에 의해 표현되는 다양한 물리적 효과를 GPU 기반으로 표현할 수 있는 통합형 프레임워크, 2) SPH 기반으로 노드의 포화도를 효율적으로 계산하고 이를 주변 Porous 입자들로 전달하는 방법, 3) 유체 흡수 및 방출 방향을 안정적으로 계산하기 위해 다르시 법칙(Darcy's law)을 기반으로 안정성을 개선시키는 방법, 4) Porous 입자들로 흡수되는 과정에서 유체의 흐름 방향에 따라 흡수되는 양을 조절하는 방법, 마지막으로 5) SPH 입자의 최대 질량이 넘지 않도록 방출할 수 있는 방법을 제시한다. 제안하는 방식의 가장 큰 장점은 모든 연산이 GPU에서 계산되고 동작하기 때문에 빠르게 옷감과 유체의 상호작용으로 표현되는 다공성 재질, 다공성 흐름, 흡수, 반사, 확산 등을 모델링할 수 있다.

• 한국유체기계학회 논문집
• /
• 제15권3호
• /
• pp.12-18
• /
• 2012

The Lattice Boltzmann Method has developed for solving the Boltzmann equation in Cartesian domains containing immersed boundaries of arbitrary geometrical complexity moving with prescribed kinematics. When a immersed boundaries are sweeping the fixed fluid node, refilling the node information in a vicinity of fluid nodes is one of the important issues in Lattice Boltzmann Method. In this study, we propose a simple refill algorithm for the particle distribution function based on a proper velocity, density and strain rate to enhance accuracy and stability of the method. The refill scheme based on a asymptotic analysis of LBGK model has improved accuracy than interpolation schemes. The proposed scheme in this study is validated by the simulations of an impulsively started rotating circular cylinder to investigate adaptability for fluid-structure interaction (FSI) problem. This refill scheme has improved stability and accuracy especially at high Reynolds number region.

• Korea-Australia Rheology Journal
• /
• 제19권4호
• /
• pp.233-242
• /
• 2007

We performed a numerical investigation to find out the optimal choice of the spatial discretization in the distributed-Lagrangian-multiplier/fictitious-domain (DLM/FD) method for the solid/fluid interaction problem. The elastic solid bar attached on the bottom in a pressure-driven channel flow of a Newtonian fluid was selected as a model problem. Our formulation is based on the scheme of Yu (2005) for the interaction between flexible bodies and fluid. A fixed regular rectangular discretization was applied for the description of solid and fluid domain by using the fictitious domain concept. The hydrodynamic interaction between solid and fluid was treated implicitly by the distributed Lagrangian multiplier method. Considering a simplified problem of the Stokes flow and the linearized elasticity, two numerical factors were investigated to clarify their effects and to find the optimum condition: the distribution of Lagrangian multipliers and the solid/fluid interfacial condition. The robustness of this method was verified through the mesh convergence and a pseudo-time step test. We found that the fluid stress in a fictitious solid domain can be neglected and that the Lagrangian multipliers are better to be applied on the entire solid domain. These results will be used to extend our study to systems of elastic particle in the Stokes flow, and of particles in the viscoelastic fluid.

• Nuclear Engineering and Technology
• /
• 제48권1호
• /
• pp.246-258
• /
• 2016

The shield building of AP1000 was designed to protect the steel containment vessel of the nuclear reactor. Therefore, the safety and integrity must be ensured during the plant life in any conditions such as an earthquake. The aim of this paper is to study the effect of water in the water tank on the response of the AP1000 shield building when subjected to three-dimensional seismic ground acceleration. The smoothed particle hydrodynamics method (SPH) and finite element method (FEM) coupling method is used to numerically simulate the fluid and structure interaction (FSI) between water in the water tank and the AP1000 shield building. Then the grid convergence of FEM and SPH for the AP1000 shield building is analyzed. Next the modal analysis of the AP1000 shield building with various water levels (WLs) in the water tank is taken. Meanwhile, the pressure due to sloshing and oscillation of the water in the gravity drain water tank is studied. The influences of the height of water in the water tank on the time history of acceleration of the AP1000 shield building are discussed, as well as the distributions of amplification, acceleration, displacement, and stresses of the AP1000 shield building. Research on the relationship between the WLs in the water tank and the response spectrums of the structure are also taken. The results show that the high WL in the water tank can limit the vibration of the AP1000 shield building and can more efficiently dissipate the kinetic energy of the AP1000 shield building by fluid-structure interaction.