Recurrent flow patterns as a basis for two-dimensional turbulence: Predicting statistics from structures.

A dynamical systems approach to turbulence envisions the flow as a trajectory through a high-dimensional state space [Hopf, Commun. Appl. Maths 1, 303 (1948)]. The chaotic dynamics are shaped by the unstable simple invariant solutions populating the inertial manifold. The hope has been to turn this picture into a predictive framework where the statistics of the flow follow from a weighted sum of the statistics of each simple invariant solution. Two outstanding obstacles have prevented this goal from being achieved: 1) paucity of known solutions and 2) the lack of a rational theory for predicting the required weights. Here, we describe a method to substantially solve these problems, and thereby provide compelling evidence that the probability density functions (PDFs) of a fully developed turbulent flow can be reconstructed with a set of unstable periodic orbits. Our method for finding solutions uses automatic differentiation, with high-quality guesses constructed by minimizing a trajectory-dependent loss function. We use this approach to find hundreds of solutions in turbulent, two-dimensional Kolmogorov flow. Robust statistical predictions are then computed by learning weights after converting a turbulent trajectory into a Markov chain for which the states are individual solutions, and the nearest solution to a given snapshot is determined using a deep convolutional autoencoder. In this study, the PDFs of a spatiotemporally chaotic system have been successfully reproduced with a set of simple invariant states, and we provide a fascinating connection between self-sustaining dynamical processes and the more well-known statistical properties of turbulence.

Wall shear stress and pressure patterns in aortic stenosis patients with and without aortic dilation captured by high-performance image-based computational fluid dynamics.

Spatial patterns of elevated wall shear stress and pressure due to blood flow past aortic stenosis (AS) are studied using GPU-accelerated patient-specific computational fluid dynamics. Three cases of moderate to severe AS, one with a dilated ascending aorta and two within the normal range (root diameter less than 4cm) are simulated for physiological waveforms obtained from echocardiography. The computational framework is built based on sharp-interface Immersed Boundary Method, where aortic geometries segmented from CT angiograms are integrated into a high-order incompressible Navier-Stokes solver. The key question addressed here is, given the presence of turbulence due to AS which increases wall shear stress (WSS) levels, why some AS patients undergo much less aortic dilation. Recent case studies of AS have linked the existence of an elevated WSS hotspot (due to impingement of AS on the aortic wall) to the dilation process. Herein we further investigate the WSS distribution for cases with and without dilation to understand the possible hemodynamics which may impact the dilation process. We show that the spatial distribution of elevated WSS is significantly more focused for the case with dilation than those without dilation. We further show that this focal area accommodates a persistent pocket of high pressure, which may have contributed to the dilation process through an increased wall-normal forcing. The cases without dilation, on the contrary, showed a rather oscillatory pressure behaviour, with no persistent pressure "buildup" effect. We further argue that a more proximal branching of the aortic arch could explain the lack of a focal area of elevated WSS and pressure, because it interferes with the impingement process due to fluid suction effects. These phenomena are further illustrated using an idealized aortic geometry. We finally show that a restored inflow eliminates the focal area of elevated WSS and pressure zone from the ascending aorta.

Heat transport in Rayleigh-Bénard convection with linear marginality.

Recent direct numerical simulations (DNS) and computations of exact steady solutions suggest that the heat transport in Rayleigh-Bénard convection (RBC) exhibits the classical [Formula: see text] scaling as the Rayleigh number [Formula: see text] with Prandtl number unity, consistent with Malkus-Howard's marginally stable boundary layer theory. Here, we construct conditional upper and lower bounds for heat transport in two-dimensional RBC subject to a physically motivated marginal linear-stability constraint. The upper estimate is derived using the Constantin-Doering-Hopf (CDH) variational framework for RBC with stress-free boundary conditions, while the lower estimate is developed for both stress-free and no-slip boundary conditions. The resulting optimization problems are solved numerically using a time-stepping algorithm. Our results indicate that the upper heat-flux estimate follows the same [Formula: see text] scaling as the rigorous CDH upper bound for the two-dimensional stress-free case, indicating that the linear-stability constraint fails to modify the boundary-layer thickness of the mean temperature profile. By contrast, the lower estimate successfully captures the [Formula: see text] scaling for both the stress-free and no-slip cases. These estimates are tested using marginally-stable equilibrium solutions obtained under the quasi-linear approximation, steady roll solutions and DNS data. This article is part of the theme issue 'Mathematical problems in physical fluid dynamics (part 1)'.

Exact Traveling Wave Solutions in Viscoelastic Channel Flow.

