Search

Particle tracking in turbulent flows is fundamental to the study of the transport of tracers, inertial particles or even active objects in space and time, i.e. the Lagrangian frame of reference. It provides experimental tests of theoretical predictions (e.g. for the statistics of fluid accelerations and particle dispersion) and helps to understand important natural processes where particle inertia is important (e.g. cloud microphysics). While the spatial (Eulerian) properties of turbulent flows have been studied for high, atmospheric Reynolds numbers ($R_\lambda > 10^4$), the profound difficulties in accurately tracking particles in turbulent flows have limited the Reynolds numbers in the Lagrangian reference frame to the Taylor scale Reynolds numbers $R_\lambda \lesssim 10^3$. Here we describe a setup that allowed Lagrangian particle tracking at $R_\lambda$ between 100 and 6000 in the Max Planck Variable Density Turbulence Tunnel (VDTT). We describe the imaging setup within the pressurised facility, the laser illumination, the particles and the particle dispersion mechanism. We verify that the KOBO Cellulobeads D-10 particles are suitable tracers. They carry negligible charge and their Stokes number is small over the full range of experimental conditions. We present typical data from the experiment and discuss the challenges and constraints of the setup.

Cells often migrate on curved surfaces inside the body, such as curved tissues, blood vessels or highly curved protrusions of other cells. Recent \textit{in-vitro} experiments provide clear evidence that motile cells are affected by the curvature of the substrate on which they migrate, preferring certain curvatures to others, termed ``curvotaxis". The origin and underlying mechanism that gives rise to this curvature sensitivity are not well understood. Here, we employ a ``minimal cell" model which is composed of a vesicle that contains curved membrane protein complexes, that exert protrusive forces on the membrane (representing the pressure due to actin polymerization). This minimal-cell model gives rise to spontaneous emergence of a motile phenotype, driven by a lamellipodia-like leading edge. By systematically screening the behaviour of this model on different types of curved substrates (sinusoidal, cylinder and tube), we show that minimal ingredients and energy terms capture the experimental data. The model recovers the observed migration on the sinusoidal substrate, where cells move along the grooves (minima), while avoiding motion along the ridges. In addition, the model predicts the tendency of cells to migrate circumferentially on convex substrates and axially on concave ones. Both of these predictions are verified experimentally, on several cell types. Altogether, our results identify the minimization of membrane-substrate adhesion energy and binding energy between the membrane protein complexes as key players of curvotaxis in cell migration.

Left ventricular (LV) stiffness and contractility, characterized by the end-diastolic and end-systolic pressure-volume relationships (EDPVR & ESPVR), are two important indicators of the performance of the human heart. Although much research has been conducted on EDPVR and ESPVR, no model with physically interpretable parameters combining both relationships has been presented, thereby impairing the understanding of cardiac physiology and pathology. Here, we present a model that evaluates both EDPVR and ESPVR with physical interpretations of the parameters in a unified framework. Our physics-based model fits the available experimental data and in silico results very well and outperforms existing models. With prescribed parameters, the new model is used to predict the pressure-volume relationships of the left ventricle. Our model provides a deeper understanding of cardiac mechanics and thus will have applications in cardiac research and clinical medicine., Comment: 14 pages, 8 figures

In windy conditions, the air is turbulent. The strong and intermittent velocity variations of turbulence are invisible to flying animals. Nevertheless, flying animals, not much larger than the smallest scales of turbulence, manage to maneuver these highly fluctuating conditions quite well. Here we quantify honeybee flight with time-resolved three-dimensional tracking in calm conditions and controlled turbulent winds. We find that honeybee mean speed and acceleration are only weakly correlated with the strength of turbulence. In flight, honeybees accelerate slowly and decelerate rapidly, i.e., they break suddenly during turns and then accelerate again. While this behavior is observed in both calm and turbulent conditions, it is increasingly dominant under turbulent conditions where short straight trajectories are broken by turns and increased maneuvering. This flight-crash behavior is reminiscent of turbulence itself. Our observations may help the development of flight strategies for miniature flying robotics under turbulent conditions., Comment: 18 pages, 7 figures

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 $\mu$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., Comment: 46 pages, 13 figures, 1 PDF with Supplementary Information

