scholarly journals Mean dynamics of transitional channel flow

2011 ◽  
Vol 678 ◽  
pp. 451-481 ◽  
Author(s):  
J. ELSNAB ◽  
J. KLEWICKI ◽  
D. MAYNES ◽  
T. AMEEL
Keyword(s):  
Reynolds Number ◽  
Reynolds Stress ◽  
Channel Flow ◽  
Inflection Point ◽  
Stress Gradient ◽  
Development Stage ◽  
Stress Profile ◽  
Turbulent Regime ◽  
Nonlinear Development ◽  
The Mean

The redistribution of mean momentum and vorticity, along with the mechanisms underlying these redistribution processes, is explored for post-laminar flow in fully developed, pressure driven, channel flow. These flows, generically referred to as transitional, include an instability stage and a nonlinear development stage. The central focus is on the nonlinear development stage. The present analyses use existing direct numerical simulation data sets, as well as recently reported high-resolution molecular tagging velocimetry measurements. Primary considerations stem from the emergence of the effects of turbulent inertia as represented by the Reynolds stress gradient in the mean differential statement of dynamics. The results describe the flow evolution following the formation of a non-zero Reynolds stress peak that is known to first arise near the critical layer of the most unstable disturbance. The positive and negative peaks in the Reynolds stress gradient profile are observed to undergo a relative movement toward both the wall and centreline for subsequent increases in Reynolds number. The Reynolds stress profiles are shown to almost immediately exhibit the same sequence of curvatures that exists in the fully turbulent regime. In the transitional regime, the outer inflection point in this profile physically indicates a localized zone within which the mean dynamics are dominated by inertia. These observations connect to recent theoretical findings for the fully turbulent regime, e.g. as described by Fife, Klewicki & Wei (J. Discrete Continuous Dyn. Syst., vol. 24, 2009, p. 781) and Klewicki, Fife & Wei (J. Fluid Mech., vol. 638, 2009, p. 73). In accord with momentum equation analyses at higher Reynolds number, the present observations provide evidence that a logarithmic mean velocity profile is most rapidly approximated on a sub-domain located between the zero in the Reynolds stress gradient (maximum in the Reynolds stress) and the outer region location of the maximal Reynolds stress gradient (inflection point in the Reynolds stress profile). Overall, the present findings provide evidence that the dynamical processes during the post-laminar regime and those operative in the high Reynolds number regime are connected and describable within a single theoretical framework.

Mean dynamics and transition to turbulence in oscillatory channel flow

2019 ◽  
Vol 880 ◽  
pp. 864-889
Author(s):  
Alireza Ebadi ◽  
Christopher M. White ◽  
Ian Pond ◽  
Yves Dubief
Keyword(s):  
Reynolds Stress ◽  
Channel Flow ◽  
Layer Structure ◽  
Development Stage ◽  
Transition To Turbulence ◽  
Nonlinear Development ◽  
The Mean ◽  
Balance Of Forces ◽  
Turbulent Wall ◽  
Near Wall

The mean dynamics in oscillatory channel flow is examined to investigate the dynamical mechanisms underlying the transition to turbulence in oscillatory wall-bounded flow. The analyses employ direct numerical simulation data acquired at three Stokes Reynolds numbers: $Re_{s}=648$, 801 and 1009, where the lower $Re_{s}$ flow is transitional over the entire cycle and the two higher $Re_{s}$ flows exhibit flow characteristics similar to steady turbulent wall-bounded flow during part of the cycle. The flow evolution over a half-period of oscillation for all three $Re_{s}$ is as follows: near-wall streamwise velocity streaks develop during the early accelerating portion of the cycle; then at some later point in the cycle that depends on $Re_{s}$, the near-wall streaks breakdown (demarking the onset of the nonlinear development stage), and the near-wall Reynolds stress grows explosively; the Reynolds stress remains elevated for part of the cycle before diminishing (yet remaining finite) during the late decelerating portion of the cycle. This process is then repeated indefinitely. The present findings demonstrate that transition to turbulence occurs when the nonlinear development stage begins during the accelerating portion of the cycle. This crucially leads to the diminishing importance of the centreline momentum source, the emergence of a locally accelerating/decelerating internal layer centred about the edge of the Stokes layer and the wall-normal rearrangement of the mean forces prior to the start of the decelerating portion of the cycle. The rearrangement of mean forces culminates in a four layer structure in the mean balance of forces. This is significant on a number of accounts since empirical and theoretical evidence suggests that the formation of a four layer structure is an important characteristic of a self-similar hierarchal structure that underlies logarithmic dependence of the mean velocity profile in steady turbulent wall-bounded flows (Klewicki et al., J. Fluid Mech., vol. 638, 2009, pp. 73–93). When the nonlinear development stage begins during the decelerating portion of the cycle (i.e. at $Re_{s}=648$), a four layer structure is not observed in the mean balance of forces and the flow remains weakly transitional over the entire cycle.

