Search

Dean's approximation for curved pipe flow, valid under loose coiling and high Reynolds numbers, is extended to study three-dimensional travelling waves. Two distinct types of solutions bifurcate from the Dean's classic two-vortex solution. The first type arises through a supercritical bifurcation from inviscid linear instability, and the corresponding self-consistent asymptotic structure aligns with the vortex-wave interaction theory. The second type emerges from a subcritical bifurcation by curvature-induced instabilities and satisfies the boundary region equations. Despite the subcritical nature of the second type of solutions, it is not possible to connect their solution branches to the zero-curvature limit of the pipe. However, by continuing from known self-sustained exact coherent structures in the straight pipe flow problem, another family of three-dimensional travelling waves can be shown to exist across all Dean numbers. The self-sustained solutions also possess the two high-Reynolds-number limits. While the vortex-wave interaction type of solutions can be computed at large Dean numbers, their branch remains unconnected to the Dean vortex solution branch., Comment: 24 pages, 14 figures

It has been observed that flattening the mean velocity profile of pipe flow by body force can laminarise turbulence, a promising means to reduce frictional drag substantially. To explore whether there is a more efficient body force to eliminate turbulence, we consider time-independent active body forces with varying spatial dependencies. Results confirm that when using an active force, a flattened forced laminar profile is needed to eliminate turbulence, and that a purely streamwise body force is best for laminarisation. While these results required an expensive nonlinear optimisation, it was also observed that the optimal forced profile exhibits reduced linear transient growth (TG). To determine whether the reduction of linear TG alone is a sufficient target for the laminarisation of turbulence, a linear Lagrange Multiplier technique is used to minimise TG of perturbations to the forced laminar profile. The optimal velocity profiles reveal that TG for each azimuthal wavenumber strongly depends on the radial velocity gradient at some specific radial interval. The optimal velocity profile obtained by minimising TG of perturbations of azimuthal wavenumber m = 1 is shown to be able to eliminate turbulence at Re = 2400, but a more effective reduction of TG, including perturbations of higher m, is needed for laminarisation at higher Reynolds number Re = 3000. The streaks formed with the flattened forced laminar profile reveal the mechanism of laminarisation: it is found that the lift-up mechanism in the more flattened forced laminar profile creates a lower streak (i.e. closer to the wall) that is more stable, so that turbulence is harder to maintain through streak instability. Further numerical experiments by disturbing a periodic orbit validated that the breakdown of turbulence self-sustaining mechanisms is mainly caused by suppression of the formation of streaks.

