Experimental investigation of transitional flow in a toroidal pipe

The flow instability and further transition to turbulence in a toroidal pipe (torus) with curvature (tube-to-coiling diameter) 0.049 is investigated experimentally. The flow inside the toroidal pipe is driven by a steel sphere fitted to the inner pipe diameter. The sphere is moved with constant azimuthal velocity from outside the torus by a moving magnet. The experiment is designed to investigate curved pipe flow by optical measurement techniques. Using stereoscopic particle image velocimetry, laser Doppler velocimetry and pressure drop measurements, the flow is measured for Reynolds numbers ranging from 1000 to 15000. Time- and space-resolved velocity fields are obtained and analysed. The steady axisymmetric basic flow is strongly influenced by centrifugal effects. On an increase of the Reynolds number we find a sequence of bifurcations. For Re=4075 a supercritical bifurcation to an oscillatory flow is found in which waves travel in the streamwise direction with a phase velocity slightly faster than the mean flow. The oscillatory flow is superseded by a presumably quasi-periodic flow at a further increase of the Reynolds number before turbulence sets in. The results are found to be compatible, in general, with earlier experimental and numerical investigations on transition to turbulence in helical and curved pipes. However, important aspects of the bifurcation scenario differ considerably

Viscous tilting and production of vorticity in homogeneous turbulence

Viscous depletion of vorticity is an essential and well known property of turbulent flows, balancing, in the mean, the net vorticity production associated with the vortex stretching mechanism. In this letter, we, however, demonstrate that viscous effects are not restricted to a mere destruction process, but play a more complex role in vorticity dynamics that is as important as vortex stretching. Based on the results from three dimensional particle tracking velocimetry experiments and direct numerical simulation of homogeneous and quasi-isotropic turbulence, we show that the viscous term in the vorticity equation can also locally induce production of vorticity and changes of the orientation of the vorticity vector (viscous tilting)

On turbulent entrainment and dissipation in dilute polymer solutions

We present a comparative experimental study of a turbulent flow developing in clear water and dilute polymer solutions (25 and 50 wppm polyethylene oxide). The flow is forced by a planar grid that oscillates vertically with stroke S and frequency f in a square container of initially still fluid. Two-component velocity fields are measured in a vertical plane passing through the center of the tank by using time resolved particle image velocimetry. After the forcing is initiated, a turbulent layer develops that is separated from the initially irrotational fluid by a sharp interface, the so-called turbulent/nonturbulent interface (TNTI). The turbulent region grows in time through entrainment of surrounding fluid until the fluid in the whole container is in turbulent motion. From the comparison of the experiments in clear water and polymer solutions we conclude: (i) Polymer additives modify the large scale shape of the TNTI. (ii) Both, in water and in the polymer solution the mean depth of the turbulent layer, H(t), follows the theoretical prediction for Newtonian fluids H(t)∞√Kt, where K∞S^2f is the “grid action.” (iii) We find a larger grid action for dilute polymer solutions than for water. As a consequence, the turbulent kinetic energy of the flow increases and the rate of energy input becomes higher. (iv) The entrainment rate β=v_e/v_(rms) (where v_e=dH/dt is the interface propagation velocity and v_(rms) is the root mean square of the vertical velocity) is lower for polymers (β_p≈0.7) than for water (β_w≈0.8). The measured values for β are in good agreement with similarity arguments, from which we estimate that in our experiment about 28% of the input energy is dissipated by polymers

Acceleration, pressure and related quantities in the proximity of the turbulent/non-turbulent interface

This paper presents an analysis of flow properties in the proximity of the turbulent/non-turbulent interface (TNTI), with particular focus on the acceleration of fluid particles, pressure and related small scale quantities such as enstrophy, ω2 = ωiωi, and strain, s2 = sijsij. The emphasis is on the qualitative differences between turbulent, intermediate and non-turbulent flow regions, emanating from the solenoidal nature of the turbulent region, the irrotational character of the non-turbulent region and the mixed nature of the intermediate region in between. The results are obtained from a particle tracking experiment and direct numerical simulations (DNS) of a temporally developing flow without mean shear. The analysis reveals that turbulence influences its neighbouring ambient flow in three different ways depending on the distance to the TNTI: (i) pressure has the longest range of influence into the ambient region and in the far region non-local effects dominate. This is felt on the level of velocity as irrotational fluctuations, on the level of acceleration as local change of velocity due to pressure gradients, Du/Dt ∂u/∂t − p/ρ, and, finally, on the level of strain due to pressure-Hessian/strain interaction, (D/Dt)(s2/2) (∂/∂t)(s2/2) −sijp,ij > 0; (ii) at intermediate distances convective terms (both for acceleration and strain) as well as strain production −sijsjkski > 0 start to set in. Comparison of the results at Taylor-based Reynolds numbers Reλ = 50 and Reλ = 110 suggests that the distances to the far or intermediate regions scale with the Taylor microscale λ or the Kolmogorov length scale η of the flow, rather than with an integral length scale; (iii) in the close proximity of the TNTI the velocity field loses its purely irrotational character as viscous effects lead to a sharp increase of enstrophy and enstrophy-related terms. Convective terms show a positive peak reflecting previous findings that in the laboratory frame of reference the interface moves locally with a velocity comparable to the fluid velocity fluctuation