Elasto-inertial turbulence (EIT) is a new, two-dimensional chaotic flow state observed in polymer solutions with possible connections to inertialess elastic turbulence and drag-reduced Newtonian turbulence. In this Letter, we argue that the origins of EIT are fundamentally different from Newtonian turbulence by finding a dynamical connection between EIT and an elasto-inertial linear instability recently found at high Weissenberg numbers [Garg et al., Phys. Rev. Lett. 121, 024502 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.024502]. This link is established by isolating the first known exact coherent structures in viscoelastic parallel flows-nonlinear elasto-inertial traveling waves (TWs)-borne at the linear instability and tracking them down to substantially lower Weissenberg numbers where EIT exists. These TWs have a distinctive "arrowhead" structure in the polymer stretch field and can be clearly recognized albeit transiently in EIT as well as being attractors for EIT dynamics if the Weissenberg number is sufficiently large. Our findings suggest that the dynamical systems picture in which Newtonian turbulence is built around the coexistence of many (unstable) simple invariant solutions populating phase space carries over to EIT, though these solutions rely on elasticity to exist.

Using nonlinear transient growth to construct the minimal seed for shear flow turbulence.

Linear transient growth analysis is commonly used to suggest the structure of disturbances which are particularly efficient in triggering transition to turbulence in shear flows. We demonstrate that the addition of nonlinearity to the analysis can substantially change the prediction made in pipe flow from simple two-dimensional streamwise rolls to a spanwise and cross-stream localized three-dimensional state. This new nonlinear optimal is demonstrably more efficient in triggering turbulence than the linear optimal indicating that there are better ways to design perturbations to achieve transition.

Connection between nonlinear energy optimization and instantons.

How systems transit between different stable states under external perturbation is an important practical issue. We discuss here how a recently developed energy optimization method for identifying the minimal disturbance necessary to reach the basin boundary of a stable state is connected to the instanton trajectory from large deviation theory of noisy systems. In the context of the one-dimensional Swift-Hohenberg equation, which has multiple stable equilibria, we first show how the energy optimization method can be straightforwardly used to identify minimal disturbances-minimal seeds-for transition to specific attractors from the ground state. Then, after generalizing the technique to consider multiple, equally spaced-in-time perturbations, it is shown that the instanton trajectory is indeed the solution of the energy optimization method in the limit of infinitely many perturbations provided a specific norm is used to measure the set of discrete perturbations. Importantly, we find that the key features of the instanton can be captured by a low number of discrete perturbations (typically one perturbation per basin of attraction crossed). This suggests a promising new diagnostic for systems for which it may be impractical to calculate the instanton.

Connections between centrifugal, stratorotational, and radiative instabilities in viscous Taylor-Couette flow.

The "Rayleigh line" µ=η^{2}, where µ=Ω_{o}/Ω_{i} and η=r_{i}/r_{o} are respectively the rotation and radius ratios between inner (subscript i) and outer (subscript o) cylinders, is regarded as marking the limit of centrifugal instability (CI) in unstratified inviscid Taylor-Couette flow, for both axisymmetric and nonaxisymmetric modes. Nonaxisymmetric stratorotational instability (SRI) is known to set in for anticyclonic rotation ratios beyond that line, i.e., η^{2}<µ<1 for axially stably stratified Taylor-Couette flow, but the competition between CI and SRI in the range µ<η^{2} has not yet been addressed. In this paper, we establish continuous connections between the two instabilities at finite Reynolds number Re, as previously suggested by Le Bars and Le Gal [Phys. Rev. Lett. 99, 064502 (2007)PRLTAO0031-900710.1103/PhysRevLett.99.064502], making them indistinguishable at onset. Both instabilities are also continuously connected to the radiative instability at finite Re. These results demonstrate the complex impact viscosity has on the linear stability properties of this flow. Several other qualitative differences with inviscid theory were found, among which are the instability of a nonaxisymmetric mode localized at the outer cylinder without stratification and the instability of a mode propagating against the inner cylinder rotation with stratification. The combination of viscosity and stratification can also lead to a "collision" between (axisymmetric) Taylor vortex branches, causing the axisymmetric oscillatory state already observed in past experiments. Perhaps more surprising is the instability of a centrifugal-like helical mode beyond the Rayleigh line, caused by the joint effects of stratification and viscosity. The threshold µ=η^{2} seems to remain, however, an impassable instability limit for axisymmetric modes, regardless of stratification, viscosity, and even disturbance amplitude.

Time-stepping approach for solving upper-bound problems: Application to two-dimensional Rayleigh-Bénard convection.

An alternative computational procedure for numerically solving a class of variational problems arising from rigorous upper-bound analysis of forced-dissipative infinite-dimensional nonlinear dynamical systems, including the Navier-Stokes and Oberbeck-Boussinesq equations, is analyzed and applied to Rayleigh-Bénard convection. A proof that the only steady state to which this numerical algorithm can converge is the required global optimal of the relevant variational problem is given for three canonical flow configurations. In contrast with most other numerical schemes for computing the optimal bounds on transported quantities (e.g., heat or momentum) within the "background field" variational framework, which employ variants of Newton's method and hence require very accurate initial iterates, the new computational method is easy to implement and, crucially, does not require numerical continuation. The algorithm is used to determine the optimal background-method bound on the heat transport enhancement factor, i.e., the Nusselt number (Nu), as a function of the Rayleigh number (Ra), Prandtl number (Pr), and domain aspect ratio L in two-dimensional Rayleigh-Bénard convection between stress-free isothermal boundaries (Rayleigh's original 1916 model of convection). The result of the computation is significant because analyses, laboratory experiments, and numerical simulations have suggested a range of exponents α and ß in the presumed Nu∼Pr(α)Ra(ß) scaling relation. The computations clearly show that for Ra≤10(10) at fixed L=2√[2],Nu≤0.106Pr(0)Ra(5/12), which indicates that molecular transport cannot generally be neglected in the "ultimate" high-Ra regime.

