Search

One of the outstanding problems in the dynamical systems approach to turbulence is to find a sufficient number of invariant solutions to characterize the underlying dynamics of turbulence [Annu. Rev. Fluid Mech. 44, 203 (2012)10.1146/annurev-fluid-120710-101228]. As a practical matter, the solutions can be difficult to find. To improve this situation, we show how to find periodic orbits and equilibria in plane Couette flow by projecting pseudorecurrent segments of turbulent trajectories onto the left-singular vectors of the Navier-Stokes equations linearized about the relevant mean flow (resolvent modes). The projections are, subsequently, used to initiate Newton-Krylov-hookstep searches, and new (relative) periodic orbits and equilibria are discovered. We call the process project-then-search and validate the process by first applying it to previously known fixed point and periodic solutions. Along the way, we find new branches of equilibria, which include bifurcations from previously known branches, and new periodic orbits that closely shadow turbulent trajectories in state space.

The effect of streamwise-varying steady transpiration on turbulent pipe flow is examined using direct numerical simulation at fixed friction Reynolds number $\mathit{Re}_{{\it\tau}}=314$. The streamwise momentum equation reveals three physical mechanisms caused by transpiration acting in the flow: modification of Reynolds shear stress, steady streaming and generation of non-zero mean streamwise gradients. The influence of these mechanisms has been examined by means of a parameter sweep involving transpiration amplitude and wavelength. The observed trends have permitted identification of wall transpiration configurations able to reduce or increase the overall flow rate $-36.1\,\%$ and $19.3\,\%$, respectively. Energetics associated with these modifications are presented. A novel resolvent formulation has been developed to investigate the dynamics of pipe flows with a constant cross-section but with time-mean spatial periodicity induced by changes in boundary conditions. This formulation, based on a triple decomposition, paves the way for understanding turbulence in such flows using only the mean velocity profile. Resolvent analysis based on the time-mean flow and dynamic mode decomposition based on simulation data snapshots have both been used to obtain a description of the reorganization of the flow structures caused by the transpiration. We show that the pipe flows dynamics are dominated by a critical-layer mechanism and the waviness induced in the flow structures plays a role on the streamwise momentum balance by generating additional terms.

We generate predictions for the fluctuating pressure field in turbulent pipe flow by reformulating the resolvent analysis of McKeon and Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382) in terms of the so-called primitive variables. Under this analysis, the nonlinear convective terms in the Fourier-transformed Navier–Stokes equations (NSE) are treated as a forcing that is mapped to a velocity and pressure response by the resolvent of the linearized Navier–Stokes operator. At each wavenumber–frequency combination, the turbulent velocity and pressure field are represented by the most-amplified (rank-1) response modes, identified via a singular value decomposition of the resolvent. We show that these rank-1 response modes reconcile many of the key relationships among the velocity field, coherent structure (i.e. hairpin vortices), and the high-amplitude wall-pressure events observed in previous experiments and direct numerical simulations (DNS). A Green’s function representation shows that the pressure fields obtained under this analysis correspond primarily to the fast pressure contribution arising from the linear interaction between the mean shear and the turbulent wall-normal velocity. Recovering the slow pressure requires an explicit treatment of the nonlinear interactions between the Fourier response modes. By considering the velocity and pressure fields associated with the triadically consistent mode combination studied by Sharma and McKeon (J. Fluid Mech., vol. 728, 2013, pp. 196–238), we identify the possibility of an apparent amplitude modulation effect in the pressure field, similar to that observed for the streamwise velocity field. However, unlike the streamwise velocity, for which the large scales of the flow are in phase with the envelope of the small-scale activity close to the wall, we expect there to be a $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\pi /2$ phase difference between the large-scale wall-pressure and the envelope of the small-scale activity. Finally, we generate spectral predictions based on a rank-1 model assuming broadband forcing across all wavenumber–frequency combinations. Despite the significant simplifying assumptions, this approach reproduces trends observed in previous DNS for the wavenumber spectra of velocity and pressure, and for the scale-dependence of wall-pressure propagation speed.

