Growth of vortical disturbances entrained in the entrance region of a circular pipe

2021 ◽  
Vol 932 ◽  
Author(s):  
Pierre Ricco ◽  
Claudia Alvarenga

The development and growth of unsteady three-dimensional vortical disturbances entrained in the entry region of a circular pipe is investigated by asymptotic and numerical methods for Reynolds numbers between $1000$ and $10\,000$ , based on the pipe radius and the bulk velocity. Near the pipe mouth, composite asymptotic solutions describe the dynamics of the oncoming disturbances, revealing how these disturbances are altered by the viscous layer attached to the pipe wall. The perturbation velocity profiles near the pipe mouth are employed as rigorous initial conditions for the boundary-region equations, which describe the flow in the limit of low frequency and large Reynolds number. The disturbance flow is initially primarily present within the base-flow boundary layer in the form of streamwise-elongated vortical structures, i.e. the streamwise velocity component displays an intense algebraic growth, while the cross-flow velocity components decay. Farther downstream the disturbance flow occupies the whole pipe, although the base flow is mostly inviscid in the core. The transient growth and subsequent viscous decay are confined in the entrance region, i.e. where the base flow has not reached the fully developed Poiseuille profile. Increasing the Reynolds number and decreasing the frequency causes more intense perturbations, whereas small azimuthal wavelengths and radial characteristic length scales intensify the viscous dissipation of the disturbance. The azimuthal wavelength that causes the maximum growth is found. The velocity profiles are compared successfully with available experimental data and the theoretical results are helpful to interpret the only direct numerical dataset of a disturbed pipe-entry flow.

1973 ◽  
Vol 59 (3) ◽  
pp. 537-559 ◽  
Author(s):  
C. W. Van Atta ◽  
T. T. Yeh

Some of the statistical characteristics of the breakdown coefficient, defined as the ratio of averages over different spatial regions of positive variables characterizing the fine structure and internal intermittency in high Reynolds number turbulence, have been investigated using experimental data for the streamwise velocity derivative ∂u/∂tmeasured in an atmospheric boundary layer. The assumptions and predictions of the hypothesis of scale similarity developed by Novikov and by Gurvich & Yaglom do not adequately describe or predict the statistical characteristics of the breakdown coefficientqr,lof the square of the streamwise velocity derivative. Systematic variations in the measured probability densities and consistent variations in the measured moments show that the assumption that the probability density of the breakdown coefficient is a function only of the scale ratio is not satisfied. The small positive correlation between adjoint values ofqr,land measurements of higher moments indicate that the assumption that the probability densities for adjoint values ofqr,lare statistically independent is also not satisfied. The moments ofqr,ldo not have the simple power-law character that is a consequence of scale similarity.As the scale ratiol/rchanges, the probability density ofqr,levolves from a sharply peaked, highly negatively skewed density for large values of the scale ratio to a very symmetrical distribution when the scale ratio is equal to two, and then to a highly positively skewed density as the scale ratio approaches one. There is a considerable effect of heterogeneity on the values of the higher moments, and a small but measurable effect on the mean value. The moments are roughly symmetrical functions of the displacement of the shorter segment from the centre of the larger one, with a minimum value when the shorter segment is centrally located within the larger one.


Author(s):  
Nariman Ashrafi

The nonlinear stability and bifurcation of the one-dimensional channel (Poiseuille) flow is examined for a Johnson-Segalman fluid. The velocity and stress are represented by orthonormal functions in the transverse direction to the flow. The flow field is obtained from the conservation and constitutive equations using the Galerkin projection method. Both inertia and normal stress effects are included. The stability picture is dramatically influenced by the viscosity ratio. The range of shear rate or Weissenberg number for which the base flow is unstable increases from zero as the fluid deviates from the Newtonian limit as decreases. Typically, two turning points are observed near the critical Weissenberg numbers. The transient response is heavily influenced by the level of inertia. It is found that the flow responds oscillatorily. When the Reynolds number is small, and monotonically at large Reynolds number when elastic effects are dominated by inertia.


