• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 12th Asian Congress of Fluid Mechanics
• /
• pp.69.3-69.3
• /
• 2008

• 한국해양공학회지
• /
• 제25권4호
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• 대한기계학회논문집B
• /
• 제40권7호
• /
• pp.441-448
• /
• 2016

프로펠러에 의한 추력은 유체의 유입 속도와 익의 회전속도에 의해 생성되며 그 성능을 전진비, 추력계수, 동력계수와 같은 무차원수로 나타내고 있다. 이 연구에서 회전체의 성능을 분석하기 위한 수치적 방법으로 STL형식의 회전체 형상을 인식할 수 있는 가상경계법을 적용한 격자볼쯔만법을 제안한다. 이 가상경계법으로 프로펠러의 회전에 의한 유동을 구현하기 위해서 프로펠러의 표면 격자점에서 속도와 유동장의 격자점에서 유속의 차를 이용하여 계산한 체적력을 볼쯔만방정식의 외력항으로 적용하게 된다. 제안한 방법을 검증하기 위하여 4개의 익을 가지고 있는 프로펠러를 이용해 레이놀즈수가 100, 500, 1000이고 전진비가 0.2~1.4일 때 유동해석을 수행하였으며 그 결과로 부터 전형적인 프로펠러의 성능특성을 얻을 수 있었다. 높은 레이놀즈수와 전진비를 갖는 유동에서 해석 안정성을 확보하기 위해서는 익의 표면에 구성한 최대 격자의 크기와 유동장에 구성한 격자 크기의 비가 3 이하로 유지해야 하며 충분히 긴 후류영역을 확보할 필요가 있다.

• Geomechanics and Engineering
• /
• 제15권6호
• /
• pp.1161-1171
• /
• 2018

The interaction between particles and fluid was investigated by IB-SEM numerical method which is a combination of combing the spectral/hp element method and the rigid immersed boundary method. The accuracy of this numerical method was verified based on the computed results with the traditional body-fitted mesh in numerical simulation of the flow through the cylinder. Then the governing equations of particles motion and contact in fluid are constructed. The movement of the particles and the interaction between the fluid and the particles are investigated. This method avoided the problem of low computational efficiency and error caused by the re-division of the grid when the solids moved. Finally, the movement simulation of multi particles in the fluid was carried out, which can provide a completely new numerical simulation method.

• 한국전산유체공학회지
• /
• 제20권2호
• /
• pp.103-108
• /
• 2015

Fluid-body interaction analysis of floating body with six degree-of-freedom motion is presented. In this study, three-dimensional incompressible Navier-Stokes equations are employed as a governing equation. The numerical method is based on a finite-volume approach on a cartesian grid together with a fractional-step method. To represent the body motion, the immersed boundary method for direct forcing is employed. In order to simulate the coupled six degree-of-freedom motion, Euler's equations based on rigid body dynamics are utilized. To represent the complex body shape, level-set based algorithm is utilized. In order to describe the free surface motion, the volume of fluid method utilizing the tangent of hyperbola for interface capturing scheme is employed. This study showed three different continuums(air, water and body) are simultaneously simulated by newly developed code. To demonstrate the applicability of the current approach, two different problems(dam-breaking with stationary obstacle and water entry) are simulated and all results are validated.

• 설비공학논문집
• /
• 제30권2호
• /
• pp.58-67
• /
• 2018

This study investigated the effect of variation in the angle of the elliptic cylinder as well as the presence of circular cylinder on natural convection inside a square enclosure. The Rayleigh number was varied between $10^3$ and $10^6$, and the Prandtl number was fixed to 0.7. In the present study, the angle of the elliptic cylinder was changed from $0^{\circ}$ to $90^{\circ}$, and the perimeter of the elliptic cylinder was same as that of the circular cylinder. The immersed boundary method was used to capture the virtual wall boundary of the inner cylinder. With the increasing angle of the elliptic cylinder, the surface-averaged Nusselt numbers on the cylinder and the enclosure increased. In the Rayleigh number range considered in the present study, the surface-averaged Nusselt number on the elliptic cylinder over = $45^{\circ}$ was higher than that of the circular cylinder. The effect of elliptic cylinder's angle on natural convection in the enclosure was analyzed according to the flow and thermal fields, and the distributions of the Nusselt number.

• 한국전산유체공학회지
• /
• 제12권2호
• /
• pp.37-45
• /
• 2007

The hybrid Cartesian/immersed boundary method is applied to simulate fluid-structure interaction of a two-dimensional orbiting flexible foil. The elastic deformation of the flexible foil is modelled based on the dynamic equation of a thin-plate. At each time step, the locations and velocities of the Lagrangian control points on the flexible foil are used to reconstruct the boundary conditions for the flow solver based on the hybrid staggered/non-staggered grid. To test the developed code, the flow fields around a flapping elliptical wing are calculated. The time history of the vertical force component and the evolution of the vorticity fields are compared with recent other computations and good agreement is achieved. For the orbiting flexible foil, the vorticity fields are compared with those of the case without the deformation. The combined effects of the angle of attack and the orbit on the deformation are investigated. The grid independency study is carried out for the computed time history of the deformation at the tip.

• 설비공학논문집
• /
• 제22권11호
• /
• pp.727-734
• /
• 2010

The present study numerically investigates the motion of a solid body suspended in the square enclosure with natural convection. A two-dimensional circular cylinder levitated thermally has been simulated by using thermal lattice Boltzmann method(TLBM) with the direct-forcing immersed boundary method. To deal with the ascending, falling or levitation of a circular cylinder in natural convection, the immersed boundary method is expanded and coupled with the TLBM. The circular cylinder is located at the bottom of a square enclosure with no restriction on the motion and freely migrates due to the Boussinesq approximation which is employed for the coupling between the flow and temperature fields. For different density ratio between the cylinder and the fluid, the motion characteristics of the circular cylinder for various Grashof numbers have been carried out. The Prandtl number is fixed as 0.7.

• 대한기계학회논문집B
• /
• 제33권2호
• /
• pp.107-115
• /
• 2009

A Computational study was carried out in order to investigate the aerodynamic characteristics of circular cylinder moving near the wavy wall at a low Reynolds number of 50. Lattice Boltzmann method was used to simulate the flow field and immersed boundary method was combined to represent the moving cylinder and wavy wall regardless of the constructed grid in the domain. The aerodynamics characteristics of the cylinder moving near the wavy wall were represented by the comparing the lifting coefficients with various altitudes (H/D) and wave length and amplitudes of wavy wall. It indicated that the twice of increasing-decreasing variations of lifting coefficient are obtained while the cylinder moves near the wavy wall. The first variation is obtained where the cylinder locates near the peak of the wavy wall. Another variation occurs when the distance to the wavy wall becomes longer after passing the peak. It was also classified that three different patterns of relation between the lifting and drag coefficient of the cylinder. However, the classification is limited to the case of the same order of altitude, amplitude and wave length of the wavy wall.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제9권6호
• /
• pp.598-612
• /
• 2017

The recently developed Harmonic Polynomial Cell (HPC) method has been proved to be a promising choice for solving potential-flow Boundary Value Problem (BVP). In this paper, a flux method is proposed to consistently deal with the Neumann boundary condition of the original HPC method and enhance the accuracy. Moreover, fixed mesh algorithm with free surface immersed is developed to improve the computational efficiency. Finally, a two dimensional (2D) multi-block strategy coupling boundary-fitted mesh and fixed mesh is proposed. It limits the computational costs and preserves the accuracy. A fully nonlinear 2D numerical wave tank is developed using the improved HPC method as a verification.