@article{5943, abstract = {The hairpin instability of a jet in a crossflow (JICF) for a low jet-to-crossflow velocity ratio is investigated experimentally for a velocity ratio range of R ∈ (0.14, 0.75) and crossflow Reynolds numbers ReD ∈ (260, 640). From spectral analysis we characterize the Strouhal number and amplitude of the hairpin instability as a function of R and ReD. We demonstrate that the dynamics of the hairpins is well described by the Landau model, and, hence, that the instability occurs through Hopf bifurcation, similarly to other hydrodynamical oscillators such as wake behind different bluff bodies. Using the Landau model, we determine the precise threshold values of hairpin shedding. We also study the spatial dependence of this hydrodynamical instability, which shows a global behaviour.}, author = {Klotz, Lukasz and Gumowski, Konrad and Wesfreid, José Eduardo}, journal = {Journal of Fluid Mechanics}, pages = {386--406}, publisher = {Cambridge University Press}, title = {{Experiments on a jet in a crossflow in the low-velocity-ratio regime}}, doi = {10.1017/jfm.2018.974}, volume = {863}, year = {2019}, } @article{5878, abstract = {We consider the motion of a droplet bouncing on a vibrating bath of the same fluid in the presence of a central potential. We formulate a rotation symmetry-reduced description of this system, which allows for the straightforward application of dynamical systems theory tools. As an illustration of the utility of the symmetry reduction, we apply it to a model of the pilot-wave system with a central harmonic force. We begin our analysis by identifying local bifurcations and the onset of chaos. We then describe the emergence of chaotic regions and their merging bifurcations, which lead to the formation of a global attractor. In this final regime, the droplet’s angular momentum spontaneously changes its sign as observed in the experiments of Perrard et al.}, author = {Budanur, Nazmi B and Fleury, Marc}, issn = {1089-7682}, journal = {Chaos: An Interdisciplinary Journal of Nonlinear Science}, number = {1}, publisher = {AIP Publishing}, title = {{State space geometry of the chaotic pilot-wave hydrodynamics}}, doi = {10.1063/1.5058279}, volume = {29}, year = {2019}, } @article{6413, abstract = {Phase-field methods have long been used to model the flow of immiscible fluids. Their ability to naturally capture interface topological changes is widely recognized, but their accuracy in simulating flows of real fluids in practical geometries is not established. We here quantitatively investigate the convergence of the phase-field method to the sharp-interface limit with simulations of two-phase pipe flow. We focus on core-annular flows, in which a highly viscous fluid is lubricated by a less viscous fluid, and validate our simulations with an analytic laminar solution, a formal linear stability analysis and also in the fully nonlinear regime. We demonstrate the ability of the phase-field method to accurately deal with non-rectangular geometry, strong advection, unsteady fluctuations and large viscosity contrast. We argue that phase-field methods are very promising for quantitatively studying moderately turbulent flows, especially at high concentrations of the disperse phase.}, author = {Song, Baofang and Plana, Carlos and Lopez Alonso, Jose M and Avila, Marc}, issn = {03019322}, journal = {International Journal of Multiphase Flow}, pages = {14--24}, publisher = {Elsevier}, title = {{Phase-field simulation of core-annular pipe flow}}, doi = {10.1016/j.ijmultiphaseflow.2019.04.027}, volume = {117}, year = {2019}, } @article{6978, abstract = {In pipes and channels, the onset of turbulence is initially dominated by localizedtransients, which lead to sustained turbulence through their collective dynamics. In thepresent work, we study numerically the localized turbulence in pipe flow and elucidate astate space structure that gives rise to transient chaos. Starting from the basin boundaryseparating laminar and turbulent flow, we identify transverse homoclinic orbits, thepresence of which necessitates a homoclinic tangle and chaos. A direct consequence ofthe homoclinic tangle is the fractal nature of the laminar-turbulent boundary, which wasconjectured in various earlier studies. By mapping the transverse intersections between thestable and unstable manifold of a periodic orbit, we identify the gateways that promote anescape from turbulence.}, author = {Budanur, Nazmi B and Dogra, Akshunna and Hof, Björn}, journal = {Physical Review Fluids}, number = {10}, pages = {102401}, publisher = {American Physical Society}, title = {{Geometry of transient chaos in streamwise-localized pipe flow turbulence}}, doi = {10.1103/PhysRevFluids.4.102401}, volume = {4}, year = {2019}, } @article{7397, abstract = {Polymer additives can substantially reduce the drag of turbulent flows and the upperlimit, the so called “maximum drag reduction” (MDR) asymptote is universal, i.e. inde-pendent of the type of polymer and solvent used. Until recently, the consensus was that,in this limit, flows are in a marginal state where only a minimal level of turbulence activ-ity persists. Observations in direct numerical simulations using minimal sized channelsappeared to support this view and reported long “hibernation” periods where turbu-lence is marginalized. In simulations of pipe flow we find that, indeed, with increasingWeissenberg number (Wi), turbulence expresses long periods of hibernation if the domainsize is small. However, with increasing pipe length, the temporal hibernation continuouslyalters to spatio-temporal intermittency and here the flow consists of turbulent puffs sur-rounded by laminar flow. Moreover, upon an increase in Wi, the flow fully relaminarises,in agreement with recent experiments. At even larger Wi, a different instability is en-countered causing a drag increase towards MDR. Our findings hence link earlier minimalflow unit simulations with recent experiments and confirm that the addition of polymersinitially suppresses Newtonian turbulence and leads to a reverse transition. The MDRstate on the other hand results from a separate instability and the underlying dynamicscorresponds to the recently proposed state of elasto-inertial-turbulence (EIT).}, author = {Lopez Alonso, Jose M and Choueiri, George H and Hof, Björn}, issn = {1469-7645}, journal = {Journal of Fluid Mechanics}, pages = {699--719}, publisher = {CUP}, title = {{Dynamics of viscoelastic pipe flow at low Reynolds numbers in the maximum drag reduction limit}}, doi = {10.1017/jfm.2019.486}, volume = {874}, year = {2019}, } @phdthesis{6957, abstract = {In many shear flows like pipe flow, plane Couette flow, plane Poiseuille flow, etc. turbulence emerges subcritically. Here, when subjected to strong enough perturbations, the flow becomes turbulent in spite of the laminar base flow being linearly stable. The nature of this instability has puzzled the scientific community for decades. At onset, turbulence appears in localized patches and flows are spatio-temporally intermittent. In pipe flow the localized turbulent structures are referred to as puffs and in planar flows like plane Couette and channel flow, patches arise in the form of localized oblique bands. In this thesis, we study the onset of turbulence in channel flow in direct numerical simulations from a dynamical system theory perspective, as well as by performing experiments in a large aspect ratio channel. The aim of the experimental work is to determine the critical Reynolds number where turbulence first becomes sustained. Recently, the onset of turbulence has been described in analogy to absorbing state phase transition (i.e. directed percolation). In particular, it has been shown that the critical point can be estimated from the competition between spreading and decay processes. Here, by performing experiments, we identify the mechanisms underlying turbulence proliferation in channel flow and find the critical Reynolds number, above which turbulence becomes sustained. Above the critical point, the continuous growth at the tip of the stripes outweighs the stochastic shedding of turbulent patches at the tail and the stripes expand. For growing stripes, the probability to decay decreases while the probability of stripe splitting increases. Consequently, and unlike for the puffs in pipe flow, neither of these two processes is time-independent i.e. memoryless. Coupling between stripe expansion and creation of new stripes via splitting leads to a significantly lower critical point ($Re_c=670+/-10$) than most earlier studies suggest. While the above approach sheds light on how turbulence first becomes sustained, it provides no insight into the origin of the stripes themselves. In the numerical part of the thesis we investigate how turbulent stripes form from invariant solutions of the Navier-Stokes equations. The origin of these turbulent stripes can be identified by applying concepts from the dynamical system theory. In doing so, we identify the exact coherent structures underlying stripes and their bifurcations and how they give rise to the turbulent attractor in phase space. We first report a family of localized nonlinear traveling wave solutions of the Navier-Stokes equations in channel flow. These solutions show structural similarities with turbulent stripes in experiments like obliqueness, quasi-streamwise streaks and vortices, etc. A parametric study of these traveling wave solution is performed, with parameters like Reynolds number, stripe tilt angle and domain size, including the stability of the solutions. These solutions emerge through saddle-node bifurcations and form a phase space skeleton for the turbulent stripes observed in the experiments. The lower branches of these TW solutions at different tilt angles undergo Hopf bifurcation and new solutions branches of relative periodic orbits emerge. These RPO solutions do not belong to the same family and therefore the routes to chaos for different angles are different. In shear flows, turbulence at onset is transient in nature. Consequently,turbulence can not be tracked to lower Reynolds numbers, where the dynamics may simplify. Before this happens, turbulence becomes short-lived and laminarizes. In the last part of the thesis, we show that using numerical simulations we can continue turbulent stripes in channel flow past the 'relaminarization barrier' all the way to their origin. Here, turbulent stripe dynamics simplifies and the fluctuations are no longer stochastic and the stripe settles down to a relative periodic orbit. This relative periodic orbit originates from the aforementioned traveling wave solutions. Starting from the relative periodic orbit, a small increase in speed i.e. Reynolds number gives rise to chaos and the attractor dimension sharply increases in contrast to the classical transition scenario where the instabilities affect the flow globally and give rise to much more gradual route to turbulence.}, author = {Paranjape, Chaitanya S}, issn = {2663-337X}, keywords = {Instabilities, Turbulence, Nonlinear dynamics}, pages = {138}, publisher = {Institute of Science and Technology Austria}, title = {{Onset of turbulence in plane Poiseuille flow}}, doi = {10.15479/AT:ISTA:6957}, year = {2019}, } @article{7197, abstract = {During bacterial cell division, the tubulin-homolog FtsZ forms a ring-like structure at the center of the cell. This Z-ring not only organizes the division machinery, but treadmilling of FtsZ filaments was also found to play a key role in distributing proteins at the division site. What regulates the architecture, dynamics and stability of the Z-ring is currently unknown, but FtsZ-associated proteins are known to play an important role. Here, using an in vitro reconstitution approach, we studied how the well-conserved protein ZapA affects FtsZ treadmilling and filament organization into large-scale patterns. Using high-resolution fluorescence microscopy and quantitative image analysis, we found that ZapA cooperatively increases the spatial order of the filament network, but binds only transiently to FtsZ filaments and has no effect on filament length and treadmilling velocity. Together, our data provides a model for how FtsZ-associated proteins can increase the precision and stability of the bacterial cell division machinery in a switch-like manner.}, author = {Dos Santos Caldas, Paulo R and Lopez Pelegrin, Maria D and Pearce, Daniel J. G. and Budanur, Nazmi B and Brugués, Jan and Loose, Martin}, issn = {2041-1723}, journal = {Nature Communications}, publisher = {Springer Nature}, title = {{Cooperative ordering of treadmilling filaments in cytoskeletal networks of FtsZ and its crosslinker ZapA}}, doi = {10.1038/s41467-019-13702-4}, volume = {10}, year = {2019}, } @article{6069, abstract = {Electron transport in two-dimensional conducting materials such as graphene, with dominant electron–electron interaction, exhibits unusual vortex flow that leads to a nonlocal current-field relation (negative resistance), distinct from the classical Ohm’s law. The transport behavior of these materials is best described by low Reynolds number hydrodynamics, where the constitutive pressure–speed relation is Stoke’s law. Here we report evidence of such vortices observed in a viscous flow of Newtonian fluid in a microfluidic device consisting of a rectangular cavity—analogous to the electronic system. We extend our experimental observations to elliptic cavities of different eccentricities, and validate them by numerically solving bi-harmonic equation obtained for the viscous flow with no-slip boundary conditions. We verify the existence of a predicted threshold at which vortices appear. Strikingly, we find that a two-dimensional theoretical model captures the essential features of three-dimensional Stokes flow in experiments.}, author = {Mayzel, Jonathan and Steinberg, Victor and Varshney, Atul}, issn = {2041-1723}, journal = {Nature Communications}, publisher = {Springer Nature}, title = {{Stokes flow analogous to viscous electron current in graphene}}, doi = {10.1038/s41467-019-08916-5}, volume = {10}, year = {2019}, } @article{6014, abstract = {Speed of sound waves in gases and liquids are governed by the compressibility of the medium. There exists another type of non-dispersive wave where the wave speed depends on stress instead of elasticity of the medium. A well-known example is the Alfven wave, which propagates through plasma permeated by a magnetic field with the speed determined by magnetic tension. An elastic analogue of Alfven waves has been predicted in a flow of dilute polymer solution where the elastic stress of the stretching polymers determines the elastic wave speed. Here we present quantitative evidence of elastic Alfven waves in elastic turbulence of a viscoelastic creeping flow between two obstacles in channel flow. The key finding in the experimental proof is a nonlinear dependence of the elastic wave speed cel on the Weissenberg number Wi, which deviates from predictions based on a model of linear polymer elasticity.}, author = {Varshney, Atul and Steinberg, Victor}, issn = {2041-1723}, journal = {Nature Communications}, publisher = {Springer Nature}, title = {{Elastic alfven waves in elastic turbulence}}, doi = {10.1038/s41467-019-08551-0}, volume = {10}, year = {2019}, } @article{6779, abstract = {Recent studies suggest that unstable recurrent solutions of the Navier-Stokes equation provide new insights into dynamics of turbulent flows. In this study, we compute an extensive network of dynamical connections between such solutions in a weakly turbulent quasi-two-dimensional Kolmogorov flow that lies in the inversion symmetric subspace. In particular, we find numerous isolated heteroclinic connections between different types of solutions—equilibria, periodic, and quasiperiodic orbits—as well as continua of connections forming higher-dimensional connecting manifolds. We also compute a homoclinic connection of a periodic orbit and provide strong evidence that the associated homoclinic tangle forms the chaotic repeller that underpins transient turbulence in the symmetric subspace.}, author = {Suri, Balachandra and Pallantla, Ravi Kumar and Schatz, Michael F. and Grigoriev, Roman O.}, issn = {2470-0053}, journal = {Physical Review E}, number = {1}, publisher = {American Physical Society}, title = {{Heteroclinic and homoclinic connections in a Kolmogorov-like flow}}, doi = {10.1103/physreve.100.013112}, volume = {100}, year = {2019}, } @article{6486, abstract = {Based on a novel control scheme, where a steady modification of the streamwise velocity profile leads to complete relaminarization of initially fully turbulent pipe flow, we investigate the applicability and usefulness of custom-shaped honeycombs for such control. The custom-shaped honeycombs are used as stationary flow management devices which generate specific modifications of the streamwise velocity profile. Stereoscopic particle image velocimetry and pressure drop measurements are used to investigate and capture the development of the relaminarizing flow downstream these devices. We compare the performance of straight (constant length across the radius of the pipe) honeycombs with custom-shaped ones (variable length across the radius) and try to determine the optimal shape for maximal relaminarization at minimal pressure loss. The optimally modified streamwise velocity profile is found to be M-shaped, and the maximum attainable Reynolds number for total relaminarization is found to be of the order of 10,000. Consequently, the respective reduction in skin friction downstream of the device is almost by a factor of 5. The break-even point, where the additional pressure drop caused by the device is balanced by the savings due to relaminarization and a net gain is obtained, corresponds to a downstream stretch of distances as low as approximately 100 pipe diameters of laminar flow.}, author = {Kühnen, Jakob and Scarselli, Davide and Hof, Björn}, issn = {1528901X}, journal = {Journal of Fluids Engineering}, number = {11}, publisher = {ASME}, title = {{Relaminarization of pipe flow by means of 3D-printed shaped honeycombs}}, doi = {10.1115/1.4043494}, volume = {141}, year = {2019}, } @article{6228, abstract = {Following the recent observation that turbulent pipe flow can be relaminarised bya relatively simple modification of the mean velocity profile, we here carry out aquantitative experimental investigation of this phenomenon. Our study confirms thata flat velocity profile leads to a collapse of turbulence and in order to achieve theblunted profile shape, we employ a moving pipe segment that is briefly and rapidlyshifted in the streamwise direction. The relaminarisation threshold and the minimumshift length and speeds are determined as a function of Reynolds number. Althoughturbulence is still active after the acceleration phase, the modulated profile possessesa severely decreased lift-up potential as measured by transient growth. As shown,this results in an exponential decay of fluctuations and the flow relaminarises. Whilethis method can be easily applied at low to moderate flow speeds, the minimumstreamwise length over which the acceleration needs to act increases linearly with theReynolds number.}, author = {Scarselli, Davide and Kühnen, Jakob and Hof, Björn}, issn = {14697645}, journal = {Journal of Fluid Mechanics}, pages = {934--948}, publisher = {Cambridge University Press}, title = {{Relaminarising pipe flow by wall movement}}, doi = {10.1017/jfm.2019.191}, volume = {867}, year = {2019}, } @article{6508, abstract = {Segregation of maternal determinants within the oocyte constitutes the first step in embryo patterning. In zebrafish oocytes, extensive ooplasmic streaming leads to the segregation of ooplasm from yolk granules along the animal-vegetal axis of the oocyte. Here, we show that this process does not rely on cortical actin reorganization, as previously thought, but instead on a cell-cycle-dependent bulk actin polymerization wave traveling from the animal to the vegetal pole of the oocyte. This wave functions in segregation by both pulling ooplasm animally and pushing yolk granules vegetally. Using biophysical experimentation and theory, we show that ooplasm pulling is mediated by bulk actin network flows exerting friction forces on the ooplasm, while yolk granule pushing is achieved by a mechanism closely resembling actin comet formation on yolk granules. Our study defines a novel role of cell-cycle-controlled bulk actin polymerization waves in oocyte polarization via ooplasmic segregation.}, author = {Shamipour, Shayan and Kardos, Roland and Xue, Shi-lei and Hof, Björn and Hannezo, Edouard B and Heisenberg, Carl-Philipp J}, issn = {10974172}, journal = {Cell}, number = {6}, pages = {1463--1479.e18}, publisher = {Elsevier}, title = {{Bulk actin dynamics drive phase segregation in zebrafish oocytes}}, doi = {10.1016/j.cell.2019.04.030}, volume = {177}, year = {2019}, } @article{7001, author = {Schwayer, Cornelia and Shamipour, Shayan and Pranjic-Ferscha, Kornelija and Schauer, Alexandra and Balda, M and Tada, M and Matter, K and Heisenberg, Carl-Philipp J}, issn = {1097-4172}, journal = {Cell}, number = {4}, pages = {937--952.e18}, publisher = {Cell Press}, title = {{Mechanosensation of tight junctions depends on ZO-1 phase separation and flow}}, doi = {10.1016/j.cell.2019.10.006}, volume = {179}, year = {2019}, } @article{6189, abstract = {Suspended particles can alter the properties of fluids and in particular also affect the transition fromlaminar to turbulent flow. An earlier study [Mataset al.,Phys. Rev. Lett.90, 014501 (2003)] reported howthe subcritical (i.e., hysteretic) transition to turbulent puffs is affected by the addition of particles. Here weshow that in addition to this known transition, with increasing concentration a supercritical (i.e.,continuous) transition to a globally fluctuating state is found. At the same time the Newtonian-typetransition to puffs is delayed to larger Reynolds numbers. At even higher concentration only the globallyfluctuating state is found. The dynamics of particle laden flows are hence determined by two competinginstabilities that give rise to three flow regimes: Newtonian-type turbulence at low, a particle inducedglobally fluctuating state at high, and a coexistence state at intermediate concentrations.}, author = {Agrawal, Nishchal and Choueiri, George H and Hof, Björn}, issn = {10797114}, journal = {Physical Review Letters}, number = {11}, publisher = {American Physical Society}, title = {{Transition to turbulence in particle laden flows}}, doi = {10.1103/PhysRevLett.122.114502}, volume = {122}, year = {2019}, } @article{291, abstract = {Over the past decade, the edge of chaos has proven to be a fruitful starting point for investigations of shear flows when the laminar base flow is linearly stable. Numerous computational studies of shear flows demonstrated the existence of states that separate laminar and turbulent regions of the state space. In addition, some studies determined invariant solutions that reside on this edge. In this paper, we study the unstable manifold of one such solution with the aid of continuous symmetry reduction, which we formulate here for the simultaneous quotiening of axial and azimuthal symmetries. Upon our investigation of the unstable manifold, we discover a previously unknown traveling-wave solution on the laminar-turbulent boundary with a relatively complex structure. By means of low-dimensional projections, we visualize different dynamical paths that connect these solutions to the turbulence. Our numerical experiments demonstrate that the laminar-turbulent boundary exhibits qualitatively different regions whose properties are influenced by the nearby invariant solutions.}, author = {Budanur, Nazmi B and Hof, Björn}, journal = {Physical Review Fluids}, number = {5}, publisher = {American Physical Society}, title = {{Complexity of the laminar-turbulent boundary in pipe flow}}, doi = {10.1103/PhysRevFluids.3.054401}, volume = {3}, year = {2018}, } @article{17, abstract = {Creeping flow of polymeric fluid without inertia exhibits elastic instabilities and elastic turbulence accompanied by drag enhancement due to elastic stress produced by flow-stretched polymers. However, in inertia-dominated flow at high Re and low fluid elasticity El, a reduction in turbulent frictional drag is caused by an intricate competition between inertial and elastic stresses. Here we explore the effect of inertia on the stability of viscoelastic flow in a broad range of control parameters El and (Re,Wi). We present the stability diagram of observed flow regimes in Wi-Re coordinates and find that the instabilities' onsets show an unexpectedly nonmonotonic dependence on El. Further, three distinct regions in the diagram are identified based on El. Strikingly, for high-elasticity fluids we discover a complete relaminarization of flow at Reynolds number in the range of 1 to 10, different from a well-known turbulent drag reduction. These counterintuitive effects may be explained by a finite polymer extensibility and a suppression of vorticity at high Wi. Our results call for further theoretical and numerical development to uncover the role of inertial effect on elastic turbulence in a viscoelastic flow.}, author = {Varshney, Atul and Steinberg, Victor}, journal = {Physical Review Fluids}, number = {10}, publisher = {American Physical Society}, title = {{Drag enhancement and drag reduction in viscoelastic flow}}, doi = {10.1103/PhysRevFluids.3.103302}, volume = {3}, year = {2018}, } @article{16, abstract = {We report quantitative evidence of mixing-layer elastic instability in a viscoelastic fluid flow between two widely spaced obstacles hindering a channel flow at Re 1 and Wi 1. Two mixing layers with nonuniform shear velocity profiles are formed in the region between the obstacles. The mixing-layer instability arises in the vicinity of an inflection point on the shear velocity profile with a steep variation in the elastic stress. The instability results in an intermittent appearance of small vortices in the mixing layers and an amplification of spatiotemporal averaged vorticity in the elastic turbulence regime. The latter is characterized through scaling of friction factor with Wi and both pressure and velocity spectra. Furthermore, the observations reported provide improved understanding of the stability of the mixing layer in a viscoelastic fluid at large elasticity, i.e., Wi 1 and Re 1 and oppose the current view of suppression of vorticity solely by polymer additives.}, author = {Varshney, Atul and Steinberg, Victor}, journal = {Physical Review Fluids}, number = {10}, publisher = {American Physical Society}, title = {{Mixing layer instability and vorticity amplification in a creeping viscoelastic flow}}, doi = {10.1103/PhysRevFluids.3.103303}, volume = {3}, year = {2018}, } @article{519, abstract = {This study treats with the influence of a symmetry-breaking transversal magnetic field on the nonlinear dynamics of ferrofluidic Taylor-Couette flow – flow confined between two concentric independently rotating cylinders. We detected alternating ‘flip’ solutions which are flow states featuring typical characteristics of slow-fast-dynamics in dynamical systems. The flip corresponds to a temporal change in the axial wavenumber and we find them to appear either as pure 2-fold axisymmetric (due to the symmetry-breaking nature of the applied transversal magnetic field) or involving non-axisymmetric, helical modes in its interim solution. The latter ones show features of typical ribbon solutions. In any case the flip solutions have a preferential first axial wavenumber which corresponds to the more stable state (slow dynamics) and second axial wavenumber, corresponding to the short appearing more unstable state (fast dynamics). However, in both cases the flip time grows exponential with increasing the magnetic field strength before the flip solutions, living on 2-tori invariant manifolds, cease to exist, with lifetime going to infinity. Further we show that ferrofluidic flow turbulence differ from the classical, ordinary (usually at high Reynolds number) turbulence. The applied magnetic field hinders the free motion of ferrofluid partials and therefore smoothen typical turbulent quantities and features so that speaking of mildly chaotic dynamics seems to be a more appropriate expression for the observed motion. }, author = {Altmeyer, Sebastian}, journal = {Journal of Magnetism and Magnetic Materials}, pages = {427 -- 441}, publisher = {Elsevier}, title = {{Non-linear dynamics and alternating ‘flip’ solutions in ferrofluidic Taylor-Couette flow}}, doi = {10.1016/j.jmmm.2017.12.073}, volume = {452}, year = {2018}, } @article{5996, abstract = {In pipes, turbulence sets in despite the linear stability of the laminar Hagen–Poiseuille flow. The Reynolds number ( ) for which turbulence first appears in a given experiment – the ‘natural transition point’ – depends on imperfections of the set-up, or, more precisely, on the magnitude of finite amplitude perturbations. At onset, turbulence typically only occupies a certain fraction of the flow, and this fraction equally is found to differ from experiment to experiment. Despite these findings, Reynolds proposed that after sufficiently long times, flows may settle to steady conditions: below a critical velocity, flows should (regardless of initial conditions) always return to laminar, while above this velocity, eddying motion should persist. As will be shown, even in pipes several thousand diameters long, the spatio-temporal intermittent flow patterns observed at the end of the pipe strongly depend on the initial conditions, and there is no indication that different flow patterns would eventually settle to a (statistical) steady state. Exploiting the fact that turbulent puffs do not age (i.e. they are memoryless), we continuously recreate the puff sequence exiting the pipe at the pipe entrance, and in doing so introduce periodic boundary conditions for the puff pattern. This procedure allows us to study the evolution of the flow patterns for arbitrary long times, and we find that after times in excess of advective time units, indeed a statistical steady state is reached. Although the resulting flows remain spatio-temporally intermittent, puff splitting and decay rates eventually reach a balance, so that the turbulent fraction fluctuates around a well-defined level which only depends on . In accordance with Reynolds’ proposition, we find that at lower (here 2020), flows eventually always resume to laminar, while for higher ( ), turbulence persists. The critical point for pipe flow hence falls in the interval of $2020 , which is in very good agreement with the recently proposed value of . The latter estimate was based on single-puff statistics and entirely neglected puff interactions. Unlike in typical contact processes where such interactions strongly affect the percolation threshold, in pipe flow, the critical point is only marginally influenced. Interactions, on the other hand, are responsible for the approach to the statistical steady state. As shown, they strongly affect the resulting flow patterns, where they cause ‘puff clustering’, and these regions of large puff densities are observed to travel across the puff pattern in a wave-like fashion.}, author = {Vasudevan, Mukund and Hof, Björn}, issn = {1469-7645}, journal = {Journal of Fluid Mechanics}, pages = {76--94}, publisher = {Cambridge University Press}, title = {{The critical point of the transition to turbulence in pipe flow}}, doi = {10.1017/jfm.2017.923}, volume = {839}, year = {2018}, }