Hostname: page-component-7b7567c8c-kwt9d Total loading time: 0 Render date: 2026-11-14T12:34:23.037Z Has data issue: false hasContentIssue false

Numerical study of flow-induced vibration of a circular cylinder with attached flexible splitter plate at low $Re$

Published online by Cambridge University Press:  10 October 2019

Tulsi Ram Sahu
Affiliation:
Department of Aerospace Engineering, Indian Institute of Technology Kanpur, Kanpur, UP 208016, India
Mohd Furquan
Affiliation:
Department of Aerospace Engineering, Indian Institute of Technology Kanpur, Kanpur, UP 208016, India
Sanjay Mittal*
Affiliation:
Department of Aerospace Engineering, Indian Institute of Technology Kanpur, Kanpur, UP 208016, India
*
†Email address for correspondence: smittal@iitk.ac.in

Abstract

Flow-induced vibration (FIV) of an elastically mounted circular cylinder with an attached splitter plate in uniform flow is studied numerically via a stabilized space–time finite element method. The Reynolds number based on the cylinder diameter $D$ and the free-stream speed is restricted to 150. The ratio of the density of the body to that of the fluid, for the major part of the study, is 10. Two different reduced speeds are defined to quantify the compliance of the elastic support and flexibility of the splitter plate, respectively: $U_{s}^{\ast }$ based on the natural frequency of the spring–mass system and $U_{p}^{\ast }$ based on the fundamental natural frequency of the plate. Flow past a stationary cylinder ($U_{s}^{\ast }=0$) with a flexible splitter plate of length $3.5D$ is studied at different values of $U_{p}^{\ast }$. The vibration response of the plate exhibits lock-in with various eigenmodes of the plate in different ranges of $U_{p}^{\ast }$. The onsets of these lock-in regions are abrupt and hysteretic. The elastically mounted cylinder, without the splitter plate, undergoes large-amplitude vortex-induced vibration (VIV) for $4<U_{s}^{\ast }<7$. These large-amplitude oscillations are a consequence of synchronization, wherein the vortex shedding frequency locks in to the cylinder oscillation frequency. A rigid splitter plate attached to the cylinder reduces significantly the peak amplitude during VIV. Increasing the length of the plate from $1.5D$ to $3.5D$ only marginally affects the peak amplitude. It, however, leads to a wider range of lock-in. Unlike the case of an isolated cylinder, the lock-in and desynchronization regimes are not well demarcated in the presence of the splitter plate. Further, galloping is observed beyond a critical value of $U_{s}^{\ast }$; the amplitude of vibration increases with an increase in $U_{s}^{\ast }$ while the vibration frequency is relatively low and remains nearly constant. Increase in plate length delays, in terms of $U_{s}^{\ast }$, the onset of galloping. It is also found that the flexibility of the plate affects the maximum oscillation amplitude in the VIV regime. It also dictates the presence/absence of galloping. The system behaves similar to an isolated cylinder for a very flexible plate. The response is devoid of galloping, but relatively large amplitude of oscillation is observed during lock-in. The behaviour of the cylinder with a stiff plate is similar to that with the rigid one. The galloping instability sets in when the flexibility of the plate is less than a certain value ($U_{p}^{\ast }<4.7$, approximately for $U_{s}^{\ast }=22$). The VIV and galloping are separated by a range of $U_{s}^{\ast }$ in which the flow is either steady, for longer plates, or exhibits very weak vortex shedding. In the VIV regime, the plate tip and cylinder vibrate in phase for low $U_{p}^{\ast }$; their motion is out of phase for larger $U_{p}^{\ast }$. The change in phase is also associated with change in the frequency of vibration. At low $U_{p}^{\ast }$, the frequency of vibration is close to the first natural frequency of the system, while at high $U_{p}^{\ast }$ it becomes closer to the second natural frequency. The vibration amplitude of the cylinder is close to maximum in the VIV regime for $U_{s}^{\ast }=6$. Computations for various $U_{p}^{\ast }$, for $U_{s}^{\ast }=6$ and $22$, are utilized to determine optimal flexibility that leads to minimal FIV. The effect of the length of the flexible splitter plate, mass ratio and damping ratio is studied. A strategy is proposed to utilize the computations from various combinations of $U_{s}^{\ast }$ and $U_{p}^{\ast }$ to choose the appropriate flexibility of the attached splitter plate to minimize FIV.

