Relaminarization by Steady Modification of the Streamwise Velocity Profile in a Pipe

We show that a rather simple, steady modification of the streamwise velocity profile in a pipe can lead to a complete collapse of turbulence and the flow fully relaminarizes. Two different devices, a stationary obstacle (inset) and a device which injects fluid through an annular gap close to the wall, are used to control the flow. Both devices modify the streamwise velocity profile such that the flow in the center of the pipe is decelerated and the flow in the near wall region is accelerated. We present measurements with stereoscopic particle image velocimetry to investigate and capture the development of the relaminarizing flow downstream these devices and the specific circumstances responsible for relaminarization. We find total relaminarization up to Reynolds numbers of 6000, where the skin friction in the far downstream distance is reduced by a factor of 3.4 due to relaminarization. In a smooth straight pipe the flow remains completely laminar downstream of the control. Furthermore, we show that transient (temporary) relaminarization in a spatially confined region right downstream the devices occurs also at much higher Reynolds numbers, accompanied by a significant local skin friction drag reduction. The underlying physical mechanism of relaminarization is attributed to a weakening of the near-wall turbulence production cycle.


Introduction
The control of turbulence to reduce skin friction is of great interest in a wide variety of technological applications, as frictional drag is a heavy consumer of energy and thus a possibly avoidable source of significant operating costs and massive carbon emissions. Numerous techniques for the control of turbulent flows have been proposed over the years. However, only a few methods of control have already been successfully implemented in practical engineering devices.
Flow control is commonly classified as either active or passive depending on whether or not external energy-input is required [1]. The most sophisticated and on a theoretical basis most elegant method is active feedback (closed loop) control of turbulence (see e.g. [2][3][4][5][6][7]). However, the practical implementation is technically highly demanding and expensive due to the required sensors for realtime flow measurements and elaborate actuators to control the flow. Predetermined (open loop) active techniques are usually characterized by greater simplicity and comparative ease of implementation. Among such techniques, modifying the flow through large-scale spanwise forcing created by boundary motion (wall oscillation or transverse traveling-wave excitation) or a body force has produced promising results (see e.g. [8][9][10][11][12][13][14]).
The wide range of passive control techniques comprises diverse approaches such as surface modifications by means of e.g. riblets [15], grooves [16], shark-skin surfaces [17], hydrophobic walls [18,19] or by forcing small optimal perturbations [20]. Furthermore, modifications of the fluid employing polymer additives [21,22] and modifications of the flow field by means of honeycombs and screens [23,24].
Most of the above mentioned methods have in common that they reduce skin friction and decrease the turbulence level by some amount in the control area but can not totally extinguish turbulence and have no or little effect further downstream. However, the ultimate goal of turbulence control in terms of energy saving is relaminarization (also denoted as turbulent-laminar transition, reversion, retransition, and laminarization) of the flow, leaving aside the discussion on possible sublaminar friction in duct flows [25,26].
Occasional evidence of relaminarization not determined by dissipation and the Reynolds number has been found when a turbulent flow is subjected to effects of acceleration, suction, blowing, magnetic fields, stratification, rotation, curvature and heating [27]. For example, a peculiarity of curved pipes is that the threshold for the onset of turbulence is postponed and occurs at Reynolds numbers considerably larger than in straight pipes. As a consequence, low Reynolds number turbulent flow emerging from a straight pipe may relaminarize in a subsequent curved pipe [28,29].
In particular, the effect of a favorable pressure gradient (FPG) on a boundary layer and the related case of accelerated pipe flow have received considerable attention concerning relaminarization. In accelerated pipe flow, i.e. during and subsequent to a rapid increase of the flow rate of an initially turbulent flow, the flow has been observed to transiently first visit a quasi-laminar state and then again undergo a process of transition that resembles the laminar-turbulent transition (see e.g. [30,31]; see also [32,33] for accelerated flow in a channel). And a strong FPG imposed on a boundary layer has been found to have a damping effect on the growth of perturbations [34] and to cause a temporary state of relaminarization (see e.g. [35][36][37][38][39][40][41][42]). In experiments, the FPG is often imposed on the flow by means of various types of convergence in wind tunnels, but can also appear at a heated wall in a setting with buoyancy-aided mixed convection (see e.g. [43] and references therein). Closely related, e.g. [44] successfully employs a 'moving surface boundary-layer control' to prevent or delay the separation of the boundary layer from the wall, where the moving surface is provided by rotating cylinders. Regarding pipe flow, [45] observed temporary relaminarization by fluid injected through a porous wall segment of the pipe. All these authors showed that a wall-bounded turbulent shear flow may relaminarize or, more accurately, transiently tend to a laminar-like or quasi-laminar state under certain suitable conditions, even if the Reynolds number is above criticality. However, in all cases inevitable retransition to turbulence was found at a later stage, rendering the use of the word relaminarization somehow misleading as it characterizes an intermediate but not the final result. Although various parameters have been proposed to quantify e.g. the acceleration level needed for relaminarization (for a compilation see [36]) there is neither agreement on a precise criterion for the occurrence of relaminarization and how it can be triggered, nor on how its onset may be recognized.
Profound understanding of ways to control turbulence has emerged from studies elucidating the near-wall flow dynamics. Obviously the near-wall region is crucial to the dynamics of attached shear flows, as it is the region of the highest rate of turbulent energy production and of the maximum turbulent intensities. The dominant structures of the near-wall region are the streamwise velocity streaks and the quasi-streamwise vortices, and the dominant dynamics is a cyclic process characterized by the formation of velocity streaks from the advection of the mean profile by streamwise vortices, and the generation of vortices from the instability of the streaks, also referred to as self-sustaining process (see e.g. [46][47][48][49]). By clever use of 'wrong physics' for numerical experiments on modified turbulent channels [50] showed that a local cycle of turbulence regeneration exists in the near wall region, which is independent of the outer flow. By numerically interrupting the cycle at various places they observed a breakdown of turbulence and eventual relaminarization. Also the aforementioned riblets are believed to directly weaken the quasi-streamwise vortices of turbulence regeneration, i.e. proposing that the drag reduction by riblets is due to the weakening of these vortices by the increase in spanwise friction at the wall [51].
Turbulence arises from finite amplitude perturbations. At the lowest Reynolds numbers where turbulence can be observed (given a large enough initial perturbation), it does so in the form of localized structures called puffs which have a constant length and advect along the pipe [52,53]. Hof et al. [54] observed an immediate collapse of single puffs in the intermittent regime at relatively low Reynolds numbers if two puffs were triggered too close to each other. They related the breakdown to a flattened streamwise velocity profile induced by the trailing puff, weakening the turbulence regeneration cycle beyond recovery. Their observation is also in qualitative agreement with predictions from a model explaining the emergence of fully turbulent flow in pipes and rectangular ducts [55]. From this model it was inferred that the excitability to the turbulent state only depends on the streamwise velocity component and hence an appropriate manipulation of the streamwise velocity profile may destroy the turbulent state.
Further investigations by [56] have recently shown that a modification of the velocity profile in a pipe by several means can lead to a complete collapse of turbulence and the flow can be forced to fully relaminarize also at higher Reynolds numbers. The annihilation of turbulence was achieved by a steady active open loop manipulation of the streamwise velocity component alone, greatly simplifying control efforts. In their numerical simulations [56] added an appropriate radially dependent body force term, F (r), to the equation of motion, modifying the streamwise velocity profile to a more plug-like one. Furthermore, they presented four different experimental techniques to modify the velocity profile of turbulent flow, such that the resulting profile was (more) plug shaped and flat or even had velocity overshoots in the near wall region of the pipe. The first technique employed four rotors located inside the pipe to vigorously stir the flow. The second technique used wall-normal injection of additional fluid through 25 small holes placed consecutively in a helical fashion around the pipe. The third technique was by injecting fluid through an annular gap at the wall to accelerate the flow close to the wall. The fourth approach was by means of a movable pipe segment which was used to locally accelerate the flow at the wall. In all cases turbulence was shown to decay and the flow eventually relaminarized completely, up to Reynolds numbers of 100,000 in DNS and 40,000 in experiments. The experiments demonstrated that relaminarization occurred as a direct result of a particular velocity and shear-stress distribution especially in the wall region. The modified profile was shown to specifically suppress transient growth such that vortices do not efficiently create streaks.
Pipes and pipelines are central for the distribution of fluids throughout society, ranging from small diameter tubes in domestic settings to pipes in industrial plants and to large scale pipelines. The overall pumping costs surmount to billions of Euros per year and the frictional losses encountered in these flows are responsible for a significant part of the global energy consumption. Thus a huge amount of energy (pumping power due to skin friction drag) could be saved if flows in pipes were laminar instead of turbulent. Pipe flow has a property that makes it particularly attractive for a way of turbulence control which totally extinguishes turbulence and establishes a laminar flow: the laminar state is stable to infinitesimal perturbations at all flow speeds [57]. Consequently, once relaminarization is achieved, the flow will remain laminar provided that the pipe is straight and smooth. Turbulence will only return if a sufficiently strong disturbance is encountered. It is therefore not necessary to apply relaminarization control techniques throughout the pipe but instead it suffices to implement control stations at fixed locations (for example behind bends) to ensure laminar flow in the straight downstream pipe sections.
In the present investigation we want to further explore the effect, scope and consequences of a modified streamwise velocity profile which relaminarizes the flow. In order to do so we insert an obstacle partially blocking the pipe (method 1) and inject fluid through an annular gap at the wall (method 2). Both approaches force the streamwise velocity profile in a similar manner. The inserted obstacle in a purely passive way, while the injection device can be considered steady open loop forcing. The obstacle device is better suited for measurements of the downstream flow field at various distances, as it can be easily mounted at any arbitrary position within the pipe. The device to inject fluid through an annular gap is better suited to continuously vary the amount of acceleration close to the wall and fathom out the maximum Reynolds number where relaminarization is possible. We present measurements with stereoscopic particle image velocimetry and investigate the development of the flow downstream these devices. Furthermore we try to elucidate the specific circumstances responsible for relaminarization. Both devices are solely used for demonstration and proofof-principle purposes. Detailed considerations on a possible net energy gain achievable with such devices are beyond the scope of the present manuscript.
The outline of the paper is as follows. In the next section we describe the experimental facility and the devices used to control the flow. Selected results of our extensive investigations and measurements are presented in Section 3. In Section 4 we discuss the results and compare them to previous investigations concerned with relaminarization.