In this position paper, a large group of interdisciplinary experts outlines response strategies against the spread of SARS-CoV-2 in the winter of 2021/2022 in Germany. We review the current state of the COVID-19 pandemic, from incidence and vaccination efficacy to hospital capacity. Building on this situation assessment, we illustrate various possible scenarios for the winter, and detail the mechanisms and effectiveness of the non-pharmaceutical interventions, vaccination, and booster. With this assessment, we want to provide orientation for decision makers about the progress and mitigation of COVID-19., Comment: in German, Extensive expert assessment on COVID-19 response policies for the winter 2021/22

Arrhythmia of the heart is a dangerous and potentially fatal condition. The current widely used treatment is the implantable cardioverter defibrillator (ICD), but it is invasive and affects the patient's quality of life. The sonogenetic treatment technique proposed here focuses transthoracic ultrasound on the heart, noninvasively stimulates endogenous stretch-activated Piezo1 ion channels on the focal region's cardiomyocyte plasma membrane, and restores normal heart rhythm. In contrast to anchoring the implanted ICD lead at a fixed position in the myocardium, the size and position of the ultrasound focal region can be selected dynamically by adjusting the signal phases of every piezoelectric chip on the wearable ultrasonic phased array, and it allows novel and efficient defibrillations. Based on the developed interdisciplinary electromechanical model of sonogenetic treatment, our analysis shows that the proposed ultrasound intensity and frequency will be safe and painless for humans and well below the limits established by the U.S. Food and Drug Administration., Comment: 14 pages, 10 figures

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., Comment: 13 pages, 7 figures

A hallmark of fluid turbulence theory is the universal power law scaling of the velocity difference statistics between two points in space in the inertial range between the large energy injection scale and the small energy dissipation scale. Even at the highest Reynolds numbers available, laboratory and natural flows such universal power laws have not been convincingly demonstrated. Here we show for the decaying active grid turbulence of the Max Planck Variable Density Turbulence Tunnel that the velocity difference statistics at high Reynolds numbers do not exhibit a power law, but have a universal functional form independent of the Reynolds number. We separate this functional form from the power law exponent and discuss potential consequences for turbulence modelling., Comment: arXiv admin note: substantial text overlap with arXiv:2006.10993

We present results on the infection risk from SARS-CoV-2 under different scenarios based on measured particle size-dependent mask penetration, measured total inward leakage, measured human aerosol emission for sizes from 10nm to 1mm, and re-hydration on inhalation. Well-mixed room models significantly underestimate the risk of infection for short and direct exposure. To this end, we estimate the upper bound for infection risk with the susceptible in the infectious exhalation cloud or wearing masks by having the masked susceptible inhale the entire exhalation of a masked infectious. Social distances without a mask, even at 3m between speaking individuals results in an upper bound of 90\% for risk of infection after a few minutes. If both wear a surgical mask, the risk of infection for the person speaking remains below 26\% even after 60 minutes. When both the infectious and susceptible wear a well-fitting FFP2 mask, the upper bound for risk is reduced by a factor of 60 compared to surgical masks. In both cases, face leakage is very important. For FFP2 masks, leakage is low in the nasal region and directed upward, which can be further reduced significantly by applying double-sided medical tape there. Considering that the calculated upper bound greatly overestimates the risk of infection, and the fact that with a poorly worn mask even the upper bound we calculated is very low, we conclude that wearing a mask, even with some leakage, provides excellent third party and self-protection.

Airborne diseases can be transmitted by infectious aerosols in the near field, i.e., in close proximity, or in the far field, i.e., by infectious aerosols that are well mixed within the indoor air. Is it possible to say which mode of disease transmission is predominant in large indoor spaces? We addressed this question by measuring the transport of aerosols equivalent to the size of human respiratory particles in two large hardware stores (V>10000 m$^3$). We found that aerosol concentrations in both stores decreased rapidly and almost independently of aerosol size, despite the different ventilation systems. A persistent and directional airflow on the order of a few cm/s was observed in both stores. Consequently, aerosol dynamics in such open settings can be expected to be dominated by turbulent dispersion and sweeping, and the accumulation of infectious aerosols in the indoor air is unlikely to contribute significantly to the risk of infection as long as the occupancy of the store is not too high. Under these conditions, well-fitting face masks are an excellent means of preventing disease transmission by human aerosols.

