• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.2
• /
• pp.261-268
• /
• 2002

Turbulent flow past a square cylinder confined in a channel is numerically investigated by Large Eddy Simulation(LES). The main objectives of this study are to verify the experimental results of Nakagawa et al.[Exp. in Fluids, Vol. 27, 3, pp. 284∼294, 1999] by LES and to obtain related flow information in detail. The LES results obtained are in excellent agreement with the experiment both qualitatively and quantitatively. The passive paticles numerically released into the flow field clearly show the barman vortex street. However, the vortices shed from the cylinder are significantly affected by the presence of the plates. Futhermore, periodic and alternating vortex-rollups are observed in the vicinity of the plates. The rolled-up vortex is convected downstream together with the corresponding Karman vortex forming a counter-rotating vortex pair. It is also revealed that the cylinder greatly enhances mixing process of the flow.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.627-632
• /
• 2000

Developing turbulent flows in a rotating 90 degree bend with square cross-section were measured by a hot-wire anemometer. The six orientation hot-wire technique was applied to measured the distributions of 3 mean velocities and 6 Reynolds stress components. Effects of Coriolis and centrifugal forces caused by the curvature and rotation of bend on the mean motion and turbulence structures were experimentally investigated Productive addition of Coriolis and centrifugal forces to the outward radial direction in the entrance region of bend increases the secondary flow intensity according to the rotational speeds. However, after 45 degree of bend, centrifugal force due to the rotation of bend may promote the break down of counter rotating vortex pair into multi-cellular pattern, thereby decreasing the production rate of turbulence energy and Reynolds stresses.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.577-589
• /
• 2005

The present study has numerically investigated two-dimensional laminar flow over a steadily rotating circular cylinder with a uniform planar shear, where the free-stream velocity varies linearly across the cylinder. It aims to find the combined effect of rotation and shear on the flow. Numerical simulations using the immersed boundary method are performed for the ranges of $-2.5{\le}\alpha{\le}2.5$ and $0{\le}K{\le}0.2$ at a fixed Reynolds number of Re=100, where a and K are respectively the dimensionless rotational speed and velocity gradient. Results show that the positive shear, with the upper side having the higher free-stream velocity than the lower one, favors the effect of the counter-clockwise rotation $(\alpha<0)$ but countervails that of the clockwise rotation $(\alpha>0)$. Accordingly, the absolute critical rotational speed, below which vortex shedding occurs, decreases with increasing K for $(\alpha>0)$, but increases for $\alpha>0$. The vortex shedding frequency increases with increasing \alpha (including the negative) and the variation becomes steeper with increasing K. The mean lift slightly decreases with increasing K regardless of the rotational direction. However, the mean drag and the amplitudes of the lift- and drag-fluctuations strongly depend on the direction. They all decrease with increasing K for $\alpha>0$, but increase for $\alpha<0$. Flow statistics as well as instantaneous flow folds are presented to identify the characteristics of the flow and then to understand the underlying mechanism.

• Journal of the Korean Society of Safety
• /
• v.27 no.5
• /
• pp.35-41
• /
• 2012

The Numerical simulation was performed on the flow field around the two-dimensional rectangular bluff body in order to complement the previous experimental results of the bluff body stabilized flames [1]. For both fuel ejection configurations against an oxidizer stream, the flame stability was affected mainly by vortex structure and mixing field near bluff body. FDS(Fire Dynamic Simulator) based on the LES(Large Eddy Simulation) was employed to clarify the isothermal mixing characteristic and wake flow pattern around bluff body. The air used atmosphere and the fuel used methane. The result of counter flow configuration shows that the flow field depends on air velocity but the mixing field is influenced on the fuel velocity. At low fuel velocity the fuel mole fraction is below the flammable limit and hence the mixing is insufficient to react. Therefore, as the result, the flame formed at low fuel velocity is characterized by non-premixed flames. For the flow field of co-flow configuration, flame stability was affected by fuel velocity as well as air velocity. the vortex generated by fuel stream has counter rotating direction against the air stream. Therefore, the momentum ratio between air and fuel stream was important to decide the flame blow out limit, which is result in the characteristic of the partially premixed reacting wake near extinction.

• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.33-46
• /
• 2011

Concave surface boundary-layer flows are subjected to centrifugal instability which results in the formation of streamwise counter-rotating vortices. Such boundary layer flows have been experimentally investigated on concave surfaces of 1 m and 2 m radius of curvature. In the experiments, to obtain uniform vortex wavelengths, thin perturbation wires placed upstream and perpendicular to the concave surface leading edge, were used to pre-set the wavelengths. Velocity contours were obtained from hot-wire anemometer velocity measurements. The most amplified vortex wavelengths can be pre-set by the spanwise spacing of the thin wires and the free-stream velocity. The velocity contours on the cross-sectional planes at several streamwise locations show the growth and breakdown of the vortices. Three different vortex growth regions can be identified. The occurrence of a secondary instability mode is also shown as mushroom-like structures as a consequence of the non-linear growth of the streamwise vortices. Wall shear stress measurements on concave surface of 1 m radius of curvature reveal that the spanwise-averaged wall shear stress increases well beyond the flat plate boundary layer values. By pre-setting much larger or much smaller vortex wavelength than the most amplified one, the splitting or merging of the streamwise vortices will respectively occur.

