Non-linear instability analysis of the two-dimensional Navier-Stokes equation: The Taylor-Green vortex problem

Abstract

An enstrophy-based non-linear instability analysis of the Navier-Stokes equation for two-dimensional (2D) flows is presented here, using the Taylor-Green vortex (TGV) problem as an example. This problem admits a time-dependent analytical solution as the base flow, whose instability is traced here. The numerical study of the evolution of the Taylor-Green vortices shows that the flow becomes turbulent, but an explanation for this transition has not been advanced so far. The deviation of the numerical solution from the analytical solution is studied here using a high accuracy compact scheme on a non-uniform grid (NUC6), with the fourth-order Runge-Kutta method. The stream function-vorticity (ψ, ω) formulation of the governing equations is solved here in a periodic square domain with four vortices at t = 0. Simulations performed at different Reynolds numbers reveal that numerical errors in computations induce a breakdown of symmetry and simultaneous fragmentation of vortices. It is shown that the actual physical instability is triggered by the growth of disturbances and is explained by the evolution of disturbance mechanical energy and enstrophy. The disturbance evolution equations have been traced by looking at (a) disturbance mechanical energy of the Navier-Stokes equation, as described in the work of Sengupta et al., “Vortex-induced instability of an incompressible wall-bounded shear layer,” J. Fluid Mech. 493, 277–286 (2003), and (b) the creation of rotationality via the enstrophy transport equation in the work of Sengupta et al., “Diffusion in inhomogeneous flows: Unique equilibrium state in an internal flow,” Comput. Fluids 88, 440–451 (2013).