scholarly journals Resolvent analysis of turbulent channel flow with manipulated mean velocity profile

2020 ◽  
Vol 15 (3) ◽  
pp. JFST0014-JFST0014
Author(s):  
Riko UEKUSA ◽  
Aika KAWAGOE ◽  
Yusuke NABAE ◽  
Koji FUKAGATA
1967 ◽  
Vol 27 (2) ◽  
pp. 253-272 ◽  
Author(s):  
W. C. Reynolds ◽  
W. G. Tiederman

The Orr-Sommerfeld stability problem has been studied for velocity profiles appropriate to turbulent channel flow. The intent was to provide an evaluation of Malkus's theory that the flow assumes a state of maximum dissipation, subject to certain constraints, one of which is that the mean velocity profile is marginally stable. Dissipation rates and neutral stability curves were obtained for a representative two-parameter family of velocity profiles. Those in agreement with experimental profiles were found to be stable; the marginally stable profile of greatest dissipation was not in good agreement with experiments. An explanation for the apparent success of Malkus's theory is offered.


2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.


2019 ◽  
Vol 863 ◽  
pp. 1190-1203 ◽  
Author(s):  
Sabarish B. Vadarevu ◽  
Sean Symon ◽  
Simon J. Illingworth ◽  
Ivan Marusic

We study the evolution of velocity fluctuations due to an isolated spatio-temporal impulse using the linearized Navier–Stokes equations. The impulse is introduced as an external body force in incompressible channel flow at $Re_{\unicode[STIX]{x1D70F}}=10\,000$. Velocity fluctuations are defined about the turbulent mean velocity profile. A turbulent eddy viscosity is added to the equations to fix the mean velocity as an exact solution, which also serves to model the dissipative effects of the background turbulence on large-scale fluctuations. An impulsive body force produces flow fields that evolve into coherent structures containing long streamwise velocity streaks that are flanked by quasi-streamwise vortices; some of these impulses produce hairpin vortices. As these vortex–streak structures evolve, they grow in size to be nominally self-similar geometrically with an aspect ratio (streamwise to wall-normal) of approximately 10, while their kinetic energy density decays monotonically. The topology of the vortex–streak structures is not sensitive to the location of the impulse, but is dependent on the direction of the impulsive body force. All of these vortex–streak structures are attached to the wall, and their Reynolds stresses collapse when scaled by distance from the wall, consistent with Townsend’s attached-eddy hypothesis.


2017 ◽  
Vol 830 ◽  
pp. 300-325 ◽  
Author(s):  
Hiroyuki Abe ◽  
Robert Anthony Antonia

Integration across a fully developed turbulent channel flow of the transport equations for the mean and turbulent parts of the scalar dissipation rate yields relatively simple relations for the bulk mean scalar and wall heat transfer coefficient. These relations are tested using direct numerical simulation datasets obtained with two isothermal boundary conditions (constant heat flux and constant heating source) and a molecular Prandtl number Pr of 0.71. A logarithmic dependence on the Kármán number $h^{+}$ is established for the integrated mean scalar in the range $h^{+}\geqslant 400$ where the mean part of the total scalar dissipation exhibits near constancy, whilst the integral of the turbulent scalar dissipation rate $\overline{\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}}$ increases logarithmically with $h^{+}$. This logarithmic dependence is similar to that established in a previous paper (Abe & Antonia, J. Fluid Mech., vol. 798, 2016, pp. 140–164) for the bulk mean velocity. However, the slope (2.18) for the integrated mean scalar is smaller than that (2.54) for the bulk mean velocity. The ratio of these two slopes is 0.85, which can be identified with the value of the turbulent Prandtl number in the overlap region. It is shown that the logarithmic $h^{+}$ increase of the integrated mean scalar is intrinsically associated with the overlap region of $\overline{\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}}$, established for $h^{+}$ (${\geqslant}400$). The resulting heat transfer law also holds at a smaller $h^{+}$ (${\geqslant}200$) than that derived by assuming a log law for the mean temperature.


2009 ◽  
Vol 23 (03) ◽  
pp. 509-512 ◽  
Author(s):  
SUHUA SHEN ◽  
JIANZHONG LIN

To explore the rheological property in turbulent channel flow of fiber suspensions, the equation of probability distribution function for mean fiber orientation and the Reynolds averaged Navier-Stokes equation with the term of additional stress resulted from fibers were solved with numerical methods to get the distributions of the mean velocity and turbulent kinetic energy. The simulation results show that the effect of fibers on turbulent channel flow is equivalent to an additional viscosity. The turbulent velocity profiles of fiber suspension become gradually sharper by increasing the fiber concentration and/or decreasing the Reynolds number. The turbulent kinetic energy will increase with increasing Reynolds number and fiber concentration.


