Source author record

Eberhard Bodenschatz

Eberhard Bodenschatz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
21topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

Experimental measurement of respiratory particles dispersed by wind instruments and analysis of the associated risk of infection transmission

Activities such as singing or playing a wind instrument release respiratory particles into the air that may contain pathogens and thus pose a risk for infection transmission. Here we report measurements of the size distribution, number, and volume concentration of exhaled particles from 31 healthy musicians playing 20 types of wind instruments using aerosol spectrometry and in-line holography in a strictly controlled cleanroom environment. We find that playing wind instruments carries a lower risk of airborne disease transmission than speaking or singing. We attribute this to the fact that the resonators of wind instruments act as filters for particles >10 $μ$m in diameter. We have also measured the size-dependent filtering properties of different types of filters that can be used as instrument masks. Based on these measurements, we calculated the risk of airborne transmission of SARS-CoV-2 in different near- and far-field scenarios with and without masking and/or distancing. We conclude that in all cases where there is a possibility that the musician is infectious, the only safe measure to prevent airborne transmission of the disease is the use of well-fitting and well-filtering masks for the instrument and the susceptible person.

preprint2022arXiv

Scaling in Decaying Turbulence at High Reynolds Numbers

The way the increment statistics of turbulent velocity fluctuations scale with the increment size is a centerpiece of turbulence theories. We report data on decaying turbulence in the Max Planck Variable Density Turbulence Tunnel (VDTT), which show an approach of the inertial range statistics toward a nontrivial shape at small scales. By correcting for the contributions of energy decay to the large-scale statistics with a model, we find the scaling exponent of the second-order velocity increment statistics to be independent of the Reynolds number and equal to $0.693\pm0.003$ for $2000\lesssim R_λ \lesssim 6000$. This is evidence of a universal inertial range at high Reynolds numbers.

preprint2021arXiv

Experimental observation of the elastic range scaling in turbulent flow with polymer additives

Minute amount of long chain flexible polymer dissolved in a turbulent flow can drastically change flow properties, such as reducing the drag and enhancing mixing. One fundamental riddle is how these polymer additives interact with the eddies of different spatial scales existing in the turbulent flow and in turn alter the turbulence energy transfer. Here we show how turbulent kinetic energy is transferred through deferent scales in the presence of the polymer additives. In particular, we observed experimentally the emerging of a new scaling range, referred to as the elastic range, where increasing amount of energy is transferred by the elasticity of the polymers. In addition, the existence of the elastic range prescribes the scaling of high-order velocity statistics. Our findings have important implications to many turbulence systems such as turbulence in plasmas or superfluids where interaction between turbulent eddies and other nonlinear physical mechanisms are often involved.

preprint2021arXiv

Fabrication of free-standing Pt nanowires for use as thermal anemometry probes in turbulence measurements

We report a robust fabrication method for patterning free-standing Pt nanowires for the use as thermal anemometry probes for small-scale turbulence measurements. Using e-beam lithography, high aspect ratio Pt nanowires (~300 nm width, ~70 $μ$m length, ~100 nm thickness) were patterned on the surface of oxidized silicon (Si) wafers. Combining precise wet etching processes with dry etching processes, these Pt nanowires have been successfully released free-standing between two silicon dioxide (SiO2) beams supported on Si cantilevers. Moreover, the unique design of the bridge holding the device allowed to release the device gently without damaging the Pt nanowires. The total fabrication time was minimized by restricting the use of e-beam lithography to the patterning of the Pt nanowires while standard photolithography was employed for other parts of the devices. We demonstrate that the fabricated sensors are suitable for turbulence measurements when operated in a constant-current mode. A robust calibration between output voltage and fluid velocity was established over the velocity range from 0.5 m s-1 to 5 m s-1 in an SF6 atmosphere at a pressure of 2 bar and a temperature of 21°C. The sensing signal from the nanowires showed negligible drift over a period of several hours. Moreover, we confirmed that the nanowires are able to withstand high dynamic pressures by testing them in air at room temperature velocities up to 55 m/s.

preprint2021arXiv

Generation of intense dissipation in high Reynolds number turbulence

Intense fluctuations of energy dissipation rate in turbulent flows result from the self-amplification of strain rate via a quadratic nonlinearity, with contributions from vorticity (via the vortex stretching mechanism) and the pressure Hessian tensor, which we analyze here using direct numerical simulations of isotropic turbulence in periodic domains of up to $12288^3$ grid points, and Taylor-scale Reynolds numbers in the range $140-1300$. We extract the statistics of various terms involved in amplification of strain and additionally condition them on the magnitude of strain. We find that strain is overall self-amplified by the quadratic nonlinearity, and depleted via vortex stretching; whereas pressure Hessian acts to redistribute strain fluctuations towards the mean-field and thus depleting intense strain. Analyzing the intense fluctuations of strain in terms of its eigenvalues reveals that the net amplification is solely produced by the third eigenvalue, resulting in strong compressive action. In contrast, the self-amplification terms acts to deplete the other two eigenvalues, whereas vortex stretching acts to amplify them, both effects canceling each other almost perfectly. The effect of the pressure Hessian for each eigenvalue is qualitatively similar to that of vortex stretching, but significantly weaker in magnitude. Our results conform with the familiar notion that intense strain is organized in sheet-like structures, which are in the vicinity of, but never overlap with regions of intense vorticity due to fundamental differences in their amplifying mechanisms.

