Subharmonic instability of large-scale wavy structures in two-dimensional channels

TL;DR

This paper investigates the stability of large-scale wavy structures in two-dimensional channel turbulence, revealing a subharmonic instability at high Reynolds numbers that may generate turbulence.

Contribution

It introduces a combined DNS and Floquet analysis to identify a subharmonic instability in large-scale structures at high Reynolds numbers, advancing understanding of turbulence onset.

Findings

Large-scale wavy structures are stable at Re=3000.

A subharmonic torsional mode is unstable at Re=200000.

The unstable mode deforms and splits into multiple wave trains, indicating turbulence generation.

Abstract

A particular interest on two-dimensional turbulence is the inverse energy cascade from small to large sales, which leads to an energy condensation accompanied by the formation of large-scale vortical structures. Indeed, such a phenomenon is observed in the two-dimensional channel (2DCH) with large Reynolds numbers, where prominent large-scale wavy structures play a central role in the momentum and energy transfer across the inhomogeneous wall-normal direction \citep{Falkovich2018}. Yet, the instability of these wavy structures remains poorly understood, and it is unknown whether they have the capacity to generate turbulence. To address this, we first conduct the direct numerical simulation (DNS) of Navier-Stokes equations for 2DCH, then extract the large-scale wavy structures through the singular value decomposition, and finally perform a Floquet-based secondary instability analysis. Two bulk Reynolds numbers are examined in particular, i.e. $Re = 3000$ and $Re = 200000$, which lie on opposite sides of the transitional regime near $Re \approx 10000$ and cover the previously reported simulation domain. At $Re = 3000$, the large-scale wavy structure is found to be linearly stable, consistent with the laminar-like DNS flow field. However, at $Re = 200000$, a subharmonic torsional mode is identified, which leads to a definite growth rate ($\lambda_r = 0.18$) for the wavy structures with a half wave-length shift. Temporal reconstruction shows that this unstable mode deforms and splits into multiple wave trains and evolves in the opposite phase. Compared to the TS (Tollmien-Schlichting) wave of laminar flow, the subharmonic mode found here offers a novel understanding for the generation of turbulence in larger Reynolds number two-dimensional channels.