Numerical study of three-dimensional Kolmogorov flow at high Reynolds numbers

High-resolution numerical simulations (with up to 2563 modes) are performed for three-dimensional flow driven by the large-scale constant force fy = F cos(x) in a periodic box of size L = 2π (Kolmogorov flow). High Reynolds number is attained by solving the Navier-Stokes equations with hyperviscosity (-1)h+1Δh (h = 8). It is shown that the mean velocity profile of Kolmogorov flow is nearly independent of Reynolds number and has the ‘laminar’ form vy = V cos(x) with a nearly constant eddy viscosity. Nevertheless, the flow is highly turbulent and intermittent even at large scales. The turbulent intensities, energy dissipation rate and various terms in the energy balance equation have the simple coordinate dependence a + b cos(2x) (with a, b constants). This makes Kolmogorov flow a good model to explore the applicability of turbulence transport approximations in open time-dependent flows. It turns out that the standard expression for effective (eddy) viscosity used in K-[Escr ] transport models overpredicts the effective viscosity in regions of high shear rate and should be modified to account for the non-equilibrium character of the flow. Also at large scales the flow is anisotropic but for large Reynolds number the flow is isotropic at small scales. The important problem of local isotropy is systematically studied by measuring longitudinal and transverse components of the energy spectra and crosscorrelation spectra of velocities and velocity-pressure-gradient spectra. Cross-spectra which should vanish in the case of isotropic turbulence decay only algebraically but somewhat faster than corresponding isotropic correlations. It is verified that the pressure plays a crucial role in making the flow locally isotropic. It is demonstrated that anisotropic large-scale flow may be considered locally isotropic at scales which are approximately ten times smaller than the scale of the flow.