• International Journal of Fluid Machinery and Systems
• /
• v.9 no.4
• /
• pp.382-393
• /
• 2016

Flow over a bluff body is an attractive research field in thermal engineering. In the present study, laminar flow over a confined heated square cylinder using CuO-Water nanofluid is considered. Unsteady two-dimensional Navier-Stokes and energy equations are solved numerically using finite volume method (FVM). Recent correlations for the thermal conductivity and viscosity of nanofluids, which are function of nanoparticle volume fraction, temperature and nanoparticle diameter, have been employed. The results of numerical solution are obtained for Richardson number, nanoparticle volume fractions and nanoparticle diameters ranges of 0-1, 1-5% and 30-100 nm respectively for a fixed Reynolds number of Re = 150. At a given volume concentration, the investigations reveal that the decreasing in size of nanoparticles produces an increase in heat transfer rates from the square cylinder and a decrease in amplitude of the lift coefficient. Also, the increment of Nusselt number is more pronounced at higher concentrations and higher Richardson numbers.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1754-1759
• /
• 2003

Two-dimensional flow over a pair of circular cylinders in side-by-side arrangements at low Reynolds numbers has been numerically investigated in this study. Numerical simulations are performed, using the immersed boundary method, in the ranges of $40{\leq}Re{\leq}160$ and $g^{\ast}<5$, where Re and $g^{\ast}$ are, respectively, the Reynolds number and the spacing between the two cylinder surfaces divided by the cylinder diameter. Results show that total six kinds of wake patterns are observed over the ranges: antiphase-synchronized, inphase-synchronized, flip-flopping, single bluff-body, deflected, and steady wake patterns. It is found that the characteristics of the flow significantly depends both on the Reynolds number and gap spacing, with the latter much stronger than the former. Instantaneous flow fields, time traces, flow statistics and so on are presented to identify the wake patterns and then to understand the underlying mechanism. It is remarkable that, for the deflected wake pattern, the gap flow is deflected invariably to the cylinder of higher drag coefficient and the deflection way does not change at all. Moreover, the bifurcation phenomena where either of two wake patterns can occur are found at certain flow conditions.

• Journal of Mechanical Science and Technology
• /
• v.19 no.3
• /
• pp.877-886
• /
• 2005

Two-dimensional incompressible Navier-Stokes equations are solved using SIMPLER method in the intrinsic curvilinear coordinates system to study the unsteady viscous flow physics over two-dimensional ellipses. Unsteady viscous flows over various thickness-to-chord ratios of 0.6, 0.8, 1.0, and 1.2 elliptic cylinders are simulated at different Reynolds numbers of 200, 400, and 1,000. This study is focused on the understanding the effects of Reynolds number and elliptic cylinder thickness on the drag and lift forces. The present numerical solutions are compared with available experimental and numerical results and show a good agreement. Through this study, it is observed that the Reynolds number and the cylinder thickness affect significantly the frequencies of the force oscillations as well as the mean values and the amplitudes of the drag and lift forces.

• Wind and Structures
• /
• v.29 no.1
• /
• pp.75-84
• /
• 2019

A steady slot suction near the free-end leading edge of a finite-length square cylinder was used to control its aerodynamic forces and vortex-induced vibration (VIV). The freestream oncoming flow velocity ($U_{\infty}$) was from 3.8 m/s to 12.8 m/s. The width of the tested cylinder d = 40 mm and aspect ratio H/d = 5, where H was the height of the cylinder. The corresponding Reynolds number was from 10,400 to 35,000. The tested suction ratio Q, defined as the ratio of suction velocity ($U_s$) at the slot over the oncoming flow velocity at which the strongest VIV occurs ($U_{\nu}$), ranged from 0 to 3. It was found that the free-end slot suction can effectively attenuate the VIV of a cantilevered square cylinder. In the experiments, the RMS value of the VIV amplitude reduced quickly with Q increasing from 0 to 1, then kept approximately constant for $Q{\geq}1$. The maximum reduction of the VIV occurs at Q = 1, with the vibration amplitude reduced by 92%, relative to the uncontrolled case. Moreover, the overall fluctuation lift of the finite-length square cylinder was also suppressed with the maximum reduction of 87%, which occurred at Q = 1. It was interesting to discover that the free-end shear flow was sensitive to the slot suction near the leading edge. The turbulent kinetic energy (TKE) of the flow over the free end was the highest at Q = 1, which may result in the strongest mixing between the high momentum free-end shear flow and the near wake.

