An official website of the United States government
The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($$Re$$)$$2000 \le Re \le 5000$$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$$Re^\gamma$$with the critical exponent$$\gamma \approx -1.25$$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$$\gamma \approx -1.25$$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.
more »
« less
Vorticity cascade and turbulent drag in wall-bounded flows: plane Poiseuille flow
Drag for wall-bounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wall-normal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number$$Re_\tau =1000$$. The net transfer from the wall of spanwise vorticity created by downstream pressure drop is due to two large opposing fluxes, one which is ‘down-gradient’ or outward from the wall, where most vorticity concentrates, and the other which is ‘up-gradient’ or toward the wall and acting against strong viscous diffusion in the near-wall region. We present evidence that the up-gradient/down-gradient transport occurs by a mechanism of correlated inflow/outflow and spanwise vortex stretching/contraction that was proposed by Lighthill. This mechanism is essentially Lagrangian, but we explicate its relation to the Eulerian anti-symmetric vorticity flux tensor. As evidence for the mechanism, we study (i) statistical correlations of the wall-normal velocity and of wall-normal flux of spanwise vorticity, (ii) vorticity flux cospectra identifying eddies involved in nonlinear vorticity transport in the two opposing directions and (iii) visualizations of coherent vortex structures which contribute to the transport. The ‘D-type’ vortices contributing to down-gradient transport in the log layer are found to be attached, hairpin-type vortices. However, the ‘U-type’ vortices contributing to up-gradient transport are detached, wall-parallel, pancake-shaped vortices with strong spanwise vorticity, as expected by Lighthill's mechanism. We discuss modifications to the attached eddy model and implications for turbulent drag reduction.
more »
« less
- Award ID(s):
- 2103874
- PAR ID:
- 10511081
- Publisher / Repository:
- CUP
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 974
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Direct numerical simulations of spanwise-rotating turbulent channel flow with a parabolic bump on the bottom wall are employed to investigate the effects of rotation on flow separation. Four rotation rates,$$Ro_b := 2\varOmega H/U_b = \pm 0.42$$,$$\pm$$1.0, are compared with the non-rotating scenario. The mild adverse pressure gradient induced by the lee side of the bump allows for a variable pressure-induced separation. The separation region is reduced (increased) when the bump is on the anti-cyclonic (cyclonic) side of the channel, compared with the non-rotating separation. The total drag is reduced in all rotating cases. Through several mechanisms, rotation alters the onset of separation, reattachment and wake recovery. The mean momentum deficit is found to be the key. A physical interpretation of the ratio between the system rotation and mean shear vorticity,$$S:=\varOmega /\varOmega _s$$, provides the mechanisms regarding stability thresholds$$S=-0.5$$and$$-$$1. The rotation effects are explained accordingly, with reference to the dynamics of several flow structures. For anti-cyclonic separation, particularly, the interaction between the Taylor–Görtler vortices and hairpin vortices of wall-bounded turbulence is proven to be responsible for the breakdown of the separating shear layer. A generalized argument is made regarding the essential role of near-wall deceleration and resultant ejection of enhanced hairpin vortices in destabilizing an anti-cyclonic flow. This mechanism is anticipated to have broad impacts on other applications in analogy to rotating shear flows, such as thermal convection and boundary layers over concave walls.more » « less
-
Two common definitions of the spatially local rate of kinetic energy cascade at some scale$$\ell$$in turbulent flows are (i) the cubic velocity difference term appearing in the ‘scale-integrated local Kolmogorov–Hill’ equation (structure-function approach), and (ii) the subfilter-scale energy flux term in the transport equation for subgrid-scale kinetic energy (filtering approach). We perform a comparative study of both quantities based on direct numerical simulation data of isotropic turbulence at Taylor-scale Reynolds number 1250. While in the past observations of negative subfilter-scale energy flux (backscatter) have led to debates regarding interpretation and relevance of such observations, we argue that the interpretation of the local structure-function-based cascade rate definition is unambiguous since it arises from a divergence term in scale space. Conditional averaging is used to explore the relationship between the local cascade rate and the local filtered viscous dissipation rate as well as filtered velocity gradient tensor properties such as its invariants. We find statistically robust evidence of inverse cascade when both the large-scale rotation rate is strong and the large-scale strain rate is weak. Even stronger net inverse cascading is observed in the ‘vortex compression’$$R>0$$,$$Q>0$$quadrant, where$$R$$and$$Q$$are velocity gradient invariants. Qualitatively similar but quantitatively much weaker trends are observed for the conditionally averaged subfilter-scale energy flux. Flow visualizations show consistent trends, namely that spatially, the inverse cascade events appear to be located within large-scale vortices, specifically in subregions when$$R$$is large.more » « less
-
Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative$$\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$$and residual$$\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$$components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in$$\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$$, whereas$$\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$$offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.more » « less
-
AbstractPrevious theoretical and simulation results indicate that anisotropic porous materials have the potential to reduce turbulent skin friction in wall-bounded flows. This study experimentally investigates the influence of anisotropy on the drag response of porous substrates. A family of anisotropic periodic lattices was manufactured using 3D printing. Rod spacing in different directions was varied systematically to achieve different ratios of streamwise, wall-normal, and spanwise bulk permeabilities ($$\kappa _{xx}$$ ,$$\kappa _{yy}$$ , and$$\kappa _{zz}$$ ). The 3D printed materials were flush-mounted in a benchtop water channel. Pressure drop measurements were taken in the fully developed region of the flow to systematically characterize drag for materials with anisotropy ratios$$\frac{\kappa _{xx}}{\kappa _{yy}} \in [0.035,28.6]$$ . Results show that all materials lead to an increase in drag compared to the reference smooth wall case over the range of bulk Reynolds numbers tested ($$\hbox {Re}_b \in [500,4000]$$ ). However, the relative increase in drag is lower for streamwise-preferential materials. We estimate that the wall-normal permeability for all tested cases exceeded the threshold identified in previous literature ($$\sqrt{\kappa _{yy}}^+> 0.4$$ ) for the emergence of energetic spanwise rollers similar to Kelvin–Helmholtz vortices, which can increase drag. The results also indicate that porous walls exhibit a departure from laminar behavior at different values for bulk Reynolds numbers depending on the geometry. Graphical abstractmore » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2023.609