preprint2020arXiv

Population Distribution in the Wake of a Sphere

The fluid physics of the heat and mass transfer from an object in its wake has much importance for natural phenomena as well as for many engineering applications. Here, we report numerical results on the population density of the spatial distribution of fluid velocity, pressure, scalar concentration and scalar fluxes of a wake flow past a sphere in the steady wake regime (Reynolds number 25 to 285). We find the population density to be well described by a Lorentzian distribution. We observe this apparently universal form both in the symmetric wake regime and in the more complex three dimensional wake structure of the steady oblique regime with Reynolds number larger than 225. The population density distribution identifies the increase in dimensionless kinetic energy and scalar fluxes with the increase in Reynolds number, whereas the dimensionless scalar population density shows negligible variation with the Reynolds number.

preprint2019arXiv

Boundary Zonal Flow in Rotating Turbulent Rayleigh-Bénard Convection

For rapidly rotating turbulent Rayleigh--Bénard convection in a slender cylindrical cell, experiments and direct numerical simulations reveal a boundary zonal flow (BZF) that replaces the classical large-scale circulation. The BZF is located near the vertical side wall and enables enhanced heat transport there. Although the azimuthal velocity of the BZF is cyclonic (in the rotating frame), the temperature is an anticyclonic traveling wave of mode one whose signature is a bimodal temperature distribution near the radial boundary. The BZF width is found to scale like $Ra^{1/4}Ek^{2/3}$ where the Ekman number $Ek$ decreases with increasing rotation rate.

preprint2019arXiv

Extreme velocity gradients in turbulent flows

Fully turbulent flows are characterized by intermittent formation of very localized and intense velocity gradients. These gradients can be orders of magnitude larger than their typical value and lead to many unique properties of turbulence. Using direct numerical simulations of the Navier-Stokes equations with unprecedented small-scale resolution, we characterize such extreme events over a significant range of turbulence intensities, parameterized by the Taylor-scale Reynolds number ($R_λ$). Remarkably, we find the strongest velocity gradients to empirically scale as $τ_K^{-1} R_λ^β$, with $β\approx 0.775 \pm 0.025$, where $τ_K$ is the Kolmogorov time scale (with its inverse, $τ_K^{-1}$, being the {r.m.s.} of velocity gradient fluctuations). Additionally, we observe velocity increments across very small distances $r \le η$, where $η$ is the Kolmogorov length scale, to be as large as the {r.m.s.} of the velocity fluctuations. Both observations suggest that the smallest length scale in the flow behaves as $ηR_λ^{-α}$, with $α= β- \frac{1}{2}$, which is at odds with predictions from existing phenomenological theories. We find that extreme gradients are arranged in vortex tubes, such that strain conditioned on vorticity grows on average slower than vorticity, approximately as a power law with an exponent $γ< 1$, which weakly increases with $R_λ$. Using scaling arguments, we get $β=(2-γ)^{-1}$, which suggests that $β$ would also slowly increase with $R_λ$. We conjecture that approaching the limit of infinite $R_λ$, the flow is overall smooth, with intense velocity gradients over scale $ ηR_λ^{-1/2}$, corresponding to $β= 1$.

preprint2019arXiv

Flagella-like beating of a single microtubule

Kinesin motors can induce a buckling instability in a microtubule with a fixed minus end. Here we show that by modifying the surface with a protein-repellent functionalization and using clusters of kinesin motors, the microtubule can exhibit persistent oscillatory motion, resembling the beating of sperm flagella. The observed period is of the order of 1 min. From the experimental images we theoretically determine a distribution of motor forces that explains the observed shapes using a maximum likelihood approach. A good agreement is achieved with a small number of motor clusters acting simultaneously on a microtubule. The tangential forces exerted by a cluster are mostly in the range 0 - 8 pN towards the microtubule minus end, indicating the action of 1 or 2 kinesin motors. The lateral forces are distributed symmetrically and mainly below 10 pN, while the lateral velocity has a strong peak around zero. Unlike well-known models for flapping filaments, kinesins are found to have a strong "pinning" effect on the beating filaments. Our results suggest new strategies to utilize molecular motors in dynamic roles that depend sensitively on the stress built-up in the system.

preprint2016arXiv

Azimuthal diffusion of the large-scale-circulation plane, and absence of significant non-Boussinesq effects, in turbulent convection near the ultimate-state transition