Experimental Facility
The experimental setup consists of a basic pipe flow facility constituting fully turbulent flow in a straight long glass-pipe and two different kinds of interchangeable stationary flow management devices (FMDs) to control the flow. In the following section the facility and the different FMDs are described in detail.

Facility
The setup consists mainly of a glass-pipe with inner diameter D = 30 ± 0.01 mm and a total length of 9 m (300 D) made of 1 meter sections. Figure 1 shows a sketch of the facility. Gravity driven water enters the pipe from a reservoir. The flow rate Q m is measured by means of an electromagnetic flowmeter (fm 1) in the supply pipe. The Reynolds number (Re = UD/ν, where U is the mean velocity and ν the kinematic viscosity of the fluid) can be adjusted by means of a valve (v 1) in the supply pipe.
The velocity vector field is measured ∼ 250 D downstream from the inlet at the position of the lightsheet. The measurement plane is perpendicular to the streamwise flow direction (pipe z−axis). All three velocity components within the plane are recorded using a highspeed stereo PIV system (Lavision GmbH) consisting of a continuous laser (Fingco 532H-2W) and two Phantom V10 high-speed cameras with a full resolution of 2400 × 1900 px. Spherical glass particles (sphericel, mean diameter 13 μm, Potter Industries) are used to seed the flow. Around the measurement plane the pipe is encased by a water filled prism such that the optical axes of the cameras are perpendicular to the air-water interface to reduce refraction and distortion of the images.
A differential pressure sensor (DP 103, Validyne) is used to measure the pressure drop p between two pressure taps. It is used for a straightforward detection of the flow state based on the large difference between turbulent and laminar friction factors (f T and f L ). For this purpose the sensor is placed downstream the FMD as indicated in Fig. 1, and the pressure taps are separated by 23 D in the streamwise direction.
For reference a non-dimensional Cartesian coordinate-system (x, y, z) = (x/D,ỹ/D,z/D) is used as shown in Fig. 1. The origin of the coordinate system is always located at the downstream end of the FMD. The respective Cartesian velocity components (ũ,ṽ,w) are made non-dimensional either with the mean velocity U , yielding (u, v, w), or with the friction velocity w τ , yielding y + and w + . To calculate w τ we use w τ = U(f/8) 0.5 [58], where the friction factor f is obtained from pressure drop measurements. Velocity fluctuations are denoted by , <> denotes cross sectional averaging and a bar averaging over time.