1999 ◽  
Author(s):  
Hidesada Kanda

Abstract A conceptual model was constructed for the problem of determining in circular pipes the conditions under which the transition from laminar to turbulent flow occurs, so that it becomes possible to calculate the minimum critical Reynolds number. Up until now this problem has been investigated by stability theory with disturbances at the pipe inlet. However, the minimum critical Reynolds number has not yet been obtained theoretically. Hence, the author took up the problem directly from many previous experimental investigations and found that (i) plots of the transition length versus the Reynolds number show that the transition occurs in the entrance region under the condition of a natural disturbance, and (ii) plots of the critical Reynolds number versus the ratio of bellmouth diameter to the pipe diamter show that with larger shapes of bellmouths, laminar flow will persist to higher Reynolds numbers. The problem is thus defined clearly as: Under the condition of an ordinary disturbance, the transition from laminar to turbulent flow occurs in the entrance region of a straight circular pipe, then the Reynolds number takes a minimum value of about 2000.


2021 ◽  
Vol 927 ◽  
Author(s):  
Pierre Ricco ◽  
Claudia Alvarenga

The entrainment of free-stream unsteady three-dimensional vortical disturbances in the entry region of a channel is studied via matched asymptotic expansions and by numerical means. The interest is in flows at Reynolds numbers where experimental studies have documented the occurrence of intense transient growth, despite the flow being stable according to classical stability analysis. The analytical description of the vortical perturbations at the channel mouth reveals how the oncoming disturbances penetrate into the wall-attached shear layers and amplify downstream. The effects of the channel confinement, the streamwise pressure gradient and the viscous/inviscid interplay between the oncoming disturbances and the boundary-layer perturbations are discussed. The composite perturbation velocity profiles are employed as initial conditions for the unsteady boundary-region perturbation equations. At a short distance from the channel mouth, the disturbance flow is mostly confined within the shear layers and assumes the form of streamwise-elongated streaks, while farther downstream the viscous disturbances permeate the whole channel although the base flow is still mostly inviscid in the core. Symmetrical disturbances exhibit a more significant growth than anti-symmetrical disturbances, the latter maintaining a nearly constant amplitude for several channel heights downstream before growing transiently, a unique feature not reported in open boundary layers. The disturbances are more intense as the frequency decreases or the bulk Reynolds number increases. We compute the spanwise wavelengths that cause the most intense downstream growth and the threshold wall-normal wavelengths below which the perturbations are damped through viscous dissipation.


2019 ◽  
Vol 864 ◽  
pp. 5-44 ◽  
Author(s):  
David Fabre ◽  
Raffaele Longobardi ◽  
Paul Bonnefis ◽  
Paolo Luchini

The unsteady axisymmetric flow through a circular aperture in a thin plate subjected to harmonic forcing (for instance under the effect of an incident acoustic wave) is a classical problem first considered by Howe (Proc. R. Soc. Lond. A, vol. 366, 1979, pp. 205–223), using an inviscid model. The purpose of this work is to reconsider this problem through a numerical resolution of the incompressible linearized Navier–Stokes equations (LNSE) in the laminar regime, corresponding to $Re=[500,5000]$. We first compute a steady base flow which allows us to describe the vena contracta phenomenon in agreement with experiments. We then solve a linear problem allowing us to characterize both the spatial amplification of the perturbations and the impedance (or equivalently the Rayleigh conductivity), which is a key quantity to investigate the response of the jet to acoustic forcing. Since the linear perturbation is characterized by a strong spatial amplification, the numerical resolution requires the use of a complex mapping of the axial coordinate in order to enlarge the range of Reynolds number investigated. The results show that the impedances computed with $Re\gtrsim 1500$ collapse onto a single curve, indicating that a large Reynolds number asymptotic regime is effectively reached. However, expressing the results in terms of conductivity leads to substantial deviation with respect to Howe’s model. Finally, we investigate the case of finite-amplitude perturbations through direct numerical simulations (DNS). We show that the impedance predicted by the linear approach remains valid for amplitudes up to order $10^{-1}$, despite the fact that the spatial evolution of the perturbations in the jet is strongly nonlinear.


