Search

The article "Dynamics of turbulent energy and dissipation in channel flow" delves into the intricate relationship between turbulent kinetic energy (TKE), turbulence dissipation rate (TDR), and turbulence production rate (TPR) in fully developed turbulent channel flow. The research uncovers time lags between these variables, with correlations that vary based on Reynolds number. The findings emphasize the complex nature of energy transfer and dissipation within different regions of the channel, showcasing correlations that involve both time and wall-distance lags. This study offers valuable insights into the self-sustaining process of turbulent flows, with implications for understanding turbulence dynamics across various contexts. [Extracted from the article]

In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ($Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components. [ABSTRACT FROM AUTHOR]

In the present study, the shape of a two-dimensional cylinder is optimised to minimise the mean drag in laminar unsteady flow under a noisy environment. A small inline stochastic oscillation in the free-stream velocity, which follows the Ornstein–Uhlenbeck process, is considered for the noise. The small noise is found to yield a large random fluctuation in instantaneous drag of the cylinder due to the effect of added mass. Subject to the strong random fluctuation of drag, the shape optimisation is performed using an ensemble-variation-based method (EnVar), as the conventional adjoint-based optimisation is not applicable to such a flow environment with unknown free-stream noise. The optimised cylinder geometry is found to be a nearly-symmetric slender oval at a low Reynolds number. As the Reynolds number is increased, two optimal shapes emerge: one is identical to the oval obtained at the low Reynolds number, and the other is an asymmetric oval, the rear side of which is more slender than the front side, reminiscent of an aerofoil. Despite the large random fluctuation in the instantaneous drag, the optimal cylinder shapes obtained for different levels of the upstream noise are found to be almost identical. It is shown that the robust nature of the optimal cylinder shape originates from the limited influence of the small upstream noise on the mean flow properties of the cylinder wake. Finally, the optimised cylinder primarily reduces the pressure component of the drag, associated mainly with vortex shedding in the wake, and this is achieved by marginally increasing the viscous drag through the shape change. [ABSTRACT FROM AUTHOR]

A set of several exact coherent states in plane Couette flow is computed under spanwise wall oscillation control, with a range of wall oscillation amplitudes and periods $({A_w}, T)$. It is found that the wall oscillation generally stabilises the upper branch of the equilibrium solutions and achieves the corresponding drag reduction, while it influences modestly the lower branch. The stabilisation effect is found to increase with the oscillation amplitude with an optimal time period around ${T^{+}} \approx 100$. The exact coherent states reproduce some key dynamical behaviours of streaks observed in previous studies, while exhibiting the rich coherent structure dynamics that cannot be extracted from a phase average of turbulent states. Visualisation of state portraits shows that the size of the state space supporting turbulent solution is reduced by the spanwise wall oscillation, and the upper-branch equilibrium solutions become less repelling, with many of their unstable manifolds being stabilised. This change of the state space dynamics leads to a significant reduction in lifetime of turbulence. Finally, the main stabilisation mechanism of the exact coherent states is found to be the suppression of the lift-up effect of streaks, explaining why previous linear analyses have been so successful for turbulence stabilisation modelling and the resulting drag reduction. [ABSTRACT FROM AUTHOR]

The computation of invariant solutions and the visualisation of the associated state space have played a key role in the understanding of transition and the self-sustaining process in wall-bounded shear flows. In this study, an extension of this approach is sought for a turbulent flow which explicitly exhibits multi-scale behaviour. The minimal unit of multi-scale near-wall turbulence, which resolves two adjacent spanwise integral length scales of motion, is considered using a shear stress-driven flow model (Doohan, Willis & Hwang J. Fluid Mech., vol. 913, 2021, A8). The edge state, 26 travelling waves and two periodic orbits are computed, which represent either the large- or small-scale self-sustaining processes. Given that the spanwise length scales are not widely separated here, it could be envisaged that turbulent trajectories visit these solutions in the state space. Considering the intra- and inter-scale dynamics of the flow, numerous phase portraits are examined, but the turbulent state is not found to approach any of these solutions. A detailed analysis reveals that this is due to the lack of scale interaction processes captured by the invariant solutions, including the mean–fluctuation interaction, the energy cascade in the streamwise wavenumber space and the cascade-driven energy production discovered recently. There is a single solution that resembles turbulence much more than the others, which captures two-scale energetics and a scale interaction process involving energy feeding from small to large spanwise scales through the subharmonic sinuous streak instability mode. Based on these observations, it is conjectured that the state space view of turbulent trajectories wandering between solutions would need suitable refinement to model multi-scale turbulence, when each solution does not represent multi-scale processes of turbulence. In particular, invariant solutions that are inherently multi-scale would be required. [ABSTRACT FROM AUTHOR]

Recent numerical experiments have shown that the temporal dynamics of isolated energy-containing eddies in the hierarchy of wall-bounded turbulence are governed by the self-sustaining process (SSP). However, high-Reynolds-number turbulence is a multi-scale phenomenon and exhibits interaction between the structures of different scales, but the dynamics of such multi-scale flows are poorly understood. In this study, the temporal dynamics of near-wall turbulent flow with two integral length scales of motion are investigated using a shear stress-driven flow model (Doohan et al., J. Fluid Mech., vol. 874, 2019, pp. 606–638), with a focus on identifying scale interaction processes through the governing equations and relating these to the SSPs at each scale. It is observed that the dynamics of the energy cascade from large to small scales is entirely determined by the large-scale SSP and the timing of the corresponding inter-scale turbulent transport coincides with the large-scale streak breakdown stage. Furthermore, the characteristic time scales of the resulting small-scale dissipation match those of the large-scale SSP, indicative of non-equilibrium turbulent dissipation dynamics. A new scale interaction process is identified, namely that the transfer of wall-normal energy from large to small scales drives small-scale turbulent production via the Orr mechanism. While the main outcome of this driving process appears to be the transient amplification of localised small-scale velocity structures and their subsequent dissipation, it also has an energising effect on the small-scale SSP. Finally, the feeding of energy from small to large scales is impelled by the small-scale SSP and coincides with the small-scale streak instability stage. The streamwise feeding process seems to be related to the subharmonic sinuous streak instability mode in particular and leads to the formation of the wall-reaching part of high-speed large-scale streaks. [ABSTRACT FROM AUTHOR]

The spectral energetics of a quasilinear (QL) model is studied in uniform shear turbulence. For the QL approximation, the velocity is decomposed into a mean averaged in the streamwise direction and the remaining fluctuation. The equations for the mean are fully considered, while the equations for the fluctuation are linearised around the mean. The QL model exhibits an energy cascade in the spanwise direction, but this is mediated by highly anisotropic small-scale motions unlike that in direct numerical simulation mediated by isotropic motions. In the streamwise direction the energy cascade is shown to be completely inhibited in the QL model, resulting in highly elevated spectral energy intensity residing only at the streamwise integral length scales. It is also found that the streamwise wavenumber spectra of turbulent transport, obtained with the classical Reynolds decomposition, statistically characterizes the instability of the linearised fluctuation equations. Further supporting evidence of this claim is presented by carrying out a numerical experiment, in which the QL model with a single streamwise Fourier mode is found to generate the strongest turbulence for $L_x/L_z=1\sim 3$ , consistent with previous findings ($L_x$ and $L_z$ are the streamwise and spanwise computational domains, respectively). Finally, the QL model is shown to completely ignore the role of slow pressure in the fluctuations, resulting in a significant damage of pressure-strain transport at all length scales. This explains the anisotropic turbulence of the QL model throughout the entire wavenumber space as well as the inhibited nonlinear regeneration of streamwise vortices in the self-sustaining process. [ABSTRACT FROM AUTHOR]

Kessler (Nature, vol. 313, 1985, pp. 218–220) first showed that plume-like structures spontaneously appear from both stationary and flowing suspensions of gyrotactic microswimmers in a vertical pipe. Recently, it has been shown that there exist multiple steady, axisymmetric and axially uniform solutions to such a system (Bees & Croze, Proc. R. Soc. A, vol. 466, 2010, pp. 2057–2077). In the present study, we generalise this finding by reporting that a countably infinite number of such solutions emerge as the Richardson number increases. Linear stability, weakly nonlinear and fully nonlinear analyses are performed, revealing that each of the solutions arises from the destabilisation of a uniform suspension. The countability of the solutions is due to the finite flow domain, while the transcritical nature of the bifurcation is because of the cylindrical geometry, which breaks the horizontal symmetry of the system. It is further shown that there exists a maximum threshold of achievable downward flow rate for each solution if the flow is to remain steady, as varying the pressure gradient can no longer increase the flow rate from the solution. All of the solutions found are unstable, except for the one arising at the lowest Richardson number, implying that they would play a role in the transient dynamics in the route from a uniform suspension to the fully developed gyrotactic pattern. [ABSTRACT FROM AUTHOR]

Townsend's model of attached eddies for boundary layers is revisited within a quasi-linear approximation. The velocity field is decomposed into a mean profile and fluctuations. While the mean is obtained from the nonlinear equations, the fluctuations are modelled by replacing the nonlinear self-interaction terms with an eddy-viscosity-based turbulent diffusion and stochastic forcing. Under this particular approximation, the resulting fluctuation equations remain linear, enabling solutions to be superposed, the same theoretical idea used in the original attached eddy model. By leveraging this feature, the stochastic forcing is determined self-consistently by solving an optimisation problem which minimises the difference between the Reynolds shear stresses from the mean and fluctuation equations, subject to a constraint that the averaged Reynolds shear-stress spectrum is sufficiently smooth in the spatial wavenumber space. The proposed quasi-linear approximation is subsequently applied to channel flow for Reynolds number $Re_{\unicode[STIX]{x1D70F}}$ ranging from 500 to 20 000. The best result is obtained when the Reynolds stress is calculated by retaining only the two leading proper orthogonal decomposition modes, which further filters out the modelling artefact caused by the unphysical stochastic forcing. In this case, the resulting turbulence intensity profile and energy spectra exhibit the same qualitative behaviour as direct numerical simulation (DNS) data throughout the entire wall-normal domain, while reproducing the early theoretical predictions of the original attached eddy model within a controlled approximation to the Navier–Stokes equations. Finally, the proposed quasi-linear approximation reveals that the peak streamwise and spanwise turbulence intensities may deviate slightly from the logarithmic scaling with the Reynolds number for $Re_{\unicode[STIX]{x1D70F}}\gtrsim 10\,000$ , and the supporting evidence is presented using the existing DNS data. [ABSTRACT FROM AUTHOR]

An inner-scaled, shear stress-driven flow is considered as a model of independent near-wall turbulence as the friction Reynolds number $Re_{\unicode[STIX]{x1D70F}}\rightarrow \infty$. In this limit, the model is applicable to the near-wall region and the lower part of the logarithmic layer of various parallel shear flows, including turbulent Couette flow, Poiseuille flow and Hagen–Poiseuille flow. The model is validated against damped Couette flow and there is excellent agreement between the velocity statistics and spectra for the wall-normal height  $y^{+}. A near-wall flow domain of similar size to the minimal unit is analysed from a dynamical systems perspective. The edge and fifteen invariant solutions are computed, the first discovered for this flow configuration. Through continuation in the spanwise width $L_{z}^{+}$ , the bifurcation behaviour of the solutions over the domain size is investigated. The physical properties of the solutions are explored through phase portraits, including the energy input and dissipation plane, and streak, roll and wave energy space. Finally, a Reynolds number is defined in outer units and the high- $Re$ asymptotic behaviour of the equilibria is studied. Three lower branch solutions are found to scale consistently with vortex–wave interaction (VWI) theory, with wave forcing localising around the critical layer. [ABSTRACT FROM AUTHOR]

A new set of exact coherent states in the form of a travelling wave is reported in plane channel flow. They are continued over a range in $Re$ from approximately $2600$ up to $30\,000$ , an order of magnitude higher than those discovered in the transitional regime. This particular type of exact coherent states is found to be gradually more localised in the near-wall region on increasing the Reynolds number. As larger spanwise sizes $L_{z}^{+}$ are considered, these exact coherent states appear via a saddle-node bifurcation with a spanwise size of $L_{z}^{+}\simeq 50$ and their phase speed is found to be $c^{+}\simeq 11$ at all the Reynolds numbers considered. Computation of the eigenspectra shows that the time scale of the exact coherent states is given by $h/U_{cl}$ in channel flow at all Reynolds numbers, and it becomes equivalent to the viscous inner time scale for the exact coherent states in the limit of $Re\rightarrow \infty$. The exact coherent states at several different spanwise sizes are further continued to a higher Reynolds number, $Re=55\,000$ , using the eddy-viscosity approach (Hwang & Cossu, Phys. Rev. Lett. , vol. 105, 2010, 044505). It is found that the continued exact coherent states at different sizes are self-similar at the given Reynolds number. These observations suggest that, on increasing Reynolds number, new sets of self-sustaining coherent structures are born in the near-wall region. Near this onset, these structures scale in inner units, forming the near-wall self-sustaining structures. With further increase of Reynolds number, the structures that emerged at lower Reynolds numbers subsequently evolve into the self-sustaining structures in the logarithmic region at different length scales, forming a hierarchy of self-similar coherent structures as hypothesised by Townsend (i.e. attached eddy hypothesis). Finally, the energetics of turbulent flow is discussed for a consistent extension of these dynamical systems notions to high Reynolds numbers. [ABSTRACT FROM AUTHOR]

We investigate the nonlinear phase-space dynamics of plane Couette flow and plane Poiseuille flow under the action of opposition control at low Reynolds numbers in domains close to the minimal unit. In Couette flow, the effect of the control is analysed by focussing on a pair of non-trivial equilibrium solutions. It is found that the control only slightly modifies the statistics, turbulent skin friction and phase-space projection of the lower-branch equilibrium solution, which, in this case, is in fact identical to the edge state. On the other hand, the upper-branch equilibrium solution and mean turbulent state are modified considerably when the control is applied. In phase space, they gradually approach the lower-branch equilibrium solution on increasing the control amplitude, and this results in an elevation of the critical Reynolds number at which the equilibrium solutions first occur via a saddle-node bifurcation. It is also found that the upper-branch equilibrium solution is stabilised by the control. In Poiseuille flow, we study an unstable periodic orbit on the edge state and find that it, too, is modified very little by opposition control. We again observe that the turbulent state gradually approaches the edge state in phase space as the control amplitude is increased. In both flows, we find that the control significantly reduces the fluctuating strength of the turbulent state in phase space. However, the reduced distance between the turbulent trajectory and the edge state yields a significant reduction in turbulence lifetimes for both Couette and Poiseuille flow. This demonstrates that opposition control greatly increases the probability of the trajectory escaping from the turbulent state, which takes the form of a chaotic saddle. [ABSTRACT FROM AUTHOR]

Large-scale motions in wall-bounded turbulent flows are frequently interpreted as resulting from an aggregation process of smaller-scale structures. Here, we explore the alternative possibility that such large-scale motions are themselves self-sustained and do not draw their energy from smaller-scale turbulent motions activated in buffer layers. To this end, it is first shown that large-scale motions in turbulent Couette flow at Re=2150 self-sustain, even when active processes at smaller scales are artificially quenched by increasing the Smagorinsky constant Cs in large-eddy simulations (LES). These results are in agreement with earlier results on pressure-driven turbulent channel flows. We further investigate the nature of the large-scale coherent motions by computing upper- and lower-branch nonlinear steady solutions of the filtered (LES) equations with a Newton-Krylov solver, and find that they are connected by a saddle-node bifurcation at large values of Cs. Upper-branch solutions for the filtered large-scale motions are computed for Reynolds numbers up to Re =2187 using specific paths in the Re-Cs parameter plane and compared to large-scale coherent motions. Continuation to Cs= 0 reveals that these large-scale steady solutions of the filtered equations are connected to the Nagata-Clever-Busse-Waleffe branch of steady solutions of the Navier-Stokes equations. In contrast, we find it impossible to connect the latter to buffer-layer motions through a continuation to higher Reynolds numbers in minimal flow units. [ABSTRACT FROM AUTHOR]

The linear growth of the spanwise correlation length scale with the distance from the wall in the logarithmic region of wall-bounded turbulent flows has been understood as a reflection of Townsend’s attached eddies. Based on this observation, in the present study, we perform a numerical experiment, which simulates energy-containing motions only at a given spanwise length scale in the logarithmic region, using their self-sustaining nature found recently. The self-sustaining energy-containing motions at each of the spanwise length scales are found to be self-similar with respect to the given spanwise length. Furthermore, their statistical structures are consistent with those of the attached eddies in the original theory, providing direct evidence on the existence of Townsend’s attached eddies. It is shown that a single self-sustaining attached eddy is composed of two distinct elements, one of which is a long streaky motion reaching the near-wall region, and the other is a relatively short vortical structure carrying all the velocity components. For the given spanwise length ${\it\lambda}_{z}$ between ${\it\lambda}_{z}^{+}=100$ and ${\it\lambda}_{z}\simeq 1.5h$, where $h$ is half the height of the channel, the former is found to be self-similar along $y\simeq 0.1{\it\lambda}_{z}$ and ${\it\lambda}_{x}\simeq 10{\it\lambda}_{z}$, while the latter is self-similar along $y\simeq 0.5{\it\lambda}_{z}\sim 0.7{\it\lambda}_{z}$ and ${\it\lambda}_{x}\simeq 2{\it\lambda}_{z}\sim 3{\it\lambda}_{z}$ where $y$ is the wall-normal direction. The scaling suggests that the smallest attached eddy would be a near-wall coherent motion in the form of a streak and quasi-streamwise vortices aligned to that, whereas the largest one would be an outer motion with a very-large-scale motion (VLSM) and large-scale motions (LSMs) aligned to that. The attached eddies in between, the size of which is proportional to their distance from the wall, contribute to the logarithmic region and fill the space caused by the length scale separation. The scaling is also found to yield behaviour consistent with the emergence of $k_{x}^{-1}$ spectra in a number of previous studies. Finally, a further discussion is provided, in particular on Townsend’s inactive motion and several recent theoretical findings. [ABSTRACT FROM PUBLISHER]

Numerical experiments that remove turbulent motions wider than ${ \lambda }_{z}^{+ } \simeq 100$ are carried out up to ${\mathit{Re}}_{\tau } = 660$ in a turbulent channel. The artificial removal of the wide outer turbulence is conducted with spanwise minimal computational domains and an explicit filter that effectively removes spanwise uniform eddies. The mean velocity profile of the remaining motions shows very good agreement with that of the full simulation below ${y}^{+ } \simeq 40$, and the near-wall peaks of the streamwise velocity fluctuation scale very well in the inner units and remain almost constant at all the Reynolds numbers considered. The self-sustaining motions narrower than ${ \lambda }_{z}^{+ } \simeq 100$ generate smaller turbulent skin friction than full turbulent motions, and their contribution to turbulent skin friction gradually decays with the Reynolds number. This finding suggests that the role of the removed outer structures becomes increasingly important with the Reynolds number; thus one should aim to control the large scales for turbulent drag reduction at high Reynolds numbers. In the near-wall region, the streamwise and spanwise velocity fluctuations of the motions of ${ \lambda }_{z}^{+ } \leq 100$ reveal significant lack of energy at long streamwise lengths compared to those of the full simulation. In contrast, the losses of the wall-normal velocity and the Reynolds stress are not as large as those of these two variables. This implies that the streamwise and spanwise velocities of the removed motions penetrate deep into the near-wall region, while the wall-normal velocity and the Reynolds stress do not. [ABSTRACT FROM PUBLISHER]

The linear response to stochastic and optimal harmonic forcing of small coherent perturbations to the turbulent channel mean flow is computed for Reynolds numbers ranging from Reτ = 500 to 20000. Even though the turbulent mean flow is linearly stable, it is nevertheless able to sustain large amplifications by the forcing. The most amplified structures consist of streamwise-elongated streaks that are optimally forced by streamwise-elongated vortices. For streamwise-elongated structures, the mean energy amplification of the stochastic forcing is found to be, to a first approximation, inversely proportional to the forced spanwise wavenumber while it is inversely proportional to its square for optimal harmonic forcing in an intermediate spanwise wavenumber range. This scaling can be explicitly derived from the linearized equations under the assumptions of geometric similarity of the coherent perturbations and of logarithmic base flow. Deviations from this approximate power-law regime are apparent in the pre-multiplied energy amplification curves that reveal a strong influence of two different peaks. The dominant peak scales in outer units with the most amplified spanwise wavelength of λz ≈ 3.5h, while the secondary peak scales in wall units with the most amplified λz+ ≈ 80. The associated optimal perturbations are almost independent of the Reynolds number when, respectively, scaled in outer and inner units. In the intermediate wavenumber range, the optimal perturbations are approximatively geometrically similar. Furthermore, the shape of the optimal perturbations issued from the initial value, the harmonic forcing and the stochastic forcing analyses are almost indistinguishable. The optimal streaks corresponding to the large-scale peak strongly penetrate into the inner layer, where their amplitude is proportional to the mean-flow profile. At the wavenumbers corresponding to the large-scale peak, the optimal amplifications of harmonic forcing are at least two orders of magnitude larger than the amplifications of the variance of stochastic forcing and both increase with the Reynolds number. This confirms the potential of the artificial forcing of optimal large-scale streaks for the flow control of wall-bounded turbulent flows. [ABSTRACT FROM AUTHOR]

Understanding the origin of large-scale structures in high-Reynolds-number wall turbulence has been a central issue over a number of years. Recently, Rawat et al. (J. Fluid Mech., vol. 782, 2015, pp. 515–540) have computed invariant solutions for the large-scale structures in turbulent Couette flow at $Re_{\unicode[STIX]{x1D70F}}\simeq 128$ using an overdamped large-eddy simulation with the Smagorinsky model to account for the effect of the surrounding small-scale motions. Here, we extend this approach to Reynolds numbers an order of magnitude higher in turbulent channel flow, towards the regime where the large-scale structures in the form of very-large-scale motions (long streaky motions) and large-scale motions (short vortical structures) emerge energetically. We demonstrate that a set of invariant solutions can be computed from simulations of the self-sustaining large-scale structures in the minimal unit (domain of size $L_{x}=3.0h$ streamwise and $L_{z}=1.5h$ spanwise) with midplane reflection symmetry at least up to $Re_{\unicode[STIX]{x1D70F}}\simeq 1000$. By approximating the surrounding small scales with an artificially elevated Smagorinsky constant, a set of equilibrium states are found, labelled upper- and lower-branch according to their associated drag. It is shown that the upper-branch equilibrium state is a reasonable proxy for the spatial structure and the turbulent statistics of the self-sustaining large-scale structures. [ABSTRACT FROM PUBLISHER]