Orifice plate obstacle (method 1)
The first device to control the flow is a stationary obstacle or inset which can be mounted within the pipe. Figure 2 shows a sketch of the device. It is made of a custom made The device can be mounted concentrically within the pipe at any axial position by means of three small streamwise ribs (interference fit). As a result 78.6% of the pipe are blocked at the upstream end of the device. The flow is divided into two separate parts. One part of the flow has to go through the 1 mm annular gap along the pipe wall, while the other part of the flow goes through the perforated plate in the bulk. The device thus forms a deliberate obstacle or obstruction acting as a spatially extended volume forcing on the flow. Its sole purpose is to tailor the velocity distribution of the flow at the downstream end of the device (in the plane z = 0) in a controlled way. For reference in the following text the device is referred to as obstacle-FMD.
Devices with several different blockage ratios, gap widths, more or less holes with different diameters have also been tested. Although there is certainly still room for further optimization, it was found that with reasonable effort the specific dimensions mentioned above seem to work best in terms of relaminarization capability. The device is quite sensitive to changes. E.g., already if the diameter d 3 of the holes is ≤ 2.8 mm or ≥ 4 mm instead of 3.3 mm and all other dimensions are left unchanged, the device is not able to relaminarize the flow at any Reynolds number.

Annular gap injection nozzle (method 2)
To complement and extend the purely passive control mechanism of the obstacle-FMD we designed a device which allows to inject fluid into the main pipe through a small annular concentric gap close to the pipe wall as shown in Fig. 3. The modification of the velocity profile resulting from this active open-loop control of the flow is similar to the stationary obstacle, but the amount of injected fluid and hence the level of acceleration close to the pipe wall can be controlled and continuously adjusted via a valve in the feeding line. In the reported setup the feeding line is implemented as a bypass that takes fluid from upstream the main pipe and re-injects it through the device such that no net mass flux is added to the system. However, the injection device could just as well be used without the bypass by injecting fluid from an external reservoir via a pump. This method was found to be equally effective and relaminarization up to similar Re values could be achieved. Fig. 3 Sketch of the control device acting as injection nozzle through an annular gap (1 mm-FMD and 2 mm-FMD). Fluid is taken from the main pipe via a bypass and then re-injected through a concentric gap close to the wall. The device is mounted within two pipe segments. The bypass is equipped with a valve (v 2) and a flowmeter (fm 2). The flow direction is from left to right. Drawing not to scale We tested two different devices by systematically varying the injected flow rate at different Reynolds numbers, one device with an annular gap of 1 mm and one with an annular gap of 2 mm. In the following these devices are referred to as 1 mm-FMD and 2 mm-FMD.
As sketched in Fig. 3 the main pipe is slightly narrowed in a short range just upstream the injection point (1 mm-FMD: . At a small backward facing step (z = 0, the plane of confluence) the fluid coming from the bypass is axially injected into the main pipe through the annular gap close to the wall. Based on the measured bypass flow rate Q bp and the open gap area the mean velocity U inj of the injected flow in the plane of confluence can be calculated. As the total flow rate Q m is measured in the supply pipe (see Section 2.1), the specified Reynolds number in the main pipe applies to the range upstream the bypass and downstream the confluence at z = 0. The Reynolds number Re n and the mean velocity U n in the slightly narrowed part of the device is calculated based on d 1 and Q m − Q bp .
In all measurements reported the flow rate was chosen such that Re n 3000 to ensure fully turbulent flow in the slightly narrowed section (as for Re n 2800 the main flow may become (intermittently) laminar already without any forcing [59]). For the 2 mm-FMD this criterion was met for Re 4000.

Results
In this section we present the results of the measurements split into two subsections for each kind of device. Representative sets of stereoscopic PIV-measurements are shown for Reynolds numbers where the flow relaminarizes and where it does not. The FMDs were first investigated at different Reynolds numbers by means of mere visualization. Neutrally buoyant anisotropic particles [60] were added as tracer particles and the flow in the pipe was illuminated by means of LED string lights along the whole length of the pipe to be able to observe the development of the flow field downstream of the FMDs. The observations were also recorded with a video camera. A selected example of a relaminarizing flow is made available in the Online Supporting Material.
From visual observations in combination with pressure drop measurements we obtained reliable information whether or not the flow completely relaminarized at a certain Reynolds number. Although relaminarization is a gradual process [27] and the flow tends only asymptotically to the fully developed laminar parabolic Hagen-Poiseuille profile [61], it is accompanied by drastic changes in the structure and dynamical behavior of the flow. The mean velocity profile departs from the well-known law of the wall, the friction factor exhibits a substantial decrease and Reynolds stresses become negligible, the turbulence intensity drastically reduces, etc. Visually most striking and conspicuous is the total decay of any eddying motion, as all velocities normal to the streamwise direction fade away.
More difficult to capture and quantify is the transient state of ongoing relaminarization, also referred to as laminarescent. It signifies the earlier stages of relaminarization -loosely, a precursor to relaminarization -in which large departures occur from the turbulent state [27]. However, remaining perturbations in the flow may trigger retransition to turbulence. It is only when the flow has reached a sufficiently laminar state that it remains laminar for the remainder of the pipe. As will be shown further down, monitoring the flow at z ≈ 100−150 provides a clear and unambiguous indication whether or not the flow has relaminarized completely.

Measurements of uncontrolled turbulent flow for reference
To be able to compare the controlled, relaminarizing flow to uncontrolled fully turbulent flow we first measured the flow without FMDs. The respective measurements are going to be used for reference in the remainder of the present work. The velocity profiles measured at Re = 3800, Re = 4500 and Re = 5000 normalized by the inner variables are plotted in Fig. 4 in the usual semi-logarithmic manner. For comparison, the law of the wall is plotted with the Kármán constant k = 0.4 and the additive constant A = 5.5 (to compensate the low-Reynolds number effect, see [62]). Furthermore, DNS results by [63] for pipe flow at Re = 5300 are shown. Eggels et al. [63] demonstrated that at low Reynolds numbers the mean velocity profile in a pipe fails to conform to the law of the wall, in contrast to channel flow.  [62], w + = 2.5 ln y + + 5.5) and DNS in a pipe from [63]. The viscous sublayer is indicated by w + = y + This is reaffirmed by our measurements. While for Re = 3800 our measured velocity profile is slightly above the DNS data of [63], their data coincide very well with our results for Re = 4500 and Re = 5000. Eggels et al. [63] stated a best-fit of k = 0.35 with an additive constant A = 4.8. For Re = 5000 we find a best-fit of k = 0.33 and A = 4.2 for 40 < y + < 120 or k = 0.34 and A = 4.6 for 40 < y + < 100. The near wall region is only partly resolved in the measured velocity profiles. The smallest values of y + are 7.9, 9.1 and 10 respectively, corresponding to ≈ 0.08 mm.
The rms values of the streamwise fluctuating velocity w rms = (w 2 ), normalized by the friction velocity, are plotted in Fig. 5 for Re = 5000. For comparison, DNS results by [63] for pipe flow at Re = 5300 are shown. The peak value (2.6) of the measured rms is in excellent agreement and the position of the peak (y + = 18.2) as well as the overall trend of the data is in good accordance with previous investigations (see also [64]).