Knowing the physicochemical properties of exhaled droplets and aerosol particles is a prerequisite for a detailed mechanistic understanding and effective prevention of the airborne transmission of infectious human diseases. This article provides a critical review and synthesis of scientific knowledge on the number concentrations, size distributions, composition, mixing state, and related properties of respiratory particles emitted upon breathing, speaking, singing, coughing, and sneezing. We derive and present a parameterization of respiratory particle size distributions based on five lognormal modes related to different origins in the respiratory tract, which can be used to trace and localize the sources of infectious particles. This approach may support the medical treatment as well as the risk assessment for aerosol and droplet transmission of infectious diseases. It was applied to analyze which respiratory activities may drive the spread of specific pathogens, such as Mycobacterium tuberculosis, influenza viruses, and SARS-CoV-2 viruses. The results confirm the high relevance of vocalization for the transmission of SARS-CoV-2 as well as the usefulness of physical distancing, face masks, room ventilation, and air filtration as preventive measures against COVID-19 and other airborne infectious diseases.

Propelling microorganisms through fluids and moving fluids along cellular surfaces are essential biological functions accomplished by long, thin structures called motile cilia and flagella, whose regular, oscillatory beating breaks the time-reversal symmetry required for transport. Although top-down experimental approaches and theoretical models have allowed us to broadly characterize such organelles and propose mechanisms underlying their complex dynamics, constructing minimal systems capable of mimicking ciliary beating and identifying the role of each component remains a challenge. Here we report the bottom-up assembly of a minimal synthetic axoneme, which we call a synthoneme, using biological building blocks from natural organisms, namely pairs of microtubules and cooperatively associated axonemal dynein motors. We show that upon provision of energy by ATP, microtubules undergo rhythmic bending by cyclic association-dissociation of dyneins. Our simple and unique beating minimal synthoneme represents a self-organized nanoscale biomolecular machine that can also help understand the mechanisms underlying ciliary beating., Comment: 14 pages, 8 figures

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., Comment: 24pages,5 figures

In the case of airborne diseases, pathogen copies are transmitted by droplets of respiratory tract fluid that are exhaled by the infectious and, after partial or full drying, inhaled as aerosols by the susceptible. The risk of infection in indoor environments is typically modelled using the Wells-Riley model or a Wells-Riley-like formulation, usually assuming the pathogen dose follows a Poisson distribution (mono-pathogen assumption). Aerosols that hold more than one pathogen copy, i.e. poly-pathogen aerosols, break this assumption even if the aerosol dose itself follows a Poisson distribution. For the largest aerosols where the number of pathogen in each aerosol can sometimes be several hundred or several thousand, the effect is non-negligible, especially in diseases where the risk of infection per pathogen is high. Here we report on a generalization of the Wells-Riley model and dose-response models for poly-pathogen aerosols by separately modeling each number of pathogen copies per aerosol, while the aerosol dose itself follows a Poisson distribution. This results in a model for computational risk assessment suitable for mono-/poly-pathogen aerosols. We show that the mono-pathogen assumption significantly overestimates the risk of infection for high pathogen concentrations in the respiratory tract fluid. The model also includes the aerosol removal due to filtering by the individuals which becomes significant for poorly ventilated environments with a high density of individuals, and systematically includes the effects of facemasks in the infectious aerosol source and sink terms and dose calculations., Comment: updated file with link to software on GitHub

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 $\mu$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{\deg}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., Comment: 25 pages, 8 figures, Revised Version

Interruptions in nonlinear wave propagation, commonly referred to as wave breaks, are typical of many complex excitable systems. In the heart they lead to fatal rhythm disorders, the so-called arrhythmias, which are one of the main causes of sudden death in the industrialized world. Progress in the treatment and therapy of cardiac arrhythmias requires a detailed understanding of the triggers and dynamics of these wave breaks. In particular, two very important questions are: 1) What determines the potential of a wave break to initiate re-entry? and 2) How do these breaks evolve such that the system is able to maintain spatiotemporally chaotic electrical activity? Here we approach these questions numerically using optogenetics in an in silico model of human atrial tissue that has undergone chronic atrial fibrillation (cAF) remodelling. In the lesser known sub-threshold illumination r\'egime, we discover a new mechanism of wave break initiation in cardiac tissue that occurs for gentle slopes of the restitution characteristics. This mechanism involves "conditioning" or reshaping the wave profile from front to back, such that, removal of the external light source causes rapid recovery of cells at the waveback, leading to the creation of vulnerable windows for sustained re-entry in spatially extended systems.