We present measurements of the orientation $θ_0$ and temperature amplitude $δ$ of the large-scale circulation in a cylindrical sample of turbulent Rayleigh-Benard convection (RBC) with aspect ratio $Γ\equiv D/L = 1.00$ ($D$ and $L$ are the diameter and height respectively) and for the Prandtl number $Pr \simeq 0.8$. Results for $θ_0$ revealed a preferred orientation with upflow in the West, consistent with a broken azimuthal invariance due to Earth's Coriolis force [see \cite{BA06b}]. They yielded the azimuthal diffusivity $D_θ$ and a corresponding Reynolds number $Re_θ$ for Rayleigh numbers over the range $2\times 10^{12} < Ra < 1.5\times 10^{14}$. In the classical state ($Ra < 2\times 10^{13}$) the results were consistent with the measurements by \cite{BA06a} for $Ra < 10^{11}$ and $Pr = 4.38$ which gave $Re_θ \propto Ra^{0.28}$, and with the Prandtl-number dependence $Re_θ \propto Pr^{-1.2}$ as found previously also for the velocity-fluctuation Reynolds number $Re_V$ \cite[]{HGBA15b}. At larger $Ra$ the data for $Re_θ(Ra)$ revealed a transition to a new state, known as the "ultimate" state, which was first seen in the Nusselt number $Nu(Ra)$ and in $Re_V(Ra)$ at $Ra^*_1 \simeq 2\times 10^{13}$ and $Ra^*_2 \simeq 8\times 10^{13}$. In the ultimate state we found $Re_θ \propto Ra^{0.40\pm 0.03}$. Recently \cite{SU15} claimed that non-Oberbeck-Boussinesq effects on the Nusselt and Reynolds numbers of turbulent RBC may have been interpreted erroneously as a transition to a new state. We demonstrate that their reasoning is incorrect and that the transition observed in the Göttingen experiments and discussed in the present paper is indeed to a new state of RBC referred to as "ultimate".

preprint2016arXiv

Extreme fluctuations of the relative velocities between droplets in turbulent airflow

We compare experiments and direct numerical simulations to evaluate the accuracy of the Stokes-drag model, which is used widely in studies of inertial particles in turbulence. We focus on statistics at the dissipation scale and on extreme values of relative particle velocities for moderately inertial particles (St < 1). The probability distributions of relative velocities in the simulations were qualitatively similar to those in the experiments. The agreement improved with increasing Stokes number and decreasing relative velocity. Simulations underestimated the probability of extreme events, which suggests that the Stokes drag model misses important dynamics. Nevertheless, the scaling behavior of the extreme events in both the experiments and the simulations can be captured by the same multi-fractal model.

preprint2016arXiv

Spatio-temporal Patterns in Inclined Layer Convection

This paper reports on a theoretical analysis of the rich variety of spatio-temporal patterns observed recently in inclined layer convection at medium Prandtl number when varying the inclination angle $γ$ and the Rayleigh number $R$. The present numerical investigation of the inclined layer convection system is based on the standard Oberbeck-Boussinesq equations. The patterns are shown to originate from a complicated competition of buoyancy-driven and shear-flow driven pattern forming mechanisms. The former are expressed as \rm{longitudinal} convection rolls with their axes oriented parallel to the incline, the latter as perpendicular \rm{transverse} rolls. Along with conventional methods to study roll patterns and their stability, we employ direct numerical simulations in large spatial domains, comparable with the experimental ones. As a result, we determine the phase diagram of the characteristic complex 3D convection patterns above onset of convection in the $γ-R$ plane, and find that it compares very well with the experiments. In particular we demonstrate that interactions of specific Fourier modes, characterized by a resonant interaction of their wavevectors in the layer plane, are key to understanding the pattern morphologies.

preprint2016arXiv

Spontaneous concentrations of solids through two-way drag forces between gas and sedimenting particles

The behaviour of sedimenting particles depends on the dust-to-gas ratio of the fluid. Linear stability analysis shows that solids settling in the Epstein drag regime would remain homogeneously distributed in non-rotating incompressible fluids, even when dust-to-gas ratios reach unity. However, the non-linear evolution has not been probed before. Here, we present numerical calculations indicating that in a particle-dense mixture solids spontaneously mix out of the fluid and form swarms overdense in particles by at least a factor 10. The instability is caused by mass-loaded regions locally breaking the equilibrium background stratification. The driving mechanism depends on non-linear perturbations of the background flow and shares some similarity to the streaming instability in accretion discs. The resulting particle-rich swarms may stimulate particle growth by coagulation. In the context of protoplanetary discs, the instability could be relevant for aiding small particles to settle to the midplane in the outer disc. Inside the gas envelopes of protoplanets, enhanced settling may lead to a reduced dust opacity, which facilitates the contraction of the envelope. We show that the relevant physical set up can be recreated in a laboratory setting. This will allow our numerical calculations to be investigated experimentally in the future.

preprint2015arXiv

Irreversibility and small-scale generation in 3D turbulent flows

In three-dimensional turbulent flows energy is supplied at large scales and cascades down to the smallest scales where viscosity dominates. The flux of energy through scales implies the generation of small scales from larger ones, which is the fundamental reason for the irreversibility of the dynamics of turbulent flows. As we showed recently, this irreversibility manifests itself by an asymmetry of the probability distribution of the instantaneous power $p$ of the forces acting on fluid elements. In particular, the third moment of $p$ was found to be negative. Yet, a physical connection between the irreversibility manifested in the distribution of $p$ and the energy flux or small-scale generation in turbulence has not been established. Here, with analytical calculations and support from numerical simulations of fully developed turbulence, we connect the asymmetry in the power distribution, {\it i.e.}, the negative value of $\langle p^3 \rangle$, to the generation of small scales, or more precisely, to the amplification (stretching) of vorticity in turbulent flows. Our result is the first step towards a quantitative understanding of the origin of the irreversibility observed at the level of individual Lagrangian trajectories in turbulent flows.

preprint2015arXiv