Orifice plate obstacle (inset)
For the visual observation of the flow field we placed the obstacle-FMD 80 D downstream of the pipe entrance where the turbulent flow can be considered fully developed. We varied the Reynolds number between 2500 and 5000 and found a continuously laminar flow field at z = 150 up to Re ≈ 4200. At slightly higher Reynolds numbers we observed laminar-turbulent intermittency. For Re 4400 we found only turbulent flow at z = 150. Inspection of the flow field in the closer downstream vicinity (0 < z 10) of the obstacle-FMD shows a rather turbulent zone right downstream the device where strong cross-stream motion is notably present. But for Re 4200 all perturbations in the flow field decay until the flow appears clearly laminar at z ≈ 50. This evolution is well observable in Fig. 6,  [65]), where the camera follows the flow at approximately the mean velocity on its journey from upstream of the device until it has relaminarized downstream. A laminarescent part of the flow as in Fig. 6 at z = 15, where the turbulence intensity is obviously reduced by some amount in the downstream vicinity of the device, could be observed up to Re ≈ 4800.
To characterize the flow downstream the device we took PIV-measurements for two representative cases at Re = 3800 and Re = 4500 and 14 axial stations downstream the FMD: at z = 2, 3, 4, 6,9,12,15,18,24,30,40,50,70 and 100. The development of the time averaged streamwise velocityw in the plane y = 0 is shown in Fig. 7a and b. Each profile is calculated from 200 independent vector fields and azimuthally and temporally averaged. For reference, the most left profile (ref) is showing the respective measurement of uncontrolled turbulent flow. Figure 7c provides the respective levels of the centerline velocity for quantification. As can be seen from the figure, at Re = 3800 the flow fully relaminarizes, while at Re = 4500 retransition to turbulence takes place.
A particularly noteworthy effect of the device on the time-averaged velocity profiles is the increase in the velocity gradient at the wall. At z = 2 the profiles at both Reynolds numbers exhibit characteristic overshoots of fast fluid close to the wall (peak at x ± 0.44) and a clearly decreased centerline velocity (∼ 1 and ∼ 0.9 respectively) causing a M-shaped At this Reynolds number it takes yet another ∼ 120 D [61] to reach the fully developed state under perfect conditions (the laminar profile is sensitive to small disturbances, e.g. due to thermal convection or minor misalignments of the pipe segments [66]). But the centerline velocity of the flow at Re = 4500 keeps increasing only until z ≈ 30. The mean turbulent centerline velocity is clearly surpassed at z 20. However, in the range 30 z 40 the centerline velocity starts to fall back toward the turbulent mean.
The scenario of complete relaminarization at Re = 3800 and of final retransition to turbulence at Re = 4500 is also apparent from the axial development of the streamwise and inplane turbulence intensity ti w = < w rms U > and ti u,v = < u rms + v rms 2U > as plotted in Fig. 8a and b. Right downstream the device the levels of turbulence intensity are slightly increased above the turbulent reference flow. It should be noted however that in the bulk region of the cross-section right downstream of the obstacle-FMD remaining flow perturbations induced by the perforated plate at the upstream end of the obstacle-FMD may still be significant, as L tot = 200 mm (see Fig. 2) is not sufficiently long to regain streamwise invariance (see e.g. [67,68]). The increase in turbulence intensity can hence not be clearly attributed to the accelerated part of the flow coming through the gap. Already at z 3 (Re = 4500) ti w drops below the level of uncontrolled turbulent flow. The same applies regarding ti uv at z = 4 (Re = 3800) and z = 6 (Re = 4500) respectively. While for Re = 3800 ti uv decreases exponentially until z ∼ 30 where the noise level of the sensor is reached, a qualitatively similar decrease in ti uv at Re = 4500 is perceived only until z = 15. For 15 z 30 ti uv is almost constant, followed by a steep increase which even leads to an overshoot above the mean turbulent level for 40 z < 70. Further downstream the flow seems to resemble an uncontrolled turbulent flow again. ti w exhibits a qualitatively very similar development, but the increase back to the turbulent level at Re = 4500 starts a little earlier (around z = 15) and the overshoot above the turbulent level takes place more rapidly and pronouncedly (+ 30% at z = 40). At Re = 3800 ti w seems to reach the laminar (noise) level in the range z = 10 − 15. The remaining velocity fluctuations can be attributed to measurement noise.
To elucidate the downstream development of turbulence fluctuations within the crosssection of the pipe, the squared streamwise velocity fluctuations (w 2 rms ) at selected axial stations are depicted in Fig. 9. The reference measurement of uncontrolled turbulent flow (see also Fig. 5) at the respective Reynolds number is plotted as a dashed line. In the core region of the cross-section right downstream of the obstacle-FMD the flow is only slightly changed for both Reynolds numbers compared to the uncontrolled flow, recognizable by the increased rms level for x 0.2 at z = 2. This is apparently due to the perturbations caused by the perforated plate at the upstream end of the obstacle-FMD. Closer to the wall a more prominent deviation from the uncontrolled flow is visible, depicting also a clear cut qualitative difference between the measurements at Re = 3800 and Re = 4500. For Re = 3800 the maximum of streamwise fluctuations close to the wall has decreased below the level of uncontrolled flow already at z = 2. Additionally, the maximum has significantly moved toward the centerline (apparent until z = 4). As z increases there is a substantial reduction in turbulence level within the whole cross-section. Initially the turbulent fluctuations decay very quickly close to the wall and more slowly in the bulk region. For z > 30 the flow can be considered laminar and the overall turbulence intensity does not change significantly anymore (see also Fig. 8).
For Re = 4500 on the other hand the maximum of streamwise fluctuations at z = 2 has not dropped significantly and is comparable to uncontrolled turbulence. A double-humped, much thickened shape of the near-wall peak can be observed instead (visible up to z = 4). However, for z = 4 − 24 the peak close to the wall clearly decreases beyond the level of uncontrolled flow and the overall turbulence intensity is substantially reduced. Yet, at z = 50 the near wall region looks similar to the reference flow again. In the core region the rms level is even increased beyond the reference flow, yielding the increased level in turbulence intensity as depicted in Fig. 8. Finally, at z = 100, the flow has completely returned to the uncontrolled level.