Turbulent fluid flows are ubiquitous in nature and technology, and are mathematically described by the incompressible Navier-Stokes equations (INSE). A hallmark of turbulence is spontaneous generation of intense whirls, resulting from amplification of the fluid rotation-rate (vorticity) by its deformation-rate (strain). This interaction, encoded in the non-linearity of INSE, is non-local, i.e., depends on the entire state of the flow, constituting a serious hindrance in turbulence theory and in establishing regularity of INSE. Here, we unveil a novel aspect of this interaction, by separating strain into local and non-local contributions utilizing the Biot-Savart integral of vorticity in a sphere of radius R. Analyzing highly-resolved numerical turbulent solutions to INSE, we find that when vorticity becomes very large, the local strain over small R surprisingly counteracts further amplification. This uncovered self-attenuation mechanism is further shown to be connected to local Beltramization of the flow, and could provide a direction in establishing the regularity of INSE., Comment: 8 pages, 3 figures

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.

The secondary activation of aerosols impacts the life cycle of a cloud. A detailed understanding is necessary for reliable climate prediction. Recent laboratory experiments demonstrate that aerosols can be activated in the wake of precipitating hydrometeors. However, many quantitative aspects of this wake-induced activation remain unclear. Here, we report a detailed numerical investigation of the activation potential of wake-induced supersaturation. By Lagrangian tracking of aerosols we show that a significant fraction of aerosols are activated in the supersaturated wake. These 'lucky aerosols' are entrained in the wake's vortices and reside in the supersaturated environment sufficiently long to be activated. Our analyses show that wake-induced activation can contribute at a level similar to other well known secondary production processes.

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_{\lambda} \lesssim 6000$. This is evidence of a universal inertial range at high Reynolds numbers., Comment: 9 pages, 5 figures updated estimate for integral length scale, updated Fig. 4, updated error on scaling exponent

An essential ingredient of turbulent flows is the vortex stretching mechanism, which emanates from the non-linear interaction of vorticity and strain-rate tensor and leads to formation of extreme events. We analyze the statistical correlations between vorticity and strain rate by using a massive database generated from very well resolved direct numerical simulations of forced isotropic turbulence in periodic domains. The grid resolution is up to $12288^3$, and the Taylor-scale Reynolds number is in the range $140-1300$. In order to understand the formation and structure of extreme vorticity fluctuations, we obtain statistics conditioned on enstrophy (vorticity-squared). The magnitude of strain, as well as its eigenvalues, is approximately constant when conditioned on weak enstrophy; whereas they grow approximately as power laws for strong enstrophy, which become steeper with increasing $R_\lambda$. We find that the well-known preferential alignment between vorticity and the intermediate eigenvector of strain tensor is even stronger for large enstrophy, whereas vorticity shows a tendency to be weakly orthogonal to the most extensive eigenvector (for large enstrophy). Yet the dominant contribution to the production of large enstrophy events arises from the most extensive eigendirection, the more so as $R_\lambda$ increases. Nevertheless, the stretching in intense vorticity regions is significantly depleted, consistent with the kinematic properties of weakly-curved tubes in which they are organized. Further analysis reveals that intense enstrophy is primarily depleted via viscous diffusion, though viscous dissipation is also significant. Implications for modeling are nominally addressed as appropriate., Comment: 17 pages, 8 figures

Cells have the ability to detect electric fields and respond to them with directed migratory movement. Investigations identified genes and proteins that play important roles in defining the migration efficiency. Nevertheless, the sensing and transduction mechanisms underlying directed cell migration are still under discussion. We use Dictyostelium discoideum cells as model system for studying eukaryotic cell migration in DC electric fields. We have defined the temporal electric persistence to characterize the memory that cells have in a varying electric field. In addition to imposing a directional bias, we observed that the electric field influences the cellular kinematics by accelerating the movement of cells along their paths. Moreover, the study of vegetative and briefly starved cells provided insight into the electrical sensing of cells. We found evidence that conditioned medium of starved cells was able to trigger the electrical sensing of vegetative cells that would otherwise not orient themselves in the electric field. This observation may be explained by the presence of the conditioned medium factor (CMF), a protein secreted by the cells, when they begin to starve. The results of this study give new insights into understanding the mechanism that triggers the electrical sensing and transduces the external stimulus into directed cell migration. Finally, the observed increased mobility of cells over time in an electric field could offer a novel perspective towards wound healing assays.