1976 ◽  
Vol 77 (1) ◽  
pp. 153-174 ◽  
Author(s):  
Dietrich K. Lezius ◽  
James P. Johnston

The stability of laminar and turbulent channel flow is examined for cases where Coriolis forces are introduced by steady rotation about an axis perpendicular to the plane of mean flow. Linearized equations of motion are derived for small disturbances of the Taylor type. Conditions for marginal stability in laminar Couette and Poiseuille flow correspond, in part, to the analogous solutions of buoyancy-driven convection instabilities in heated fluid layers, and to those of Taylor instabilities in the flow between rotating cylinders. In plane Poiseuille flow with rotation, the critical disturbance mode occurs at a Reynolds number of Rec= 88.53 and rotation number Ro= 0.5. At higher Reynolds numbers, unstable conditions canexist over the range of rotation numbers given by 0 < Ro< 3, provided the undisturbed flow remains laminar. A two-layer model is devised to investigate the onset of longitudinal instabilities in turbulent flow. The linear disturbance equations are solved essentially in their laminar form, whereby the velocity gradient of laminar flow is replaced by a numerically computed profile for the gradient of the turbulent mean velocity. The turbulent stress levels in the stable and unstable flow regions are represented by integrated averages of the eddy viscosity. Onset of instability for Reynolds numbers between 6000 and 35 000 is predicted to occur at Ro= 0.022, a value in remarkable agreement with the experimentally observed appearance of roll instabilities in rotating turbulent channel flow.


Author(s):  
Saeed Jafari ◽  
Mohammad Rahnama

Generalized Lattice Boltzmann Equation (GLBE) was used for computation of turbulent channel flow for which Large Eddy Simulation (LES) was employed as a turbulence model. The subgrid-Scale turbulence effects were simulated through a Shear-Improved Smagorinsky Model (SISM) which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. Results of mean velocity distribution across the channel showed good correspondence with Direct Numerical Simulation (DNS) data. Negligible discrepancies were observed for computed turbulent statistics between present computations and those reported from DNS. Three-dimensional instantaneous vorticity contours showed complex vortical structures appeared in such flow geometries. It is concluded that such framework is capable of predicting accurate results for turbulent channel flow without adding significant complication and computational cost to the standard Smagorinsky model.


2013 ◽  
Vol 721 ◽  
pp. 627-668 ◽  
Author(s):  
Rémi Zamansky ◽  
Ivana Vinkovic ◽  
Mikhael Gorokhovski

AbstractThis paper focuses on the characterization and the stochastic modelling of the fluid acceleration in turbulent channel flow. In the first part, the acceleration is studied by direct numerical simulation (DNS) at three Reynolds numbers (${\mathit{Re}}_{\ast } = {u}_{\ast } h/ \nu = 180$, 590 and 1000). It is observed that whatever the wall distance is, the norm of acceleration is log-normally distributed and that the variance of the norm is very close to its mean value. It is also observed that from the wall to the centreline of the channel, the orientation of acceleration relaxes statistically towards isotropy. On the basis of dimensional analysis, a universal scaling law for the acceleration norm is proposed. In the second part, in the framework of the norm/orientation decomposition, a stochastic model of the acceleration is introduced. The stochastic model for the norm is based on fragmentation process which evolves across the channel with the wall distance. Simultaneously the orientation is simulated by a random walk on the surface of a unit sphere. The process is generated in such a way that the mean components of the orientation vector are equal to zero, whereas with increasing wall distance, all directions become equally probable. In the third part, the models are assessed in the framework of large-eddy simulation with stochastic subgrid acceleration model (LES–SSAM), introduced recently by Sabel’nikov, Chtab-Desportes & Gorokhovski (Euro. Phys. J. B, vol. 80, 2011, p. 177–187), and designed to account for the intermittency at subgrid scales. Computations by LES–SSAM and its assessment using DNS data show that the prediction of important statistics to characterize the flow, such as the mean velocity, the energy spectra at small scales, the viscous and turbulent stresses, the distribution of the acceleration can be considerably improved in comparison with standard LES. In the last part of this paper, the advantage of LES–SSAM in accounting for the subgrid flow structure is demonstrated in simulation of particle-laden turbulent channel flows. Compared to standard LES, it is shown that for different Stokes numbers, the particle dynamics and the turbophoresis effect can be predicted significantly better when LES–SSAM is applied.


Sign in / Sign up

Export Citation Format

Share Document