skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, January 16 until 2:00 AM ET on Friday, January 17 due to maintenance. We apologize for the inconvenience.


Title: Forced oscillations of a cylinder in the flow of viscoelastic fluids

We investigate the effects of fluid elasticity on the flow forces and the wake structure when a rigid cylinder is placed in a viscoelastic flow and is forced to oscillate sinusoidally in the transverse direction. We consider a two-dimensional, uniform, incompressible flow of viscoelastic fluid at$Re=100$, and use the FENE-P model to represent the viscoelastic fluid. We study how the flow forces and the wake patterns change as the amplitude of oscillations,$A^*$, the frequency of oscillations (inversely proportional to a reduced velocity,$U^*$), the Weissenberg number,$Wi$, the square of maximum polymer extensibility,$L^2$, and the viscosity ratio,$\beta$, change individually. We calculate the lift coefficient in phase with cylinder velocity to determine the range of different system parameters where self-excited oscillations might occur if the cylinder is allowed to oscillate freely. We also study the effect of fluid elasticity on the added mass coefficient as these parameters change. The maximum elastic stress of the fluid occurs in between the vortices that are observed in the wake. We observe a new mode of shedding in the wake of the cylinder: in addition to the primary vortices that are also observed in the Newtonian flows, secondary vortices that are caused entirely by the viscoelasticity of the fluid are observed in between the primary vortices. We also show that, for a constant$Wi$, the strength of the polymeric stresses increases with increasing reduced velocity or with decreasing amplitude of oscillations.

 
more » « less
Award ID(s):
2126175
PAR ID:
10537476
Author(s) / Creator(s):
; ;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
975
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

     
    more » « less
  2. The motion of a freely rotating prolate spheroid in a simple shear flow of a dilute polymeric solution is examined in the limit of large particle aspect ratio,$\kappa$. A regular perturbation expansion in the polymer concentration,$c$, a generalized reciprocal theorem, and slender body theory to represent the velocity field of a Newtonian fluid around the spheroid are used to obtain the$O(c)$correction to the particle's orientational dynamics. The resulting dynamical system predicts a range of orientational behaviours qualitatively dependent upon$c\, De$($De$is the imposed shear rate times the polymer relaxation time) and$\kappa$and quantitatively on$c$. At a small but finite$c\, De$, the particle spirals towards a limit cycle near the vorticity axis for all initial conditions. Upon increasing$\kappa$, the limit cycle becomes smaller. Thus, ultimately the particle undergoes a periodic motion around and at a small angle from the vorticity axis. At moderate$c\, De$, a particle starting near the flow–gradient plane departs it monotonically instead of spirally, as this plane (a limit cycle at smaller$c\, De$) obtains a saddle and an unstable node. The former is close to the flow direction. Upon further increasing$c\, De$, the saddle node changes to a stable node. Therefore, depending upon the initial condition, a particle may either approach a periodic orbit near the vorticity axis or obtain a stable orientation near the flow direction. Upon further increasing$c\, De$, the limit cycle near the vorticity axis vanishes, and the particle aligns with the flow direction for all starting orientations.

     
    more » « less
  3. The wake flow past an axisymmetric body of revolution at a diameter-based Reynolds number$Re=u_{\infty }D/\nu =5000$is investigated via a direct numerical simulation. The study is focused on identification of coherent vortical motions and the dominant frequencies in this flow. Three dominant coherent motions are identified in the wake: the vortex shedding motion with the frequency of$St=fD/u_{\infty }=0.27$, the bubble pumping motion with$St=0.02$, and the very-low-frequency (VLF) motion originated in the very near wake of the body with the frequency$St=0.002$$0.005$. The vortex shedding pattern is demonstrated to follow a reflectional symmetry breaking mode, whereas the vortex loops are shed alternatingly from each side of the vortex shedding plane, but are subsequently twisted and tangled, giving the resulting wake structure a helical spiraling pattern. The bubble pumping motion is confined to the recirculation region and is a result of a Görtler instability. The VLF motion is related to a stochastic destabilisation of a steady symmetric mode in the near wake and manifests itself as a slow, precessional motion of the wake barycentre. The VLF mode with$St=0.005$is also detectable in the intermediate wake and may be associated with a low-frequency radial flapping of the shear layer.

     
    more » « less
  4. Electrophoresis is the motion of a charged colloidal particle in an electrolyte under an applied electric field. The electrophoretic velocity of a spherical particle depends on the dimensionless electric field strength$\beta =a^*e^*E_\infty ^*/k_B^*T^*$, defined as the ratio of the product of the applied electric field magnitude$E_\infty ^*$and particle radius$a^*$, to the thermal voltage$k_B^*T^*/e^*$, where$k_B^*$is Boltzmann's constant,$T^*$is the absolute temperature, and$e^*$is the charge on a proton. In this paper, we develop a spectral element algorithm to compute the electrophoretic velocity of a spherical, rigid, dielectric particle, of fixed dimensionless surface charge density$\sigma$over a wide range of$\beta$. Here,$\sigma =(e^*a^*/\epsilon ^*k_B^*T^*)\sigma ^*$, where$\sigma ^*$is the dimensional surface charge density, and$\epsilon ^*$is the permittivity of the electrolyte. For moderately charged particles ($\sigma ={O}(1)$), the electrophoretic velocity is linear in$\beta$when$\beta \ll 1$, and its dependence on the ratio of the Debye length ($1/\kappa ^*$) to particle radius (denoted by$\delta =1/(\kappa ^*a^*)$) agrees with Henry's formula. As$\beta$increases, the nonlinear contribution to the electrophoretic velocity becomes prominent, and the onset of this behaviour is$\delta$-dependent. For$\beta \gg 1$, the electrophoretic velocity again becomes linear in field strength, approaching the Hückel limit of electrophoresis in a dielectric medium, for all$\delta$. For highly charged particles ($\sigma \gg 1$) in the thin-Debye-layer limit ($\delta \ll 1$), our computations are in good agreement with recent experimental and asymptotic results.

     
    more » « less
  5. This experimental study investigates the fluid–structure–surface interactions of a flexibly mounted rigid plate in axial flow, focusing on flow-induced vibration (FIV) response and vortex dynamics of the system within a reduced velocity range of$U^*=0.29\unicode{x2013}8.73$, corresponding to a Reynolds number range of$Re=518\unicode{x2013}15\,331$. The plate, with one and two degrees of freedom (DoFs) for pitching and plunging oscillations, is examined at various submerged heights near the free surface. Results show that the plate exhibits divergence instability at low reduced velocities in both 1DoF and 2DoF systems. As the flow velocity surpasses a critical reduced velocity, periodic limit-cycle oscillations (LCOs) occur, increasing in amplitude until a second critical reduced velocity is reached. Beyond this point, LCOs are suppressed, and the plate experiences an increased static divergence angle with further flow velocity increase. The proximity to the free surface significantly influences the FIV response, with decreasing submerged heights leading to reduced LCO amplitudes and a shift of instabilities to higher reduced velocities. Vortex dynamics are analysed using time-resolved volumetric particle tracking velocimetry and hydrogen bubble flow visualisation. The analysis reveals disruptions in the symmetric flow field near the free surface, causing elongation and fragmentation of vortices in the wake of the plate, as well as vortex coupling. Proper orthogonal decomposition (POD) identifies dominant coherent structures, including leading-edge and trailing-edge vortices, captured in the first and second paired modes. On the other hand, higher POD modes capture the interaction of vortices in the wake and near the free surface.

     
    more » « less