scholarly journals Mean dynamics of transitional boundary-layer flow

2011 ◽  
Vol 682 ◽  
pp. 617-651 ◽  
Author(s):  
J. KLEWICKI ◽  
R. EBNER ◽  
X. WU
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Reynolds Stress ◽  
Boundary Layer Flow ◽  
Layer Structure ◽  
Momentum Equation ◽  
Stress Profile ◽  
Transitional Regime ◽  
The Mean ◽  
Layer Flow

The dynamical mechanisms underlying the redistribution of mean momentum and vorticity are explored for transitional two-dimensional boundary-layer flow at nominally zero pressure gradient. The analyses primarily employ the direct numerical simulation database of Wu & Moin (J. Fluid Mech., vol. 630, 2009, p. 5), but are supplemented with verifications utilizing subsequent similar simulations. The transitional regime is taken to include both an instability stage, which effectively generates a finite Reynolds stress profile, −ρuv(y), and a nonlinear development stage, which progresses until the terms in the mean momentum equation attain the magnitude ordering of the four-layer structure revealed by Wei et al. (J. Fluid Mech., vol. 522, 2005, p. 303). Self-consistently applied criteria reveal that the third layer of this structure forms first, followed by layers IV and then II and I. For the present flows, the four-layer structure is estimated to be first realized at a momentum thickness Reynolds number Rθ = U∞ θ/ν ≃ 780. The first-principles-based theory of Fife et al. (J. Disc. Cont. Dyn. Syst. A, vol. 24, 2009, p. 781) is used to describe the mean dynamics in the laminar, transitional and four-layer regimes. As in channel flow, the transitional regime is marked by a non-negligible influence of all three terms in the mean momentum equation at essentially all positions in the boundary layer. During the transitional regime, the action of the Reynolds stress gradient rearranges the mean viscous force and mean advection profiles. This culminates with the segregation of forces characteristic of the four-layer regime. Empirical and theoretical evidence suggests that the formation of the four-layer structure also underlies the emergence of the mean dynamical properties characteristic of the high-Reynolds-number flow. These pertain to why and where the mean velocity profile increasingly exhibits logarithmic behaviour, and how and why the Reynolds stress distribution develops such that the inner normalized position of its peak value, ym+, exhibits a Reynolds number dependence according to $y_m^+ {\,\simeq\,} 1.9 \sqrt{\delta^+}$.

scholarly journals Passive scalars in turbulent channel flow at high Reynolds number

2016 ◽  
Vol 788 ◽  
pp. 614-639 ◽  
Author(s):  
Sergio Pirozzoli ◽  
Matteo Bernardini ◽  
Paolo Orlandi
Keyword(s):  
Reynolds Number ◽  
Channel Flow ◽  
High Reynolds Number ◽  
Correction Term ◽  
Scalar Dissipation ◽  
Channel Height ◽  
Passive Scalars ◽  
The Core ◽  
The Mean ◽  
High Reynolds