Lagrangian view of time irreversibility of fluid turbulence

A turbulent flow is maintained by an external supply of kinetic energy, which is eventually dissipated into heat at steep velocity gradients. The scale at which energy is supplied greatly differs from the scale at which energy is dissipated, the more so as the turbulent intensity (the Reynolds number) is larger. The resulting energy flux over the range of scales, intermediate between energy injection and dissipation, acts as a source of time irreversibility. As it is now possible to follow accurately fluid particles in a turbulent flow field, both from laboratory experiments and from numerical simulations, a natural question arises: how do we detect time irreversibility from these Lagrangian data? Here we discuss recent results concerning this problem. For Lagrangian statistics involving more than one fluid particle, the distance between fluid particles introduces an intrinsic length scale into the problem. The evolution of quantities dependent on the relative motion between these fluid particles, including the kinetic energy in the relative motion, or the configuration of an initially isotropic structure can be related to the equal-time correlation functions of the velocity field, and is therefore sensitive to the energy flux through scales, hence to the irreversibility of the flow. In contrast, for single-particle Lagrangian statistics, the most often studied velocity structure functions cannot distinguish the "arrow of time." Recent observations from experimental and numerical simulation data, however, show that the change of kinetic energy following the particle motion, is sensitive to time-reversal. We end the survey with a brief discussion of the implication of this line of work.

preprint2015arXiv

Turbulence attenuation by large neutrally buoyant particles

Turbulence modulation by inertial-range-size, neutrally-buoyant particles is investigated experimentally in a von Kármán flow. Increasing the particle volume fraction $Φ_\mathrm{v}$, maintaining constant impellers Reynolds number attenuates the fluid turbulence. The inertial-range energy transfer rate decreases as $\proptoΦ_\mathrm{v}^{2/3}$, suggesting that only particles located on a surface affect the flow. Small-scale turbulent properties, such as structure functions or acceleration distribution, are unchanged. Finally, measurements hint at the existence of a transition between two different regimes occurring when the average distance between large particles is of the order of the thickness of their boundary layers.

preprint2015arXiv

Two-Particle Dispersion in Weakly Turbulent Thermal Convection

We present results from a numerical study of particle dispersion in the weakly nonlinear regime of Rayleigh-Bénard convection of a fluid with Prandtl number around unity, where bi-stability between ideal straight convection rolls and weak turbulence in the form of Spiral Defect Chaos exists. While Lagrangian pair statistics has become a common tool for studying fully developed turbulent flows at high Reynolds numbers, we show that key characteristics of mass transport can also be found in convection systems that show no or weak turbulence. Specifically, for short times, we find Richardson-like $t^3$ scaling of pair dispersion. We explain our findings quantitatively with a model that captures the interplay of advection and diffusion. For long times we observe diffusion-like dispersion of particles that becomes independent of the individual particles' stochastic movements. The spreading rate is found to depend on the degree of spatio-temporal chaos.

preprint2014arXiv

Logarithmic Spatial Variations and Universal $f^{-1}$ Power Spectra of Temperature Fluctuations in Turbulent Rayleigh-Bénard Convection

We report measurements of the temperature variance $σ^2(z,r)$ and frequency power spectrum $P(f,z,r)$ ($z$ is the distance from the sample bottom and $r$ the radial coordinate) in turbulent Rayleigh-Bénard convection (RBC) for Rayleigh numbers $\textrm{Ra} = 1.6\times10^{13}$ and $1.1\times10^{15}$ and for a Prandtl number $\textrm{Pr} \simeq 0.8$ for a sample with a height $L = 224$ cm and aspect ratio $D/L = 0.50$ ($D$ is the diameter). For $z/L$ less than or similar to $0.1$ $σ^2(z,r)$ was consistent with a logarithmic dependence on $z$, and there was a universal (independent of $\textrm{Ra}$, $r$, and $z$) normalized spectrum which, for $0.02$ less than or similar to $fτ_0$ less than or similar to $0.2$, had the form $P(fτ_0) = P_0 (fτ_0)^{-1}$ with $P_0 =0.208 \pm 0.008$ a universal constant. Here $τ_0 = \sqrt{2R}$ where $R$ is the radius of curvature of the temperature autocorrelation function $C(τ)$ at $τ= 0$. For $z/L \simeq 0.5$ the measurements yielded $P(fτ_0) \sim (fτ_0)^{-α}$ with $α$ in the range from 3/2 to 5/3. All the results are similar to those for velocity fluctuations in shear flows at sufficiently large Reynolds numbers, suggesting the possibility of an analogy between the flows that is yet to be determined in detail.

preprint2014arXiv

Logarithmic temperature profiles of turbulent Rayleigh-Bénard convection in the classical and ultimate state for a Prandtl number of 0.8

We report on experimental determinations of the temperature field in the interior (bulk) of turbulent Rayleigh-Benard convection for a cylindrical sample with aspect ratio (diameter over height) of 0.50, both in the classical and in the ultimate state. The Prandtl number was close to 0.8. We find a "logarithmic layer" in which the temperature varies as A*ln(z/L) + B with the distance z from the bottom plate of the sample. The amplitude A varies with radial position r. In the classical state these results are in good agreement with direct numerical simulations (DNS); in the ultimate state there are as yet no DNS. A close analogy between the temperature field in the classical state and the "Law of the Wall" for the time-averaged down-stream velocity in shear flow is discussed.

