Experiments and direct numerical simulations of two-dimensional turbulence

Asymptotic turbulent friction in 2D rough-walled flows

The friction f is the property of wall-bounded flows that sets the pumping cost of a pipeline, the draining capacity of a river, and other variables of practical relevance. For highly turbulent rough-walled pipe flows, f depends solely on the roughness length scale r, and the f - r relation may be expressed by the Strickler empirical scaling f ∝ r ^1/3 Here, we show experimentally that for soap film flows that are the two-dimensional (2D) equivalent of highly turbulent rough-walled pipe flows, f ∝ r and the f - r relation is not the same in 2D as in 3D. Our findings are beyond the purview of the standard theory of friction but consistent with a competing theory in which f is linked to the turbulent spectrum via the spectral exponent α: In 3D, α = 5/3 and the theory yields f ∝ r ^1/3; in 2D, α = 3 and the theory yields f ∝ r.

Stability of a directional Marangoni flow

Marangoni flows result from surface-tension gradients, and these flows occur over finite distances on the surface, but the subsequent secondary flows can be observed on much larger lengthscales. These flows play major roles in various phenomena, from foam dynamics to microswimmer propulsion. We show here that if a Marangoni flow of soluble surfactants is confined laterally, the flow forms an inertial surface jet. A full picture of the flows on the surface is exhibited, and the velocity profile of the jet is predicted analytically, and is successfully compared with the experimental measurements. Moreover, this straight jet eventually destabilizes into meanders. A quantitative comparison between the theory and our experimental observations yields a very good agreement in terms of critical wavelengths. The characterization and understanding of the 2D flows generated by confined Marangoni spreading is a first step to understand the role of inertial effects in the Marangoni flows with and without confinement.

Instability onset and scaling laws of an auto-oscillating turbulent flow in a complex plasma

We study a complex plasma under microgravity conditions that is first stabilized with an oscillating electric field. Once the stabilization is stopped, the so-called heartbeat instability develops. We study how the kinetic energy spectrum changes during and after the onset of the instability and compare with the double cascade predicted by Kraichnan and Leith for two-dimensional turbulence. The onset of the instability manifests clearly in the ratio of the reduced rates of cascade of energy and enstrophy and in the power-law exponents of the energy spectra.

Reynolds-number dependence of the longitudinal dispersion in turbulent pipe flow

In Taylor's theory, the longitudinal dispersion in turbulent pipe flows approaches, on long time scales, a diffusive behavior with a constant diffusivity K_{L}, which depends empirically on the Reynolds number Re. We show that the dependence on Re can be determined from the turbulent energy spectrum. By using the intimate connection between the friction factor and the longitudinal dispersion in wall-bounded turbulence, we predict different asymptotic scaling laws of K_{L}(Re) depending on the different turbulent cascades in two-dimensional turbulence. We also explore numerically the K_{L}(Re) dependence in turbulent channel flows with smooth and rough walls using a lattice Boltzmann method.

Lattice Boltzmann simulation for forced two-dimensional turbulence

The direct numerical simulations of forced two-dimensional turbulent flow are presented by using the lattice Boltzmann method. The development of an energy-enstrophy double cascade is investigated in the two cases of external force of two-dimensional turbulence, Gaussian force and Kolmogorov force. It is found that the friction force is a necessary condition of the occurrence of a double cascade. The energy spectrum k(-3) in the enstrophy inertial range is in accord with the classical Kraichnan theory for both external forces. The energy spectrum of the Gaussian force case in an inverse cascade is k(-2); however, the Kolmogorov force drives the k(-5/3) energy in a backscatter cascade. The result agrees with Scott's standpoint, which describes nonrobustness of the two-dimensional turbulent inverse cascade. Also, intermittency is found for the enstrophy cascade in two cases of the external force form. Intermittency refers to the nonuniform distribution of saddle points in the two-dimensional turbulent flow.

Scaling of near-wall flows in quasi-two-dimensional turbulent channels

