• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.18-27
• /
• 2002

In order to improve the accuracy of the solution near the boundary in an analysis of viscous flow around an arbitrary boundary which move and be deformed using an Eulerian concept, a level-set based grid deformation method is introduced to concentrate grid points near the boundary. This paper presents a new monitor function which can easily control the level of the concentration of grid points along the boundary. Computations for steady flow around a semi-circular cylinder mounted on the bottom of the flow domain were carried out to check the improvement of the solution using the adaptive grid system with an immersed boundary method. The present numerical results show a good agreement with the solutions obtained by a body fitted grid system and more accurate solutions than those computed with non-adaptive grid system. For the validation of mechanical usefulness of the present method, an expanded analysis of flow around multi-body fixed in the flow domain was carried out. Finally, the present moving adaptive grid method was applied to a two-dimensional bubble rise problem. The computed results show well adapted grid points around the boundary of the bubble at every time and a good agreement with the result calculated by fixed grid system.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.11 s.242
• /
• pp.1199-1208
• /
• 2005

The flow fields including velocities, turbulence intensities, Reynolds shear stress and turbulent kinetic energy were investigated using particle image velocimetry(PIV) to study the flow characteristics around two square cylinders in a tandem arrangement. The experiments were carried out in the range of the spacing from 1.0 to 4.0 widths of cylinder, Reynolds number of 5.3$\times$10$^{3}$ and 1.6$\times$10$^{4}$ respectively. Discontinuous jumping at the drag coefficient variation was found for two cylinders simultaneously when the spacing between two cylinders is varied. This phenomenon is attributed to a sudden change of the flow pattern which depends on the reattachment of the shear layer separated from the upstream cylinder. Near such a critical spacing, the changes of the flow fields as well as the effect of Reynolds number were studied in detail.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.7
• /
• pp.629-636
• /
• 2013

Pneumatic levitation is based upon Bernoulli's principle. However, this method is known to require a large gas flow rate that can lead to an increase in the cost of products. In this case, the gas flow rate should be increased, and the compressible effects of the gas may be of practical importance. In the present study, a computational fluid dynamics method has been used to obtain insights into Bernoulli levitation flows. Three-dimensional compressible Navier-Stokes equations in combination with the SST k-${\omega}$ turbulence model were solved using a fully implicit finite volume scheme. The gas flow rate, workpiece diameter,and clearance gap between the workpiece and the circular cylinder were varied to investigate the flow characteristics inside. It is known that there is an optimal clearance gap for the lifting force and that increasing the supply gas flow rate results in a larger lifting force.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.3
• /
• pp.169-176
• /
• 2019

In the present study, a CFD/DSMC hybrid method performed by a coupled analysis between the CFD method and the DSMC method was developed to obtain the flow information on the rarefied gas flows effectively. Flow simulations around the high speed vehicles on the transition flow regimes were conducted by using the developed method. The FRESH-FX vehicle made of cone and cylinder shapes was considered for the simulations. The results of the hybrid method were compared with the results of the pure CFD and the pure DSMC method to confirm the reliability and efficiency of the hybrid method. It was found that the gradient and the intensity of the shock waves were weakened due to the relatively low density on the transition flow regime. It was confirmed that the results of the hybrid analysis were different to those of the pure CFD analysis and almost identical to those of the pure DSMC analysis. In addition, the computational time of the hybrid method was reduced than that of the pure DSMC method. As a result, it was obtained that the validity and the efficiency of the CFD/DSMC hybrid method.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.3
• /
• pp.213-221
• /
• 2012

The present study numerically investigates the interaction between a free-surface and flow around a circular cylinder over a moving wall. In order to simulate the flow past the circular cylinder over a moving wall near a free-surface, this study has adopted the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. Numerical simulation is performed for a Reynolds numbers of 100 in the range of $0.25{\leq}g/D{\leq}2.00$ and $0.5{\leq}h/D{\leq}2.00$, where g/D and h/D are the gaps between the cylinder and a moving wall and the cylinder and a free-surface normalized by cylinder diameter D, respectively. According to g/D and h/D, the vortex structures have been classified into three patterns of the two-row, one-row, steady elongation. In general, both of g/D and h/D have the large values which mean the cylinder is far away from the wall and the free-surface, two-row vortex structure forms in the wake. When g/D decreases, the two-row vortex structure gradually transfers into the one-row vortex structure. When the g/D reveals the critical value below which the flow becomes steady state, resulting in the steady elongation vortex.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2637-2649
• /
• 1996