A new theory of coherent structure in wall turbulence is presented. The theory is the first to predict packets of hairpin vortices and other structure in turbulence, and their dynamics, based on an analysis of the Navier-Stokes equations, under an assumption of a turbulent mean profile. The assumption of the turbulent mean acts as a restriction on the class of possible structures. It is shown that the coherent structure is a manifestation of essentially low-dimensional flow dynamics, arising from a critical layer mechanism. Using the decomposition presented in McKeon & Sharma (J. Fluid Mech, 658, 2010), complex coherent structure is recreated from minimal superpositions of response modes predicted by the analysis, which take the form of radially-varying travelling waves. By way of example, simple combinations of these modes are offered that predicts hairpins and modulated hairpin packets. The phase interaction also predicts important skewness and correlation results known in the literature. It is also shown that the very large scale motions act to organise hairpin-like structures such that they co-locate with areas of low streamwise momentum, by a mechanism of locally varying the shear profile. The relationship between Taylor's hypothesis and coherence is discussed and both are shown to be the consequence of the localisation of the response modes around the critical layer. A pleasing link is made to the classical laminar inviscid theory, whereby the essential mechanism underlying the hairpin vortex is captured by two obliquely interacting Kelvin-Stuart (cat's eye) vortices. Evidence for the theory is presented based on comparison to observations of structure reported in the experimental, transitional flow and turbulent flow numerical simulation literature., Accepted for publication in the Journal of Fluid Mechanics. 41 pages

Previous work has established the usefulness of the resolvent operator that maps the terms nonlinear in the turbulent fluctuations to the fluctuations themselves. Further work has described the self-similarity of the resolvent arising from that of the mean velocity profile. The orthogonal modes provided by the resolvent analysis describe the wall-normal coherence of the motions and inherit that self-similarity. In this contribution, we present the implications of this similarity for the nonlinear interaction between modes with different scales and wall-normal locations. By considering the nonlinear interactions between modes, it is shown that much of the turbulence scaling behaviour in the logarithmic region can be determined from a single arbitrarily chosen reference plane. Thus, the geometric scaling of the modes is impressed upon the nonlinear interaction between modes. Implications of these observations on the self-sustaining mechanisms of wall turbulence, modelling and simulation are outlined., Comment: Accepted for publication in Philosophical Transactions of the Royal Society A. arXiv admin note: substantial text overlap with arXiv:1409.6047

A novel method to estimate unsteady aerodynamic force coefficients from pointwise velocity measurements is presented. The methodology is based on a resolvent-based reduced-order model which requires the mean flow to obtain physical flow structures and pointwise measurement to calibrate their amplitudes. A computationally-affordable time-stepping methodology to obtain resolvent modes in non-trivial flow domains is introduced and compared to previous existing matrix-free and matrix-forming strategies. The technique is applied to the unsteady flow around an inclined square cylinder at low Reynolds number. The potential of the methodology is demonstrated through good agreement between the fluctuating pressure distribution on the cylinder and the temporal evolution of the unsteady lift and drag coefficients predicted by the model and those computed by direct numerical simulation., In review

The relationship between Koopman mode decomposition, resolvent mode decomposition and exact invariant solutions of the Navier-Stokes equations is clarified. The correspondence rests upon the invariance of the system operators under symmetry operations such as spatial translation. The usual interpretation of the Koopman operator is generalised to permit combinations of such operations, in addition to translation in time. This invariance is related to the spectrum of a spatio-temporal Koopman operator, which has a travelling wave interpretation. The relationship leads to a generalisation of dynamic mode decomposition, in which symmetry operations are applied to restrict the dynamic modes to span a subspace subject to those symmetries. The resolvent is interpreted as the mapping between the Koopman modes of the Reynolds stress divergence and the velocity field. It is shown that the singular vectors of the resolvent (the resolvent modes) are the optimal basis in which to express the velocity field Koopman modes where the latter are not a priori known.

