Vortex statistics for turbulence in a container with rigid boundaries

Effective merging dynamics of two and three fluid vortices: application to two-dimensional decaying turbulence

We present a kinetic theory of two-dimensional decaying turbulence in the context of two-body and three-body vortex merging processes. By introducing the equations of motion for two or three vortices in the effective noise due to all the other vortices, we demonstrate analytically that a two-body mechanism becomes inefficient at low vortex density n<<1. When the more efficient three-body vortex mergings are considered (involving vortices of different signs), we show that n~t(-ξ), with ξ=1. We generalize this argument to three-dimensional geostrophic turbulence, finding ξ=5/4, in excellent agreement with direct Navier-Stokes simulations [McWilliams et al., J. Fluid Mech. 401, 1 (1999)].

Regional statistics in confined two-dimensional decaying turbulence

Two-dimensional decaying turbulence in a square container has been simulated using the lattice Boltzmann method. The probability density function (PDF) of the vorticity and the particle distribution functions have been determined at various regions of the domain. It is shown that, after the initial stage of decay, the regional area averaged enstrophy fluctuates strongly around a mean value in time. The ratio of the regional mean and the overall enstrophies increases monotonously with increasing distance from the wall. This function shows a similar shape to the axial mean velocity profile of turbulent channel flows. The PDF of the vorticity peaks at zero and is nearly symmetric considering the statistics in the overall domain. Approaching the wall, the PDFs become skewed owing to the boundary layer.

Quasi-two-dimensional turbulence in shallow fluid layers: the role of bottom friction and fluid layer depth

The role of bottom friction and the fluid layer depth in numerical simulations and experiments of freely decaying quasi-two-dimensional turbulence in shallow fluid layers has been investigated. In particular, the power-law behavior of the compensated kinetic energy E0(t)=E(t)e(2lambda t), with E(t) the total kinetic energy of the flow and lambda the bottom-drag coefficient, and the compensated enstrophy Omega(0)(t)=Omega(t)e(2lambda t), with Omega(t) the total enstrophy of the flow, have been studied. We also report on the scaling exponents of the ratio Omega(t)/E(t), which is considered as a measure of the characteristic length scale in the flow, for different values of lambda. The numerical simulations on square bounded domains with no-slip boundaries revealed bottom-friction independent power-law exponents for E0(t), Omega(0)(t), and Omega(t)/E(t). By applying a discrete wavelet packet transform technique to the numerical data, we have been able to compute the power-law exponents of the average number density of vortices rho(t), the average vortex radius a(t), the mean vortex separation r(t), and the averaged normalized vorticity extremum omega(ext)(t)/square root E(t). These decay exponents proved to be independent of the bottom friction as well. In the experiments we have varied the fluid layer depth, and it was found that the decay exponents of E0(t), Omega(0)(t), Omega(t)/E(t), and omega(ext)(t)/square root E(t) are virtually independent of the fluid layer depth. The experimental data for rho(t) and a(t) are less conclusive; power-law exponents obtained for small fluid layer depths agree with those from previously reported experiments, but significantly larger power-law exponents are found for experiments with larger fluid layer depths.

Dissipation of kinetic energy in two-dimensional bounded flows

The role of no-slip boundaries as an enstrophy source in two-dimensional (2D) flows has been investigated for high Reynolds numbers. Numerical simulations of normal and oblique dipole-wall collisions are performed to investigate the dissipation of the kinetic energy E(t), and the evolution of the enstrophy Omega(t) and the palinstrophy P(t). It is shown for large Reynolds numbers that dE(t)/dt=-2Omega(t)/Re proportional, variant 1/sqrt[Re] instead of the familiar relation dE(t)/dt proportional, variant 1/Re as found for 2D unbounded flows.