The law of the wall and the log law rule the near-wall mean velocity profile of three-dimensional turbulent flows. These well-known laws, which are validated by legions of experiments and simulations, may be universal. Here, using a soap-film channel, we report the first experimental test of these laws in quasi-two-dimensional turbulent channel flows under two disparate turbulent spectra. We find that despite the differences with three-dimensional flows, the laws prevail, albeit with notable distinctions: the two parameters of the log law are markedly distinct from their three-dimensional counterpart; further, one parameter (the von Kármán constant) is independent of the spectrum whereas the other (the offset of the log law) depends on the spectrum. Our results suggest that the classical theory of scaling in wall-bounded turbulence is incomplete wherein a key missing element is the link with the turbulent spectrum.

Transition of the scaling law in inverse energy cascade range caused by a nonlocal excitation of coherent structures observed in two-dimensional turbulent fields

We numerically investigate the inverse energy cascade range of two-dimensional Navier-Stokes turbulence. Our focus is on the universality of the Kolmogorov's phenomenology. In our direct numerical simulations, two types of forcing processes, the random forcing and the deterministic forcing, are employed besides the systematically varied numerical parameters. We first calculate the two-dimensional Navier-Stokes equations and confirm that results in the quasi steady state are consistent with the classical phenomenology for both types of forcing processes. It is also found that the difference in forcing process appears after the inverse energy cascade range reaches the system size; the dipole coherent vortices emerge and grow only when the random forcing is adopted. Then we add a large-scale drag term to the Navier-Stokes equations to obtain the statistically stationary state. When the random forcing is used, the scaling exponent of the energy spectrum in the stationary state starts to differ from the predicted -5/3 in the inverse energy cascade range as the infrared Reynolds number Re(d) increases, where Re(d) is defined as k(f)/k(d) with the forcing wave number k(f) and the large-scale drag wave number k(d). That can be interpreted as a transition phenomenon in which the local maximum vorticity grows like an order parameter caused by excitation of strong coherent vortices. Strong coherent vortices emerge and grow after the quasi steady state and destroy the scaling law when Re(d) is over a critical value. These coherent vortices are not due to the finite-size effect, unlike the dipole coherent vortices. On the other hand, when the deterministic forcing is adopted, strong coherent vortices are hardly seen and the -5/3 scaling law holds independently of Re(d). We examine the cases of the combination of both types of forcing processes and find that formation of such coherent vortices is sensitive to the mechanism of the external forcing process as well as the numerical parameters. Several types of large-scale drag terms are also tested and their insignificant influence on these qualitative properties is revealed.

Testing a missing spectral link in turbulence

Although the cardinal attribute of turbulence is the velocity fluctuations, these fluctuations have been ignored in theories of the frictional drag of turbulent flows. Our goal is to test a new theory that links the frictional drag to the spectral exponent α, a property of the velocity fluctuations in a flow. We use a soap-film channel wherein for the first time the value of α can be switched between 3 and 5/3, the two theoretically possible values in soap-film flows. To induce turbulence with α = 5/3, we make one of the edges of the soap-film channel serrated. Remarkably, the new theory of the frictional drag holds in both soap-film flows (for either value of the spectral exponent α) and ordinary pipe flows (where α = 5/3), even though these types of flow are governed by different equations.

Evidence for the double cascade scenario in two-dimensional turbulence

Statistical features of homogeneous, isotropic, two-dimensional turbulence is discussed on the basis of a set of direct numerical simulations up to the unprecedented resolution 32768(2). By forcing the system at intermediate scales, narrow but clear inertial ranges develop both for the inverse and for direct cascades where the two Kolmogorov laws for structure functions are simultaneously observed. The inverse cascade spectrum is found to be consistent with Kolmogorov-Kraichnan prediction and is robust with respect the presence of an enstrophy flux. The direct cascade is found to be more sensible to finite size effects: the exponent of the spectrum has a correction with respect theoretical prediction which vanishes by increasing the resolution.

Spontaneous angular momentum generation of two-dimensional fluid flow in an elliptic geometry