Information

Type
JFM Papers
Copyright
© 2019 Cambridge University Press 
Figure 0

Figure 1.Schematic of the computational set-up (not to scale) for a spring-mounted circular cylinder with flexible splitter plate. Also shown are the boundary conditions.

Figure 1

Table 1.Description of some of the key symbols. The scales used for non-dimensionalization of the quantities are placed in square brackets.

Figure 2

Table 2.The various configurations of cylinder and splitter plate studied in this work.

Figure 3

Table 3.Flow past a spring-mounted circular cylinder with flexible splitter plate for$Re=150$, $m^{\ast }=10$, $U_{s}^{\ast }=6$, $U_{p}^{\ast }=12.8$, $L=3.5D$ and $B=2\,\%$: effect of mesh resolution and time step size ($\unicode[STIX]{x0394}t$).

Figure 4

Figure 2.Free vibration of a flexible splitter plate attached to a fixed cylinder (${\mathcal{C}}_{stat}{\mathcal{P}}_{vib}$): variation of non-dimensionalized natural frequencies, $F_{np}^{k}$, with reduced velocity, $U_{p}^{\ast }$, for the first four modes of vibration. The corresponding mode shapes of the plate are shown for $L=3.5D$.

Figure 5

Figure 3.Free vibration of an elastically supported circular cylinder with flexible plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of non-dimensionalized natural frequency, $F_{nsp}^{k}$, with the reduced velocity, $U_{p}^{\ast }$, for various values of $U_{s}^{\ast }$. The natural frequencies associated with different eigenmodes lie in non-overlapping regions, shown with different shades. These regions are bounded by the $F_{np}^{k}$ versus $U_{p}^{\ast }$ curves, drawn with solid lines, which are also shown in figure 2.

Figure 6

Figure 4.Free vibration of an elastically supported circular cylinder with flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): figure 3 that shows variation of natural frequency with $U_{p}^{\ast }$ and $U_{s}^{\ast }$ is re-plotted in terms of similarity variables $U_{s}^{\ast }F_{nsp}^{k}$ and $U_{p}^{\ast }F_{nsp}^{k}$. Curves I, II, III and IV, shown in different colours, correspond to first four natural frequencies. Also shown are the different mode shapes for $L=3.5D$.

Figure 7

Figure 5.Flow past a stationary cylinder with rigid splitter plate of length $L$ (${\mathcal{C}}_{stat}{\mathcal{P}}_{stat}$): variation of (a) root mean square of the lift coefficient, $C_{L_{rms}}$, and (b) time-averaged drag coefficient, $\bar{C}_{D}$, with $Re$. The data points are shown for only those $Re$ for which the flow is unsteady.

Figure 8

Figure 6.The $Re=150$ flow past a stationary cylinder (${\mathcal{C}}_{stat}$) and with an attached rigid splitter plate of length $L$ (${\mathcal{C}}_{stat}{\mathcal{P}}_{stat}$): time variation of (a) lift and (b) drag coefficient. Here $L=0$ corresponds to the ${\mathcal{C}}_{stat}$ case.

Figure 9

Figure 7.The $Re=150$ flow past a stationary cylinder (${\mathcal{C}}_{stat}$) and with an attached rigid splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{stat}$): vorticity field for (a) an isolated cylinder ($L=0$); and a cylinder with a rigid plate of length (b) $L=1.5D$, (c) $L=2.5D$ and (d) $L=3.5D$. All panels correspond to a time instant when the $C_{L}$ is at its peak value for the fully developed unsteady flow. Also shown are the time-averaged streamlines.

