Search

Although the roll-streak (R-S) is fundamentally involved in the dynamics of wall-turbulence, the physical mechanism responsible for its formation and maintenance remains controversial. In this work we investigate the dynamics maintaining the R-S in turbulent Poiseuille flow at R=1650. Spanwise collocation is used to remove spanwise displacement of the streaks and associated flow components, which isolates the streamwise-mean flow R-S component and the second-order statistics of the streamwise-varying fluctuations that are collocated with the R-S. This streamwise-mean/fluctuation partition of the dynamics facilitates exploiting insights gained from the analytic characterization of turbulence in the second-order statistical state dynamics (SSD), referred to as S3T, and its closely associated restricted nonlinear dynamics (RNL) approximation. Symmetry of the statistics about the streak centerline permits separation of the fluctuations into sinuous and varicose components. The Reynolds stress forcing induced by the sinuous and varicose fluctuations acting on the R-S is shown to reinforce low- and high-speed streaks respectively. This targeted reinforcement of streaks by the Reynolds stresses occurs continuously as the fluctuation field is strained by the streamwise-mean streak and not intermittently as would be associated with streak-breakdown events. The Reynolds stresses maintaining the streamwise-mean roll arise primarily from the dominant POD modes of the fluctuations, which can be identified with the time average structure of optimal perturbations growing on the streak. These results are consistent with a universal process of R-S growth and maintenance in turbulent shear flow arising from roll forcing generated by straining turbulent fluctuations, which was identified using the S3T SSD., Comment: 33 pages, 26 figures, submitted for publication to the Journal of Fluid Mechanics

Turbulence in wall-bounded flows is characterized by stable statistics. Although, in many turbulent systems, this stable statistical state corresponds to a stable fixed point of an associated statistical state dynamics (SSD) closed at second order, referred to as S3T, this is not the case for wall turbulence. In wall-turbulence the trajectory of the statistical state is on a transient chaotic attractor in the S3T statistical state space and the time-mean statistical state is neither a stable fixed point of this SSD nor, if it is maintained as an equilibrium, is it stable. Nevertheless, sufficiently small perturbations from the ensemble/time-mean state relax back to the mean statistical state following an effective linear dynamics. In this work, the dynamics of spanwise uniform perturbations to the time-mean flow are studied using a linear inverse model (LIM) to identify the linear operator governing the ensemble stability of the ensemble/time-mean state by obtaining the time-mean stability properties over the transient attractor of the turbulence identified by the S3T SSD. The ensemble/time-mean stability of an unstable equilibrium can be understood by noting that even when every member of an ensemble is unstable the ensemble mean may be stable with perturbations following stable dynamics. While simplifying insight into turbulent flows has been obtained by identifying and studying ensemble mean statistical states, less attention has been accorded to identifying and studying the ensemble mean dynamics. We show that in the case of wall turbulence, even though stable fixed point SSD equilibria are not available to allow the application of traditional perturbation analysis methods to identify the perturbation stability of the mean state, an effective linear stability analysis can be obtained to identify the perturbation dynamics of the ensemble/time-mean statistical state., Comment: 42 pages, 12 figures, submitted to Phys. Rev. Fluids

Turbulence in the restricted nonlinear (RNL) dynamics is analyzed and compared with DNS of Poiseuille turbulence at $R=1650$. The structures are obtained by POD analysis of the two components of the flow partition used in RNL dynamics: the streamwise-mean flow and fluctuations. POD analysis of the streamwise-mean flow indicates that the dominant POD modes, in both DNS and RNL, are roll-streaks harmonic in the spanwise. However, we conclude that these POD modes do not occur in isolation but rather are Fourier components of a coherent roll-streak structure. POD analysis of the fluctuations in DNS and RNL reveals similar complex structures consisting in part of oblique waves collocated with the streak. The origin of these structures is identified by their correspondence to POD modes predicted using a stochastic turbulence model (STM). These predicted POD modes are dominated by the optimally growing structures on the streak, which the STM predicts correctly to be of sinuous oblique wave structure. This close correspondence between the roll-streak structure and the associated fluctuations in DNS, RNL and the STM implies that the self-sustaining mechanism operating in DNS is essentially the same as that in RNL, which has been previously associated with optimal perturbation growth on the streak., Comment: 31 pages, accepted by the Journal of Fluid Mechanics

