Research Papers: Fundamental Issues and Canonical Flows

Flow Unsteadiness and Stability Characteristics of Low-Re Flow Past an Inclined Triangular Cylinder

[+] Author and Article Information
Wei Zhang

Faculty of Mechanical Engineering
and Automation,
Zhejiang Sci-Tech University,
Hangzhou, Zhejiang 310018, China
e-mail: zhangwei@zstu.edu.cn

Hui Yang

Faculty of Mechanical Engineering
and Automation,
Zhejiang Sci-Tech University,
Hangzhou, Zhejiang 310018, China
e-mail: yanghui@zstu.edu.cn

Hua-Shu Dou

Faculty of Mechanical Engineering
and Automation,
Zhejiang Sci-Tech University,
Hangzhou, Zhejiang 310018, China
e-mail: huashudou@yahoo.com

Zuchao Zhu

Faculty of Mechanical Engineering
and Automation,
Zhejiang Sci-Tech University,
Hangzhou, Zhejiang 310018, China
e-mail: zhuzuchao@zstu.edu.cn

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received February 16, 2017; final manuscript received June 16, 2017; published online September 1, 2017. Assoc. Editor: Praveen Ramaprabhu.

J. Fluids Eng 139(12), 121203 (Sep 01, 2017) (14 pages) Paper No: FE-17-1099; doi: 10.1115/1.4037277 History: Received February 16, 2017; Revised June 16, 2017

The present study investigates the two-dimensional flow past an inclined triangular cylinder at Re = 100. Numerical simulation is performed to explore the effect of cylinder inclination on the aerodynamic quantities, unsteady flow patterns, time-averaged flow characteristics, and flow unsteadiness. We also provide the first global linear stability analysis and sensitivity analysis on the targeted physical problem for the potential application of flow control. The objective of this work is to quantitatively identify the effect of cylinder inclination on the characteristic quantities and unsteady flow patterns, with emphasis on the flow unsteadiness and instability. Numerical results reveal that the flow unsteadiness is generally more pronounced for the base-facing-like cylinders (α → 60 deg) where separation occurs at the front corners. The inclined cylinder reduces the velocity deficiency in the near-wake, and the reduction in far-wake is the most notable for the α = 30 deg cylinder. The transverse distributions of several quantities are shifted toward the negative y-direction, such as the maximum velocity deficiency and maximum/minimum velocity fluctuation. Finally, the global stability and sensitivity analysis show that the spatial structures of perturbed velocities are quite similar for α ≤ 30 deg and the temporal growth rate of perturbation is sensitive to the near-wake flow, while for α ≥ 40 deg there are remarkable transverse expansion and streamwise elongation of the perturbed velocities, and the growth rate is sensitive to the far-wake flow.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Chatterjee, D. , and Gupta, S. K. , 2015, “ Numerical Study of the Laminar Flow Past a Rotating Square Cylinder at Low Spinning Rates,” ASME J. Fluids Eng., 137(2), p. 021204. [CrossRef]
Supradeepan, K. , and Roy, A. , 2015, “ Low Reynolds Number Flow Characteristics for Two Side by Side Rotating Cylinders,” ASME J. Fluids Eng., 137(10), p. 101204. [CrossRef]
Chauhan, M. K. , Dutta, S. , Gandhi, B. K. , and More, B. S. , 2016, “ Experimental Investigation of Flow Over a Transversely Oscillating Square Cylinder at Intermediate Reynolds Number,” ASME J. Fluids Eng., 138(5), p. 051105. [CrossRef]
Kumar, S. , 2016, “ Effect of Channel Inlet Blockage on the Wake Structure of a Rotationally Oscillating Cylinder,” ASME J. Fluids Eng., 138(12), p. 121203. [CrossRef]
Zielinska, B. J. A. , and Wesfreid, J. E. , 1995, “ On the Spatial Structure of Global Modes in Wake Flow,” Phys. Fluids, 7(6), pp. 1418–1424. [CrossRef]
De, A. K. , and Dalal, A. , 2006, “ Numerical Simulation of Unconfined Flow Past a Triangular Cylinder,” Int. J. Numer. Methods Fluids, 52(7), pp. 801–821. [CrossRef]
Zeitoun, O. , Ali, M. , and Nuhait, A. , 2011, “ Convective Heat Transfer Around a Triangular Cylinder in an Air Cross Flow,” Int. J. Therm. Sci., 50(9), pp. 1685–1697. [CrossRef]
Marquet, O. , Sipp, D. , and Jacquin, L. , 2008, “ Sensitivity Analysis and Passive Control of Cylinder Flow,” J. Fluid Mech., 615, pp. 221–252. [CrossRef]
Zhang, W. , and Samtaney, R. , 2017, “ Effect of Corner Radius in Stabilizing the Low-Re Flow Past a Cylinder,” ASME J. Fluids Eng., (accepted).
De, A. K. , and Dalal, A. , 2007, “ Numerical Study of Laminar Forced Convection Fluid Flow and Heat Transfer From a Triangular Cylinder Placed in a Channel,” ASME J. Heat Transfer, 129(5), pp. 646–656. [CrossRef]
Bao, Y. , Zhou, D. , and Zhao, Y.-J. , 2010, “ A Two-Step Taylor-Characteristic-Based Galerkin Method for Incompressible Flows and Its Application to Flow Over Triangular Cylinder With Different Incidence Angles,” Int. J. Numer. Methods Fluids, 62(11), pp. 1181–1208.
Tu, J.-H. , Zhou, D. , Bao, Y. , Han, Z.-L. , and Li, R.-D. , 2014, “ Flow Characteristics and Flow Induced Forces of a Stationary and Rotating Triangular Cylinder With Different Incidence Angles at Low Reynolds Numbers,” J. Fluids Struct., 45, pp. 107–123. [CrossRef]
Alawadhi, E. M. , 2013, “ Numerical Simulation of Fluid Flow Past an Oscillating Triangular Cylinder in a Channel,” ASME J. Fluids Eng., 135(4), p. 041202. [CrossRef]
Zang, Y. , Street, R. L. , and Koseff, J. R. , 1994, “ A Non-Staggered Grid, Fractional Step Method for Time-Dependent Incompressible Navier–Stokes Equations in Curvilinear Coordinates,” J. Comput. Phys., 114(1), pp. 18–33. [CrossRef]
Åkervik, E. , Brandt, L. , Henningson, D. S. , Hœpffner, J. , Marxen, O. , and Schlatter, P. , 2006, “ Steady Solutions of the Navier–Stokes Equations by Selective Frequency Damping,” Phys. Fluids, 18(6), p. 068102. [CrossRef]
Zhang, W. , and Samtaney, R. , 2016, “ Numerical Simulation and Global Linear Stability Analysis of Low-Re Flow Past a Heated Circular Cylinder,” Int. J. Heat Mass Transfer, 98, pp. 584–595. [CrossRef]
Zhang, W. , and Samtaney, R. , 2016, “ BiGlobal Linear Stability Analysis on Low-Re Flow Past an Airfoil at High Angle of Attack,” Phys. Fluids, 28(4), p. 044105. [CrossRef]
Lehoucq, R. B. , Sorensen, D. C. , and Yang, C. , 1998, ARPACK Users' Guide: Solution of Largescale Eigenvalue Problems With Implicitly Restarted Arnoldi Methods, SIAM, Philadelphia, PA. [CrossRef]
Amestoy, P. R. , Duff, I. S. , L'Excellent, J. Y. , and Koster, J. , 2001, “ A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling,” SIAM J. Matrix Anal. Appl., 23(1), pp. 15–41. [CrossRef]
Amestoy, P. R. , Guermouche, A. , L'Excellent, J. Y. , and Pralet, S. , 2006, “ Hybrid Scheduling for the Parallel Solution of Linear Systems,” Parallel Comput., 32(2), pp. 136–156. [CrossRef]
Giannetti, F. , and Luchini, P. , 2007, “ Structural Sensitivity of the First Instability of the Cylinder Wake,” J. Fluid Mech., 581, pp. 167–197. [CrossRef]
Zielinska, B. J. A. , Goujon-Durand, S. , Dušek, J. , and Wesfreid, J. E. , 1997, “ Strongly Nonlinear Effect in Unstable Wakes,” Phys. Rev. Lett., 79(20), pp. 3893–3896. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic of the physical model