Figure 10

Figure 8.The $Re=150$ flow past a stationary cylinder (${\mathcal{C}}_{stat}$) and with an attached rigid splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{stat}$): fluctuation in vorticity field for (a) an isolated cylinder ($L=0$); and a cylinder with a rigid plate of length (b) $L=1.5D$, (c) $L=2.5D$ and (d) $L=3.5D$. All panels correspond to a time instant when the $C_{L}$ is at its peak value for the fully developed unsteady flow.

Figure 11

Figure 9.The $Re=150$ flow past a freely vibrating cylinder (${\mathcal{C}}_{vib}$): variation of (a) non-dimensionalized amplitude of cylinder displacement $A_{c}^{\ast }$ and (b) frequency ratio $f_{s}^{\ast }$ with reduced velocity $U_{s}^{\ast }$. Here $f_{s}^{\ast }$ is the ratio of the cylinder vibration frequency ($f_{v}$) to the natural frequency of the spring–mass system ($f_{ns}$). The lock-in regime is shaded in coral. Also marked are the desynchronization region (DS), initial branch (IB) and lower branch (LB). The insets in (a) show the instantaneous vorticity field for the fully developed unsteady flow for DS ($U_{s}^{\ast }=2$) and LB ($U_{s}^{\ast }=5$).

Figure 12

Figure 10.The $Re=150$ flow past a fixed cylinder with a flexible splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{vib}$) with $L=3.5D$: variation of non-dimensionalized (a) amplitude $A_{p}^{\ast }$ and (b) frequency $F_{v}$ of the plate-tip displacement with reduced speed $U_{p}^{\ast }$. Also shown for reference are the non-dimensional vortex shedding frequency for a stationary cylinder with rigid splitter plate, $F_{vrp_{0}}$, and fixed cylinder without splitter plate, $F_{v_{0}}$.

Figure 13

Figure 11.The $Re=150$ flow past a fixed cylinder with flexible splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{vib}$) with $L=3.5D$: shape of the flexible splitter plate at various time instants in one cycle of oscillation (top) and the corresponding time history of plate-tip displacement (bottom) for (a) $U_{p}^{\ast }=5$, (b) $U_{p}^{\ast }=60$ and (c) $U_{p}^{\ast }=100$. For $U_{p}^{\ast }=5$, the deflection of the plate is amplified by a factor of five for clarity. Images for $U_{p}^{\ast }=100$ do not correspond to a fully developed limit cycle, but during the time when the flow is still developing.

Figure 14

Figure 12.The $Re=150$ flow past a fixed cylinder with flexible splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{vib}$) with $L=3.5D$: dominant proper orthogonal decomposition modes for (a) $U_{p}^{\ast }=5$, (b) $U_{p}^{\ast }=60$ and (c) $U_{p}^{\ast }=100$. The numerals indicate ranking of the modes in decreasing order of their energy content.

Figure 15

Figure 13.The $Re=150$ flow past a fixed cylinder with flexible splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{vib}$) with $L=3.5D$: instantaneous vorticity field for (a) $U_{p}^{\ast }=5$, (b) $U_{p}^{\ast }=60$ and (c) $U_{p}^{\ast }=100$ at an instant when the displacement of the plate tip is at its maximum value. The image for $U_{p}^{\ast }=100$ does not correspond to a fully developed limit cycle, but during the time when the flow is still developing.

Figure 16

Figure 14.The $Re=150$ flow past a spring-mounted circular cylinder, with (${\mathcal{C}}_{vib}{\mathcal{P}}_{stat}$) and without (${\mathcal{C}}_{vib}$) a rigid splitter plate: variation of (a) normalized amplitude of cylinder displacement $A_{c}^{\ast }$ and (b) frequency ratio $f_{s}^{\ast }$ with reduced velocity $U_{s}^{\ast }$. Three plate lengths are considered. The various regimes of response for $L=3.5D$ are marked with different shades: VIV (yellow), steady (grey) and galloping (coral).