For rapidly rotating turbulent Rayleigh--B\'enard 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.

In nature interactions between biopolymers and motor proteins give rise to biologically essential emergent behaviours. Besides cytoskeleton mechanics, active nematics arise from such interactions. Here we present a study on 3D active nematics made of microtubules, kinesin motors and depleting agent. It shows a rich behaviour evolving from a nematically ordered space-filling distribution of microtubule bundles toward a flattened and contracted 2D ribbon that undergoes a wrinkling instability and subsequently transitions into a 3D active turbulent state. The wrinkle wavelength is independent of the ATP concentration and our theoretical model describes its relation with the appearance time. We compare the experimental results with a numerical simulation that confirms the key role of kinesin motors in cross-linking and sliding the microtubules. Our results on the active contraction of the network and the independence of wrinkle wavelength on ATP concentration are important steps forward for the understanding of these 3D systems., Comment: 12 pages, 5 figures

The formation of ice in mixed phase clouds greatly impacts Earths hydrologic cycle. The intensity, distribution and frequency of precipitation as well as radiative properties of clouds in the mid latitudes are strongly influenced by the number concentration of ice particles. A long standing riddle in mixed phase clouds is the frequent observation of measured ice particle concentrations several orders of magnitude higher than measured ice nucleating particle concentrations. Here, we report laboratory observations of copious cloud droplets and ice crystals formed in the wake of a warm, falling water drop. Aerosols were activated in the transient regions of very high supersaturation due to evaporative mixing in the wake. We extend these results to typical mixed phase atmospheric conditions, and our calculations show that the induced evaporative supersaturation may significantly enhance the activated ice nuclei concentration in the particles wake., Comment: 15 pages, 7 figures

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.

Using volumetric velocity data from a turbulent laboratory water flow and numerical simulations of homogeneous, isotropic turbulence, we present a direct experimental and numerical assessment of Kolmogorov's first refined similarity hypothesis based on three-dimensional measurements of the local energy dissipation rate $\epsilon_r$ measured at dissipative scales $r$. We focus on the properties of the stochastic variables $V_L = \Delta u(r)/(r \epsilon_r)^{1/3}$ and $V_T = \Delta v(r)/(r\epsilon_r)^{1/3}$, where $\Delta u(r)$ and $\Delta v(r)$ are longitudinal and transverse velocity increments. Over one order of magnitude of scales $r$ within the dissipative range, the distributions of $V_L$ and $V_T$ from both experiment and simulation collapse when parameterised by a suitably defined local Reynolds number, providing the first conclusive experimental evidence in support of the first refined similarity hypothesis and its universality., Comment: As published in PRFluids, 5 pages, 2 figures

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_\lambda$). Remarkably, we find the strongest velocity gradients to empirically scale as $\tau_K^{-1} R_\lambda^{\beta}$, with $\beta \approx 0.775 \pm 0.025$, where $\tau_K$ is the Kolmogorov time scale (with its inverse, $\tau_K^{-1}$, being the {r.m.s.} of velocity gradient fluctuations). Additionally, we observe velocity increments across very small distances $r \le \eta$, where $\eta$ 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 $\eta R_\lambda^{-\alpha}$, with $\alpha = \beta - \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 $\gamma < 1$, which weakly increases with $R_\lambda$. Using scaling arguments, we get $\beta=(2-\gamma)^{-1}$, which suggests that $\beta$ would also slowly increase with $R_\lambda$. We conjecture that approaching the limit of infinite $R_\lambda$, the flow is overall smooth, with intense velocity gradients over scale $ \eta R_\lambda^{-1/2}$, corresponding to $\beta = 1$., Comment: 14 pages, 8 figures