Annular gap injection nozzle
With the 1 mm-FMD and the 2 mm-FMD, the devices which allow to inject fluid into the main pipe through an annular gap close to the pipe wall, the amount of injected fluid (Q bp ) and hence the level of acceleration close to the pipe wall can be controlled and continuously adjusted via a valve in the bypass. To investigate which flow rates and mean flow velocities of injected flow result in complete relaminarization, we increased the Reynolds number in increments of 500 and varied the flow rate through the bypass, i.e. the amount of fluid injected through the gap. Up to Re ∼ = 6000 we could find injection flow rates for both gap widths which cause a continuously laminar flow field at z = 150 and the remainder of the pipe. I.e., when we measure the pressure drop sufficiently far downstream we find either turbulent or laminar pressure drop depending on the respective injection flow rate, nothing in between. Figure 10 shows the measured ratios Q bp /(Q m − Q bp ) of the flow rates of injected flow to the main flow for each gap width which lead to relaminarization. Figure 11 shows the respective velocity ratios U inj /U n of injected flow to the main flow. In both figures max and min depict the maximum and minimum flow rate (velocity) ratios that would cause total relaminarization downstream. All flow rate (velocity) ratios in between also relaminarize the flow. While at lower Reynolds numbers a relatively broad range of flow rate ratios is suitable for relaminarization, the range quickly narrows down with increasing Reynolds number. For both devices it is rather the maximum of the flow rate ratio which decreases from initially ∼ 30% (1 mm-FMD) and ∼ 35% (2 mm-FMD), while the minimum amount of injected flow stays relatively constant at 18 and 25% for the 1 mm-FMD and 2 mm-FMD respectively. Concerning the velocity ratio this implies possible injection velocities in a broad range of 1.5-3 (1 mm-FMD) and 1-1.5 (2 mm-FMD) times higher than the mean flow velocity in the pipe at Re = 4000. At Re = 6000 the injection velocities necessary for total relaminarization narrow down to 1.5 and 1 as compared to the velocity of the main flow. For both gap widths the curves intersect slightly above Re = 6000, indicating that for higher Reynolds numbers no complete relaminarization is possible with the present devices. Visual inspection of the flow field in the downstream vicinity (0 < z 50) of the FMDs at Reynolds numbers and flow rate ratios between min and max shows a scenario which is very similar to the one during relaminarization downstream the obstacle-FMD as described in Section 3.2: the flow field in the direct downstream vicinity (0 < z 10) of the device shows a turbulent flow and turbulence levels look comparable to the level upstream the device. However, for Re 6000 all visible perturbations in the flow field quickly decay downstream until the flow is clearly perceivable laminar at z ≈ 30 − 50. Furthermore, the flow stays unambiguously laminar for the remainder of the pipe.
If the ratio of injected fluid is just below the minimum or just above the maximum depicted in Fig. 10, a transiently relaminarizing (laminarescent) part of the flow can be observed as in Fig. 6 at z = 15, where the turbulence intensity is obviously reduced by some amount in the downstream vicinity of the device. I.e., incipient (transient) relaminarization at earlier stages occurs, but finally the flow always returns to a fully turbulent state at around z ≈ 30 − 50.
As the scenario for injection flow rates between the indicated minimum and maximum and Re 6000 is so similar to the one downstream the obstacle-FMD at Re = 3800 (see e.g. Figs. 6 and 7) and for flow rates closely below the minimum or above the maximum similar to to the one downstream the obstacle-FMD at Re = 4500 respectively, no further pictures or measurements of the downstream evolution are shown for the 1 mm-FMD and the 2 mm-FMD. Figure 12 presents PIV-measurements of 4 different injection flow rates with the 2 mm-FMD at a single axial station, namely z = 2.5, which is the closest distance to the device possible for our stereoscopic PIV-measurements. The four cases c1-c4 of different injection flow rates are indicated in Fig. 10b by the symbols at Re = 4000 and Re = 5000. c1 represents a case where the flow rate of injected fluid is just beneath the minimum flow rate which is necessary for full relaminarization. c2 and c3 represent cases where full relaminarization is observed. At c4 the flow rate is already somewhat above the maximum, meaning that the controlled flow exhibits features of relaminarization but finally returns to a turbulent state downstream (similar to c1). The streamwise velocity profiles in Fig. 12 clearly  Fig. 10b by a . A measured profile of uncontrolled turbulent flow is shown for reference (ref) show the increasing injection flow rate close to the wall. While for c1 only a minor hump is visible, for c2-c4 the peak close to the wall (stably at x ∼ = 0.41 for all injection flow rates and both Reynolds numbers) is manifest. Again, the appearance of the velocity profiles (at least for c4) could be described as M-shaped with a pronounced plateau in the center.
The respective streamwise velocity fluctuations at z = 2.5 are plotted in Fig. 13. For the injection case c1 the near wall peak is greatly reduced compared to the reference flow and moved from x ∼ = 0.44 to x ∼ = 0.33 for both Reynolds numbers. All streamwise fluctuations close to the wall have almost disappeared, and in the core region the rms level is already reduced too. However, the flow does not relaminarize in the long term. For injection cases c2 and c3 the reduction of the near wall peak is even more distinct, especially at Re = 5000, where no significant peak is present at all. Also the reduction in the core region is more pronounced. Although for c3 the fluctuation level close to the wall is increased as compared to c2, noticeably especially for Re = 4000, c2 and c3 show clearly reduced fluctuation levels and accordingly relaminarize completely downstream. For the injection case c4 it seems that, although the absolute level of the near wall peak is still reduced below the uncontrolled turbulent reference, a double humped, thickened shape of the peak can be observed. c4 eventually turns fully turbulent further downstream.
For Reynolds numbers > 6000 we visually also observed a transiently relaminarizing state in the section z 6 − 30 up to Re ≈ 10000, similar to the scenario described in pronounced. The higher the Reynolds number the higher the turbulent fraction, where laminar patches would quickly shrink downstream due to the faster propagation of the turbulent fronts [55].
To further investigate the evolution of the flow right downstream the injection device and detect temporally and spatially confined (partial) relaminarization we measured the mean streamwise pressure gradient between two pressure taps at z = 6 and z = 29 respectively (i.e. over a 23 D stretch). The first pressure tap is supposed to be sufficiently far away from z = 0 to avoid direct influence from the backward facing step of the injection device. At the same time the second pressure tap is supposed to be sufficiently close to the device to cover the whole area of transient relaminarization suggested by visual observations and Figs. 7 and 8.
The measured variation of p with the injection flow rate at several Reynolds numbers is plotted in Fig. 14. For all cases depicted we find significant drag reduction at certain injection flow rates, and the trend is similar for increasing the injection from zero: above an injected flowrate of ∼ 15% a clear drag reduction can be observed, suggesting the onset of relaminarization for flowrate ratios roughly between 10-15% for this specific device. The highest decrease in p is found for flowrate ratios around ∼ 20-25%. With even higher injection flow rates p returns toward the initial uncontrolled value. For Re = 4000 and Re = 5000 p decreases by up to ∼ 30% for injection flow rates of ∼ 20-25%. Interestingly, the maximum drag reduction over the 23 D stretch downstream of the injection point is not achieved for those injection flow rates necessary for complete relaminarization further downstream (indicated by filled symbols, see also Fig. 10), but for flow rates which are slightly below.
For Re = 7000 and Re = 9000 the drag reduction is even larger (down to 65% of the uncontrolled flow), despite the fact that the flows return to turbulence further downstream. At Re = 11000 and Re = 13000 the decrease in p is still significant (by 18 and 12% respectively), yet already much less. However, the measurements clearly indicate transient relaminarization also at these higher Reynolds numbers, i.e. temporary relaminarization in a spatially confined region right downstream the injection device.