A model-based description of the scaling and radial location of turbulent fluctuations in turbulent pipe flow is presented and used to illuminate the scaling behaviour of the very large scale motions. The model is derived by treating the nonlinearity in the perturbation equation (involving the Reynolds stress) as an unknown forcing, yielding a linear relationship between the velocity field response and this nonlinearity. We do not assume small perturbations. We examine propagating helical velocity response modes that are harmonic in the wall-parallel directions and in time, permitting comparison of our results to experimental data. The steady component of the velocity field that varies only in the wall-normal direction is identified as the turbulent mean profile. A singular value decomposition of the resolvent identifies the forcing shape that will lead to the largest velocity response at a given wavenumber–frequency combination. The hypothesis that these forcing shapes lead to response modes that will be dominant in turbulent pipe flow is tested by using physical arguments to constrain the range of wavenumbers and frequencies to those actually observed in experiments. An investigation of the most amplified velocity response at a given wavenumber–frequency combination reveals critical-layer-like behaviour reminiscent of the neutrally stable solutions of the Orr–Sommerfeld equation in linearly unstable flow. Two distinct regions in the flow where the influence of viscosity becomes important can be identified, namely wall layers that scale withR+1/2and critical layers where the propagation velocity is equal to the local mean velocity, one of which scales withR+2/3in pipe flow. This framework appears to be consistent with several scaling results in wall turbulence and reveals a mechanism by which the effects of viscosity can extend well beyond the immediate vicinity of the wall. The model reproduces inner scaling of the small scales near the wall and an approach to outer scaling in the flow interior. We use our analysis to make a first prediction that the appropriate scaling velocity for the very large scale motions is the centreline velocity, and show that this is in agreement with experimental results. Lastly, we interpret the wall modes as the motion required to meet the wall boundary condition, identifying the interaction between the critical and wall modes as a potential origin for an interaction between the large and small scales that has been observed in recent literature as an amplitude modulation of the near-wall turbulence by the very large scales.

The objective of this study was to assess the effect of monosodium glutamate (MSG) on broiler performance, blood chemistry and the brain tissues. A total of 120 day-old Lohman broiler chicks were randomly divided into four groups, with 3 replicates (30 birds per treatment). Four approximately isonitrogenous and isocaloric diets with different levels of MSG viz. 0, 0.25, 0.5 and 1.0% were formulated. Feed and water were provided ad libitum throughout the experimental period of 49 days. The results revealed that total feed intake was significantly higher (P

We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and channel flows are well represented by just few modes of the model of McKeon & Sharma J. Fl. Mech. 658, 356 (2010). This model provides modes that act as a basis to decompose the velocity field, ordered by their amplitude of response to forcing arising from the interaction between scales. The model was originally derived from the Navier-Stokes equations to represent turbulent flows and has been used to explain coherent structure and to predict turbulent statistics. This establishes a surprising new link between the two distinct approaches to understanding turbulence., 21 figures

This paper employs the theoretical framework developed by Luhar et al. (J. Fluid Mech., 768, 415-441) to consider the design of compliant walls for turbulent skin friction reduction. Specifically, the effects of simple spring-damper walls are contrasted with the effects of more complex walls incorporating tension, stiffness and anisotropy. In addition, varying mass ratios are tested to provide insight into differences between aerodynamic and hydrodynamic applications. Despite the differing physical responses, all the walls tested exhibit some important common features. First, the effect of the walls (positive or negative) is greatest at conditions close to resonance, with sharp transitions in performance across the resonant frequency or phase speed. Second, compliant walls are predicted to have a more pronounced effect on slower-moving structures because such structures generally have larger wall-pressure signatures. Third, two-dimensional (spanwise constant) structures are particularly susceptible to further amplification. These features are consistent with many previous experiments and simulations, suggesting that mitigating the rise of such two-dimensional structures is essential to designing performance-improving walls. For instance, it is shown that further amplification of such large-scale two-dimensional structures explains why the optimal anisotropic walls identified by Fukagata et al. via DNS (J. Turb., 9, 1-17) only led to drag reduction in very small domains. The above observations are used to develop design and methodology guidelines for future research on compliant walls.

A novel reduced-order model for time-varying nonlinear flows arising from a resolvent decomposition based on the time-mean flow is proposed. The inputs required for the model are the mean-flow field and a small set of velocity time-series data obtained at isolated measurement points, which are used to fix relevant frequencies, amplitudes and phases of a limited number of resolvent modes that, together with the mean flow, constitute the reduced-order model. The technique is applied to derive a model for the unsteady three-dimensional flow in a lid-driven cavity at a Reynolds number of 1200 that is based on the two-dimensional mean flow, three resolvent modes selected at the most active spanwise wavenumber, and either one or two velocity probe signals. The least-squares full-field error of the reconstructed velocity obtained using the model and two point velocity probes is of the order of 5 % of the lid velocity, and the dynamical behaviour of the reconstructed flow is qualitatively similar to that of the complete flow.

This paper considers the convergence upon grid refinement of control strategies derived from ODE aproxima- tions of diffusive boundary-controlled linear PDE systems. It focuses specifically on the Dirichlet boundary control of the heat equation as a canonical model for more general diffusive PDE systems. It treats two classes of problems: the controllability problem (that is, the determination of a control distribution to move a system exactly from a specified initial state to a specified terminal state in finite time) and the state feedback control problem (that is, the determination of an optimal feedback rule u = Kx which minimizes some quadratic cost function J measuring both the state of the system and the control input), in the latter problem focusing specifically, for simplicity, on the infinite-horizon (that is, constant-gain) case. Both classes of problems require special attention beyond the usual considerations commonly known for the control of low-order ODE systems. Specifically, convergence of the control strategies upon refinement of the ODE approximation used in the controller calculation is not guaranteed. Note that the present study considers sine, finite difference, and Chebyshev discretizations, all of which provide consistent results, indicating that the results reported are not a spurious artificat of any particular numerical discretization.

This paper extends the resolvent formulation proposed by McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382) to consider turbulence–compliant wall interactions. Under this formulation, the turbulent velocity field is expressed as a linear superposition of propagating modes, identified via a gain-based decomposition of the Navier–Stokes equations. Compliant surfaces, modelled as a complex wall admittance linking pressure and velocity, affect the gain and structure of these modes. With minimal computation, this framework accurately predicts the emergence of the quasi-two-dimensional propagating waves observed in recent direct numerical simulations. Further, the analysis also enables the rational design of compliant surfaces, with properties optimized to suppress flow structures energetic in wall turbulence. It is shown that walls with unphysical negative damping are required to interact favourably with modes resembling the energetic near-wall cycle, which could explain why previous studies have met with limited success. Positive-damping walls are effective for modes resembling the so-called very-large-scale motions, indicating that compliant surfaces may be better suited for application at higher Reynolds number. Unfortunately, walls that suppress structures energetic in natural turbulence are also predicted to have detrimental effects elsewhere in spectral space. Consistent with previous experiments and simulations, slow-moving spanwise-constant structures are particularly susceptible to further amplification. Mitigating these adverse effects will be central to the development of compliant coatings that have a net positive influence on the flow.