We study passive scalars in turbulent plane channels at computationally high Reynolds number, thus allowing us to observe previously unnoticed effects. The mean scalar profiles are found to obey a generalized logarithmic law which includes a linear correction term in the whole lower half-channel, and they follow a universal parabolic defect profile in the core region. This is consistent with recent findings regarding the mean velocity profiles in channel flow. The scalar variances also exhibit a near universal parabolic distribution in the core flow and hints of a sizeable log layer, unlike the velocity variances. The energy spectra highlight the formation of large scalar-bearing eddies with size proportional to the channel height which are caused by a local production excess over dissipation, and which are clearly visible in the flow visualizations. Close correspondence of the momentum and scalar eddies is observed, with the main difference being that the latter tend to form sharper gradients, which translates into higher scalar dissipation. Another notable Reynolds number effect is the decreased correlation of the passive scalar field with the vertical velocity field, which is traced to the reduced effectiveness of ejection events.

Properties of the mean recirculation region in the wakes of two-dimensional bluff bodies

1997 ◽  
Vol 351 ◽  
pp. 167-199 ◽  
Author(s):  
S. BALACHANDAR ◽  
R. MITTAL ◽  
F. M. NAJJAR
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Reynolds Stress ◽  
Vortex Shedding ◽  
Reynolds Numbers ◽  
Recirculation Region ◽  
Bluff Bodies ◽  
Two Dimensional ◽  
Stress Components ◽  
The Mean

The properties of the time- and span-averaged mean wake recirculation region are investigated in separated flows over several different two-dimensional bluff bodies. Ten different cases are considered and they divide into two groups: cylindrical geometries of circular, elliptic and square cross-sections and the normal plate. A wide Reynolds number range from 250 to 140000 is considered, but in all the cases the attached portion of the boundary layer remains laminar until separation. The lower Reynolds number data are from direct numerical simulations, while the data at the higher Reynolds number are obtained from large-eddy simulation and the experimental work of Cantwell & Coles (1983), Krothapalli (1996, personal communication), Leder (1991) and Lyn et al. (1995). Unlike supersonic and subsonic separations with a splitter plate in the wake, in all the cases considered here there is strong interaction between the shear layers resulting in Kármán vortex shedding. The impact of this fundamental difference on the distribution of Reynolds stress components and pressure in relation to the mean wake recirculation region (wake bubble) is considered. It is observed that in all cases the contribution from Reynolds normal stress to the force balance of the wake bubble is significant. In fact, in the cylinder geometries this contribution can outweigh the net force from the shear stress, so that the net pressure force tends to push the bubble away from the body. In contrast, in the case of normal plate, owing to the longer wake, the net contribution from shear stress outweighs that from the normal stress. At higher Reynolds numbers, separation of the Reynolds stress components into incoherent contributions provides more insight. The behaviour of the coherent contribution, arising from the dominant vortex shedding, is similar to that at lower Reynolds numbers. The incoherent contribution to Reynolds stress, arising from small-scale activity, is compared with that of a canonical free shear layer. Based on these observations a simple extension of the wake model (Sychev 1982; Roshko 1993a, b) is proposed.

scholarly journals Effect of turbulent fluctuations on the drag and lift forces on a towed sphere and its boundary layer

2013 ◽  
Vol 721 ◽  
pp. 155-179 ◽  
Author(s):  
Holger Homann ◽  
Jérémie Bec ◽  
Rainer Grauer
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Spherical Particle ◽  
Drag Force ◽  
Turbulent Regime ◽  
Mean Velocity ◽  
Turbulent Fluctuations ◽  
The Mean ◽  
Lift Forces ◽  
Drag And Lift

