Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Computational study of subcritical response in flow past a circular cylinder

C. D. Cantwell and D. Barkley
Phys. Rev. E 82, 026315 – Published 25 August 2010

Abstract

Flow past a circular cylinder is investigated in the subcritical regime, below the onset of Bénard-von Kármán vortex shedding at Reynolds number Rec47. The transient response of infinitesimal perturbations is computed. The domain requirements for obtaining converged results is discussed at length. It is shown that energy amplification occurs as low as Re=2.2. Throughout much of the subcritical regime the maximum energy amplification increases approximately exponentially in the square of Re reaching 6800 at Rec. The spatiotemporal structure of the optimal transient dynamics is shown to be transitory Bénard-von Kármán vortex streets. At Re42 the long-time structure switches from exponentially increasing downstream to exponentially decaying downstream. Three-dimensional computations show that two-dimensional structures dominate the energy growth except at short times.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
10 More
  • Received 19 May 2010

DOI:https://doi.org/10.1103/PhysRevE.82.026315

©2010 American Physical Society

Authors & Affiliations

C. D. Cantwell* and D. Barkley

  • Mathematics Institute and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom

  • *c.cantwell@imperial.ac.uk
  • d.barkley@warwick.ac.uk

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 82, Iss. 2 — August 2010

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1
    Diagram of the cylinder geometry (not to scale), showing the inflow, outflow and cross-stream dimensions referenced later. Also marked are the separation streamlines and the downstream stagnation point, xs.Reuse & Permissions
  • Figure 2
    Figure 2
    Representative base flow. Streamlines of the flow at Re=40. The computational domain is much larger than the region shown. The main flow field uses a contour spacing of 0.25, while a smaller spacing of 0.01 is used to highlight the structure of the recirculation bubble.Reuse & Permissions
  • Figure 3
    Figure 3
    Representative spectral-elemental mesh. Dimensions are Li=45, Lc=45, and Lo=125 (refer to Fig. 1 for definitions).Reuse & Permissions
  • Figure 4
    Figure 4
    Convergence of optimal growth calculation with mesh geometry by (a) in-flow length Li and (b) cross-stream length Lc. Points indicate the computed values. Optimal growth is for a time horizon of τ=20. Percentage errors are shown relative to the calculation using Li=65 and Lc=65, respectively.Reuse & Permissions
  • Figure 5
    Figure 5
    (Color online) Energy of the optimal initial condition for Re=45, τ=100 shown on two scales: (a) matching that of the figures shown later in the paper and, (b) showing an extended scale which highlights the upstream tail of the perturbation.Reuse & Permissions
  • Figure 6
    Figure 6
    Optimal energy growth at Reynolds numbers from Re=5 to Re=50 in increments of Re=5. Points indicate the computed values. The case Re=50 is above Rec and is shown as dashed. The horizontal line is an estimate of the maximum growth in the subcritical regime (see text).Reuse & Permissions
  • Figure 7
    Figure 7
    Contour plot of optimal energy growth in the subcritical regime, with contour levels as indicated. The thick black curve denotes the contour of no-growth, G=1.Reuse & Permissions
  • Figure 8
    Figure 8
    Maximum growth Gmax as a function of the square of Reynolds number. The vertical dotted line corresponds to Rec. Gcmax, the maximum growth at Rec, is indicated with a point and is computed differently from cases Re46 (see text). The dashed-dotted line highlights the relationship given by Eq. (10).Reuse & Permissions
  • Figure 9
    Figure 9
    Energy of evolving perturbations computed at Re=40 for three time horizons: τ=20, τ=65, and τ=95. The three curves touch the optimal growth envelope (circles) at their respective τ values. Qualitatively similar evolution is seen over a large range of optimal perturbations.Reuse & Permissions
  • Figure 10
    Figure 10
    (Color online) Contours of energy showing the linear evolution of perturbations at Re=20 (left) and Re=40 (right). The panels are snapshots at 10 time-unit intervals from t=0 (bottom) to t=50 (top).Reuse & Permissions
  • Figure 11
    Figure 11
    (Color online) Same evolution as in Fig. 10, Re=20 (left) and Re=40 (right), but viewed in terms of vorticity. The vorticity of a linear superposition of the base flow and the perturbation is shown at snapshots separated by 10 time units from t=0 (bottom) up to t=50 (top). The maximum vorticity is 1.5.Reuse & Permissions
  • Figure 12
    Figure 12
    (Color online) Plots showing the space-time evolution of energy in the perturbation u. Each tile covers x between 0 and 125 (horizontal) and t between 0 and 125 (vertical), and shows the energy of perturbations sampled on the centerline y=0. Each row corresponds to a different Reynolds number, specifically: Re=20 (a,b,c), Re=30 (d,e,f), Re=40 (g,h,i), Re=45 (j,k,l), and Re=50 (m,n,o). For each Reynolds number we show the evolution of the dominant mode (left column), the first sub-dominant mode (center column), and a combination (right column) revealing the envelope of the perturbation as explained in the text. The scale of each row of tiles is normalized by the maximum energy in the right column over the space-time domain, with white corresponding to the highest energy and black to zero energy.Reuse & Permissions
  • Figure 13
    Figure 13
    Optimal growth as function of spanwise wavenumber β at Re=40. β=0 is dominant for long time horizons, but higher β may provide slightly larger growth at short time horizons.Reuse & Permissions
  • Figure 14
    Figure 14
    Contour plot of optimal growth at Re=40. The thicker black line denotes the contour of no growth.Reuse & Permissions
  • Figure 15
    Figure 15
    Convergence of the base flow stagnation point with mesh dimensions. Points indicate the computed values. In (a) Lc is fixed at 25 and in (b) Li is fixed at 25.Reuse & Permissions
  • Figure 16
    Figure 16
    Streamwise velocity profiles of base flows at Re=5, (a) and (b), and at Re=46, (c) and (d), showing variation with Li and Lc. For variations in Li, we fix Lc=25. For variations in Lc, we fix Li=25.Reuse & Permissions
  • Figure 17
    Figure 17
    Convergence of critical Reynolds number Rec with mesh size. Points indicate the computed values.Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review E

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×