Figure 17

Figure 15.The $Re=150$ flow past a spring mounted circular cylinder with a rigid splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{stat}$): the instantaneous vorticity field for $L=3.5D$ and (a) $U_{s}^{\ast }=6$, (b) $U_{s}^{\ast }=12$, (c) $U_{s}^{\ast }=22$ and (d) $L=1.5D$ and $U_{s}^{\ast }=22$, when the cylinder displacement is maximum.

Figure 18

Figure 16.The $Re=150$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of the normalized amplitude of (a) cylinder oscillation $A_{c}^{\ast }$ and (b) plate tip vibration $A_{p}^{\ast }$ with reduced velocity $U_{s}^{\ast }$. (c) The phase, $\unicode[STIX]{x1D6F7}$, between the vibration of cylinder and plate tip.

Figure 19

Figure 17.The $Re=150$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of frequency ratio calculated using the natural frequency of the (a) first ($f_{sp1}^{\ast }$) and (b) second ($f_{sp2}^{\ast }$) mode of the system ${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$ with reduced velocity $U_{s}^{\ast }$.

Figure 20

Table 4.The $Re=150$, $L=3.5D$ flow past a spring-mounted cylinder with flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): effect of plate flexibility, $U_{p}^{\ast }$, on the extent of the various flow regimes in terms of $U_{s}^{\ast }$ and the peak response of the cylinder and plate tip.

Figure 21

Figure 18.The $Re=150$, $U_{s}^{\ast }=6$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of (a) normalized amplitude of cylinder oscillation $A_{c}^{\ast }$ and plate-tip oscillation $A_{p}^{\ast }$, (b) phase difference between the motion of cylinder and plate tip $\unicode[STIX]{x1D6F7}$ and (c) frequency ratios $f_{sp1}^{\ast }$ and $f_{sp2}^{\ast }$ calculated, respectively, using the first and second natural frequencies of the ${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$ modes with plate flexibility $U_{p}^{\ast }$. The points corresponding to the transition state are shown with solid symbols in (b,c).

Figure 22

Figure 19.The $Re=150$, $U_{s}^{\ast }=6$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): time histories of cylinder location and that of the plate tip for (a) $U_{p}^{\ast }=5.36$ (regime-i) and (b) $U_{p}^{\ast }=12.8$ (regime-ii). Vorticity field at a time instant when the cylinder displacement is maximum, and marked in the time histories, is shown below the respective time histories. For reference, the undeformed structure is also shown.

Figure 23

Figure 20.The $Re=150$, $U_{s}^{\ast }=6$, $U_{p}^{\ast }=5.68$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): (a) time histories of cylinder and plate-tip displacements and (b) vorticity field showing $P+S$ mode of vortex shedding. For comparison, the $2S$ mode of vortex shedding for $U_{p}^{\ast }=12.8$ is shown in (c).

Figure 24

Figure 21.The $Re=150$, $U_{s}^{\ast }=22$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of (a) normalized amplitude of cylinder oscillation $A_{c}^{\ast }$ and plate-tip oscillation $A_{p}^{\ast }$, (b) phase difference between the motion of cylinder and plate tip $\unicode[STIX]{x1D6F7}$ and (c) frequency ratios $f_{sp1}^{\ast }$ and $f_{sp2}^{\ast }$ calculated, respectively, using the first and second natural frequencies of the ${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$ modes with plate flexibility $U_{p}^{\ast }$.

Figure 25

Figure 22.The $Re=150$, $U_{s}^{\ast }=22$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): time histories of cylinder and plate tip positions for (a) $U_{p}^{\ast }=4.68$ (galloping) and (b) $U_{p}^{\ast }=7.4$ (VIV). Vorticity field at a time instant when the cylinder displacement is maximum is shown below the respective time histories. For reference, the undeformed structure is also shown.

