skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Structured Input-Output Analysis of Oblique Turbulent Bands in Transitional Plane Couette-Poiseuille Flow
In this work, we apply structured input-output analysis to study optimal perturbations and dominant flow patterns in transitional plane Couette-Poiseuille flow. The results demonstrate that this approach predicts the high structured gain of perturbations with wavelengths corresponding to the oblique turbulent bands observed in experiments. The inclination angles of these structures and their Reynolds number dependence are also consistent with previously observed trends. Reynolds number scalings of the maximally amplified structures for an intermediate laminar profile that is equally balanced between plane Couette and Poiseuille flow show an exponent that is at the midpoint of previously computed values for these two flows. However, the dependence of these scaling exponents on the shape of laminar flow as the relative contribution moves from predominately plane Couette to Poiseuille flow is not monotonic and our analysis indicates the emergence of different optimal perturbation structures through the parameter regime. Finally we adapt our approach to estimate the advection speeds of oblique turbulent bands in plane Couette flow and Poiseuille flow by computing their phase speed. The results show good agreement with prior predictions of the convection speeds of these structures from direct numerical simulations, which suggests that this framework has further potential in examining the dynamics of these structures.  more » « less
Award ID(s):
1652244
PAR ID:
10396185
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Structured Input-Output Analysis of Oblique Turbulent Bands in Transitional Plane Couette-Poiseuille Flow
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We employ a recently introduced structured input–output analysis (SIOA) approach to analyse streamwise and spanwise wavelengths of flow structures in stably stratified plane Couette flow. In the low-Reynolds-number ( $Re$ ) low-bulk Richardson number ( $$Ri_b$$ ) spatially intermittent regime, we demonstrate that SIOA predicts high amplification associated with wavelengths corresponding to the characteristic oblique turbulent bands in this regime. SIOA also identifies quasi-horizontal flow structures resembling the turbulent–laminar layers commonly observed in the high- $Re$ high- $$Ri_b$$ intermittent regime. An SIOA across a range of $$Ri_b$$ and $Re$ values suggests that the classical Miles–Howard stability criterion ( $$Ri_b\leq 1/4$$ ) is associated with a change in the most amplified flow structures when the Prandtl number is close to one ( $$Pr\approx 1$$ ). However, for $$Pr\ll 1$$ , the most amplified flow structures are determined by the product $$PrRi_b$$ . For $$Pr\gg 1$$ , SIOA identifies another quasi-horizontal flow structure that we show is principally associated with density perturbations. We further demonstrate the dominance of this density-associated flow structure in the high $Pr$ limit by constructing analytical scaling arguments for the amplification in terms of $Re$ and $Pr$ under the assumptions of unstratified flow (with $$Ri_b=0$$ ) and streamwise invariance. 
    more » « less
  2. Input–output analysis of transitional channel flows has proven to be a valuable analytical tool for identifying important flow structures and energetic motions. The traditional approach abstracts the nonlinear terms as forcing that is unstructured, in the sense that this forcing is not directly tied to the underlying nonlinearity in the dynamics. This paper instead employs a structured-singular-value-based approach that preserves certain input–output properties of the nonlinear forcing function in an effort to recover the larger range of key flow features identified through nonlinear analysis, experiments and direct numerical simulation (DNS) of transitional channel flows. Application of this method to transitional plane Couette and plane Poiseuille flows leads to not only the identification of the streamwise coherent structures predicted through traditional input–output approaches, but also the characterization of the oblique flow structures as those requiring the least energy to induce transition, in agreement with DNS studies, and nonlinear optimal perturbation analysis. The proposed approach also captures the recently observed oblique turbulent bands that have been linked to transition in experiments and DNS with very large channel size. The ability to identify the larger amplification of the streamwise varying structures predicted from DNS and nonlinear analysis in both flow regimes suggests that the structured approach allows one to maintain the nonlinear effects associated with weakening of the lift-up mechanism, which is known to dominate the linear operator. Capturing this key nonlinear effect enables the prediction of a wider range of known transitional flow structures within the analytical input–output modelling paradigm. 
    more » « less
  3. Inward radial Rayleigh-Be'nard-Poiseuille flow can exhibit a buoyancy-driven instability when the Rayleigh number exceeds a critical value. Furthermore, similar to plane Rayleigh-Be'nard-Poiseuille flow, a viscous Tollmien-Schlichting instability can occur when the Reynolds number is high enough. Direct numerical simulations were carried out with a compressible Navier-Stokes code in cylindrical coordinates to investigate the spatial stability of the inward radial flow inside the collector of a hypothetical solar chimney power plant. The convective terms were discretized with fifth-order-accurate upwind-biased compact finite-differences and the viscous terms were discretized with fourth-order-accurate compact finite differences. For cases with buoyancy-driven instability, steady three-dimensional waves are strongly amplified. The spatial growth rates vary significantly in the radial direction and lower azimuthal mode numbers are amplified closer to the outflow. Traveling oblique modes are amplified as well. The growth rates of the oblique modes decrease with increasing frequency. In addition to the purely radial flow, a spiral flow with swept inflow was examined. Overall lower growth rates are observed for the spiral flow compared to the radial flow. Different from the radial flow, the relative wave angles and growth rates of the left and right traveling oblique modes are not identical. A plane RBP case with viscosity-driven instability by Chung et al. was considered as well. The reported growth rate and phase speed were matched with good accuracy. 
    more » « less
  4. We present numerical simulations of laminar and turbulent channel flow of an elastoviscoplastic fluid. The non-Newtonian flow is simulated by solving the full incompressible Navier–Stokes equations coupled with the evolution equation for the elastoviscoplastic stress tensor. The laminar simulations are carried out for a wide range of Reynolds numbers, Bingham numbers and ratios of the fluid and total viscosity, while the turbulent flow simulations are performed at a fixed bulk Reynolds number equal to 2800 and weak elasticity. We show that in the laminar flow regime the friction factor increases monotonically with the Bingham number (yield stress) and decreases with the viscosity ratio, while in the turbulent regime the friction factor is almost independent of the viscosity ratio and decreases with the Bingham number, until the flow eventually returns to a fully laminar condition for large enough yield stresses. Three main regimes are found in the turbulent case, depending on the Bingham number: for low values, the friction Reynolds number and the turbulent flow statistics only slightly differ from those of a Newtonian fluid; for intermediate values of the Bingham number, the fluctuations increase and the inertial equilibrium range is lost. Finally, for higher values the flow completely laminarizes. These different behaviours are associated with a progressive increases of the volume where the fluid is not yielded, growing from the centreline towards the walls as the Bingham number increases. The unyielded region interacts with the near-wall structures, forming preferentially above the high-speed streaks. In particular, the near-wall streaks and the associated quasi-streamwise vortices are strongly enhanced in an highly elastoviscoplastic fluid and the flow becomes more correlated in the streamwise direction. 
    more » « less
  5. The heat and mass transfer characteristics of a simple shear flow over a surface covered with staggered herringbone structures are numerically investigated using the lattice Boltzmann method. Two flow motions are identified. The first is a spiral flow oscillation above the herringbone structures that advect heat and mass from the top plane to herringbone structures. The second is a flow recirculation in the grooves between the ridges that advect heat and mass from the area around the tips of the structures to their side walls and the bottom surfaces. These two basic flow motions couple together to form a complex transport mechanism. The results show that when advective heat and mass transfer takes effect at relatively large Reynolds and Schmidt numbers, the dependence of the total transfer rate on Schmidt number follows a power law, with the exponent being the same as that in the Dittus–Boelter equation for turbulent heat transfer. As the Reynolds number increases, the dependence of the total transfer rate on the Reynolds number also approaches a power law, and the exponent is close to that in the Dittus–Boelter equation. 
    more » « less