A new active-closed-loop linear time-invariant robust control strategy for reducing drag in a channel flow is presented. Using the theory of passive systems, controllers were designed so as to make the closed-loop system as close to passive as is possible. Physically, this means minimising the maximum bound on perturbation energy production within the system. Although the resulting controllers are linear, the nonlinearity, inherent in all fluidic systems, is accounted for and included as an exogenous feedback forcing. Passivity-based controllers were generated for a Re = 2230 channel flow for a range of spatial wavenumber pairs and then analysed. They were shown to be capable of greatly restricting transient energy growth.

This paper extends the resolvent analysis of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382) to consider flow control techniques that employ linear control laws, focusing on opposition control (Choi, Moin & Kim, J. Fluid Mech., vol. 262, 1994, pp. 75–110) as an example. Under this formulation, the velocity field for turbulent pipe flow is decomposed into a series of highly amplified (rank-1) response modes, identified from a gain analysis of the Fourier-transformed Navier–Stokes equations. These rank-1 velocity responses represent propagating structures of given streamwise/spanwise wavelength and temporal frequency, whose wall-normal footprint depends on the phase speed of the mode. Opposition control, introduced via the boundary condition on wall-normal velocity, affects the amplification characteristics (and wall-normal structure) of these response modes; a decrease in gain indicates mode suppression, which leads to a decrease in the drag contribution from that mode. With basic assumptions, this rank-1 model reproduces trends observed in previous direct numerical simulation and large eddy simulation, without requiring high-performance computing facilities. Further, a wavenumber–frequency breakdown of control explains the deterioration of opposition control performance with increasing sensor elevation and Reynolds number. It is shown that slower-moving modes localized near the wall (i.e. attached modes) are suppressed by opposition control. Faster-moving detached modes, which are more energetic at higher Reynolds number and more likely to be detected by sensors far from the wall, are further amplified. These faster-moving modes require a phase lag between sensor and actuator velocity for suppression. Thus, the effectiveness of opposition control is determined by a trade-off between the modes detected by the sensor. However, it may be possible to develop control strategies optimized for individual modes. A brief exploration of such mode-optimized control suggests the potential for significant performance improvement.

