• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.6
• /
• pp.591-597
• /
• 2013

Because most existing non-Newtonian models are not suitable for application to the lattice Boltzmann method, theoretical and numerical studies in this regard remain challenging. In this study, the hydrokinetic (HK) model was modified and applied to a 3D sudden contraction-expansion channel flow, and the characteristics of the HK model flow were evaluated to generate non-trivial predictions in three-dimensional strong shear flows. The HK model is very efficient for application to the lattice Boltzmann method because it utilizes the shear rate and relaxation time. However, the simulation would be unstable in a high shear flow field because the local relaxation time sharply decreases with an increase in the shear rate in a strong shear flow field. In the HK model, it may become necessary to truncate the relaxation time and non-dimensional parameter to obtain stable numerical results.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.11
• /
• pp.689-695
• /
• 2016

The performance of proton exchange membrane fuel cells (PEMFC) is strongly related to the water flow and accumulation in the gas diffusion layer (GDL) and catalyst layer. Understanding the behavior of fluid from the characteristics of the media is crucial for the improvement of the performance and design of the GDL. In this paper, a numerical method is proposed to calculate the design parameters of the GDL, i.e., permeability, tortuosity, and effective diffusivity. The fluid flow in a channel filled with randomly packed hard spheres is simulated to validate the method. The flow simulation was performed by lattice Boltzmann method with bounce back condition for the solid volume fraction in the porous media, with different values of porosities. Permeability, which affects the flow, was calculated from the average pressure drop and the velocity in the porous media. Tortuosity, calculated by the ratio the average path length of the randomly injected massless particles to the thickness of the porous media, and the resultant effective diffusivity were in good agreement with the theoretical model. The suggested method can be used to calculate the parameters of real GDL accurately without any modification.

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

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

• Journal of computational fluids engineering
• /
• v.11 no.3 s.34
• /
• pp.28-36
• /
• 2006

NSBT(No Slip Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the solid wall of obstacle in a flow field. In our research, NSBT was used to perform LBM simulation of a flow over a circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of NSBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a circular and square cylinder in a channel and the results were compared against the solution of Navier-Stokes equation. The simulations were performed in a moderate range of Reynolds number at each cylinder position to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is in the range of 200$\sim$250. For the gap parameter $\gamma$ = 2 cases at Re > 240, the vortex shedding were symmetric and it resembled the Karmann vortex. As the cylinder approached to one wall, the vorticity significantly reduced in length while the vorticity on the other side elongated and the vorticity combined with the wall boundary-layer vorticity. The resultant $C_d$ by LBM concurred with the results of DNS simulation performed by previous researchers.

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

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.61-67
• /
• 2014

Since the most of the existing non-Newtonian models are not adequate to apply to the lattmce Boltzmann method, it is a challenging task from both the theoretical and the numerical points of view. In this research the hydro-kinetic model was modified and applied to the 3-D moving sphere in the circular channel flow and the characteristics of the shear thinning effect by the HK-model was evaluated and the condition of ${\Gamma}$ in the model was suggested for the stable simulation to generate non-trivial prediction in three dimension strong shear flows. On the wall boundaries of circular channel the curved wall surface treatment with constant velocity condition was applied and the bounceback condition was applied on the sphere wall to simulate the relative motion of the sphere. The condition is adequate at the less blockage than 0.7 but It may need to apply a multi-scale concept of grid refinement at the narrow flow region. to obtain the stable numerical results.