Flow characteristics and aerodynamic forces acting on an elliptic cylinder placed in a plane boundary layer were investigated experimentally. Four cylinder models with axis ratio(major axis to minor axis, AR=A/B) of 1, 2, 3, and 4 having the same equivalent diameter were used in this experiment. The Reynolds number based on the equivalent diameter $D_e$(=20mm) was 13,000. In the case of circular cylinder, regular vortex shedding occurs for the cylinder gaps larger than G/B=0.3 and is not almost related to the boundary layer thickness. But, for the elliptic cylinders, the vortex shedding frequency is increased with increasing the gap ratio (G/B) and the axis ratio (AR) of elliptic cylinders. The maximum drag coefficient acting on a circular cylinder is mainly affected by the boundary layer thickness. But, the elliptic cylinders(AR$\geq$2), except for the smaller gap G/B<0.2, show a nearly constant drag coefficient which is much smaller than that of a circular cylinder. The base pressure on the flat plate decreases with increasing the axis ratio(AR) of the elliptic cylinder. In the case of a circular cylinder, the base pressure has the minimum value at the gap ratio G/B=0.4, but it occurs at G/D=2 for elliptic cylinders. The mean velocity of the cylinder wake is quickly recovered at a small cylinder height ratio(H/$\delta$), but the turbulent intensity is rapidly recovered at a large cylinder height ratio(H/$\delta$). The effective wake region in the plane boundary layer is shrinkaged with increasing the axis ratio(AR) of elliptic cylinder. And the drag coefficient and streamwise turbulent intensity of the elliptic cylinder with AR=4 are less than half of those for the circular cylinder(AR=1).

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.368-378
• /
• 2005

The effects of rotary oscillation on the unsteady laminar flow past a circular cylinder. are numerically investigated in the present study. The numerical solutions for the 20 Wavier-Stokes equation are obtained using a finite volume method Tn the framework of an overlapping grid system. The vortex formation behind a circular cylinder and the hydrodynamics of wake flows for different rotary oscillation conditions are analyzed from the results of numerical simulation. The lock-on region is defined as the region that the natural shedding frequency due to the Karmann Vortex shedding and the forcing frequency due to the forced oscillating a cylinder are nearly same, and the quasi-periodic states are observed around that region. At the intersection between lock-on and non-lock-on region the shedding frequency is bifurcated. After the bifurcation, one frequency fellows the forcing frequency($S_f$) and the other returns to the natural shedding frequency($St_0$). in the quasi-periodic states, the variation of magnitudes and relevant phase changes of $C_L$ with forcing phase are examined.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.1
• /
• pp.1-9
• /
• 2009

Three-dimensional characteristics of fluid flow and heat transfer around a wavy circular cylinder having sinusoidal variation in cross sectional area along the spanwise direction are numerically investigated using the immersed boundary method. The three different wavelengths of ${\pi}4$, ${\pi}3$ and ${\pi}2$ at the fixed wavy amplitude of 0.1 have been considered to investigate the effects of waviness especially on the forced convection heat transfer around a wavy cylinder when the Reynolds and Prandtl numbers are 300 and 0.71, respectively. The present computational results for a wavy cylinder are compared with those for a smooth cylinder. The time- and total surface-averaged Nusselt number for a wavy cylinder with ${\lambda}={\pi}/2$ is larger than that for a smooth cylinder, whereas that with ${\lambda}={\pi}/4$ and ${\pi}/3$ is smaller than that for a smooth cylinder. However, because the surface area exposed to heat transfer for a wavy cylinder is larger than that for a smooth cylinder, the total heat transfer rate for a wavy cylinder with different wavelengths of ${\lambda}={\pi}/4$, ${\pi}/3$ and ${\pi}/2$ is larger than that for a smooth cylinder.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.17 no.6
• /
• pp.861-867
• /
• 2014

In the present study, the flow around a circular cylinder at $Re=3.6{\time}10^6$ is numerically simulated using URANS approach. The objective of this study is to evaluate the turbulence models(Realizable k-${\varepsilon}$, RNG k-${\varepsilon}$) through the prediction of the unsteady flow characteristics around the cylinder. The time-averaged drag coefficients and vortex shedding phenomenon in the wake region are compared to available experimental data and other numerical results. The simulation with Realizable k-${\varepsilon}$ model is found to be more dissipative due to large eddy viscosity predicted in the wake region while the simulation with RNG k-${\varepsilon}$ model predicts a complex vortex shedding phenomenon with more coherent structures realistically.