We combine resolvent-mode decomposition with techniques from convex optimization to optimally approximate velocity spectra in a turbulent channel. The velocity is expressed as a weighted sum of resolvent modes that are dynamically significant, non-empirical, and scalable with Reynolds number. To optimally represent direct numerical simulations (DNS) data at friction Reynolds number 2003, we determine the weights of resolvent modes as the solution of a convex optimization problem. Using only 12 modes per wall-parallel wavenumber pair and temporal frequency, we obtain close agreement with DNS-spectra, reducing the wall-normal and temporal resolutions used in the simulation by three orders of magnitude.

Compressive sampling is well-known to be a useful tool used to resolve the energetic content of signals that admit a sparse representation. The broadband temporal spectrum acquired from point measurements in wall-bounded turbulence has precluded the prior use of compressive sampling in this kind of flow, however it is shown here that the frequency content of flow fields that have been Fourier transformed in the homogeneous spatial (wall-parallel) directions is approximately sparse, giving rise to a compact representation of the velocity field. As such, compressive sampling is an ideal tool for reducing the amount of information required to approximate the velocity field. Further, success of the compressive sampling approach provides strong evidence that this representation is both physically meaningful and indicative of special properties of wall turbulence. Another advantage of compressive sampling over periodic sampling becomes evident at high Reynolds numbers, since the number of samples required to resolve a given bandwidth with compressive sampling scales as the logarithm of the dynamically significant bandwidth instead of linearly for periodic sampling. The combination of the Fourier decomposition in the wall-parallel directions, the approximate sparsity in frequency, and empirical bounds on the convection velocity leads to a compact representation of an otherwise broadband distribution of energy in the space defined by streamwise and spanwise wavenumber, frequency, and wall-normal location. The data storage requirements for reconstruction of the full field using compressive sampling are shown to be significantly less than for periodic sampling, in which the Nyquist criterion limits the maximum frequency that can be resolved. Conversely, compressive sampling maximizes the frequency range that can be recovered if the number of samples is limited, resolving frequencies up to several times higher than the mean sampling rate. It is proposed that the approximate sparsity in frequency and the corresponding structure in the spatial domain can be exploited to design simulation schemes for canonical wall turbulence with significantly reduced computational expense compared with current techniques.

We study the Reynolds number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (2010), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier-Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behavior with Reynolds number on the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it was shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. We identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_��\approx 10^3 - 10^{10}$)., The paper is to appear in the Journal of Fluid Mechanics

We evaluate the efficacy of a gain-based rank-1 model, developed by McKeon & Sharma (J. Fluid Mech., 2010), for representing the energy spectra and the streamwise/wall-normal co-spectrum in a turbulent channel. This is motivated by our previous observation that the streamwise turbulent energy intensity is well approximated by the rank-1 model subject to a broadband forcing in the wall-parallel directions and a properly selected temporal intensity. In the present study, the evaluation is based on finding the optimal forcing spectrum that minimizes the deviation between the two-dimensional velocity spectra at different wall-normal locations obtained from direct numerical simulations at friction Reynolds number 2003 (Hoyas & Jiminenez, Phys. Fluids, 2006) and from the rank-1 model at equal Reynolds number. It is shown that the optimally forced rank-1 model captures the streamwise energy spectrum for streamwise wavelengths smaller than approximately 1000 viscous units throughout the channel. For larger wavelengths, the streamwise spectrum is matched in the outer region of the channel, i.e. wall-normal distances larger than approximately 0.15 times the channel half-height, and the mismatch close to the wall results in less than 5 percent error in the inner-scaled peak of the streamwise energy intensity. In addition, we show that the rank-1 model with optimal forcing captures the essential features of the wall-normal and spanwise spectra and the streamwise/wall-normal co-spectrum. We observe that the predicted magnitudes of the latter three spectra are smaller in the rank-1 model compared to the simulation results suggesting that a higher-order or different rank-1 model may be necessary for accurate representation of these spectra.

