• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.3 s.258
• /
• pp.273-282
• /
• 2007

Numerical calculations are carried out for the natural convection induced by temperature difference between a cold outer square cylinder and a hot inner circular cylinder. A two-dimensional solution for unsteady natural convection is obtained, using the immersed boundary method (IBM) to model an inner circular cylinder based on finite volume method, for different Rayleigh numbers varying over the range of $10^4\;to\;10^6$. The study goes further to investigate the effect of an inner cylinder location on the heat transfer and fluid flow. The location of inner circular cylinder is changed vertically along the center-line of square enclosure. The number, size and formation of cell strongly depend on Rayleigh number and the position of inner circular cylinder. The changes in heat transfer quantities have been presented.

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

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

• Wind and Structures
• /
• v.30 no.5
• /
• pp.533-546
• /
• 2020

Large Eddy Simulation (LES) is used to study the effects of steady slot suction on the aerodynamic forces of and flow around a wall-mounted finite-length square cylinder. The aspect ratio H/d of the tested cylinder is 5, where H and d are the cylinder height and width, respectively. The Reynolds number based on free-stream oncoming flow velocity U_∞ and d is 2.78×10⁴. The suction slot locates near the leading edge of the free end, with a width of 0.025d and a length of 0.9d. The suction coefficient Q (= U_s/U_∞) is varied as Q = 0, 1 and 3, where U_s is the velocity at the entrance of the suction slot. It is found that the free-end steady slot suction can effectively suppress the aerodynamic forces of the model. The maximum reduction of aerodynamic forces occurs at Q = 1, with the time-mean drag, fluctuating drag, and fluctuating lift reduced by 3.75%, 19.08%, 40.91%, respectively. For Q = 3, all aerodynamic forces are still smaller than those for Q = 0 (uncontrolled case), but obviously higher than those for Q = 1. The involved control mechanism is successfully revealed, based on the comparison of the flow around cylinder free end and the near wake for the three tested Q values.

• Wind and Structures
• /
• v.10 no.6
• /
• pp.511-532
• /
• 2007

The present study aims to generate turbulent inflow data to more accurately represent the turbulent flow around a square cylinder when the inflow turbulence level is significant. The modified random flow generation (RFG) technique in conjunction with a previously developed LES code is successfully adopted into a finite element based fluid flow solver to generate the required inflow turbulence boundary conditions for the three-dimensional (3-D) LES computations of transitional turbulent flow around a square cylinder at Reynolds number of 22,000. The near wall region is modelled without using wall approximate conditions and a wall damping coefficient is introduced into the calculation of sub-grid length scale in the boundary layer of the cylinder wall. The numerical results obtained from simulations are compared with each other and with the experimental data for different inflow turbulence boundary conditions in order to discuss the issues such as the synthetic inflow turbulence effects on the 3-D transitional flow behaviour in the near wake and the free shear layer, the basic mechanism by which stream turbulence interacts with the mean flow over the cylinder body and the prediction of integral flow parameters. The comparison among the LES results with and without inflow turbulence and the experimental data emphasizes that the turbulent inflow data generated by the present RFG technique for the LES computation can be a viable approach in accurately predicting the effects of inflow turbulence on the near wake turbulent flow characteristics around a bluff body.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.804-809
• /
• 2003

In this study, 2D computations of the Aeolian tones for some obstacles (circular cylinder, square cylinder and NACA0012 airfoil) are simulated. First of all, we calculate the flow noise generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method (FDLBM). The third-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuation with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. To investigate the effect of the lattice dependence, furthermore, simulations of the Aeolian tones at the low Reynolds number radiated by a square cylinder and a NACA0012 airfoil with a blunt trailing edge at high incidence are also investigated.

• Wind and Structures
• /
• v.12 no.1
• /
• pp.21-47
• /
• 2009

The logarithmic form for turbulent flow analysis guarantees the positivity of the turbulence variables as ${\kappa}$ and ${\varepsilon}$ of the ${\kappa}-{\varepsilon}$ model by using the natural logarithm of these variables. In the present study, the logarithmic form is incorporated into the finite element solution procedure for the unsteady turbulent flow analysis. A backward facing step flow using the standard ${\kappa}-{\varepsilon}$ model and a flow around a 2D square cylinder using the modified ${\kappa}-{\varepsilon}$ model (the Kato-Launder model) are simulated. These results show that the logarithmic form effectively keeps adequate balance of turbulence variables and makes the analysis stable during transient or unsteady processes.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Journal of the Korean Society of Visualization
• /
• v.22 no.2
• /
• pp.44-54
• /
• 2024

This study investigates the effectiveness of using a plasma actuator for active control of turbulent flow around a finite square cylinder. The primary objective is to analyze the impact of plasma actuators on flow separation and wake region characteristics, which are critical for reducing drag and suppressing vortex-induced vibrations. Direct Numerical Simulation (DNS) was employed to explore the flow dynamics at various operational parameters, including different actuation frequencies and voltages. The proposed methodology employs a neural network trained using the Proximal Policy Optimization (PPO) algorithm to determine optimal control policies for plasma actuators. This network is integrated with a computational fluid dynamics (CFD) solver for real-time control. Results indicate that this deep reinforcement learning (DRL)-based strategy outperforms existing methods in controlling flow, demonstrating robustness and adaptability across various flow conditions, which highlights its potential for practical applications.

• Wind and Structures
• /
• v.5 no.1
• /
• pp.15-24
• /
• 2002

This paper reports the numerical calculations of uniform turbulent shear flow around a square cylinder. The predictions are obtained by solving the two-dimensional unsteady Navier-Stokes equations in a finite volume technique. The turbulent fluctuations are simulated by the standard $k-{\varepsilon}$ model and one of its variant which takes care of the realizability constraint in order to suppress the excessive generation of turbulence in a stagnation condition. It has been found that the Strouhal number and the mean drag coefficient are almost unaffected by the shear parameter but the mean lift coefficient is increased. The present predictions are compared with available experimental data.