Axisymmetric granular collapse: a transient 3D flow test of viscoplasticity.

A viscoplastic continuum theory has recently been proposed to model dense, cohesionless granular flows [P. Jop, Nature (London) 441, 727 (2006)10.1038/nature04801]. We confront this theory for the first time with a transient, three-dimensional flow situation--the simple collapse of a cylinder of granular matter onto a horizontal plane--by extracting stress and strain rate tensors directly from soft particle simulations. These simulations faithfully reproduce the different flow regimes and capture the observed scaling laws for the final deposit. Remarkably, the theoretical hypothesis that there is a simple stress-strain rate tensorial relationship does seem to hold across the whole flow even close to the rough boundary provided the flow is dense enough. These encouraging results suggest viscoplastic theory is more generally applicable to transient, multidirectional, dense flows and open the way for quantitative predictions in real applications.

Highly symmetric travelling waves in pipe flow.

The recent theoretical discovery of finite-amplitude travelling waves (TWs) in pipe flow has reignited interest in the transitional phenomena that Osborne Reynolds studied 125 years ago. Despite all being unstable, these waves are providing fresh insight into the flow dynamics. We describe two new classes of TWs, which, while possessing more restrictive symmetries than previously found TWs of Faisst & Eckhardt (2003 Phys. Rev. Lett. 91, 224502) and Wedin & Kerswell (2004 J. Fluid Mech. 508, 333-371), seem to be more fundamental to the hierarchy of exact solutions. They exhibit much higher wall shear stresses and appear at notably lower Reynolds numbers. The first M-class comprises the various discrete rotationally symmetric analogues of the mirror-symmetric wave found in Pringle & Kerswell (2007 Phys. Rev. Lett. 99, 074502), and have a distinctive double-layered structure of fast and slow streaks across the pipe radius. The second N-class has the more familiar separation of fast streaks to the exterior and slow streaks to the interior and looks like the precursor to the class of non-mirror-symmetric waves already known.

Coherent structures in localized and global pipe turbulence.

The recent discovery of unstable traveling waves (TWs) in pipe flow has been hailed as a significant breakthrough with the hope that they populate the turbulent attractor. We confirm the existence of coherent states with internal fast and slow streaks commensurate in both structure and energy with known TWs using numerical simulations in a long pipe. These only occur, however, within less energetic regions of (localized) "puff" turbulence at low Reynolds numbers (Re=2000-2400), and not at all in (homogeneous) "slug" turbulence at Re=2800. This strongly suggests that all currently known TWs sit in an intermediate region of phase space between the laminar and turbulent states rather than being embedded within the turbulent attractor itself. New coherent fast streak states with strongly decelerated cores appear to populate the turbulent attractor instead.

Critical behavior in the relaminarization of localized turbulence in pipe flow.

The statistics of the relaminarization of localized turbulence in a pipe are examined by direct numerical simulation. As in recent experimental data [J. Peixinho and T. Mullin, Phys. Rev. Lett. 96, 094501 (2006)10.1103/PhysRevLett.96.094501], the half-life for the decaying turbulence is consistent with the scaling (Rec-Re) -1, indicating a boundary crisis of the localized turbulent state familiar in low-dimensional dynamical systems. The crisis Reynolds number is estimated as Rec=1870, a value within 7% of the experimental value 1750. We argue that the frequently asked question, of which initial disturbances at a given Re trigger sustained turbulence in a pipe, is really two separate questions: the "local phase space" question (local to the laminar state) of what threshold disturbance at a given Re is needed to initially trigger turbulence, followed by the "global phase space" question of whether Re exceeds Rec at which point the turbulent state becomes an attractor.

Asymmetric, helical, and mirror-symmetric traveling waves in pipe flow.

New families of three-dimensional nonlinear traveling waves are discovered in pipe flow. In contrast with known waves [H. Faisst and B. Eckhardt, Phys. Rev. Lett. 91, 224502 (2003); H. Wedin and R. R. Kerswell, J. Fluid Mech. 508, 333 (2004), they possess no discrete rotational symmetry and exist at a significantly lower Reynolds numbers (Re). First to appear is a mirror-symmetric traveling wave which is born in a saddle node bifurcation at Re=773. As Re increases, "asymmetric" modes arise through a symmetry-breaking bifurcation. These look to be a minimal coherent unit consisting of one slow streak sandwiched between two fast streaks located preferentially to one side of the pipe. Helical and nonhelical rotating waves are also found, emphasizing the richness of phase space even at these very low Reynolds numbers.