While linear non-normality underlies the mechanism of energy transfer from the externally driven flow to the perturbation field that sustains turbulence, nonlinearity is also known to play an essential role. The goal of this study is to better understand the role of nonlinearity in sustaining turbulence. The method used in this study is implementation in Couette flow of a statistical state dynamics (SSD) closure at second order in a cumulant expansion of the Navier-Stokes equations in which the averaging operator is the streamwise mean. The perturbations are the deviations from the streamwise mean and two mechanisms potentially contributing to maintaining these perturbations are identified. These are parametric perturbation growth arising from interaction of the perturbations with the fluctuating mean flow and transient growth of perturbations arising from nonlinear interaction between components of the perturbation field. By the method of comparing the turbulence maintained in the SSD and in the associated direct numerical simulation (DNS) in which these mechanisms have been selectively included and excluded, parametric growth is found to maintain the perturbation field of the turbulence while the more commonly invoked mechanism of transient growth of perturbations arising from scattering by nonlinear interaction is found to suppress perturbation growth. In addition to verifying that the parametric mechanism maintains the perturbations in DNS it is also verified that the Lyapunov vectors are the structures that dominate the perturbation energy and energetics in DNS. It is further verified that these vectors are responsible for maintaining the roll circulation that underlies the self-sustaining process (SSP) and in particular the maintenance of the fluctuating streak that supports the parametric perturbation growth., Comment: arXiv admin note: substantial text overlap with arXiv:1801.02093

The no-slip boundary condition results in a velocity shear forming in fluid flow near a solid surface. This shear flow supports the turbulence characteristic of fluid flow near boundaries at Reynolds numbers above $\approx1000$ by making available to perturbations the kinetic energy of the externally forced flow. Understanding the physical mechanism underlying this energy transfer poses a fundamental question. Although qualitative understanding that this transfer involves nonlinear destabilization of the roll-streak coherent structure has been established, identification of this instability has resisted analysis. The reason this instability has resisted analysis is that its analytic expression lies in the Navier-Stokes equations (NS) expressed using statistical rather than state variables. Expressing NS as a statistical state dynamics (SSD) at second order in a cumulant expansion suffices to allow analytical identification of the nonlinear roll-streak instability underlying turbulence in wall-bounded shear flow. In this nonlinear instability the turbulent perturbation field is identified by the SSD with the Lyapunov vectors of the linear operator governing perturbation evolution about the time-dependent streamwise mean flow. In this work, the implications of the predictions of SSD analysis (that this parametric instability underlies the dynamics of turbulence in Couette flow and that the perturbation structures are the associated Lyapunov vectors) are interpreted to imply new conceptual approaches to controlling turbulence. It is shown that the perturbation component of turbulence is supported on the streamwise mean flow, which implies optimal control should be formulated to suppress perturbations from the streamwise mean. It is also shown that suppressing only the top few Lyapunov vectors on the streamwise mean vectors results in laminarization. These results are verified using DNS.

Vertically banded zonal jets are frequently observed in weakly or non-rotating stratified turbulence, with the quasi-biennial oscillation in the equatorial stratosphere and the ocean's equatorial deep jets being two examples. Explaining the formation of jets in stratified turbulence is a fundamental problem in geophysical fluid dynamics. Statistical state dynamics (SSD) provides powerful methods for analyzing turbulent systems exhibiting emergent organization, such as banded jets. In SSD, dynamical equations are written directly for the evolution of the turbulence statistics, enabling direct analysis of the statistical interactions between the incoherent component of the turbulence and the coherent large-scale structure component that underlie jet formation. A second-order closure of SSD, known as S3T, has previously been applied to show that meridionally banded jets emerge in barotropic beta-plane turbulence via a statistical instability referred to as the zonostrophic instability. Two-dimensional Boussinesq turbulence provides a simple model of non-rotating stratified turbulence analogous to the beta-plane model of planetary turbulence. Jets known as vertically sheared horizontal flows (VSHFs) often emerge in simulations of Boussinesq turbulence, but their dynamics is not yet clearly understood. In this work S3T analysis of the zonostrophic instability is extended to study VSHF emergence in two-dimensional Boussinesq turbulence using an analytical formulation of S3T amenable to perturbation stability analysis. VSHFs are shown to form via an instability that is analogous in stratified turbulence to the zonostrophic instability in beta-plane turbulence. This instability is shown to be strikingly similar to the zonostrophic instability, suggesting that jet emergence in both geostrophic and non-rotating stratified turbulence may be understood as instances of the same generic phenomenon., Comment: Re-submitted to Journal of the Atmospheric Sciences

Turbulence in wall-bounded shear flow results from a synergistic interaction between linear non-normality and nonlinearity in which non-normal growth of a subset of perturbations configured to transfer energy from the externally forced component of the turbulent state to the perturbation component maintains the perturbation energy, while the subset of energy-transferring perturbations is replenished by nonlinearity. Although it is accepted that both linear non-normality mediated energy transfer from the forced component of the mean flow and nonlinear interactions among perturbations are required to maintain the turbulent state, the detailed physical mechanism by which these processes interact in maintaining turbulence has not been determined. In this work a statistical state dynamics based analysis is performed on turbulent Couette flow at $R=600$ and a comparison to DNS is used to demonstrate that the perturbation component in Couette flow turbulence is replenished by a non-normality mediated parametric growth process in which the fluctuating streamwise mean flow has been adjusted to marginal Lyapunov stability. It is further shown that the alternative mechanism in which the subspace of non-normally growing perturbations is maintained directly by perturbation-perturbation nonlinearity does not contribute to maintaining the turbulent state. This work identifies parametric interaction between the fluctuating streamwise mean flow and the streamwise varying perturbations to be the mechanism of the nonlinear interaction maintaining the perturbation component of the turbulent state, and identifies the associated Lyapunov vectors with positive energetics as the structures of the perturbation subspace supporting the turbulence., Comment: work done in the COTURB Summer Workshop 2017, submitted to J. Phys.: Conf. Series