preprint2014arXiv

Observation of the sling effect

When cloud particles are small enough, they move with the turbulent air in the cloud. On the other hand, as particles become larger their inertia affects their motions, and they move differently than the air. These inertial dynamics impact cloud evolution and ultimately climate prediction, since clouds govern the earth's energy balances. Yet we lack a simple description of the dynamics. Falkovich et al. describes theoretically a new dynamical mechanism called the "sling effect" by which extreme events in the turbulent air cause idealized inertial cloud particles to break free from the airflow (Falkovich G, Fouxon A, Stepanov MG 2002 Nature 419, 151). The sling effect thereafter causes particle trajectories to cross each other within isolated pockets in the flow, which increases the chance of collisions that form larger particles. We combined experimental techniques that allow for precise control of a turbulent flow with three-dimensional tracking of multiple particles at unprecedented resolution. In this way, we could observe both the sling effect and crossing trajectories between real particles. We isolated the inertial sling dynamics from those caused by turbulent advection by conditionally averaging the data. We found the dynamics to be universal in terms of a local Stokes number that quantifies the local particle velocity gradients. We measured the probability density of this quantity, which shows that sharp gradients become more frequent as the global Stokes number increases. We observed that sharp compressive gradients in the airflow initiated the sling effect, and that thereafter gradients in the particle flow ran away and steepened in a way that produced singularities in the flow in finite time. During this process both the fluid motions and gravity became unimportant. The results underpin a framework for describing a crucial aspect of inertial particle dynamics.

preprint2014arXiv

On the Relaxation of Turbulence at High Reynolds Numbers

Turbulent motions in a fluid relax at a certain rate once stirring has stopped. The role of the most basic parameter in fluid mechanics, the Reynolds number, in setting the relaxation rate is not generally known. This paper concerns the high-Reynolds-number limit of the process. In a classical grid-turbulence wind-tunnel experiment that both reached higher Reynolds numbers than ever before and covered a wide range of them ($10^4 < Re = UM/ν< 5\times10^6$), we measured the relaxation rate with the unprecedented precision of about 2\%. Here $U$ is the mean speed of the flow, $M$ the forcing scale, and $ν$ the kinematic viscosity of the fluid. We observed that the relaxation rate was Reynolds-number independent, which contradicts some models and supports others.

preprint2014arXiv

Power Fluctuations and Irreversibility in Turbulence

The breaking of detailed balance, the symmetry between forward and backward probability transition between two states, is crucial to understand irreversible systems. In hydrodynamic turbulence, a far-from equilibrium system, we observe a strong manifestation of the breaking of detailed balance by following the evolution of the kinetic energy of individual fluid elements. We found in all the flows that we have investigated that fluid elements decelerate faster than they accelerate, giving rise to negative third moment of energy increments, independently of space dimensionality. The exchange of energy between fluid elements however is fundamentally different in two and three dimensions. While pressure forces do not provide net energy change to slow or fast particles in two dimensions, they tend to transfer energy from {\it slow} to {\it fast} particles in three dimensions, possibly implying a runaway of energy.

preprint2014arXiv

Time-symmetry breaking in turbulence

In three-dimensional turbulent flows, the flux of energy from large to small scales breaks time symmetry. We show here that this irreversibility can be quantified by following the relative motion of several Lagrangian tracers. We find by analytical calculation, numerical analysis and experimental observation that the existence of the energy flux implies that, at short times, two particles separate temporally slower forwards than backwards, and the difference between forward and backward dispersion grows as $t^3$. We also find the geometric deformation of material volumes, surrogated by four points spanning an initially regular tetrahedron, to show sensitivity to the time-reversal with an effect growing linearly in $t$. We associate this with the structure of the strain rate in the flow.

preprint2014arXiv

Variable Density Turbulence Tunnel Facility

The Variable Density Turbulence Tunnel (VDTT) at the Max Planck Institute for Dynamics and Self-Organization in Göttingen, Germany produces very high turbulence levels at moderate flow velocities, low power consumption and adjustable kinematic viscosity between $10^{-4} m^2/s$ and $10^{-7} m^2/s$. The Reynolds number can be varied by changing the pressure or flow rate of the gas or by using different non-flammable gases including air. The highest kinematic viscosities, and hence lowest Reynolds numbers, are reached with air or nitrogen at 0.1 bar. To reach the highest Reynolds numbers the tunnel is pressurized to 15 bar with the dense gas sulfur hexafluoride (SF$_6$). Turbulence is generated at the upstream ends of two measurement sections with grids, and the evolution of this turbulence is observed as it moves down the length of the sections. We describe the instrumentation presently in operation, which consists of the tunnel itself, classical grid turbulence generators, and state-of-the-art nano-fabricated hot-wire anemometers provided by Princeton University [Vallikivi et al. (2011) Exp. Fluids 51, 1521]. We report measurements of the characteristic scales of the flow and of turbulent spectra up to Taylor Reynolds number $R_λ\approx 1600$, higher than any other grid-turbulence experiment. We also describe instrumentation under development, which includes an active grid and a Lagrangian particle tracking system that moves down the length of the tunnel with the mean flow. In this configuration, the properties of the turbulence are adjustable and its structure is resolvable up to $R_λ\approx 8000$.