• Wind and Structures
• /
• v.26 no.5
• /
• pp.331-341
• /
• 2018

A numerical study on the flow over a square cylinder in the vicinity of a wall is conducted for different Couette-Poiseuille-based non-uniform flow with the non-dimensional pressure gradient P varying from 0 to 5. The non-dimensional gap ratio L (=$H^{\ast}/a^{\ast}$) is changed from 0.1 to 2, where $H^{\ast}$ is gap height between the cylinder and wall, and $a^{\ast}$ is the cylinder width. The governing equations are solved numerically through finite volume method based on SIMPLE algorithm on a staggered grid system. Both P and L have a substantial influence on the flow structure, time-mean drag coefficient ${\bar{C}}_D$, fluctuating (rms) lift coefficient ($C_L{^{\prime}}$), and Strouhal number St. The changes in P and L leads to four distinct flow regimes (I, II, III and IV). Following the flow structure change, the ${\bar{C}}_D$, $C_L{^{\prime}}$, and St all vary greatly with the change in L and/or P. The ${\bar{C}}_D$ and $C_L{^{\prime}}$ both grow with increasing P and/or L. The St increases with P for a given L, being less sensitive to L for a smaller P (< 2) and more sensitive to L for a larger P (> 2). A strong relationship is observed between the flow regimes and the values of ${\bar{C}}_D$, $C_L{^{\prime}}$ and St. An increase in P affects the pressure distribution more on the top surface than on bottom surface while an increase in L does the opposite.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.167-175
• /
• 2006

A parametric study has been accomplished to figure out the effects of elliptic cylinder thickness, angle of attack, and Reynolds number on the unsteady lift and drag forces exerted on the elliptic cylinder. A two-dimensional incompressible Navier-Stokes flow solver is developed based on the SIMPLER method in the body-intrinsic coordinates system to analyze the unsteady viscous flow over elliptic cylinder. Thickness-to-chord ratios of 0.2, 0.4, and 0.6 elliptic cylinders are simulated at different Reynolds numbers of 400 and 600, and angles of attack of $10^{\circ},\;20^{\circ},\;and\;30^{\circ}$. Through this study, it is observed that the elliptic cylinder thickness, angle of attack, and Reynolds number are very important parameters to decide the lift and drag forces. All these parameters also affect significantly the frequencies of the unsteady force oscillations.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.459-462
• /
• 2002

Many active and passive flow control methods have been studied since decades, but there are only few works about flow control methods using ion wind. This paper presents an experimental study on the wake control behind a circular cylinder using ion wind, a bulk motion of neutral molecules driven by locally ionized air of corona discharge. Experiments are done f3r different electrohydrodynamic numbers - the ratio of an electrical body farce to a fluid Inertial force - from 0 to 2 and for the Reynolds number ranging from $4{\times}10^3\;to\;8{\times}10^3$. Pressure distributions over a cylinder surface are measured and flow visualizations are carried out by smoke wire method. Flow visualizations confirm that ion wind affects significantly the wake structure behind a circular cylinder and pressure drag could be dramatically reduced by the superimposing ion wind.

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

Many researches were carried out to estimate heat transfer rate on a circular cylinder in a uniform flow. Various empirical correlations were suggested in the past through experimental studies, however there are considerable discrepancies in the estimated values of heat transfer coefficient. The effect of fluid physical properties on the forced convective heat transfer between a circular cylinder and the external flow was numerically investigated in the present study, The flow and temperature fields were solved using a Finite Volume Method over a wider range of Prandtl number(0.7-40,000) than existing correlations. The cold as well as the hot cylinders in the uniform liquid flow of constant temperature were investigated. A unified correlation was obtainde for both cases.

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

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.30-37
• /
• 2005

To evaluate LBM we performed the simulation of the unsteady two dimensional flow over a square cylinder in a channel in moderate Reynolds number range, $100\~500$ by using LBM and Fractional-Step method. Frist of all we compared LBM solution of Poiseuille flow applied Farout and periodic boundary conditions with the analytical solution to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured 500x100 grids. Prescribed maximum velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback boundary condition was applied to the channel and the cylinder waifs. The flow patterns and vortex shedding strouhal numbers were compared with previous research results. The flow patterns by LBM were in agreement with the flow pattern by fractional step method. Furthermore the strouhal number computed by LBM simulation result was more accurate than that of fractional step method through the comparison of the previous research results.