The simulation of turbulent flow requires many degrees of freedom to resolve all the relevant times and length scales. However, due to the dissipative nature of the Navier-Stokes equations, the long-term dynamics are expected to lie on a finite-dimensional invariant manifold with fewer degrees of freedom. In this study, we build low-dimensional data-driven models of pressure-driven flow through a circular pipe. We impose the `shift-and-reflect' symmetry to study the system in a minimal computational cell (e.g., smallest domain size that sustains turbulence) at a Reynolds number of 2500. We build these models by using autoencoders to parametrize the manifold coordinates and neural ODEs to describe their time evolution. Direct numerical simulations (DNS) typically require on the order of O(105) degrees of freedom, while our data-driven framework enables the construction of models with fewer than 20 degrees of freedom. Remarkably, these reduced order models effectively capture crucial features of the flow, including the streak breakdown. In short-time tracking, these models accurately track the true trajectory for one Lyapunov time, while at long-times, they successfully capture key aspects of the dynamics such as Reynolds stresses and energy balance. Additionally, we report a library of exact coherent states (ECS) found in the DNS with the aid of these low-dimensional models. This approach leads to the discovery of seventeen previously unknown solutions within the turbulent pipe flow system, notably featuring relative periodic orbits characterized by the longest reported periods for such flow conditions.

Direct numerical simulations (DNS) of rotating pipe flows up to $Re_{\tau} \approx 3000$ are carried out to investigate drag reduction effects associated with axial rotation, extending previous studies carried out at a modest Reynolds number (Orlandi & Fatica 1997; Orlandi & Ebstein 2000). The results show that the drag reduction, which we theoretically show to be equivalent to net power saving assuming no mechanical losses, monotonically increases as either the Reynolds number or the rotation number increases, proportionally to the inner-scaled rotational speed. Net drag reduction up to about $70\%$ is observed, while being far from flow relaminarisation. Scaling laws for the mean axial and azimuthal velocity are proposed, from which a predictive formula for the friction factor is derived. The formula can correctly represent the dependency of the friction factor on the Reynolds and rotation numbers, maintaining good accuracy for low-to-moderate rotation numbers. Examination of the turbulent structures highlights the role of rotation in widening and elongating the small-scale streaks, with subsequent suppression of sweeps and ejections. In the core part of the flow, clear weakening of large-scale turbulent motions is observed at high Reynolds numbers, with subsequent suppression of the outer-layer peak in the pre-multiplied spectra. The Fukagata-Iwamoto-Kasagi decomposition indicates that, consistent with a theoretically derived formula, the outer layer yields the largest contribution to drag reduction at increasingly high Reynolds numbers. In contrast, both the inner and the outer layers contribute to drag reduction as the rotation number increases.

We demonstrate the first successful non-invasive stabilisation of nonlinear travelling waves in a straight circular pipe using time-delayed feedback control (TDF). The main novelty in this work is in the application of a "multiple time-delayed feedback" (MTDF) approach, where several control terms are required to attenuate a broad range of unstable eigenfrequencies. We implement a gradient descent method to dynamically adjust the gain functions, $G_i(t),$ in order to reduce the need for tuning a high dimensional parameter space. We consider the effect of the control terms via an approximate linear analysis and a frequency-domain analysis, which justifies the choice of MTDF parameters. By using an adaptive approach to setting the control gains, we advocate for a broad range of MTDF terms which are then free to find their own stabilising amplitudes for a given target solution, thereby removing trial and error when the instability of the target solution is unknown a'priori. This also enables travelling waves to be stabilised from generic turbulent states and with speculative starting parameters.

We consider the influence of a transverse magnetic field on the transient growth of perturbations in a liquid-metal circular pipe flow with an electrically insulating or conducting wall. In this configuration, the mean flow profile and the amplification of perturbations are strongly affected by the applied magnetic field, leading to a rich dynamical landscape depending on its intensity. The analysis is performed for Reynolds numbers 5000 and 10 000, close to the transitional regime for moderate values of the Hartmann number, a non-dimensional parameter proportional to the applied magnetic field's intensity. Aside from a slight modification of hydrodynamic optimal perturbations at very small Hartmann numbers, we observe three other characteristic topologies of optimal perturbations depending on the intensity of the magnetic field. Their growth mechanisms differ, with the lift-up effect dominating at low Hartmann numbers and the Orr-mechanism becoming increasingly important as the magnetic field intensity is increased. In particular, we show in the intermediate regime of Hartmann numbers how transient growth occurs in two stages, with initial growth through the lift-up effect followed by a further increase of energy through the Orr-mechanism. We also conduct three-dimensional nonlinear simulations to track the time evolution of optimal perturbations, illustrating their nonlinear growth and eventual breakdown to a sustained turbulent state or the return of the system to a laminar state.

We simulated a turbulent pipe flow within the Lattice Boltzmann Method using a multiple-relaxation-time collision operator with Maxwell-Boltzmann equilibrium distribution expanded, for the sake of a more accurate description, up to the sixth order in Hermite polynomials. The moderately turbulent flow ($Re_{\tau} \approx 181.3$) is able to reproduce up to the fourth statistical moment with great accuracy, compared with other numerical schemes and with experimental data. A coherent structure identification was performed based on the most energetic streamwise turbulent mode, which revealed a surprising memory effect related to the large scale forcing scheme used to trigger the turbulent state in the pipe. We observe that the existence of large scale motions which are out of the pipe's stationary regime do not affect the detailed single-point statistical features of the flow. Furthermore, the transitions between the coherent structures of different topological modes were analyzed as a stochastic process. We find that for finely resolved data the transitions are effectively Markovian, but for larger decimation time lags, due to topological mode degeneracy, non-Markovian behavior emerges, in agreement with previous experimental studies.

In this work, we simulate the transition to turbulence in the pipe flow based on the modified NS theory incorporating the viscous fluid strength in the constitutive equations. The latter concept enriches theory by allowing for material instabilities in addition to the kinematic ones. We present results of comparative numerical simulations based on the classical NS model and the NS model enhanced with the finite viscous strength. As expected, simulations based on the classical NS model exhibit stable laminar flow in contrast to experimental observations. Conversely, simulations based on the modified NS model with viscous strength exhibit instabilities and transition to turbulence per experimental observations. The transition to turbulence is triggered by the growing material instabilities.

It is well known that buoyancy suppresses, and can even laminarise turbulence in upward heated pipe flow. Heat transfer seriously deteriorates in this case. Through a new DNS model, we confirm that the deteriorated heat transfer within convective turbulence is related to a lack of near-wall rolls, which leads to a weak mixing between the flow near the wall and centre of pipe. Having surveyed the fundamental properties of the system, we perform a nonlinear nonmodal stability analysis. it is found that, the minimal seed becomes thinner and closer to the wall, with increase of buoyancy number C. Most importantly, we show that the critical initial energy required to trigger shear-driven turbulence keeps increasing, implying that attempts to artificially trigger it may not be an efficient means to improve heat transfer at larger C. The new minimal seed, found at C=6, is localised in streamwise direction and is active in the centre of pipe. To find this branch of optimal, we took advantage of a window of linear stability. While the nonlinear optimal causes transition to convective turbulence directly at this and larger C, transition via the linear instability passes via a travelling wave or periodic orbit solutions. Detailed analysis of the periodic solution reveals three stages: growth of the unstable eigenfunction, the formation of streaks, and the decay of streaks due to suppression of the instability. Flow visualization at C up to 10 also show similar features, suggesting that convective turbulence is sustained by these three typical processes.

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.

Large-scale coherent structures are identified in turbulent pipe flow at $Re_\tau=181$ by having long lifetimes, living on large scales and travelling with a certain group velocity. A Characteristic Dynamic Mode Decomposition (CDMD) is used to detect events which meet these criteria. To this end, a temporal sequence of state vectors from Direct Numerical Simulations are rotated in space-time such that persistent dynamical modes on a hyper-surface are found travelling along its normal in space-time, which serves as the new time-like coordinate. Reconstruction of the candidate modes in physical space gives the low rank model of the flow. The modes within this subspace are highly aligned, but are separated from the remaining modes by larger angles. We are able to capture the essential features of the flow like the spectral energy distribution and Reynolds stresses with a subspace consisting of about 10 modes. The remaining modes are collected in two further subspaces, which distinguish themselves by their axial length scale and degree of isotropy.

Turbulent concentric coaxial pipe flows are numerically investigated as canonical problem addressing spanwise curvature effects on heat and momentum transfer that are encountered in various engineering applications. It is demonstrated that the wall-adapting local eddy-viscosity (WALE) model within a large-eddy simulation (LES) framework, without model parameter recalibration, has limited predictive capabilities as signalized by poor representation of wall curvature effects and notable grid dependence. The identified lack in the modeling of radial transport processes is therefore addressed here by utilizing a stochastic one-dimensional turbulence (ODT) model. A standalone ODT formulation for cylindrical geometry is used in order to asses to which extent the predictability can be expected to improve by utilizing an advanced wall-modeling modeling strategy. It is shown that ODT is capable of capturing spanwise curvature and finite Reynolds number effects for fixed adjustable ODT model parameters. Based on the analogy of heat and mass transfer, present results yield new opportunities for modeling turbulent transfer process in chemical, process, and thermal engineering., Comment: In: New Results in Numerical and Experimental Fluid Mechanics XIV -- Contributions to the 23rd STAB/DGLR Symposium Berlin, Germany, 2022, edited by Andreas Dillmann, Gerd Heller, Ewald Kr\"amer, Claus Wagner, and Julien Weiss

The linear stability of a fully-developed liquid-metal MHD pipe flow subject to a transverse magnetic field is studied numerically. Because of the lack of axial symmetry in the mean velocity profile, we need to perform a BiGlobal stability analysis. For that purpose, we develop a two-dimensional complex eigenvalue solver relying on a Chebyshev-Fourier collocation method in physical space. By performing an extensive parametric study, we show that in contrast to Hagen-Poiseuille flow known to be linearly stable for all Reynolds numbers, the MHD pipe flow with transverse magnetic field is unstable to three-dimensional disturbances at sufficiently high values of the Hartmann number and wall conductance ratio. The instability observed in this regime is attributed to the presence of velocity overspeeds in the so-called Roberts layers and the corresponding inflection points in the mean velocity profile. The nature and characteristics of the most unstable modes are investigated, and we show that they vary significantly depending on the wall conductance ratio. A major result of this paper is that the global critical Reynolds number for the MHD pipe flow with transverse magnetic field is $Re=45230$, and it occurs for a perfectly conducting pipe wall and the Hartmann number $Ha=19.7$.

The world is full of fluids that flow. The fluid nature of flowing skyrmionic quasiparticles is of fundamental physical interest and plays an essential role in the transport of many skyrmions. Here, we report the laminar and transiently disordered dynamic behaviors of many magnetic skyrmions flowing in a pipe channel. The skyrmion flow driven by a uniform current may show a lattice structural transition. The skyrmion flow driven by a non-uniform current shows a dynamically varying lattice structure. A large uniform current could result in the compression of skyrmions toward the channel edge, leading to the transition of the skyrmion pipe flow into an open-channel flow with a free surface. Namely, the width of the skyrmion flow could be adjusted by the driving current. Skyrmions on the free surface may form a single shear layer adjacent to the main skyrmion flow. In addition, although we focus on the skyrmion flow dynamics in a clean pipe channel without any pinning or defect effect, we also show that a variation of magnetic anisotropy in the pipe channel could lead to more complicated skyrmion flow dynamics and pathlines. Our results reveal the fluid nature of skyrmionic quasiparticles that may motivate future research on the complex flow physics of magnetic textures., Comment: 16 pages, 17 figures

This study aims to extract and characterize structures in fully developed pipe flow at a friction Reynolds number of $\text{Re}_\tau = 12\,400$. To do so, we employ data-driven wavelet decomposition (DDWD) [D.~Floryan and M.~D.~Graham, PNAS 118, e2021299118 (2021)], a method that combines features of proper orthogonal decomposition and wavelet analysis in order to extract energetic and spatially localized structures from data. We apply DDWD to streamwise velocity signals measured separately via a thermal anemometer at 40 wall-normal positions. The resulting localized velocity structures, which we interpret as being reflective of underlying eddies, are self-similar across streamwise extents of 40 wall units to one pipe radius, and across wall-normal positions from $y^+=350$ to $y/R=1$. Notably, the structures are similar in shape to Meyer wavelets. Projections of the data onto the DDWD wavelet subspaces are found to be self-similar as well, but in Fourier space; the bounds of self-similarity are the same as before, except streamwise self-similarity starts at a larger length scale of $450$ wall units. The evidence of self-similarity provided in this study lends further support to Townsend's attached eddy hypothesis, although we note that the self-similar structures are detected beyond the log layer and extend to large length scales.

A resolvent-based methodology is employed to obtain spatio--temporal estimates of turbulent pipe flow from probe measurements of wall shear-stress fluctuations. Direct numerical simulations (DNS) and large-eddy simulations (LES) of turbulent pipe flow at friction Reynolds number of 550 are used as databases. We consider a DNS database as the true spatio--temporal flow field, from which wall shear-stress fluctuations are extracted and considered as measurements. A resolvent-based estimator is built following our earlier work (Amaral et al. J. Fluid Mech., vol. 927, 2021, p. A17), requiring a model for the nonlinear (or forcing) terms of the Navier-Stokes equations system, which are obtained from another DNS database, as in our earlier work, and from a series of computationally cheaper LES databases with coarser grids; the underlying idea is that LES may provide accurate statistics of non-linear terms related to large-scale structures at a low computational cost. Comparisons between the DNS and the estimates indicate that sufficiently accurate results can be achieved with estimators built with statistics from LES with an order of magnitude less grid points than the DNS, with estimates closely matching the reference DNS results up to the buffer layer and reasonable agreement up to the beginning of the log layer., Comment: 18 pages, 12 figures

Turbulent pipe flow is still an essentially open area of research, boosted in the last two decades by considerable progress achieved both on the experimental and numerical frontiers, mainly related to the identification and characterization of coherent structures as basic building blocks of turbulence. It has been a challenging task, however, to detect and visualize these coherent states. We address, by means of stereoscopic particle image velocimetry, that issue with the help of a large diameter (6 inches) pipe loop, which allowed us to probe for coherent states at various moderate Reynolds numbers (5300 < Re < 29000)). Although these states have been observed at flow regimes around laminar-turbulent transition (Re $\approx$ 2300) and also at high Reynolds number pipe flow (Re $\approx$ 35000), at moderate Reynolds numbers their existence had not been observed yet by experiment. By conditionally averaging the flow fields with respect to their dominant azimuthal wavenumber of streamwise velocity streaks, we have been able to uncover the existence of ten well-defined coherent flow patterns. It turns out, as a remarkable phenomenon, that their occurrence probabilities and the total number of dominant modes do not essentially change as the Reynolds number is varied. Their occurrence probabilities are noted to be reasonably well described by a Poisson distribution, which suggests that low-speed streaks are created as a Poisson process on the pipe circular geometry., Comment: 9 pages, 14 Figures

Turbulent pipe flows exhibit organizational states (OSs) that are labelled by discrete azimuthal wavenumber modes and are reminiscent of the traveling wave solutions of low Reynolds number regimes. The discretized time evolution of the OSs, obtained through stereoscopic particle image velocimetry, is shown to be non-Markovian for data acquisition carried out at a structure-resolved sampling rate. In particular, properly defined time-correlation functions for the OS transitions are observed to decay as intriguing power laws, up to a large-eddy time horizon, beyond which they decorrelate at much faster rates. We are able to establish, relying upon a probabilistic description of the creation and annihilation of streamwise streaks, a lower-level {\it{Markovian}} model for the OS transitions, which reproduces their time-correlated behavior with meaningful accuracy. These findings indicate that the OSs are distributed along the pipe as statistically correlated packets of quasi-streamwise vortical structures., Comment: 6 pages, 7 figures