preprint2013arXiv

Tetrahedron deformation and alignment of perceived vorticity and strain in a turbulent flow

We describe the structure and dynamics of turbulence by the scale-dependent perceived velocity gradient tensor as supported by following four tracers, i.e. fluid particles, that initially form a regular tetrahedron. We report results from experiments in a von Kármán swirling water flow and from numerical simulations of the incompressible Navier-Stokes equation. We analyze the statistics and the dynamics of the perceived rate of strain tensor and vorticity for initially regular tetrahedron of size $r_0$ from the dissipative to the integral scale. Just as for the true velocity gradient, at any instant, the perceived vorticity is also preferentially aligned with the intermediate eigenvector of the perceived rate of strain. However, in the perceived rate of strain eigenframe fixed at a given time $t=0$, the perceived vorticity evolves in time such as to align with the strongest eigendirection at $t=0$. This also applies to the true velocity gradient. The experimental data at the higher Reynolds number suggests the existence of a self-similar regime in the inertial range. In particular, the dynamics of alignment of the perceived vorticity and strain can be rescaled by $t_0$, the turbulence time scale of the flow when the scale $r_0$ is in the inertial range. For smaller Reynolds numbers we found the dynamics to be scale dependent.

preprint2013arXiv

Turbulence of Dilute Polymer Solution

In fully developed three dimensional fluid turbulence the fluctuating energy is supplied at large scales, cascades through intermediate scales, and dissipates at small scales. It is the hallmark of turbulence that for intermediate scales, in the so called inertial range, the average energy flux is constant and independent of viscosity [1-3]. One very important question is how this range is altered, when an additional agent that can also transport energy is added to the fluid. Long-chain polymers dissolved at very small concentrations in the fluid are such an agent [4,5]. Based on prior work by de Gennes and Tabor [6,7] we introduce a theory that balances the energy flux through the turbulent cascade with that of the energy flux into the elastic degrees of freedom of the dilute long-chain polymer solution. We propose a refined elastic length scale, $r_\varepsilon$, which describes the effect of polymer elasticity on the turbulence energy cascade. Our experimental results agree excellently with this new energy flux balance theory.

preprint2012arXiv

Heat transport by turbulent Rayleigh-Bénard convection for $\Pra\ \simeq 0.8$ and $3\times 10^{12} \alt \Ra\ \alt 10^{15}$: Aspect ratio $Γ= 0.50$

We report experimental results for heat-transport measurements, in the form of the Nusselt number \Nu, by turbulent Rayleigh-Bénard convection in a cylindrical sample of aspect ratio $Γ\equiv D/L = 0.50$ ($D = 1.12$ m is the diameter and $L = 2.24$ m the height). The measurements were made using sulfur hexafluoride at pressures up to 19 bars as the fluid. They are for the Rayleigh-number range $3\times 10^{12} \alt \Ra \alt 10^{15}$ and for Prandtl numbers \Pra\ between 0.79 and 0.86. For $\Ra < \Ra^*_1 \simeq 1.4\times 10^{13}$ we find $\Nu = N_0 \Ra^{γ_{eff}}$ with $γ_{eff} = 0.312 \pm 0.002$, consistent with classical turbulent Rayleigh-Bénard convection in a system with laminar boundary layers below the top and above the bottom plate. For $\Ra^*_1 < \Ra < \Ra^*_2$ (with $\Ra^*_2 \simeq 5\times 10^{14}$) $γ_{eff}$ gradually increases up to $0.37\pm 0.01$. We argue that above $\Ra^*_2$ the system is in the ultimate state of convection where the boundary layers, both thermal and kinetic, are also turbulent. Several previous measurements for $Γ= 0.50$ are re-examined and compared with the present results.

preprint2012arXiv

Heat transport by turbulent Rayleigh-Bénard convection for $\Pra\ \simeq 0.8$ and $4\times 10^{11} \alt \Ra\ \alt 2\times10^{14}$: Ultimate-state transition for aspect ratio $Γ= 1.00$

We report experimental results for heat-transport measurements by turbulent Rayleigh-Bénard convection in a cylindrical sample of aspect ratio $Γ\equiv D/L = 1.00$ ($D = 1.12$ m is the diameter and $L = 1.12$ m the height). They are for the Rayleigh-number range $4\times10^{11} \alt \Ra \alt 2\times10^{14}$ and for Prandtl numbers \Pra\ between 0.79 and 0.86. For $\Ra < \Ra^*_1 \simeq 2\times 10^{13}$ we find $\Nu = N_0 \Ra^{γ_{eff}}$ with $γ_{eff} = 0.321 \pm 0.002$ and $N_0 = 0.0776$, consistent with classical turbulent Rayleigh-Bénard convection in a system with laminar boundary layers below the top and above the bottom plate and with the prediction of Grossmann and Lohse. For $\Ra > \Ra_1^*$ the data rise above the classical-state power-law and show greater scatter. In analogy to similar behavior observed for $Γ= 0.50$, we interpret this observation as the onset of the transition to the ultimate state. Within our resolution this onset occurs at nearly the same value of $\Ra_1^*$ as it does for $Γ= 0.50$. This differs from an earlier estimate by Roche {\it et al.} which yielded a transition at $\Ra_U \simeq 1.3\times 10^{11} Γ^{-2.5\pm 0.5}$. A $Γ$-independent $\Ra^*_1$ would suggest that the boundary-layer shear transition is induced by fluctuations on a scale less than the sample dimensions rather than by a global $Γ$-dependent flow mode. Within the resolution of the measurements the heat transport above $\Ra_1^*$ is equal for the two $Γ$ values, suggesting a universal aspect of the ultimate-state transition and properties. The enhanced scatter of \Nu\ in the transition region, which exceeds the experimental resolution, indicates an intrinsic irreproducibility of the state of the system.

