Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA statistical state dynamics-based study of the structure and mechanism of large-scale motions in plane Poiseuille flow
153
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The perspective of statistical state dynamics (SSD) has recently been applied to the study of mechanisms underlying turbulence in various physical systems. An example implementation of SSD is the second order closure referred to as stochastic structural stability theory (S3T), which has provided insight into the dynamics of wall turbulence and specifically the emergence and maintenance of the roll/streak structure. This closure eliminates nonlinear interactions among the perturbations has been removed, restricting nonlinearity in the dynamics to that of the mean equation and the interaction between the mean and perturbation covariance. Here, this quasi-linear restriction of the dynamics is used to study the structure and dynamics of turbulence in plane Poiseuille flow at moderately high Reynolds numbers in a closely related dynamical system, referred to as the restricted nonlinear (RNL) system. RNL simulations reveal that the essential features of wall-turbulence dynamics are retained. Remarkably, the RNL system spontaneously limits the support of its turbulence to a small set of streamwise Fourier components giving rise to a naturally minimal representation of its turbulence dynamics. Although greatly simplified, this RNL turbulence exhibits natural-looking structures and statistics. Surprisingly, even when further truncation of the perturbation support to a single streamwise component is imposed the RNL system continues to produce self-sustaining turbulent structure and dynamics. RNL turbulence at the Reynolds numbers studied is dominated by the roll/streak structure in the buffer layer and similar very-large-scale structure (VLSM) in the outer layer. Diagnostics of the structure, spectrum and energetics of RNL and DNS turbulence are used to demonstrate that the roll/streak dynamics supporting the turbulence in the buffer and logarithmic layer is essentially similar in RNL and DNS.
rate research
Read More
Turbulence structure in the quasi-linear restricted nonlinear (RNL) model is analyzed and compared with DNS of turbulent Poiseuille flow at Reynolds number R=1650. The turbulence structure is obtained by POD analysis of the two components of the flow partition used in formulating the RNL model: the streamwise-mean flow and the associated perturbations. The dominant structures are found to be similar in RNL simulations and DNS despite the neglect of perturbation-perturbation nonlinearity in the RNL formulation. POD analysis of the streamwise-mean flow indicates that the dominant structure in both RNL and DNS is a coherent roll-streak structure in which the roll is collocated with the streak in a manner configured to reinforce the streak by the lift-up process. This mechanism of roll-streak maintenance accords with analytical predictions made using the second order statistical state dynamics (SSD) model, referred to as S3T, which shares with RNL the dynamical restriction of neglecting the perturbation-perturbation nonlinearity. POD analysis of perturbations from the streamwise-mean streak reveals that similar structures characterize these perturbations in both RNL and DNS. The perturbation to the low-speed streak POD are shown to have the form of oblique waves collocated with the streak that can be identified with optimally growing structures on the streak. Given that the mechanism sustaining turbulence in RNL has been analytically characterized, this close correspondence between the streamwise-mean and perturbation structures in RNL and DNS supports the conclusion that the self-sustaining mechanism in DNS is the same as that in RNL.
The convection velocity of localized wave packet in plane-Poiseuille flow is found to be determined by a solitary wave at the centerline of a downstream vortex dipole in its mean field after deducting the basic flow. The fluctuation component following the vortex dipole oscillates with a global frequency selected by the upstream marginal absolute instability, and propagates obeying the local dispersion relation of the mean flow. By applying localized initial disturbances, a nonzero wave-packet density is achieved at the threshold state, suggesting a first order transition.
The ultimate goal of a sound theory of turbulence in fluids is to close in a rational way the Reynolds equations, namely to express the time averaged turbulent stress tensor as a function of the time averaged velocity field. This closure problem is a deep and unsolved problem of statistical physics whose solution requires to go beyond the assumption of a homogeneous and isotropic state, as fluctuations in turbulent flows are strongly related to the geometry of this flow. This links the dissipation to the space dependence of the average velocity field. Based on the idea that dissipation in fully developed turbulence is by singular events resulting from an evolution described by the Euler equations, it has been recently observed that the closure problem is strongly restricted, and that it implies that the turbulent stress is a non local function (in space) of the average velocity field, an extension of classical Boussinesq theory of turbulent viscosity. The resulting equations for the turbulent stress are derived here in one of the simplest possible physical situation, the turbulent Poiseuille flow between two parallel plates. In this case the integral kernel giving the turbulent stress, as function of the averaged velocity field, takes a simple form leading to a full analysis of the averaged turbulent flow in the limit of a very large Reynolds number. In this limit one has to match a viscous boundary layer, near the walls bounding the flow, and an outer solution in the bulk of the flow. This asymptotic analysis is non trivial because one has to match solution with logarithms. A non trivial and somewhat unexpected feature of this solution is that, besides the boundary layers close to the walls, there is another inner boundary layer near the center plane of the flow.
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.
A modal stability analysis shows that plane Poiseuille flow of an Oldroyd-B fluid becomes unstable to a `center mode with phase speed close to the maximum base-flow velocity, $U_{max}$. The governing dimensionless groups are the Reynolds number $Re = rho U_{max} H/eta$, the elasticity number $E = lambda eta/(H^2rho)$, and the ratio of solvent to solution viscosity $eta_s/eta$; here, $lambda$ is the polymer relaxation time, $H$ is the channel half-width, and $rho$ is the fluid density. For experimentally relevant values (e.g., $E sim 0.1$ and $beta sim 0.9$), the predicted critical Reynolds number, $Re_c$, for the center-mode instability is around $200$, with the associated eigenmodes being spread out across the channel. In the asymptotic limit of $E(1 -beta) ll 1$, with $E$ fixed, corresponding to strongly elastic dilute polymer solutions, $Re_c propto (E(1-beta))^{-frac{3}{2}}$ and the critical wavenumber $k_c propto (E(1-beta))^{-frac{1}{2}}$. The unstable eigenmode in this limit is confined in a thin layer near the channel centerline. The above features are largely analogous to the center-mode instability in viscoelastic pipe flow (Garg et al., Phys. Rev. Lett., 121, 024502 (2018)), and suggest a universal linear mechanism underlying the onset of turbulence in both channel and pipe flows of suffciently elastic dilute polymer solutions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