• Journal of the Korean Society of Visualization
• /
• v.15 no.1
• /
• pp.25-31
• /
• 2017

This study aims to visualize the Marangoni flow inside a droplet placed on heated hydrophobic surface and to measure its internal velocity field. The experimental result shows that the internal velocity increases with the increase of the plate temperature. In addition, the temperature difference induces the initial flow and drives the Marangoni circulation inside the droplet as soon as the evaporation starts (i.e. the thermal Marangoni flow). The fluorescence particles in the droplet trace two large-scale counter-rotating vortex pairs yielding the downwards flow along the vertical central axis. These vortex pairs gradually become small and move towards the contact line as time goes by, and this Marangoni flow sustains only for a half of the total evaporation time.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.274-277
• /
• 2012

Large eddy simulation is applied to the turbulent flow and heat transfer in straight cooling passages with varying aspect ratio. The turbulent statistics of the flow and thermal quantities are calculated and the characteristics of Nusselt number are investigated. To scrutinize near-wall streamwise vortices, a conditional sampling technique is adopted. Clockwise and counter-clockwise rotating streamwise vortices are sampled and the probability density function of the vortex circulation Reynolds number and wall Nusselt number are calculated.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.14-19
• /
• 2008

Effects of the Reynolds and Knudsen numbers on a micro-viscous pump are studied by using a Navier-Stokes code based on a finite volume method. The micro viscous pump consists of a circular rotor and a two-dimensional channel. The channel walls are treated by using a slip velocity model. The Reynolds number is studied in the range of $0.1{\sim}50$. The Knudsen number varies from 0.01 to 0.1. Numerical solutions show that the pump works efficiently when two counter rotating vortices formed on both sides of the rotor have the same size and intensity. As the Reynolds number increases, the size and intensity of the vortex on the inlet side of the pump decrease. It disappears when the Reynolds number is larger than about Re=20. The characteristics of the performance of the pump is shown to deteriorate, in terms of mean velocity and pressure rise, as the Reynolds number increases. The Knudsen number shows a different effect on the characteristics of the pump. As it increases, the mean velocity and pressure rise decrease but the characteristics of the vortex flow remains unchanged, unlike the effect of Reynolds number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.834-841
• /
• 2004

Experimental study on the three-dimensional topology of hairpin packet structures in turbulent boundary layers were carried out. Two different Reynolds number based on momentum thickness, Re$\sub$$\theta$/=514 and 934 were generated in a blowing type wind tunnel under the condition of zero pressure gradient. Simultaneous measurements of velocity fields at a wall-normal plane and wall-parallel plane by a plane PIV and a Stereo-PIV systems. The two Nd:Yag laser systems and three CCD cameras were synchronized to obtain instantaneous velocity fields at the same time. To avoid optical noise at the crossing line by the two laser light sheets, a new optical arrangement using polarization was applied. The obtained velocity fields show the existence of hairpin packet structure vividly and the idealized hairpin vortex signature is confirmed by experiment. Two counter-rotating vortex pair which reflects the cutting plane of hairpin legs are found both side of a strong streaky structure when the wall-normal plane cuts the hairpin head.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1265-1273
• /
• 1999

PIV(Particle Image Velocimetry) velocity field measurements inside the snout of a1/10 scale model of the Zn plating process were carried out at the strip speed $V_s=1.5m/s$. Aluminum powder particles ($1{\mu}m$) and atomized olive oil ($3{\mu}m$) were used as seeding particles to simulate the molten Zinc flow and deoxidization gas flow, respectively. A pulsed Nd:Yag laser and a $2K{\times}2K$ high-resolution CCD camera were synchronized for the PIV velocity field measurement. From flow visualization study, it is found that the liquid flow in the Zn pot is dominantly governed by the uprising flow caused by the rotating sink roll, with its effect on the steel strip inside the snout largely diminished by installing of the snout. The deoxidization gas flow in front of the strip inside the snout can be characterized by a large-scale vortex rotating clockwise direction formed by the moving strip. In the rear side of the strip, a counter-clockwise vortex is formed and some of the flow entrained by the moving strip impinges on the free surface of molten zinc. The liquid flow in front of the strip is governed by the flow entering the snout, caused by the spinning sink roll. Just below the free surface a counter-clockwise vortex is formed near the snout wall. The moving strip affects dominantly the flow behind the strip inside the snout, and large amount of the liquid flow follows the moving strip toward the sink roll. The thickness of the flow following the strip is very thin in the front side due to the uprising flow, however thick boundary layer is formed in the rear side of the strip. Its thickness is increased as moving downstream toward the sink roll. Inside the snout, the deoxidization gas flow above the free surface is much faster than the liquid flow in the zinc pot. Due to the larger influx of the flow following the moving strip in the rear side of the strip, higher percentage of imperfection can be anticipated on the rear surface of the strip.