Jeffrey Nichols,Robert D. Guy,Becca Thomases
Abstract:Polymer solutions can develop chaotic flows, even at low inertia. This purely elastic turbulence is well studied, but little is known about the transition to chaos. In 2D channel flow and parallel shear flow, traveling wave solutions involving coherent structures are present for sufficiently large fluid elasticity. We numerically study 2D periodic parallel shear flow in viscoelastic fluids and show that these traveling waves become oscillatory and undergo a series of period-doubling bifurcations en-route to chaos.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how the Kolmogorov flow in elastic fluids transitions from an ordered coherent structure to a chaotic state under low - inertia conditions. Specifically, through numerical simulations, the researchers explored how, in a two - dimensional periodic parallel shear flow, as the Weissenberg number (Wi) increases, the traveling wave solutions (TWS) undergo a series of period - doubling bifurcations and finally enter the process of chaotic dynamics.
### Main research contents and methods:
1. **Research background**:
- Elastic turbulence refers to the complex flow phenomenon in polymer solutions due to fluid elasticity under low Reynolds numbers or inertia - free conditions.
- Although the research on elastic turbulence has been quite extensive, the transition mechanism from order to chaos is still unclear.
2. **Research object**:
- This paper focuses on the two - dimensional Kolmogorov flow, which is a periodic shear flow and is suitable for studying elastic turbulence under low - inertia conditions.
3. **Model and method**:
- The Stokes - Oldroyd B model is used to describe the motion of the fluid, where the fluid velocity \( \mathbf{u} \) and the polymer stress tensor \( \tau \) satisfy the following equations:
\[
- \nabla p+\eta_s \Delta \mathbf{u}+\nabla \cdot \tau+\mathbf{F} = 0, \quad \nabla \cdot \mathbf{u} = 0,
\]
\[
\tau+\lambda \nabla \tau=\eta_p(\nabla \mathbf{u}+(\nabla \mathbf{u})^T),
\]
where the upper convective derivative is defined as:
\[
\nabla \tau \equiv \frac{\partial \tau}{\partial t}+\mathbf{u} \cdot \nabla \tau-(\nabla \mathbf{u}) \tau+\tau(\nabla \mathbf{u})^T.
\]
4. **Numerical simulation**:
- The Fourier pseudo - spectral method is used for numerical simulation, and the initial condition is the analytical solution plus small random perturbations.
- The simulation results show that as the Weissenberg number Wi increases, the traveling wave solutions gradually become unstable and exhibit oscillatory behavior, and finally enter a chaotic state through a series of period - doubling bifurcations.
### Main results:
1. **Linear stability analysis**:
- When Wi < 9.61, the analytical solution is stable.
- Just beyond the critical value (9.65 ≤ Wi ≤ 11.2), standing wave solutions appear.
- As Wi continues to increase, traveling wave solutions (narwhals) appear, and at Wi ≈ 24, a sub - critical bifurcation occurs, and a new low - energy stable branch appears.
- Further increasing Wi, traveling wave solutions and oscillatory solutions coexist until Wi ≳ 27.3 when chaos is entered.
2. **Period - doubling bifurcations**:
- At Wi ≈ 27.12, the period jumps from 7.72 to 15.47.
- At Wi ≈ 27.216, the period doubles again.
- The process of period - doubling bifurcations is shown by visualizing the strain energy and the frequency spectrum.
### Conclusion:
This research systematically reveals for the first time the transition mechanism from coherent structures to chaotic dynamics in elastic fluids. Through numerical simulations and linear stability analysis, the researchers found that the traveling wave solutions are gradually...