Simulations of strongly stratified turbulence often exhibit coherent large-scale structures called vertically sheared horizontal flows (VSHFs). VSHFs emerge in both two-dimensional (2D) and three-dimensional (3D) stratified turbulence with similar vertical structure. The mechanism responsible for VSHF formation is not fully understood. In this work, the formation and equilibration of VSHFs in a 2D Boussinesq model of stratified turbulence is studied using statistical state dynamics (SSD). In SSD, equations of motion are expressed directly in the statistical variables of the turbulent state. Restriction to 2D turbulence makes available an analytically and computationally attractive implementation of SSD referred to as S3T, in which the SSD is expressed by coupling the equation for the horizontal mean structure with the equation for the ensemble mean perturbation covariance. This second order SSD produces accurate statistics, through second order, when compared with fully nonlinear simulations. In particular, S3T captures the spontaneous emergence of the VSHF and associated density layers seen in simulations of turbulence maintained by homogeneous large-scale stochastic excitation. An advantage of the S3T system is that the VSHF formation mechanism, which is wave-mean flow interaction between the emergent VSHF and the stochastically excited large-scale gravity waves, is analytically understood in the S3T system. Comparison with fully nonlinear simulations verifies that S3T solutions accurately predict the scale selection, dependence on stochastic excitation strength, and nonlinear equilibrium structure of the VSHF. These results facilitate relating VSHF theory and geophysical examples of turbulent jets such as the ocean's equatorial deep jets., Comment: Revised version under review for publication in Journal of Fluid Mechanics

This paper reviews results from the study of wall-bounded turbulent flows using statistical state dynamics (SSD) that demonstrate the benefits of adopting this perspective for understanding turbulence in wall-bounded shear flows. The SSD approach used in this work employs a second-order closure which isolates the interaction between the streamwise mean and the equivalent of the perturbation covariance. This closure restricts nonlinearity in the SSD to that explicitly retained in the streamwise constant mean together with nonlinear interactions between the mean and the perturbation covariance. This dynamical restriction, in which explicit perturbation-perturbation nonlinearity is removed from the perturbation equation, results in a simplified dynamics referred to as the restricted nonlinear (RNL) dynamics. RNL systems in which an ensemble of a finite number of realizations of the perturbation equation share the same mean flow provide tractable approximations to the equivalently infinite ensemble RNL system. The infinite ensemble system, referred to as the S3T, introduces new analysis tools for studying turbulence. The RNL with a single ensemble member can be alternatively viewed as a realization of RNL dynamics. RNL systems provide computationally efficient means to approximate the SSD, producing self-sustaining turbulence exhibiting qualitative features similar to those observed in direct numerical simulations (DNS) despite its greatly simplified dynamics. Finally, we show that RNL turbulence can be supported by as few as a single streamwise varying component interacting with the streamwise constant mean flow and that judicious selection of this truncated support, or "band-limiting", can be used to improve quantitative accuracy of RNL turbulence. The results suggest that the SSD approach provides new analytical and computational tools allowing new insights into wall-turbulence., Comment: to appear in Phil. Trans. R. Soc. A

Coherent jets with most of the kinetic energy of the flow are common in atmospheric turbulence. In the gaseous planets these jets are maintained by incoherent turbulence excited by small-scale convection. Large-scale coherent waves are sometimes observed to coexist with the jets; a prominent example is Saturn's hexagonal North polar jet (NPJ). The mechanism responsible for forming and maintaining such a turbulent state remains elusive. The coherent planetary-scale component of the turbulence arises and is maintained by interaction with the incoherent small-scale turbulence component. Theoretical understanding of the dynamics of the jet/wave/turbulence coexistence regime is gained by employing a statistical state dynamics (SSD) model. Here, a second-order closure implementation of a two-layer beta-plane SSD is used to develop a theory that accounts for the structure and dynamics of the NPJ. Asymptotic analysis of the SSD equilibrium in the weak jet damping limit predicts a universal jet structure in agreement with NPJ observations. This asymptotic theory also predicts the wavenumber (six) of the prominent jet perturbation. Analysis with this model of the jet/wave/turbulence regime dynamics reveals that jet formation is controlled by the effective value of $\beta$; the required value of this parameter for correspondence with observation is obtained. As this is a robust prediction it is taken as an indirect observation of a deep poleward sloping stable layer beneath the NPJ. The slope required is obtained from observations of NPJ structure as is the small-scale turbulence excitation required to maintain the jet. The observed jet structure is then predicted by the theory as is the wave-six disturbance. This wave, which is identified with the least stable mode of the equilibrated jet, is shown to be primarily responsible for equilibrating the jet with the observed structure and amplitude., Comment: (minor title change)