The energy spectrum of incompressible turbulence is known to reveal a pileup of energy at those high wavenumbers where viscous dissipation begins to act. It is called the bottleneck effect. Based on direct numerical simulations of the incompressible Navier-Stokes equations, results from Donzis & Sreenivasan (2010) pointed to a decrease of the strength of the bottleneck with increasing intensity of the turbulence, measured by the Taylor micro-scale Reynolds number $R_{\lambda}$. Here we report first experimental results on the dependence of the amplitude of the bottleneck as a function of $R_{\lambda}$ in a wind-tunnel flow. We used an active grid in the Variable Density Turbulence Tunnel (VDTT) (see Bodenschatz et al. (2014)) to reach $R_{\lambda}$ > 5000, which is unmatched in laboratory flows of decaying turbulence. The VDTT with the active grid permitted us to measure energy spectra from flows of different $R_{\lambda}$, with the small-scale features appearing always at the same frequencies. We relate those spectra recorded to a common reference spectrum, largely eliminating systematic errors which plague hotwire measurements at high frequencies. The data are consistent with a power law for the decrease of the bottleneck strength for the finite range of $R_{\lambda}$ in the experiment., Comment: 13 pages; 11 figures v2: clarified tick labels in Figs 6 and 7

Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and facilitate growth past barriers in the metre size range. Amongst a number of mechanisms that can amplify the local density of solids, aerodynamic streaming instability (SI) is one of the most promising. This work tests the physical assumptions of models that lead to SI in protoplanetary disks (PPDs). We conduct laboratory experiments in which we track the three-dimensional motion of spherical solid particles fluidized in a low-pressure, laminar, incompressible, gas stream. The particle sizes span the Stokes-Epstein drag regime transition and the overall dust-to-gas mass density ratio is close to unity. Lambrechts et al. (2016) established the similarity of the laboratory flow to a simplified PPD model flow. We obtain experimental results suggesting an instability due to particle-gas interaction: i) there exist variations in particle concentration in the direction of the mean drag forces; ii) the particles have a tendency to 'catch up' to one another when they are in proximity; iii) particle clumping occurs on very small scales, which implies local enhancements above the background dust-to-gas mass density ratio by factors of several tens; v) the presence of these density enhancements occurs for a mean dust-to-gas mass density ratio approaching or greater than 1; v) we find evidence for collective particle drag reduction when the local particle number density becomes high and when the background gas pressure is high so that the drag is in the continuum regime. The experiments presented here are precedent-setting for observing SI under controlled conditions and may lead to a deeper understanding of how it operates in nature., Comment: Accepted for publication in A&A; abstract abridged for arXiv

The character of turbulence depends on where it develops. Turbulence near boundaries, for instance, is different than in a free stream. To elucidate the differences between flows, it is instructive to vary the structure of turbulence systematically, but there are few ways of stirring turbulence that make this possible. In other words, an experiment typically examines either a boundary layer or a free stream, say, and the structure of the turbulence is fixed by the geometry of the experiment. We introduce a new active grid with many more degrees of freedom than previous active grids. The additional degrees of freedom make it possible to control various properties of the turbulence. We show how long-range correlations in the turbulent velocity fluctuations can be shaped by changing the way the active grid moves. Specifically, we show how not only the correlation length but also the detailed shape of the correlation function depends on the correlations imposed in the motions of the grid. Until now, large-scale structure had not been adjustable in experiments. This new capability makes possible new systematic investigations into turbulence dissipation and dispersion, for example, and perhaps in flows that mimic features of boundary layers, free streams, and flows of intermediate character., Comment: This paper has been accepted to Experiments in Fluids. 25 pages, 10 figures

We investigate sources of error in acceleration statistics from Lagrangian Particle Tracking (LPT) data and demonstrate techniques to eliminate or minimise bias errors introduced during processing. Numerical simulations of particle tracking experiments in isotropic turbulence show that the main sources of bias error arise from noise due to position uncertainty and selection biases introduced during numerical differentiation. We outline the use of independent measurements and filtering schemes to eliminate these biases. Moreover, we test the validity of our approach in estimating the statistical moments and probability densities of the Lagrangian acceleration. Finally, we apply these techniques to experimental particle tracking data and demonstrate their validity in practice with comparisons to available data from literature. The general approach, which is not limited to acceleration statistics, can be applied with as few as two cameras and permits a substantial reduction in the spatial resolution and sampling rate required to adequately measure statistics of Lagrangian acceleration.