The resolvent-based analysis of wall turbulence (McKeon & Sharma 2010) is extended with an explicit treatment of the non-linearity in the Navier-Stokes equations. The equivalent of triadic interaction in the wall-normal direction is described, and it is found that the resulting forcing has a phase shift of π/2 in terms associated with symmetric spatial directions. Explicit forms for the nonlinearity in fully-developed pipe and channel flows, i.e. in cylindrical and Cartesian coordinate systems.

A direct transient growth analysis for three-dimensional perturbations to flow past a periodic array of T-106/300 low-pressure turbine fan blades is presented. The methodology is based on a singular value decomposition of the flow evolution operator, linearised about a steady or periodic base flow. This analysis yields the optimal growth modes. Previous work on global mode stability analysis of this flow geometry showed the flow is asymptotically stable, indicating a non-modal explanation of transition may be more appropriate. The present work extends previous investigations into the transient growth around a steady base flow, to higher Reynolds numbers and periodic base flows. It is found that the notable transient growth of the optimal modes suggests a plausible route to transition in comparison to modal growth for this configuration. The spatial extent and localisation of the optimal modes is examined and possible physical triggering mechanisms are discussed. It is found that for longer times and longer spanwise wavelengths, a separation in the shear layer excites the wake mode. For shorter times and spanwise wavelengths, smaller growth associated with excitation of the near wake are observed.

Bayesian estimation strategies represent the most fundamental formulation of the state estimation problem available, and apply readily to nonlinear systems with non-Gaussian uncertainties. The present paper introduces a novel method for implementing grid-based Bayesian estimation which largely sidesteps the severe computational expense that has prevented the widespread use of such methods. The method represents the evolution of the probability density function (PDF) in phase space, $p_{\x}(\x',t)$, discretized on a fixed Cartesian grid over {\it all} of phase space, and consists of two main steps: (i) Between measurement times, $p_{\x}(\x',t)$ is evolved via numerical discretization of the Kolmogorov forward equation, using a Godunov method with second-order corner transport upwind correction and a total variation diminishing flux limiter; (ii) at measurement times, $p_{\x}(\x',t)$ is updated via Bayes' theorem. Computational economy is achieved by exploiting the localised nature of $p_{\x}(\x',t)$. An ordered list of cells with non-negligible probability, as well as their immediate neighbours, is created and updated, and the PDF evolution is tracked {\it only} on these active cells. %The grid-based discretization of $p_{\x}(\x',t)$ in this approach avoids the requirement for resampling associated with particle-based representations of the PDF.

A method for finding reduced-order approximations of turbulent flow models is presented. The method preserves bounds on the production of turbulent energy in the sense of the [Formula: see text] norm of perturbations from a notional laminar profile. This is achieved by decomposing the Navier–Stokes system into a feedback arrangement between the linearized system and the remaining, normally neglected, nonlinear part. The linear system is reduced using a method similar to balanced truncation, but preserving bounds on the supply rate. The method involves balancing two algebraic Riccati equations. The bounds are then used to derive bounds on the turbulent energy production. An example of the application of the procedure to flow through a long straight pipe is presented. Comparison shows that the new method approximates the supply rate at least as well as, or better than, canonical balanced truncation.

The problems of nonlinearity and high dimension have so far prevented a complete solution of the control of turbulent flow. Addressing the problem of nonlinearity, we propose a flow control strategy which ensures that the energy of any perturbation to the target profile decays monotonically. The controller's estimate of the flow state is similarly guaranteed to converge to the true value. We present a one-time off-line synthesis procedure, which generalises to accommodate more restrictive actuation and sensing arrangements, with conditions for existence for the controller given in this case. The control is tested in turbulent channel flow ($Re_\tau=100$) using full-domain sensing and actuation on the wall-normal velocity. Concentrated at the point of maximum inflection in the mean profile, the control directly counters the supply of turbulence energy arising from the interaction of the wall-normal perturbations with the flow shear. It is found that the control is only required for the larger-scale motions, specifically those above the scale of the mean streak spacing. Minimal control effort is required once laminar flow is achieved. The response of the near-wall flow is examined in detail, with particular emphasis on the pressure and wall-normal velocity fields, in the context of Landahl's theory of sheared turbulence.