Figure 26

Figure 23.The $Re=150$, $U_{s}^{\ast }=6$, $L=1.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of (a) normalized amplitude of cylinder oscillation $A_{c}^{\ast }$ and plate-tip oscillation $A_{p}^{\ast }$, (b) phase difference between the motion of cylinder and plate tip $\unicode[STIX]{x1D6F7}$ and (c) frequency ratios $f_{sp1}^{\ast }$ and $f_{sp2}^{\ast }$ calculated, respectively, using the first and second natural frequencies of the ${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$ modes with plate flexibility $U_{p}^{\ast }$.

Figure 27

Figure 24.The $Re=150$, $U_{s}^{\ast }=22$, $L=1.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$): variation of (a) normalized amplitude of cylinder oscillation $A_{c}^{\ast }$ and plate-tip oscillation $A_{p}^{\ast }$, (b) phase difference between the motion of cylinder and plate tip $\unicode[STIX]{x1D6F7}$ and (c) frequency ratios $f_{sp1}^{\ast }$, $f_{sp2}^{\ast }$ and $f_{sp3}^{\ast }$ calculated, respectively, using the first, second and third natural frequencies of the ${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$ modes with plate flexibility $U_{p}^{\ast }$.

Figure 28

Figure 25.The $Re=150$, $U_{p}^{\ast }=7.4$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$) for three different $m^{\ast }$: variation of the normalized amplitude of (a) cylinder oscillation $A_{c}^{\ast }$ and (b) plate-tip vibration $A_{p}^{\ast }$ with reduced velocity $U_{s}^{\ast }$. Shown in (c) is the phase difference, $\unicode[STIX]{x1D6F7}$, between the vibration of cylinder and plate-tip displacements.

Figure 29

Figure 26.The $Re=150$, $U_{p}^{\ast }=7.4$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$) for three different $m^{\ast }$: variation of frequency ratio calculated using the natural frequency of the (a) first ($f_{sp1}^{\ast }$) and (b) second ($f_{sp2}^{\ast }$) mode of the system ${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$ with reduced velocity $U_{s}^{\ast }$.

Figure 30

Figure 27.The $Re=150$, $m^{\ast }=10$, $U_{p}^{\ast }=7.4$, $L=3.5D$ flow past a spring-mounted circular cylinder with a flexible splitter plate (${\mathcal{C}}_{vib}{\mathcal{P}}_{vib}$) for three different $\unicode[STIX]{x1D701}$: variation of the normalized amplitude of (a) cylinder oscillation $A_{c}^{\ast }$ and (b) plate-tip vibration $A_{p}^{\ast }$ with reduced velocity $U_{s}^{\ast }$.

Figure 31

Figure 28.Flow past a fixed cylinder with a flexible splitter plate (${\mathcal{C}}_{stat}{\mathcal{P}}_{vib}$): variation of non-dimensionalized amplitude $A_{p}^{\ast }$ of the plate-tip displacement with (a) non-dimensional stiffness of the plate $K^{\ast }$ and (b) reduced speed $U_{p}^{\ast }$. Open symbols denote data from numerical studies while filled triangles represent data from the experimental study of Shukla et al. (2013).

Figure 32

Figure 29.The $Re=150$ flow past a spring-mounted circular cylinder with splitter plate amplitude of cylinder vibration for various $U_{p}^{\ast }$ and $U_{s}^{\ast }$ considered in this study. Coral shade marks those values of $U_{p}^{\ast }$ and $U_{s}^{\ast }$ wherein high-amplitude oscillation (${\geqslant}0.3D$) occurs at some flow speed. Designs that correspond to this ratio of $U_{p}^{\ast }/U_{s}^{\ast }$ may be avoided for applications where VIV is undesirable. However, such $U_{p}^{\ast }/U_{s}^{\ast }$ values can be used for applications like energy harvesting, where a large amplitude of vibration is desirable.