We report experimental and theoretical results on spatiotemporal pattern formation in cell populations, where the parameters vary in space and time due to mechanisms intrinsic to the system, namely Dictyostelium discoideum (D.d.) in the starvation phase. We find that different patterns are formed when the populations are initialized at different developmental stages, or, when populations at different initial developmental stages are mixed. The experimentally observed patterns can be understood with a modified Kessler-Levine model that takes into account the initial spatial heterogeneity of the cell populations and a developmental path introduced by us, i.e., the time dependence of the various biochemical parameters. The dynamics of the parameter agree with known biochemical studies. Most importantly the modified model reproduces not only our results, but also the observations of an independent experiment published earlier. This shows that pattern formation can be used to understand and quantify the temporal evolution of the system parameters., Comment: 28 pages, 12 figures

In a reaction-diffusion-advection system, with a convectively unstable regime, a perturbation creates a wave train that is advected downstream and eventually leaves the system. We show that the convective instability coexists with a local absolute instability when a fixed boundary condition upstream is imposed. This boundary induced instability acts as a continuous wave source, creating a local periodic excitation near the boundary, which initiates waves traveling both up and downstream. To confirm this, we performed analytical analysis and numerical simulations of a modified Martiel-Goldbeter reaction-diffusion model with the addition of an advection term. We provide a quantitative description of the wave packet appearing in the convectively unstable regime, which we found to be in excellent agreement with the numerical simulations. We characterize this new instability and show that in the limit of high advection speed, it is suppressed. This type of instability can be expected for reaction-diffusion systems that present both a convective instability and an excitable regime. In particular, it can be relevant to understand the signaling mechanism of the social amoeba Dictyostelium discoideum that may experience fluid flows in its natural habitat., Comment: 10 pages, 13 figures, published in Chaos: An Interdisciplinary Journal of Nonlinear Science

We present results from moist convection in a mixture of pressurized sulfur hexa-flouride (liquid and vapor) and helium (gas) to model the wet and dry components of the earth's atmosphere. To allow for homogeneous nucleation, we operate the experiment close to critical conditions. We report on the nucleation of microdroplets in the wake of large cold liquid drops falling through the supersaturated atmosphere and show that the homogeneous nucleation is caused by isobaric cooling of the saturated sulfur hexaflouride vapor. Our results carry over to atmospheric clouds: falling hail and cold rain drops may enhance the heterogeneous nucleation of microdroplets in their wake under supersaturated atmospheric conditions. We also observed that under appropriate conditions settling microdroplets form a rather stable horizontal cloud layer, which separates regions of super and sub critical saturation., Comment: 6 pages, 3 figures

Using the unique capabilities of the Variable Density Turbulence Tunnel at the Max Planck Institute for Dynamics and Self-Organization, G\"{o}ttingen, we report experimental result on classical grid turbulence that uncover fine, yet important details of the structure functions in the inertial range. This was made possible by measuring extremely long time series of up to $10^{10}$ samples of the turbulent fluctuating velocity, which corresponds to $\mathcal{O}\left(10^5\right)$ large eddy turnover times. These classical grid measurements were conducted in a well-controlled environment at a wide range of high Reynolds numbers from $R_\lambda=110$ up to $R_\lambda=1600$, using both traditional hot-wire probes as well as NSTAP probes developed at Princeton University. We found that deviations from ideal scaling are anchored to the small scales and that dissipation influences the inertial-range statistics at scales larger than the near-dissipation range., Comment: 6 pages, 5 figures

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., Comment: Accepted for publication in Astronomy and Astrophysics

We present measurements of the orientation $\theta_0$ and temperature amplitude $\delta$ of the large-scale circulation in a cylindrical sample of turbulent Rayleigh-Benard convection (RBC) with aspect ratio $\Gamma \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 $\theta_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_\theta$ and a corresponding Reynolds number $Re_{\theta}$ 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_{\theta} \propto Ra^{0.28}$, and with the Prandtl-number dependence $Re_{\theta} \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_{\theta}(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_{\theta} \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\"ottingen experiments and discussed in the present paper is indeed to a new state of RBC referred to as "ultimate"., Comment: 12 pages, 4 figures, to be pub. in JFMR

We present results from a numerical study of particle dispersion in the weakly nonlinear regime of Rayleigh-B\'enard 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., Comment: 15 pages, 6 figures

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 $\gamma$ 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 $\gamma-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., Comment: 26 pages, 20 figures, 1 table. Changes in this version: extensive restructuring according to reviewer comments and modification to order of authors

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., Comment: 13 pages, 8 figures

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., Comment: accepted for publication in Science China - Physics, Mechanics & Astronomy