AbstractThe impact of turbulent fluctuations on the forces exerted by a fluid on a towed spherical particle is investigated by means of high-resolution direct numerical simulations. The measurements are carried out using a novel scheme to integrate the two-way coupling between the particle and the incompressible surrounding fluid flow maintained in a high-Reynolds-number turbulent regime. The main idea consists of combining a Fourier pseudo-spectral method for the fluid with an immersed-boundary technique to impose the no-slip boundary condition on the surface of the particle. This scheme is shown to converge as the power $3/ 2$ of the spatial resolution. This behaviour is explained by the ${L}_{2} $ convergence of the Fourier representation of a velocity field displaying discontinuities of its derivative. Benchmarking of the code is performed by measuring the drag and lift coefficients and the torque-free rotation rate of a spherical particle in various configurations of an upstream-laminar carrier flow. Such studies show a good agreement with experimental and numerical measurements from other groups. A study of the turbulent wake downstream of the sphere is also reported. The mean velocity deficit is shown to behave as the inverse of the distance from the particle, as predicted from classical similarity analysis. This law is reinterpreted in terms of the principle of ‘permanence of large eddies’ that relates infrared asymptotic self-similarity to the law of decay of energy in homogeneous turbulence. The developed method is then used to attack the problem of an upstream flow that is in a developed turbulent regime. It is shown that the average drag force increases as a function of the turbulent intensity and the particle Reynolds number. This increase is significantly larger than predicted by standard drag correlations based on laminar upstream flows. It is found that the relevant parameter is the ratio of the viscous boundary layer thickness to the dissipation scale of the ambient turbulent flow. The drag enhancement can be motivated by the modification of the mean velocity and pressure profile around the sphere by small-scale turbulent fluctuations. It is demonstrated that the variance of the drag force fluctuations can be modelled by means of standard drag correlations. Temporal correlations of the drag and lift forces are also presented.

Stability bounds on turbulent Poiseuille flow

1988 ◽  
Vol 187 ◽  
pp. 435-449 ◽  
Author(s):  
G. R. Ierley ◽  
W. V. R. Malkus
Keyword(s):  
Reynolds Number ◽  
Reynolds Stress ◽  
Poiseuille Flow ◽  
Numerical Solutions ◽  
Numerical Study ◽  
Linear Instability ◽  
Mean Flow ◽  
Starting Point ◽  
The Mean ◽  
The Stability

For steady-state turbulent flows with unique mean properties, we determine a sense in which the mean velocity is linearly supercritical. The shear-turbulence literature on this point is ambiguous. As an example, we reassess the stability of mean profiles in turbulent Poiseuille flow. The Reynolds & Tiederman (1967) numerical study is used as a starting point. They had constructed a class of one-dimensional flows which included, within experimental error, the observed profile. Their numerical solutions of the resulting Orr-Sommerfeld problems led them to conclude that the Reynolds number for neutral infinitesimal disturbances was twenty-five times the Reynolds number characterizing the observed mean flow. They found also that the first nonlinear corrections were stabilizing. In the realized flow, this latter conclusion appears incompatible with the former. Hence, we have sought a more complete set of velocity profiles which could exhibit linear instability, retaining the requirement that the observed velocity profile is included in the set. We have added two dynamically generated modifications of the mean. The first addition is a fluctuation in the curvature of the mean flow generated by a Reynolds stress whose form is determined by the neutrally stable Orr-Sommerfeld solution. We find that this can reduce the stability of the observed flow by as much as a factor of two. The second addition is the zero-average downstream wave associated with the above Reynolds stress. The three-dimensional linear instability of this modification can even render the observed flow unstable. Those wave amplitudes that just barely will ensure instability of the observed flow are determined. The relation of these particular amplitudes to the limiting conditions admitted by an absolute stability criterion for disturbances on the mean flow is found. These quantitative results from stability theory lie in the observationally determined Reynolds-Tiederman similarity scheme, and hence are insensitive to changes in Reynolds number.

scholarly journals Turbulent channel flow of an elastoviscoplastic fluid

2018 ◽  
Vol 853 ◽  
pp. 488-514 ◽  
Author(s):  
Marco E. Rosti ◽  
Daulet Izbassarov ◽  
Outi Tammisola ◽  
Sarah Hormozi ◽  
Luca Brandt
Keyword(s):  
Reynolds Number ◽  
Turbulent Flow ◽  
Friction Factor ◽  
Channel Flow ◽  
Viscosity Ratio ◽  
Turbulent Channel Flow ◽  
Turbulent Regime ◽  
Bingham Number ◽  
Wide Range ◽  
Near Wall

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.

scholarly journals Statistics and structure of spanwise rotating turbulent channel flow at moderate Reynolds numbers

