Speed and structure of turbulent fronts in pipe flow

Using extensive direct numerical simulations, the dynamics of laminar-turbulent fronts in pipe flow is investigated for Reynolds numbers between $Re=2000$ and $5500$. We here investigate the physical distinction between the fronts of weak and strong slugs both by analysing the turbulent kinetic energy budget and by comparing the downstream front motion to the advection speed of bulk turbulent structures. Our study shows that weak downstream fronts travel slower than turbulent structures in the bulk and correspond to decaying turbulence at the front. At $Re\simeq 2900$ the downstream front speed becomes faster than the advection speed, marking the onset of strong fronts. In contrast to weak fronts, turbulent eddies are generated at strong fronts by feeding on the downstream laminar flow. Our study also suggests that temporal fluctuations of production and dissipation at the downstream laminar-turbulent front drive the dynamical switches between the two types of front observed up to $Re\simeq 3200$.Comment: 14 pages, accepted for publication in Journal of Fluid Mechanic

Scaling and mechanism of the propagation speed of the upstream turbulent front in pipe flow

Scaling and mechanism of the propagation speed of turbulent fronts in pipe flow with the Reynolds number has been a long-standing problem in the past decades. Here, we derive an explicit scaling law of the upstream front speed, which approaches to a power-law scaling at high Reynolds numbers and explain the underlying mechanism. Our data show that the average wall distance of low-speed streaks at the tip of the upstream front, where transition occurs, appears to be constant in local wall units in the wide bulk-Reynolds-number range investigated, between 5000 and 60000. By further assuming that the axial propagation of velocity fluctuations at the front tip, resulting from streak instabilities, is dominated by the advection of the local mean flow, the front speed can be derived as an explicit function of the Reynolds number. The derived formula agrees well with the measured speed by front tracking. Our finding reveals the relationship between the structure and speed of a front, which enables to obtain a close approximation of the front speed based on a single velocity field without having to track the front over time

Phase-field simulation of core-annular pipe flow

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.Comment: Paper accepted for publication in International Journal of Multiphase Flo

Transition to turbulence in pulsating pipe flow

Fluid flows in nature and applications are frequently subject to periodic velocity modulations. Surprisingly, even for the generic case of flow through a straight pipe, there is little consensus regarding the influence of pulsation on the transition threshold to turbulence: while most studies predict a monotonically increasing threshold with pulsation frequency (i.e. Womersley number, $\alpha$), others observe a decreasing threshold for identical parameters and only observe an increasing threshold at low $\alpha$. In the present study we apply recent advances in the understanding of transition in steady shear flows to pulsating pipe flow. For moderate pulsation amplitudes we find that the first instability encountered is subcritical (i.e. requiring finite amplitude disturbances) and gives rise to localized patches of turbulence ("puffs") analogous to steady pipe flow. By monitoring the impact of pulsation on the lifetime of turbulence we map the onset of turbulence in parameter space. Transition in pulsatile flow can be separated into three regimes. At small Womersley numbers the dynamics are dominated by the decay turbulence suffers during the slower part of the cycle and hence transition is delayed significantly. As shown in this regime thresholds closely agree with estimates based on a quasi steady flow assumption only taking puff decay rates into account. The transition point predicted in the zero $\alpha$ limit equals to the critical point for steady pipe flow offset by the oscillation Reynolds number. In the high frequency limit puff lifetimes are identical to those in steady pipe flow and hence the transition threshold appears to be unaffected by flow pulsation. In the intermediate frequency regime the transition threshold sharply drops (with increasing $\alpha$) from the decay dominated (quasi steady) threshold to the steady pipe flow level

The rise of fully turbulent flow

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