• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.392-404
• /
• 1997

The flow induced by the oscillatory motion of a solid body is important in a number of practical problems. As the solid boundary oscillates harmonically, there is steady streaming motion invoked by the Reynolds stresses, which could cause extensive migration of the fluid during a period of fluid motion. We here analyzed the flow in a square cavity with an oscillating top wall for the parameters which make the time derivatives and the convective terms equally important in the entire cavity flow. The full Navier-Stokes equations are solved by the second-order time accurate Momentum Coupling Method which is devised by the authors. The particular numerical scheme does not need subiteration at each time step which is usually a required process to calculate the incompressible Navier-Stokes equations. The effect of two parameters, the Reynolds number and the frequency parameter, on the oscillatory flow has been investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.27-33
• /
• 2009

In this paper, the flow past a rotary oscillating circular cylinder is simulated. The high-order and high-resolution numerical schemes with the characteristic boundary conditions are used for the compressible Navier-Stokes equation. The frequencies of rotating oscillation are $0.19\;{\leq}\;S_f\;{\leq}\;0.25$ for the maximum angular $\theta_{max}=10^{\circ}$ and $17^{\circ}$. The flow conditions are Mach number of 0.3 and Reynolds number of 1000. At Lock-on and Non-lock-on region which are defined by the relation between the vortex shedding frequency and the oscillating frequency, the drag and lift coefficient are analyzed.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.22 no.11
• /
• pp.767-772
• /
• 2010

The incompressible viscous flow past a sphere under forced oscillation is numerically investigated at a Reynolds number of 300. The immersed boundary method is used to handle the sphere oscillating vertically to the streamwise direction. There are two important variables to characterize the oscillating state of a sphere. One is an oscillating amplitude normalized by the sphere diameter is set as a fixed number of 0.2. Another is the frequency ratio which is defined by $f_e/f_o$, where fe and fo are the excited frequency and the natural frequency of vortex shedding for the stationary sphere. In this study, three different frequency ratios of 0.8, 1.0 and 1.2 are considered. The results show a periodic flow with hairpin vortices shedding from upper and lower positions as well as vortical legs obliquely extended by oscillating motion of sphere. The enveloping vortical structure experience rupture twice in one period of oscillation. As the frequency of oscillation is increased, the vortical legs are getting shorter and eventually the hairpin vortices are much closer to the adjacent one.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.1
• /
• pp.88-98
• /
• 1996

The present study is concerned with the flow friction and heat transfer characteristics of the combination of various regenerator materials, using the different Darcy number and porosity, which is filled uniformly and partially in a tube under oscillating flow condition. The poros medium is adopted as Brinkmann-Forschheimer extended Darcy model. Numerical results are obtained or the flow and temperature fields and described the effect of the combination of various regenerator materials and Womersley number on the pressure drop, the heat transfer and the regenerator efficiency. The results obtained indicate that not only heat transfer between the tube wall and oscillating flow but also the pressure drop at both ends of the regenerator are increased, while the regenerator efficiency is decreased in the increase of womersley number. It is also found that the friction factor is increased as Reynolds number is increased. The comparison between the combination of the various regenerator materials and the homogeneous regenerator material shows that the regenerator efficiency can be enhanced with the proper combination of various regenerator materials even though the averaged porosity of the regenerator is same.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.6
• /
• pp.421-428
• /
• 2008

In this paper, we investigated the flow of an incompressible viscous fluid past a sphere which is oscillated one-dimensionally over flow regimes including laminar flow at Reynolds number of 100, 200 and Strouhal number of up to 5000. In order to analyze flow and estimate critical Strouhal number, we introduce three-dimensional vortex element method. With this method, separation only appears in decreasing velocity region during the high Strouhal numbers. We find out that vorticity distribution around sphere is proportionl to the Strouhal number. And we can decide that low Strouhal number is below 100, high Strouhal number is above 500 from many results. Thus the critical Strouhal number(St) effected to the flow field is expected to be 100

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.387-390
• /
• 2006

Numerical study has been conducted in two dimensions about a NACA0012 airfoil with an oscillating microflap on the surface. We show that this microflap is effective in controlling the unsteady stall at high angles of attack. We solve the compressible Navier-Stokes equations for the Reynolds numbers with an extensible chimera grid fitted to the oscillatory microflap. For turbulent calculation, we adopt the SST $k-{\omega}$ model. We investigate the parametric effect of angle of attacks, Reynolds number, and the location where the microflap is installed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.5 s.248
• /
• pp.430-437
• /
• 2006

The mixing effect is studied by comparing rotating and oscillating stirrers in the micro channel. The cases of Re=10 to 80 with various stirring speeds are considered to analysis the effect of Re and stirrer speed for the mixing. Under Re=20, the oscillating stirrer represents better mixing rate than the rotating stirrer up to the critical stirrer speed which has a maximum efficiency. Over Re=30, the results of oscillating and rotating stirrer show that the faster the stirrer speed, the higher the mixing effect within the concerned stirrer speed range and the oscillating stirrer keeps the higher mixing rate. It was found that the mixing effect is a function which has an optimum of the Reynolds number and the stirrer speed. The D2Q9 Lattice Boltzmann Method is used due to the merits of calculation for the unsteady flow with moving boundary.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2294-2303
• /
• 1993

This paper investigated the characteristics of the turbulent incompressible flow past the orifice ring in an axi-symmetric pipe. The flow field was the turbulent pulsatile flow for Reynolds number of $2{\times}10^{5}$ which was defined based on the maximum velocity and the pipe diameter at the inlet, with oscillating frequence $(f_{os})=1/4{\pi}$ which was considered as quasi-steady state frequence. In the present investigation, finite analytic method was used to solve the governing equations in Navier Stokes and turbulent transport formulations. Particularly at high Reynolds number and low oscillation frequency, the effects of orifice ring on the flow were numerically investigated. The separation zone behind the orifice ring during the acceleration phase was found to be decreased. However, during the deceleration phase, the separation behind the orifice ring for pulsatile flow continuously grow to a size even larger than that in steady flow. The pressure drop in steady flow was found to be constant and always positive while for pulsatile flow the pressure drop change with time. And large turbulent kinetic energy, dissipation rate were found to be located in the region where the flow passes through the orifics ring. The maximum turbulent kinetic energy, generally occurs along the shear layer where the velocity gradient is large.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT 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 CMBT on the complex shape of the obstacle, 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 Navier-Stokes equation with deforming mesh technique. 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 coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.