skip to main content

Attention:

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


Title: Coherent motions in a turbulent wake of an axisymmetric bluff body

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
Award ID(s):
1944568 1707075
PAR ID:
10483348
Author(s) / Creator(s):
;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
962
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Dynamic stall at low Reynolds numbers,$Re \sim O(10^4)$, exhibits complex flow physics with co-existing laminar, transitional and turbulent flow regions. Current state-of-the-art stall onset criteria use parameters that rely on flow properties integrated around the leading edge. These include the leading edge suction parameter or$LESP$(Rameshet al.,J. Fluid Mech., vol. 751, 2014, pp. 500–538) and boundary enstrophy flux or$BEF$(Sudharsanet al.,J. Fluid Mech., vol. 935, 2022, A10), which have been found to be effective for predicting stall onset at moderate to high$Re$. However, low-$Re$flows feature strong vortex-shedding events occurring across the entire airfoil surface, including regions away from the leading edge, altering the flow field and influencing the onset of stall. In the present work, the ability of these stall criteria to effectively capture and localize these vortex-shedding events in space and time is investigated. High-resolution large-eddy simulations for an SD7003 airfoil undergoing a constant-rate, pitch-up motion at two$Re$(10 000 and 60 000) and two pitch rates reveal a rich variety of unsteady flow phenomena, including instabilities, transition, vortex formation, merging and shedding, which are described in detail. While stall onset is reflected in both$LESP$and$BEF$, local vortex-shedding events are identified only by the$BEF$. Therefore,$BEF$can be used to identify both dynamic stall onset and local vortex-shedding events in space and time.

     
    more » « less
  2. 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
  3. 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
  4. Abstract

    Define theCollatz map${\operatorname {Col}} \colon \mathbb {N}+1 \to \mathbb {N}+1$on the positive integers$\mathbb {N}+1 = \{1,2,3,\dots \}$by setting${\operatorname {Col}}(N)$equal to$3N+1$whenNis odd and$N/2$whenNis even, and let${\operatorname {Col}}_{\min }(N) := \inf _{n \in \mathbb {N}} {\operatorname {Col}}^n(N)$denote the minimal element of the Collatz orbit$N, {\operatorname {Col}}(N), {\operatorname {Col}}^2(N), \dots $. The infamousCollatz conjectureasserts that${\operatorname {Col}}_{\min }(N)=1$for all$N \in \mathbb {N}+1$. Previously, it was shown by Korec that for any$\theta> \frac {\log 3}{\log 4} \approx 0.7924$, one has${\operatorname {Col}}_{\min }(N) \leq N^\theta $for almost all$N \in \mathbb {N}+1$(in the sense of natural density). In this paper, we show that foranyfunction$f \colon \mathbb {N}+1 \to \mathbb {R}$with$\lim _{N \to \infty } f(N)=+\infty $, one has${\operatorname {Col}}_{\min }(N) \leq f(N)$for almost all$N \in \mathbb {N}+1$(in the sense of logarithmic density). Our proof proceeds by establishing a stabilisation property for a certain first passage random variable associated with the Collatz iteration (or more precisely, the closely related Syracuse iteration), which in turn follows from estimation of the characteristic function of a certain skew random walk on a$3$-adic cyclic group$\mathbb {Z}/3^n\mathbb {Z}$at high frequencies. This estimation is achieved by studying how a certain two-dimensional renewal process interacts with a union of triangles associated to a given frequency.

     
    more » « less
  5. Abstract

    LetKbe an imaginary quadratic field and$p\geq 5$a rational prime inert inK. For a$\mathbb {Q}$-curveEwith complex multiplication by$\mathcal {O}_K$and good reduction atp, K. Rubin introduced ap-adicL-function$\mathscr {L}_{E}$which interpolates special values ofL-functions ofEtwisted by anticyclotomic characters ofK. In this paper, we prove a formula which links certain values of$\mathscr {L}_{E}$outside its defining range of interpolation with rational points onE. Arithmetic consequences includep-converse to the Gross–Zagier and Kolyvagin theorem forE.

    A key tool of the proof is the recent resolution of Rubin’s conjecture on the structure of local units in the anticyclotomic${\mathbb {Z}}_p$-extension$\Psi _\infty $of the unramified quadratic extension of${\mathbb {Q}}_p$. Along the way, we present a theory of local points over$\Psi _\infty $of the Lubin–Tate formal group of height$2$for the uniformizing parameter$-p$.

     
    more » « less