Spontaneous spin-up, i.e., the significant increase of the total angular momentum of a flow that initially has no net angular momentum, is very characteristic for decaying two-dimensional turbulence in square domains bounded by rigid no-slip walls. In contrast, spontaneous spin-up is virtually absent for such flows in a circular domain with a no-slip boundary. In order to acquire an understanding of this strikingly different behavior observed on the square and the circle, we consider a set of elliptic geometries with a gradual increase of the eccentricity. It is shown that a variation of the eccentricity can be used as a control parameter to tune the relative contribution of the pressure and viscous stresses in the angular momentum balance. Direct numerical simulations demonstrate that the magnitude of the torque can be related to the relative contribution of the pressure. As a consequence, the number of spin-up events in an ensemble of slightly different initial conditions depends strongly on the eccentricity. For small eccentricities, strong and rapid spin-up events are observed occasionally, whereas the majority of the runs do not show significant spin-up. Small differences in the initial condition can result in a completely different evolution of the flow and an appearance of the end state of the decay process. For sufficiently large eccentricities, all the runs in the ensemble demonstrate strong and rapid spin-up, which is consistent with the flow development on the square. It is verified that the number of spin-up events for a given eccentricity does not depend on the Reynolds number of the flow. This observation is consistent with the conjecture that it is the pressure on the domain boundaries that drives the spin-up processes.

Use of an ultrasound blood-mimicking fluid for Doppler investigations of turbulence in vitro

Turbulence is an important factor in the assessment of stenotic disease and a possible causative mechanism for thromboembolism. Previous Doppler studies of turbulence have typically used whole-blood preparations or suspensions of erythrocytes. Recently, a water-glycerol based blood-mimicking fluid (BMF) has been developed for use in Doppler ultrasound studies. This fluid has desirable ultrasound properties but it has not previously been described during in vitro investigations of turbulence intensity. We report on investigations of grid-generated and constrained-jet turbulence in an in vitro test system. The BMF was found to generate significant levels of turbulence during steady flow at physiological flow rates, producing turbulent patterns in the distal region that were consistent with previous studies. Turbulence intensity increased significantly with flow rate (p < 0.005) for both the constrained jet and the constrained grid. Based on our observations, we conclude that a water-glycerol based BMF provides a suitable working fluid during in vitro investigations of turbulence using Doppler ultrasound.

Craig's XY distribution and the statistics of Lagrangian power in two-dimensional turbulence

We examine the probability distribution function (PDF) of the energy injection rate (power) in numerical simulations of stationary two-dimensional (2D) turbulence in the Lagrangian frame. The simulation is designed to mimic an electromagnetically driven fluid layer, a well-documented system for generating 2D turbulence in the laboratory. In our simulations, the forcing and velocity fields are close to Gaussian. On the other hand, the measured PDF of injected power is very sharply peaked at zero, suggestive of a singularity there, with tails which are exponential but asymmetric. Large positive fluctuations are more probable than large negative fluctuations. It is this asymmetry of the tails which leads to a net positive mean value for the energy input despite the most probable value being zero. The main features of the power distribution are well described by Craig's XY distribution for the PDF of the product of two correlated normal variables. We show that the power distribution should exhibit a logarithmic singularity at zero and decay exponentially for large absolute values of the power. We calculate the asymptotic behavior and express the asymmetry of the tails in terms of the correlation coefficient of the force and velocity. We compare the measured PDFs with the theoretical calculations and briefly discuss how the power PDF might change with other forcing mechanisms.

Dynamics of energy condensation in two-dimensional turbulence

We report a numerical study, supplemented by phenomenological explanations, of "energy condensation" in forced 2D turbulence in a biperiodic box. Condensation is a finite size effect which occurs after the standard inverse cascade reaches the size of the system. It leads to the emergence of a coherent vortex dipole. We show that the time growth of the dipole is self-similar, and it contains most of the injected energy, thus resulting in an energy spectrum which is markedly steeper than the standard k{-5/3} one. Once the coherent component is subtracted, however, the remaining fluctuations have a spectrum close to k{-1}. The fluctuations decay slowly as the coherent part grows.