2017 ◽  
Vol 828 ◽  
pp. 424-458 ◽  
Author(s):  
Geert Brethouwer
Keyword(s):  
Reynolds Number ◽  
Channel Flow ◽  
Large Scale ◽  
Linear Part ◽  
Turbulent Channel Flow ◽  
Streamwise Velocity ◽  
Reynolds Stresses ◽  
System Rotation ◽  
The Mean ◽  
Spanwise Rotation

A study of fully developed plane turbulent channel flow subject to spanwise system rotation through direct numerical simulations is presented. In order to study both the influence of the Reynolds number and spanwise rotation on channel flow, the Reynolds number $Re=U_{b}h/\unicode[STIX]{x1D708}$ is varied from a low 3000 to a moderate 31 600 and the rotation number $Ro=2\unicode[STIX]{x1D6FA}h/U_{b}$ is varied from 0 to 2.7, where $U_{b}$ is the mean bulk velocity, $h$ the channel half-gap, $\unicode[STIX]{x1D708}$ the viscosity and $\unicode[STIX]{x1D6FA}$ the system rotation rate. The mean streamwise velocity profile displays also at higher $Re$ a characteristic linear part with a slope near to $2\unicode[STIX]{x1D6FA}$, and a corresponding linear part in the profiles of the production and dissipation rate of turbulent kinetic energy appears. With increasing $Ro$, a distinct unstable side with large spanwise and wall-normal Reynolds stresses and a stable side with much weaker turbulence develops in the channel. The flow starts to relaminarize on the stable side of the channel and persisting turbulent–laminar patterns appear at higher $Re$. If $Ro$ is further increased, the flow on the stable side becomes laminar-like while at yet higher $Ro$ the whole flow relaminarizes, although the calm periods might be disrupted by repeating bursts of turbulence, as explained by Brethouwer (Phys. Rev. Fluids, vol. 1, 2016, 054404). The influence of the Reynolds number is considerable, in particular on the stable side of the channel where velocity fluctuations are stronger and the flow relaminarizes less quickly at higher $Re$. Visualizations and statistics show that, at $Ro=0.15$ and 0.45, large-scale structures and large counter-rotating streamwise roll cells develop on the unstable side. These become less noticeable and eventually vanish when $Ro$ rises, especially at higher $Re$. At high $Ro$, the largest energetic structures are larger at lower $Re$.

Explicit algebraic subgrid stress models with application to rotating channel flow

2009 ◽  
Vol 639 ◽  
pp. 403-432 ◽  
Author(s):  
LINUS MARSTORP ◽  
GEERT BRETHOUWER ◽  
OLOF GRUNDESTAM ◽  
ARNE V. JOHANSSON
Keyword(s):  
Dynamic Model ◽  
Reynolds Stress ◽  
Channel Flow ◽  
Stress Model ◽  
Smagorinsky Model ◽  
Mean Velocity ◽  
The Mean ◽  
Stress Models ◽  
Subgrid Stress ◽  
Rotating Channel

New explicit subgrid stress models are proposed involving the strain rate and rotation rate tensor, which can account for rotation in a natural way. The new models are based on the same methodology that leads to the explicit algebraic Reynolds stress model formulation for Reynolds-averaged Navier–Stokes simulations. One dynamic model and one non-dynamic model are proposed. The non-dynamic model represents a computationally efficient subgrid scale (SGS) stress model, whereas the dynamic model is the most accurate. The models are validated through large eddy simulations (LESs) of spanwise and streamwise rotating channel flow and are compared with the standard and dynamic Smagorinsky models. The proposed explicit dependence on the system rotation improves the description of the mean velocity profiles and the turbulent kinetic energy at high rotation rates. Comparison with the dynamic Smagorinsky model shows that not using the eddy-viscosity assumption improves the description of both the Reynolds stress anisotropy and the SGS stress anisotropy. LESs of rotating channel flow at Reτ = 950 have been carried out as well. These reveal some significant Reynolds number influences on the turbulence statistics. LESs of non-rotating turbulent channel flow at Reτ = 950 show that the new explicit model especially at coarse resolutions significantly better predicts the mean velocity, wall shear and Reynolds stresses than the dynamic Smagorinsky model, which is probably the result of a better prediction of the anisotropy of the subgrid dissipation.

Sign in / Sign up

Export Citation Format

Share Document