2014 ◽  
Vol 747 ◽  
pp. 518-544 ◽  
Author(s):  
Jan Östh ◽  
Bernd R. Noack ◽  
Siniša Krajnović ◽  
Diogo Barros ◽  
Jacques Borée

AbstractWe investigate a hierarchy of eddy-viscosity terms in proper orthogonal decomposition (POD) Galerkin models to account for a large fraction of unresolved fluctuation energy. These Galerkin methods are applied to large eddy simulation (LES) data for a flow around a vehicle-like bluff body called an Ahmed body. This flow has three challenges for any reduced-order model: a high Reynolds number, coherent structures with broadband frequency dynamics, and meta-stable asymmetric base flow states. The Galerkin models are found to be most accurate with modal eddy viscosities as proposed by Rempfer & Fasel (J. Fluid Mech., vol. 260, 1994a, pp. 351–375; J. Fluid Mech. vol. 275, 1994b, pp. 257–283). Robustness of the model solution with respect to initial conditions, eddy-viscosity values and model order is achieved only for state-dependent eddy viscosities as proposed by Noack, Morzyński & Tadmor (Reduced-Order Modelling for Flow Control, CISM Courses and Lectures, vol. 528, 2011). Only the POD system with state-dependent modal eddy viscosities can address all challenges of the flow characteristics. All parameters are analytically derived from the Navier–Stokes-based balance equations with the available data. We arrive at simple general guidelines for robust and accurate POD models which can be expected to hold for a large class of turbulent flows.


2015 ◽  
Vol 783 ◽  
pp. 379-411 ◽  
Author(s):  
I. Marusic ◽  
K. A. Chauhan ◽  
V. Kulandaivelu ◽  
N. Hutchins

In this paper we study the spatial evolution of zero-pressure-gradient (ZPG) turbulent boundary layers from their origin to a canonical high-Reynolds-number state. A prime motivation is to better understand under what conditions reliable scaling behaviour comparisons can be made between different experimental studies at matched local Reynolds numbers. This is achieved here through detailed streamwise velocity measurements using hot wires in the large University of Melbourne wind tunnel. By keeping the unit Reynolds number constant, the flow conditioning, contraction and trip can be considered unaltered for a given boundary layer’s development and hence its evolution can be studied in isolation from the influence of inflow conditions by moving to different streamwise locations. Careful attention was given to the experimental design in order to make comparisons between flows with three different trips while keeping all other parameters nominally constant, including keeping the measurement sensor size nominally fixed in viscous wall units. The three trips consist of a standard trip and two deliberately ‘over-tripped’ cases, where the initial boundary layers are over-stimulated with additional large-scale energy. Comparisons of the mean flow, normal Reynolds stress, spectra and higher-order turbulence statistics reveal that the effects of the trip are seen to be significant, with the remnants of the ‘over-tripped’ conditions persisting at least until streamwise stations corresponding to $Re_{x}=1.7\times 10^{7}$ and $x=O(2000)$ trip heights are reached (which is specific to the trips used here), at which position the non-canonical boundary layers exhibit a weak memory of their initial conditions at the largest scales $O(10{\it\delta})$, where ${\it\delta}$ is the boundary layer thickness. At closer streamwise stations, no one-to-one correspondence is observed between the local Reynolds numbers ($Re_{{\it\tau}}$, $Re_{{\it\theta}}$ or $Re_{x}$ etc.), and these differences are likely to be the cause of disparities between previous studies where a given Reynolds number is matched but without account of the trip conditions and the actual evolution of the boundary layer. In previous literature such variations have commonly been referred to as low-Reynolds-number effects, while here we show that it is more likely that these differences are due to an evolution effect resulting from the initial conditions set up by the trip and/or the initial inflow conditions. Generally, the mean velocity profiles were found to approach a constant wake parameter ${\it\Pi}$ as the three boundary layers developed along the test section, and agreement of the mean flow parameters was found to coincide with the location where other statistics also converged, including higher-order moments up to tenth order. This result therefore implies that it may be sufficient to document the mean flow parameters alone in order to ascertain whether the ZPG flow, as described by the streamwise velocity statistics, has reached a canonical state, and a computational approach is outlined to do this. The computational scheme is shown to agree well with available experimental data.