Grahic Jump Location
Fig. 5

Instantaneous streamline and vorticity (ωd/U0) fields at α = 20 deg

Grahic Jump Location
Fig. 2

Multiblock grid for the α = 10 deg configuration plotted at every eighth gridline in each direction for clarity

Grahic Jump Location
Fig. 3

Time histories of the lift and drag coefficients

Grahic Jump Location
Fig. 4

Instantaneous streamline and vorticity (ωd/U0) fields at α = 0 deg

Grahic Jump Location
Fig. 6

Instantaneous streamline and vorticity (ωd/U0) fields at α = 40 deg

Grahic Jump Location
Fig. 7

Instantaneous streamline and vorticity (ωd/U0) fields at α = 60 deg

Grahic Jump Location
Fig. 8

Azimuthal distribution of time-averaged (Cp,avg) and maximum fluctuating (ΔCp) pressure coefficient along the cylinder surface

Grahic Jump Location
Fig. 14

Growth rate sensitivity to the base flow modification ∇Uσ of the unstable mode: (a)–(g) α = 0 deg (10 deg) 60 deg. The magnitude of ∇Uσ is visualized by the contour and its orientation by arrows.

Grahic Jump Location
Fig. 13

Contour of the eigenfunction: (a1)–(g1)  Re(û) and (a2)–(g2)  Re(v̂). The subfigures are plotted at α = 0 deg (10 deg) 60 deg from top to bottom. The dashed line denotes the wake centerline.

Grahic Jump Location
Fig. 9

Streamwise distribution of time-averaged velocities along the wake centerline

Grahic Jump Location
Fig. 10

Transverse distribution of time-averaged streamwise velocity (uavg − U0)/U0 at several streamwise stations

Grahic Jump Location
Fig. 11

Streamwise distribution of maximum fluctuating velocities along the wake centerline

Grahic Jump Location
Fig. 12

Transverse distribution of maximum fluctuating streamwise velocity Δu/U0 at several streamwise stations



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In