Discussion
Our results demonstrate that by increasing the flow velocity close to the pipe wall in the proposed way turbulence can either be completely annihilated or at least temporally (in a spatially confined region downstream the control) weakened considerably. A chief precondition to the relaminarization process seems to be a steady, homogeneous local acceleration of the region of the viscous sublayer and the buffer layer above the level of uncontrolled flow, while the velocity is decreased accordingly in the log law region. Once this profile modification is realized, we observe an immediate collapse of turbulence production at the wall. Thereafter turbulence intensity decays exponentially with z. The most drastic reduction in turbulence levels is found in the region near the pipe wall. If the injection (forcing of the velocity profile) is too weak or too strong or Re 6000 several effects like remaining perturbations in the flow or secondary circulation caused by the injection trigger turbulence again and inevitably lead to a turbulent flow further downstream.
The results are in good accordance with the principal process of relaminarization in previous investigations concerning the natural decay of turbulence at low (subcritical) Reynolds numbers. E.g. [69], who investigated relaminarization by reducing the initial Reynolds number, found exponential decay of the turbulence intensity in the streamwise direction too. [70] also noted that the rate of decay of turbulence fluctuations appears to be more rapid near the wall and in the middle region of the pipe rather than at intermediate positions.
Pennell et al. [45], who similarly observed temporal relaminarization triggered by fluid injected through a porous-walled pipe at Reynolds numbers above criticality, noticed a thickening of the viscous and buffer layers due to the start of injection. This caused a significant yet transient reduction of the turbulence level. In the later stages of the re-transition to full turbulence (already at 6 and 10 D downstream the injection) their rms-profiles exhibit a very similar double-humped shape as we observe in Fig. 9b at z = 2 and z = 3 and Fig. 13 for case 4. Note that the cases in which double peaks are present in the rms-profiles coincide with the existence of a strongly M-shaped velocity profile with strong velocity overshoots at the wall, and that all these profiles belong to cases which return to a full turbulent state further downstream. I.e., one peak seems to be related to the near-wall turbulence generation and the other to the inner-flank of the M-shape. This is consistent with [45] who suggested that the double-humped rms-profile is a definite characteristic of the re-transition process (to turbulence) in a pipe, and it seems also consistent with the observations of e.g. [71].
Furthermore, the qualitative trend of our data is very similar to observations by e.g. [35,[37][38][39][40][41][42] and [36] concerning the effect of a favorable pressure gradient (FPG) on a turbulent boundary layer. Under the (accelerating) influence of a FPG the viscous sublayer of a turbulent boundary layer is known to increase in thickness and the velocity in the outer region decreases. Measurements and simulations of boundary layers under a FPG have been characterized to exhibit reversion from the turbulent to the laminar state. It is also the FPG and its effect on the wall region of the flow which has been found to be primarily responsible for departures from the inner law (and, by inference, for triggering relaminarization).
Various parameters have been suggested to quantify the acceleration level in spatially accelerating flows and define the onset of relaminarization, see [36] for a compilation of 'acceleration driven laminarization parameters'. While there is not a consensus in criteria of laminarization, one of the most widely used is the acceleration parameter where z is the streamwise direction, U ∞ is the free-stream velocity in this direction and ν the kinematic viscosity. Several authors have reported that turbulence is not sustained if K is higher than a critical value in the range 2.5-3.6 × 10 −6 [35,41,72,73]. At the same time most researchers have also criticized the use of K for being based on bulk flow parameters while the laminarization phenomenon is assumed to be a boundary layer event, necessarily having its basis in boundary layer considerations. Based on Eq. 1 we calculated the acceleration parameter K = ν w 2 w z for different injection flow rates by FMD 2 at Re = 5000 with z = 2.5, i.e. we compare the (local) acceleration level at z = 2.5 to the uncontrolled turbulent flow. Figure 15 shows the result for selected flowrate ratios (0.12-0.33 in increments of 0.03) as a function of the radial direction. Furthermore, the aforementioned (global) level of K ≈ 3 × 10 −6 is indicated for reference. Note that at arbitrarily small z close to the backward-facing step of the injection device the values for K in the near-wall region would be even higher. Figure 15 suggests that in the near wall region the calculated values for K roughly match with the critical values proposed in the literature, in particular when taking into account the findings provided by Fig. 14, indicating the onset of relaminarization for flowrate ratios between 10 and 15% for FMD 2. It would need measurements much closer to z = 0 for a more accurate assessment.
However, the major drawback of K is in any case the empirical nature and the lack of a mechanistic explanation. Warnack and Fernholz [42] also pointed out that no single criterion could be used to predict the beginning or end of the breakdown of the law of the wall. They found that for flow with low acceleration where laminarization did not occur, there still was a breakdown of the law of the wall and a slight reduction in turbulent intensities. Relaminarization in boundary layers was observed as a gradual change of the turbulence properties and not catastrophic. Retransition, however, is a fast process due to the remaining turbulence structure and may be compared with bypass transition [42]. An interesting approach to the question how the specific conditions can cause the flow to laminarize on a local scale is provided by [50]. They have shown not only that the near-wall turbulence regeneration cycle is autonomous, i.e. self-sustaining independently of whether or not fluctuations are present in the region y + > 60, but also that it can be interrupted numerically at various places, leading to the decay of turbulence and to eventual laminarization. Using a filter to make the streamwise velocity (more) uniform they identified the minimum streak length needed to sustain the cycle to be between 300 and 400 wall units. Their proposed control strategy is hence to weaken or decorrelate the streaks in the region below y + = 60 and above y + = 20.
Instead of using a numerical filter on the near-wall region, we locally replace and substitute the flow by injection through the gap. The injected flow and the resulting, modified flow in the near-wall region is supposedly more homogeneous and 'streak-free'. In support of this, Fig. 16a shows several distinct (high-speed) streaks in the near wall region (the region roughly below y + = 60 marked by the dashed circle) of an instantaneous snapshot of turbulent flow at Re = 4000. However, in (b), depicting the flow field right downstream (at z = 2.5) of the 2 mm-FMD for the case c 2 (i.e. relatively small injection just sufficient for full relaminarization), the streaky structures in the near-wall region are replaced or superimposed by the injected flow coming from the gap at the wall.
Kühnen et al. [56] designate the relaminarization observed due to a modified profile also to a weakening of the near-wall cycle, yet argue from a slightly different perspective focusing on the efficiency of the 'lift up mechanism' [46]. To obtain a measure for the amplification mechanism of the regeneration cycle they consider the linearized Navier-Stokes equations and perform a transient growth (TG) analysis (following the algorithm given by [74]). The velocity profiles of all successfully modified, i.e. relaminarizing flows considered by [56] are shown to exhibit a substantially reduced transient growth. Generally the flatter the velocity profile the more the streak vortex interaction is suppressed and in the limiting case of a uniformly flat profile the lift up mechanism breaks down entirely.
We applied the same procedure here to (azimuthally and temporally averaged) velocity profiles measured at z = 2.5 downstream of the 2 mm-FMD for different injection flow rates. Assuming that the profile is fixed under the influence of the perturbation we conducted a TG analysis around the modified profiles. As shown in Fig. 17 the profiles indeed show a considerably decreasing TG with increasing injection flowrate suggesting that   Fig. 14) the TG exhibits a clear minimum amplification at an optimal flowrate ratio before it starts to rise again toward the value of unmodified flow. Interestingly, the optimal is found at slightly higher injection flowrates than those for the smallest pressure drop in Fig. 14, but z = 2.5 is not necessarily the axial location of the minimum transient growth. More refined measurements would be necessary to find the axial location of the minimum TG along z. In any case, the flowrate ratios with the least amplified perturbations are more or less exactly those which relaminarize downstream at Re = 4000 and Re = 5000 (indicated with filled symbols). The data suggest that the turbulence regeneration cycle is weakened indeed and breaks down due to a reduced efficiency of the lift up process. For a more detailed discussion of the reduced TG we refer to [56]. The present study elucidates the method and the general process of relaminarization caused by a modified streamwise velocity profile by simple means. Both the stationary obstacle-FMD and the injection nozzles (1 mm-FMD and 2 mm-FMD) can cause total relaminarization of the flow downstream and by that reduce skin friction by large amounts. As laminar pipe flow is linearly stable [57] it will stay laminar for the remainder of a smooth straight pipe. Although both devices cause an initial cost in terms of an increased pressure drag either by the pipe blockage or by the required injection for full relaminarization, it is only a question of how long the unperturbed (straight and smooth) downstream pipe section is until a net energy gain can be realized. In the present context both devices are used solely for demonstration and proof-of-principle purposes. A lot of optimization regarding the added pressure drag of the methods could be done. Detailed considerations on energy savings and a possible net energy gain achievable with the present devices are beyond the scope of the present manuscript.

