• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.7
• /
• pp.441-448
• /
• 2016

The thrust force created by a propeller depends on the incoming flow velocity and the rotational velocity of the propeller. The performance of the propeller can be described by dimensionless variables, advanced ratio, thrust coefficient, and power coefficient. This study included the application of the immersed boundary lattice Boltzmann method (IBLBM) with the stereo lithography (STL) file of the rotating object for performance analysis. The immersed boundary method included the addition of the external force term to the LB equation defined by the velocity difference between the lattice points of the propeller and the grid points in the domain. The flow by rotating a 4-blade propeller was simulated with various Reynolds numbers (Re) (including 100, 500 and 1000), with advanced ratios in the range of 0.2~1.4 to verify the suggested method. The typical tendency of the thrust efficiency of the propeller was obtained from the simulation results of different advanced ratios. It was also necessary to keep the maximum mesh size ratio of the propeller surface to a grid size below 3. Additionally, a sufficient length of the downstream region in the domain was maintained to ensure the numerical stability of the higher Re and advanced ratio flow.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.21-27
• /
• 2010

Immersed boundary lattice Boltzmann method (IBLBM) has been applied to simulate a turbulent flow over circular cylinder in a flow field effectively. Although IBLBM is very effective method to simulate the flow over a complex shape of obstacle in the flow field regardless of the constructed grids in the calculation domain, the results, however, become numerically unstable in high reynolds number flow. The most effective suggestion to archive the numerical stability in high Reynolds number flow is applying the multiple relaxation time (MRT) model instead of single relaxation time(SRT) model in the collision term of lattice Boltzmann equation. In the research MRT model for IBLBM was introduced and comparing the numerical results obtained by applying SRT and MRT. The hydraulic characteristic of cylinder in a flow field between two parallel plate at the range of $Re{\leqq}2000$represented and it is also compared the drag and lifting coefficients of the cylinder calculated by IBLBM with SRT and MRT model.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.24-29
• /
• 2011

Immersed boundary lattice Boltzmann method has been applied to analyze the characteristics of the self-propelled fish motion swimming robot. The airfoil NACA0012 with caudal fin stroke model was considered to examine the characteristics. The foil in steady forward motion and a combination of steady-state harmonic deformation produces thrust through the formation of a flow downstream from the trailing edge. The harmonic motion of the foil causes unsteady shedding of vorticity from the trailing edge, while forming the vortices at the leading edge as well. The resultant thrust is developed by the pressure difference formed on the upper and lower surface of the airfoil. and the time averaged thrust coefficient increases as Re increase in the region of $Re{\leqq}700$. The suggested numerical method is suitable to develop the fish-motion model to control the swimming robot, however It would need to extend in 3D analysis to examine the higher Re and to determine the more detail mechanism of thrust production.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.12
• /
• pp.924-930
• /
• 2008

Immersed boundary method (IBM) is the most effective method to overcome the disadvantage of LBM (Lattice Boltzmann Method) related to the limitation of the grid shape. IBM also make LBM possible to simulate flow over complex shape of obstacle without any treatment on the curved boundary. In the research, IBLBM was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of IBLBM on the moving obstacle near the wall, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of moving cylinder in the channel using IBLBM. The simulations were performed in a moderate range of Reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag and lift coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical Reynolds number for vortex shedding is Re=50 and the result is the same as the case of fixed cylinder. As the cylinder approaching to a wall (${\gamma}<2.5$), the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. When the cylinder is very closed to the wall, ${\gamma}<0.6$, the cylinder acts like blockage to block the flow between the cylinder and wall so that the vortex developed on the upper cylinder elongated and time averaged lifting and drag coefficients abruptly increase.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.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.