preprint2012arXiv

Logarithmic temperature profiles in turbulent Rayleigh-Bénard convection

We report results for the temperature profiles of turbulent Rayleigh-Bénard convection (RBC) in the interior of a cylindrical sample of aspect ratio $Γ\equiv D/L = 0.50$ ($D$ and $L$ are the diameter and height respectively). Results from experiment over the Rayleigh number range $4\times 10^{12} \alt Ra \alt 10^{15}$ for a Prandtl number $\Pra \simeq 0.8$ and from direct numerical simulation (DNS) at $Ra = 2 \times 10^{12}$ for $\Pra = 0.7$ are presented. We find that the temperature varies as $A*ln(z/L) + B$ where $z$ is the distance from the bottom or top plate. This is the case in the classical as well as in the ultimate state of RBC. From DNS we find that $A$ in the classical state decreases in the radial direction as the distance from the side wall increases and becomes small near the sample center.

preprint2012arXiv

Simultaneous 3D measurement of the translation and rotation of finite size particles and the flow field in a fully developed turbulent water flow

We report a novel experimental technique that measures simultaneously in three dimensions the trajectories, the translation, and the rotation of finite size inertial particles together with the turbulent flow. The flow field is analyzed by tracking the temporal evolution of small fluorescent tracer particles. The inertial particles consist of a super-absorbent polymer that renders them index and density matched with water and thus invisible. The particles are marked by inserting at various locations tracer particles into the polymer. Translation and rotation, as well as the flow field around the particle are recovered dynamically from the analysis of the marker and tracer particle trajectories. We apply this technique to study the dynamics of inertial particles much larger in size (Rp/η \approx 100) than the Kolmogorov length scale η in a von Kármán swirling water flow (Rλ \approx 400). We show, using the mixed (particle/fluid) Eulerian second order velocity structure function, that the interaction zone between the particle and the flow develops in a spherical shell of width 2Rp around the particle of radius Rp. This we interpret as an indication of a wake induced by the particle. This measurement technique has many additional advantages that will make it useful to address other problems such as particle collisions, dynamics of non-spherical solid objects, or even of wet granular matter.

preprint2012arXiv

Where do small, weakly inertial particles go in a turbulent flow?

We report experimental results on the dynamics of heavy particles of the size of the Kolmogorov-scale in a fully developed turbulent flow. The mixed Eulerian structure function of two-particle velocity and acceleration difference vectors <δv\cdotδa_p> was observed to increase significantly with particle inertia for identical flow conditions. We show that this increase is related to a preferential alignment between these dynamical quantities. With increasing particle density the probability for those two vectors to be collinear was observed to grow. We show that these results are consistent with the preferential sampling of strain-dominated regions by inertial particles.

preprint2011arXiv

Experimental study of the influence of anisotropy on the inertial scales of turbulence

We ask whether the scaling exponents or the Kolmogorov constants depend on the anisotropy of the velocity fluctuations in a turbulent flow with no shear. According to our experiment, the answer is no for the Eulerian second-order transverse velocity structure function. The experiment consisted of 32 loudspeaker-driven jets pointed toward the centre of a spherical chamber. We generated anisotropy by controlling the strengths of the jets. We found that the form of the anisotropy of the velocity fluctuations was the same as that in the strength of the jets. We then varied the anisotropy, as measured by the ratio of axial to radial root-mean-square (RMS) velocity fluctuations, between 0.6 and 2.3. The Reynolds number was approximately constant at around $R_λ$ = 481. In a central volume with a radius of 50 mm, the turbulence was approximately homogeneous, axisymmetric, and had no shear and no mean flow. We observed that the scaling exponent of the structure function was $0.70 \pm 0.03$, independent of the anisotropy and regardless of the direction in which we measured it. The Kolmogorov constant, $C_2$, was also independent of direction and anisotropy to within the experimental error of 4%.

preprint2011arXiv

Hydrodynamics of Hypersonic Jets: Experiments and Numerical Simulations

Stars form in regions of the galaxy that are denser and cooler than the mean interstellar medium. These regions are called Giant Molecular Clouds. At the beginning of their life, up to $10^5-10^6$ years, stars accrete matter from their rich surrounding environment and are origin of a peculiar phenomenon that is the jet emission. Jets from Young Stellar Objects (YSOs) are intensively studied by the astrophysical community by observations at different wavelengths, analytical and numerical modeling and laboratory experiments. Indications about the jet propagation and its resulting morphologies are here obtained by means of a combined study of hypersonic jets carried out both in the laboratory and by numerical simulations.

preprint2011arXiv

Path lengths in turbulence

