Youngjae Kim,Vedasri Godavarthi,Laura Victoria Rolandi,Joseph T. Klamo,Kunihiko Taira
Abstract:We investigate the effect of three-dimensionality on the synchronization characteristics of the wake behind an oscillating circular cylinder at Re = 300. Cylinder oscillations in rotation, transverse translation, and streamwise translation are considered. We utilize phase-reduction analysis, which quantifies the phase-sensitivity function of periodic flows, to examine the synchronization properties. Here, we present an ensemble-based framework for phase-reduction analysis to handle three-dimensional wakes that are not perfectly time-periodic. Based on the phase-sensitivity functions, synchronizability to three types of cylinder oscillations is evaluated. In spite of similar trends, we find that phase-sensitivity functions involving three-dimensional wakes are lower in magnitude compared to those of two-dimensional wakes, which leads to narrower conditions for synchronization to weak cylinder oscillations. We unveil that the difference between the phase-sensitivity functions of two- and three-dimensional flows is strongly correlated to the amplitude variation of the three-dimensional flow by the cylinder motions. This finding reveals that the cylinder motion modifies the three-dimensionality of the wake as well as the phase of vortex shedding, which leads to reduced phase modulation. The synchronization conditions of three-dimensional wakes, predicted by phase-reduction analysis, agree with the identification by parametric studies using direct numerical simulations for forced oscillations with small amplitudes. This study presents the potential capability of phase-reduction to study synchronization characteristics of complex flows.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to explore how three - dimensional effects influence the synchronization characteristics of the oscillating cylinder wake. Specifically, the researchers focused on the influence of three types of cylinder oscillations, namely rotation, lateral translation and longitudinal translation, on wake synchronization, and used the phase - reduction analysis method to evaluate the synchronism under these oscillation types. The core problem of the paper is to reveal the differences in synchronization characteristics between two - dimensional and three - dimensional wakes. In particular, the phase - sensitive function in the three - dimensional wake is lower, resulting in more stringent synchronization conditions for weakly oscillating cylinders. In addition, the study also found that the cylinder motion not only changes the phase of vortex shedding, but also affects the three - dimensional characteristics of the wake, thereby reducing the phase modulation. The synchronization conditions predicted by the phase - reduction analysis were verified by direct numerical simulation, demonstrating the potential of this method in studying the synchronization characteristics of complex fluids.
### Formula Summary
1. **Strouhal Number**:
\[
St=\frac{fL}{U_{\infty}}
\]
where \( f \) is the maximum energy frequency of lift fluctuations, \( L \) is the cylinder diameter, and \( U_{\infty} \) is the free - stream velocity.
2. **Lift Coefficient**:
\[
C_{L}=\frac{F_{L}}{\frac{1}{2}\rho U_{\infty}^{2}DL_{z}}
\]
where \( F_{L} \) is the lift acting on the cylinder, \( \rho \) is the fluid density, \( D \) is the cylinder diameter, and \( L_{z} \) is the spanwise length of the cylinder.
3. **Drag Coefficient**:
\[
C_{D}=\frac{F_{D}}{\frac{1}{2}\rho U_{\infty}^{2}DL_{z}}
\]
where \( F_{D} \) is the drag acting on the cylinder.
4. **Definitions of Phase and Amplitude**:
\[
\theta(t)=\theta_{0}+\tan^{- 1}\left(\frac{\tilde{C}_{L}(t)}{C_{L}(t)}\right)
\]
\[
r(t)=\sqrt{\tilde{C}_{L}(t)^{2}+C_{L}(t)^{2}}
\]
where \( \tilde{C}_{L}(t) \) is the Hilbert transform of \( C_{L}(t) \).
5. **Hilbert Transform**:
\[
\tilde{C}_{L}(t)=\frac{1}{\pi}\text{p.v.}\int_{-\infty}^{\infty}\frac{C_{L}(t)}{t-\tau}d\tau
\]
6. **Instantaneous Flow Snapshots**:
\[
u[\sigma_{1},\sigma_{2}](x,y,z,t)=\int_{S}u(x,y,z,t)\left[G(x',y';x,y,\sigma_{1})-G(x',y';x,y,\sigma_{2})\right]dx'dy'
\]
\[
G(x',y';x,y,\sigma)=\frac{1}{2\pi\sigma^{2}}\exp\left(-\frac{(x - x')^{2}+(y - y')^{2}}{2\sigma^{2}}\right)
\]
7. **Vorticity Field Decomposition**:
\[
\omega(x,t)=\sum_{m = 0}^{N_{z}/2}\hat{\omega}(x,y,t,m)e^{2\pi imz/L_{z}}=\omega_{0}(x,t)+\omega_{A}(x,t)+\omega_{B}(x,t)
\]
8. **Phase - Reduction Analysis**:
\[
\dot{\theta}=\O