Near infrared few-cycle pulses for high harmonic generation

We report on the development of tunable few-cycle pulses with central wavelengths from 1.6 um to 2 um. Theses pulses were used as a proof of principle for high harmonic generation in atomic and molecular targets. In order to generate such pulses we produced a filament in a 4 bar krypton cell. Spectral broadening by a factor of 2 to 3 of a 40 fs near infrared input pulse was achieved. The spectrally broadened output pulses were then compressed by fused silica plates down to the few-cycle regime close to the Fourier limit. The auto-correlation of these pulses revealed durations of about 3 cycles for all investigated central wavelengths. Pulses with a central wavelength of 1.7 um and up to 430 uJ energy per pulse were employed to generate high order harmonics in Xe, Ar and N2. Moving to near infrared few-cycle pulses opens the possibility to operate deeply in the non-perturbative regime with a Keldysh parameter smaller than 1. Hence, this source is suitable for the study of the non-adiabatic tunneling regime in most generating systems used for high order harmonic generation and attoscience.Comment: 12 pages, 4 figure

A Lagrangian investigation of the small-scale features of turbulent entrainment through particle tracking and direct numerical simulation

We report an analysis of small-scale enstrophy ω2 and rate of strain s2 dynamics in the proximity of the turbulent/non-turbulent interface in a flow without strong mean shear. The techniques used are three-dimensional particle tracking (3D-PTV), allowing the field of velocity derivatives to be measured and followed in a Lagrangian manner, and direct numerical simulations (DNS). In both experiment and simulation the Taylor-microscale Reynolds number is Reλ = 50. The results are based on the Lagrangian viewpoint with the main focus on flow particle tracers crossing the turbulent/non-turbulent interface. This approach allowed a direct investigation of the key physical processes underlying the entrainment phenomenon and revealed the role of small-scale non-local, inviscid and viscous processes. We found that the entrainment mechanism is initiated by self-amplification of s2 through the combined effect of strain production and pressure--strain interaction. This process is followed by a sharp change of ω2 induced mostly by production due to viscous effects. The influence of inviscid production is initially small but gradually increasing, whereas viscous production changes abruptly towards the destruction of ω2. Finally, shortly after the crossing of the turbulent/non-turbulent interface, production and dissipation of both enstrophy and strain reach a balance. The characteristic time scale of the described processes is the Kolmogorov time scale, τη. Locally, the characteristic velocity of the fluid relative to the turbulent/non-turbulent interface is the Kolmogorov velocity, u

Correlation density matrices for 1- dimensional quantum chains based on the density matrix renormalization group

A useful concept for finding numerically the dominant correlations of a given ground state in an interacting quantum lattice system in an unbiased way is the correlation density matrix. For two disjoint, separated clusters, it is defined to be the density matrix of their union minus the direct product of their individual density matrices and contains all correlations between the two clusters. We show how to extract from the correlation density matrix a general overview of the correlations as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. To determine the correlation density matrix, we calculate the ground state for a class of spinless extended Hubbard models using the density matrix renormalization group. This numerical method is based on matrix product states for which the correlation density matrix can be obtained straightforwardly. In an appendix, we give a detailed tutorial introduction to our variational matrix product state approach for ground state calculations for 1- dimensional quantum chain models. We show in detail how matrix product states overcome the problem of large Hilbert space dimensions in these models and describe all techniques which are needed for handling them in practice.Comment: 50 pages, 34 figures, to be published in New Journal of Physic

Stroboscopic observation of quantum many-body dynamics

Recent experiments have demonstrated single-site resolved observation of cold atoms in optical lattices. Thus, in the future it may be possible to take repeated snapshots of an interacting quantum many-body system during the course of its evolution. Here we address the impact of the resulting quantum (anti-)Zeno physics on the many-body dynamics. We use the time-dependent density-matrix renormalization group to obtain the time evolution of the full wave function, which is then periodically projected in order to simulate realizations of stroboscopic measurements. For the example of a one-dimensional lattice of spinless fermions with nearest-neighbor interactions, we find regimes for which many-particle configurations are stabilized or destabilized, depending on the interaction strength and the time between observations

The rise of fully turbulent flow

Over a century of research into the origin of turbulence in wallbounded shear flows has resulted in a puzzling picture in which turbulence appears in a variety of different states competing with laminar background flow. At slightly higher speeds the situation changes distinctly and the entire flow is turbulent. Neither the origin of the different states encountered during transition, nor their front dynamics, let alone the transformation to full turbulence could be explained to date. Combining experiments, theory and computer simulations here we uncover the bifurcation scenario organising the route to fully turbulent pipe flow and explain the front dynamics of the different states encountered in the process. Key to resolving this problem is the interpretation of the flow as a bistable system with nonlinear propagation (advection) of turbulent fronts. These findings bridge the gap between our understanding of the onset of turbulence and fully turbulent flows.Comment: 31 pages, 9 figure