2011 ◽  
Vol 684 ◽  
pp. 204-226 ◽  
Author(s):  
Alejandro Sevilla

AbstractThe linear spatiotemporal stability properties of axisymmetric laminar capillary jets with fully developed initial velocity profiles are studied for large values of both the Reynolds number, $\mathit{Re}= Q/ (\lrm{\pi} a\nu )$, and the Froude number, $\mathit{Fr}= {Q}^{2} / ({\lrm{\pi} }^{2} g{a}^{5} )$, where $a$ is the injector radius, $Q$ the volume flow rate, $\nu $ the kinematic viscosity and $g$ the gravitational acceleration. The downstream development of the basic flow and its stability are addressed with an approximate formulation that takes advantage of the jet slenderness. The base flow is seen to depend on two parameters, namely a Stokes number, $G= \mathit{Re}/ \mathit{Fr}$, and a Weber number, $\mathit{We}= \rho {Q}^{2} / ({\lrm{\pi} }^{2} \sigma {a}^{3} )$, where $\sigma $ is the surface tension coefficient, while its linear stability depends also on the Reynolds number. When non-parallel terms are retained in the local stability problem, the analysis predicts a critical value of the Weber number, ${\mathit{We}}_{c} (G, \mathit{Re})$, below which a pocket of local absolute instability exists within the near field of the jet. The function ${\mathit{We}}_{c} (\mathit{Re})$ is computed for the buoyancy-free jet, showing marked differences with the results previously obtained with uniform velocity profiles. It is seen that, in accounting for gravity effects, it is more convenient to express the parametric dependence of the critical Weber number with use made of the Morton and Bond numbers, $\mathit{Mo}= {\nu }^{4} {\rho }^{3} g/ {\sigma }^{3} $ and $\mathit{Bo}= \rho g{a}^{2} / \sigma $, as replacements for $G$ and $\mathit{Re}$. This alternative formulation is advantageous to describe jets of a given liquid for a known value of $g$, in that the resulting Morton number becomes constant, thereby leaving $\mathit{Bo}$ as the only relevant parameter. The computed function ${\mathit{We}}_{c} (\mathit{Bo})$ for a water jet under Earth gravity is shown to be consistent with the experimental results of Clanet and Lasheras for the transition from jetting to dripping of water jets discharging into air from long injection needles, which cannot be properly described with a uniform velocity profile assumed at the jet exit.


2013 ◽  
Vol 80 (5) ◽  
Author(s):  
Zhao C. Kong ◽  
Duncan O. Eddy ◽  
Nathan K. Martin ◽  
Brent C. Houchens

The steady, axisymmetric base flow and instabilities in a rotating lid-driven cylinder are investigated experimentally via ultrasonic Doppler velocimetry and verified with computations. The flow is governed by two parameters: the Reynolds number (based on the angular velocity of the top lid, the cylinder radius, and kinematic viscosity) and the aspect ratio (cylinder height/radius). Base states and instabilities are explored using ultrasonic Doppler velocimetry in two mixtures of glycerol and water. Velocity profiles in the cylinder are constructed for aspect ratio 2.5 and Reynolds numbers between 1000 and 3000. The results are compared to computational spectral element simulations, as well as previously published findings. The base flow velocity profiles measured by ultrasonic Doppler velocimetry are in good agreement with the numerical results below the critical Reynolds number. The same is true for time-averaged results above the critical Reynolds number. Prediction of the first axisymmetric instability is demonstrated, although not always at the expected critical Reynolds number. Advantages and limitations of ultrasonic Doppler velocimetry are discussed.


Sign in / Sign up

Export Citation Format

Share Document