By tracking tracer particles at high speeds and for long times, we study the geometric statistics of Lagrangian trajectories in an intensely turbulent laboratory flow. In particular, we consider the distinction between the displacement of particles from their initial positions and the total distance they travel. The difference of these two quantities shows power-law scaling in the inertial range. By comparing them with simulations of a chaotic but non-turbulent flow and a Lagrangian Stochastic model, we suggest that our results are a signature of turbulence.

preprint2011arXiv

Prandtl and the Göttingen School

This is a historical review of "Prandtl and the Göttingen School" from his first work in 1898 to after WWII. We report the development of the boundary layer theory and the investigations that lead Ludwig Prandtl in the spring of 1945 to use a energy cascade model to derive for fully developed turbulence the Kolmogrov length scale. Working at the end of WWII he was not aware of the work by Kolmogorov. Thus Ludwig Prandtl was before Onsager, von Weizäcker, and Heisenberg in rederiving parts of K41 theory. In our historical analysis we rely heavily on exchanges between Ludwig Prandtl with Theodore von Kármán and Goeffrey Ingram Taylor and on the original notes by Ludwig Prandtl and others.

preprint2010arXiv

On the Swimming of \textit{Dictyostelium} amoebae

Traditionally, the primary mode for locomotion of amoeboid cells was thought to be crawling on a substrate. Recently, it has been experimentally shown that \textit{Dictostelium} amoeba and neutrophils can also swim in a directed fashion. The mechanisms for amoeboid crawling and swimming were hypothesized to be similar. In this letter, we show that the shape changes generated by a crawling \textit{D. discoideum} cell are consistent with swimming.

preprint2007arXiv

Competition and bistability of ordered undulations and undulation chaos in inclined layer convection

Experimental and theoretical investigations of undulation patterns in high-pressure, inclined layer gas convection at a Prandtl number near unity are reported. Particular focus is given to the competition between the spatiotemporal chaotic state of undulation chaos and stationary patterns of ordered undulations. In experiments a competition and bistability between the two states is observed, with ordered undulations most prevalent at higher Rayleigh number. The spectral pattern entropy, spatial correlation lengths, and defect statistics are used to characterize the competing states. The experiments are complemented by a theoretical analysis of the Oberbeck-Boussinesq equations. The stability region of the ordered undulation as a function of their wavevectors and the Rayleigh number is obtained with Galerkin techniques. In addition, direct numerical simulations are used to investigate the spatiotemporal dynamics. In the simulations both ordered undulations and undulation chaos were observed dependent on initial conditions. Experiment and theory are found to agree well.

preprint1997arXiv

A New Fast Method for Determining Local Properties of Striped Patterns

From the striped coats of zebras to the ripples in windblown sand, the natural world abounds with locally banded patterns. Such patterns have been of great interest throughout history, and, in the last twenty years, scientists in a wide variety of fields have been studying the patterns formed in well-controlled experiments that yield enormous quantities of high-precision data. These experiments involving phenomena as diverse as chemical reactions in shallow layers, thermal convection in horizontal fluid layers, periodically shaken layers of sand, and the growth of slime mold colonies often display patterns that appear qualitatively similar. Methods are needed to characterize in a reasonable amount of time the differences and similarities in patterns that develop in different systems, as well as in patterns formed in one system for different experimental conditions. In this Letter, we introduce a novel, fast method for determining local pattern properties such as wavenumber, orientation, and curvature as a function of position for locally striped patterns.

preprint1996arXiv

Core Dynamics of Multi-Armed Spirals in Rayleigh-Benard Convection

We report experimental observations of the core dynamics of multi-armed, rotating spirals in Rayleigh-Benard convection for a fluid with Prandtl number near one. In addition to the large-scale rotation of the spirals, we found a cyclic core motion within a central area of radius r = appr. (n lambda/2), where n is the number of spiral arms ending in the core and lambda is the wavelength of the pattern. The dynamics of the core was much faster than the large-scale spiral rotation. We observed multi-armed spirals for which the two periods were commensurate and others for which they were incommensurate.

preprint1996arXiv

The Spatio-Temporal Structure of Spiral-Defect Chaos

We present a study of the recently discovered spatially-extended chaotic state known as spiral-defect chaos, which occurs in low-Prandtl-number, large-aspect-ratio Rayleigh-Benard convection. We employ the modulus squared of the space-time Fourier transform of time series of two-dimensional shadowgraph images to construct the structure factor ${S}({\vec k},ω)$. This analysis is used to characterize the average spatial and temporal scales of the chaotic state. We find that the correlation length and time can be described by power-law dependences on the reduced Rayleigh number $ε$. These power laws have as yet no theoretical explanation.

preprint1995arXiv

Defect--Defect Correlation Functions, Generic Scale Invariance, and the Complex Ginzburg--Landau Equation

Motivated by the idea of developing a ``hydrodynamic'' description of spatiotemporal chaos, we have investigated the defect--defect correlation functions in the defect turbulence regime of the two--dimensional, anisotropic complex Ginzburg--Landau equation. We compare our results with the predictions of generic scale invariance. Using the topological nature of the defects, we prove that defect--defect correlations cannot decay as slowly as predicted by generic scale invariance. We also present results on the fluctuations of the amplitude field $A$.