Summary and Conclusion
Initially, relaminarization due to flow acceleration close to the wall may seem counterintuitive. However, in this work (see also [56]) we demonstrate that fully turbulent flows in pipes can be completely relaminarized by rather simple means. An open-loop forcing of the streamwise velocity profile, either actively or passively, is sufficient to trigger turbulence breakdown. The present experiments demonstrate that the onset of relaminarization in a fully developed turbulent pipe flow occurs as a direct result of a particular shear-stress distribution in the wall-region. The underlying physical mechanism of relaminarization is attributed to a weakening of the near-wall turbulence production cycle.
Two different devices, a stationary obstacle (inset) and a device to inject additional fluid through an annular gap close to the wall, are used to control the flow. Both devices modify the streamwise velocity profile such that the flow in the center of the pipe is decelerated and the flow in the near wall region is accelerated. Visualization, pressure drop measurements and stereoscopic PIV measurements have been employed to determine both practice-oriented and fundamental information for relaminarizing turbulent flows in a circular pipe in the presence of devices which accelerate the near wall region. Downstream of the devices a fully developed laminar flow is established. The fluctuations are greatest just beyond the devices and die away with increasing downstream distance.
High amounts of energy (pumping power due to frictional drag) can be saved if the flow is laminar instead of turbulent. At Re = 6000, the highest Reynolds number where we achieve full relaminarization with the present devices, the skin friction in the downstream distance is reduced by a factor of 3.4 due to relaminarization. At e.g. Re = 13,000, where we find transient relaminarization in a spatially confined region right downstream the devices, the measured local drag reduction is still higher than 10%.
Future research should focus on the possible net energy gain and an optimization of the specific design in terms of pressure loss, such that an energetic break even point can be reached in a close downstream distance of the devices. Furthermore, a smarter design may allow relaminarization at much higher Reynolds numbers. In order to establish a possible cost saving potential of the presented control technique it is also necessary to determine over which distances the relaminarized flow persists under less perfect pipeline conditions.