• 대한기계학회논문집B
• /
• 제37권6호
• /
• pp.591-597
• /
• 2013

대부분의 비뉴턴 유체 모델의 경우 격자볼쯔만 법을 이용한 3차원 유동해석에 효율적으로 적용하기 어려운 문제점을 가지고 있다. 이 연구에서는 이를 해결하기 위해 개발된 격자볼쯔만법 전용 비뉴턴 유체 모델인 Hydro-Kinetic 모델을 3차원 격자볼쯔만법 해석에 적용할 수 있도록 수정하고 이를 3차원 급축소 및 급확대 유동에 적용하여 얻은 결과를 통해서 강한 전단유동장에서 HK모델 유체의 거동을 분석하였다. HK모델은 변형률과 완화시간과의 관계를 나타낸 모델로 강한 전단유동장에서 국소적으로 큰 변형률이 발생하는 경우 완화시간이 급격하게 감소하여 해석이 불안정해지는 경향이 있어 격자볼쯔만법 해석안정성을 확보하기 위해서는 HK모델의 인수 ${\Gamma}$와 완화시간을 일정한 구간에서 변화하도록 제한 할 필요가 있다.

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

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

• 대한기계학회논문집B
• /
• 제40권11호
• /
• pp.689-695
• /
• 2016

연료전지의 성능에 가장 큰 영향을 주는 요소 중에 하나가 가스확산층과 촉매층에서 물의 거동이다. 따라서 가스확산층의 특성에 따른 유체의 거동의 변화를 이해하는 것은 연료전지의 성능개선과 가스확산층의 설계를 위한 필수적인 요소이다. 이 연구에서는 가스확산층의 설계요소인 기공도, 굴곡도와 유효확산계수를 수치적으로 계산할 수 있는 방법을 제안한다. 제안한 방법의 검증을 위하여 지름이 일정한 구형입자를 이용하여 기공도가 다른 다공체를 만들고 구형입자에 Bounceback 조건을 적용한 격자 볼쯔만법 유동해석을 수행하였다. 다공체 내의 유동효과를 나타내는 투과도는 다공체에 의한 압력강하와 평균유속으로 계산하고, 질량이 없는 입자의 평균 다공체 통과 거리로부터 계산한 굴곡도와 기공도를 이용하여 계산한 유효확산계수를 Neale의 이론식과 비교하여 정확하게 일치하는 것을 확인하였다. 이 방법은 실제 다공체의 이미지를 이용한 계산에도 수정없이 이용할 수 있어 연료전지의 성능향상과 설계를 위한 가스확산층의 특성분석에 활용될 수 있다.

• 대한기계학회논문집B
• /
• 제32권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.

• 한국전산유체공학회지
• /
• 제15권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.

• 한국전산유체공학회지
• /
• 제16권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.

• 한국전산유체공학회지
• /
• 제11권3호
• /
• 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.

• 대한기계학회논문집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.

• 한국전산유체공학회지
• /
• 제19권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.