scholarly journals Suspensions of finite-size neutrally buoyant spheres in turbulent duct flow

2018 ◽  
Vol 851 ◽  
pp. 148-186 ◽  
Author(s):  
Walter Fornari ◽  
Hamid Tabaei Kazerooni ◽  
Jeanette Hussong ◽  
Luca Brandt
Keyword(s):  
Reynolds Number ◽  
Slip Velocity ◽  
Volume Fraction ◽  
Secondary Flows ◽  
Duct Flow ◽  
Turbulent Kinetic Energy Budget ◽  
Square Duct ◽  
The Core ◽  
The Mean ◽  
Neutrally Buoyant

We study the turbulent square duct flow of dense suspensions of neutrally buoyant spherical particles. Direct numerical simulations (DNS) are performed in the range of volume fractions $\unicode[STIX]{x1D719}=0{-}0.2$, using the immersed boundary method (IBM) to account for the dispersed phase. Based on the hydraulic diameter a Reynolds number of 5600 is considered. We observe that for $\unicode[STIX]{x1D719}=0.05$ and 0.1, particles preferentially accumulate on the corner bisectors, close to the corners, as also observed for laminar square duct flows of the same duct-to-particle size ratio. At the highest volume fraction, particles preferentially accumulate in the core region. For plane channel flows, in the absence of lateral confinement, particles are found instead to be uniformly distributed across the channel. The intensity of the cross-stream secondary flows increases (with respect to the unladen case) with the volume fraction up to $\unicode[STIX]{x1D719}=0.1$, as a consequence of the high concentration of particles along the corner bisector. For $\unicode[STIX]{x1D719}=0.2$ the turbulence activity is reduced and the intensity of the secondary flows reduces to below that of the unladen case. The friction Reynolds number increases with $\unicode[STIX]{x1D719}$ in dilute conditions, as observed for channel flows. However, for $\unicode[STIX]{x1D719}=0.2$ the mean friction Reynolds number is similar to that for $\unicode[STIX]{x1D719}=0.1$. By performing the turbulent kinetic energy budget, we see that the turbulence production is enhanced up to $\unicode[STIX]{x1D719}=0.1$, while for $\unicode[STIX]{x1D719}=0.2$ the production decreases below the values for $\unicode[STIX]{x1D719}=0.05$. On the other hand, the dissipation and the transport monotonically increase with $\unicode[STIX]{x1D719}$. The interphase interaction term also contributes positively to the turbulent kinetic energy budget and increases monotonically with $\unicode[STIX]{x1D719}$, in a similar way as the mean transport. Finally, we show that particles move on average faster than the fluid. However, there are regions close to the walls and at the corners where they lag behind it. In particular, for $\unicode[STIX]{x1D719}=0.05,0.1$, the slip velocity distribution at the corner bisectors seems correlated to the locations of maximum concentration: the concentration is higher where the slip velocity vanishes. The wall-normal hydrodynamic and collision forces acting on the particles push them away from the corners. The combination of these forces vanishes around the locations of maximum concentration. The total mean forces are generally low along the corner bisectors and at the core, also explaining the concentration distribution for $\unicode[STIX]{x1D719}=0.2$.

scholarly journals Experimental investigation of turbulent suspensions of spherical particles in a square duct

2018 ◽  
Vol 857 ◽  
pp. 748-783 ◽  
Author(s):  
Sagar Zade ◽  
Pedro Costa ◽  
Walter Fornari ◽  
Fredrik Lundell ◽  
Luca Brandt
Keyword(s):  
Particle Size ◽  
Reynolds Number ◽  
Pressure Drop ◽  
Single Phase ◽  
Volume Fraction ◽  
Particle Diameter ◽  
Spherical Particles ◽  
Particle Sizes ◽  
Duct Flow ◽  
Square Duct

We report experimental observations of turbulent flow with spherical particles in a square duct. Three particle sizes, namely $2H/d_{p}=40$ , 16 and 9 ( $2H$ being the duct full height and $d_{p}$ being the particle diameter), are investigated. The particles are nearly neutrally buoyant with a density ratio of 1.0035 and 1.01 with respect to the suspending fluid. Refractive index matched–particle image velocimetry (RIM–PIV) is used for fluid velocity measurement even at the highest particle volume fraction (20 %) and particle tracking velocimetry (PTV) for the particle velocity statistics for the flows seeded with particles of the two largest sizes, whereas only pressure measurements are reported for the smallest particles. Settling effects are seen at the lowest bulk Reynolds number $Re_{2H}\approx$ 10 000, whereas, at the highest $Re_{2H}\approx 27\,000$ , particles are in almost full suspension. The friction factor of the suspensions is found to be significantly larger than that of single-phase duct flow at the lower $Re_{2H}$ investigated; however, the difference decreases when increasing the flow rate and the total drag approaches the values of the single-phase flow at the higher Reynolds number considered, $Re_{2H}=27\,000$ . The pressure drop is found to decrease with the particle diameter for volume fractions lower than $\unicode[STIX]{x1D719}=10\,\%$ for nearly all $Re_{2H}$ investigated. However, at the highest volume fraction $\unicode[STIX]{x1D719}=20\,\%$ , we report a peculiar non-monotonic behaviour: the pressure drop first decreases and then increases with increasing particle size. The decrease of the turbulent drag with particle size at the lowest volume fractions is related to an attenuation of the turbulence. The drag increase for the two largest particle sizes at $\unicode[STIX]{x1D719}=20\,\%$ , however, occurs despite this large reduction of the turbulent stresses, and it is therefore due to significant particle-induced stresses. At the lowest Reynolds number, the particles reside mostly in the bottom half of the duct, where the mean velocity significantly decreases; the flow is similar to that in a moving porous bed near the bottom wall and to turbulent duct flow with low particle concentration near the top wall.

Direct numerical simulation of turbulent duct flow with inclined or V-shaped ribs mounted on one wall

2021 ◽  
Vol 932 ◽  
Author(s):  
S.V. Mahmoodi-Jezeh ◽  
Bing-Chen Wang
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Secondary Flows ◽  
Wall Turbulence ◽  
Spanwise Direction ◽  
Duct Flow ◽  
Mean Flow ◽  
Square Duct ◽  
Turbulence Structures ◽  
The Mean

In this research, highly disturbed turbulent flow of distinct three-dimensional characteristics in a square duct with inclined or V-shaped ribs mounted on one wall is investigated using direct numerical simulation. The turbulence field is highly sensitive to not only the rib geometry but also the boundary layers developed over the side and top walls. In a cross-stream plane secondary flows appear as large longitudinal vortices in both inclined and V-shaped rib cases due to the confinement of four sidewalls of the square duct. However, owing to the difference in the pattern of cross-stream secondary flow motions, the flow physics is significantly different in these two ribbed duct cases. It is observed that the mean flow structures in the cross-stream directions are asymmetrical in the inclined rib case but symmetrical in the V-shaped rib case, causing substantial differences in the momentum transfer across the spanwise direction. The impacts of rib geometry on near-wall turbulence structures are investigated using vortex identifiers, joint probability density functions between the streamwise and vertical velocity fluctuations, statistical moments of different orders, spatial two-point autocorrelations and velocity spectra. It is found that near the leeward and windward rib faces, the mean inclination angle of turbulence structures in the V-shaped rib case is greater than that of the inclined rib case, which subsequently enhances momentum transport between the ribbed bottom wall and the smooth top wall.

scholarly journals Turbulence and secondary motions in square duct flow

2018 ◽  
Vol 840 ◽  
pp. 631-655 ◽  
Author(s):  
Sergio Pirozzoli ◽  
Davide Modesti ◽  
Paolo Orlandi ◽  
Francesco Grasso
Keyword(s):  
Reynolds Number ◽  
Turbulent Flows ◽  
Good Accuracy ◽  
Streamwise Velocity ◽  
Secondary Flows ◽  
Duct Flow ◽  
Streamwise Vorticity ◽  
Square Duct ◽  
Developed Turbulence ◽  
Number Variation

We study turbulent flows in pressure-driven ducts with square cross-section through direct numerical simulation in a wide enough range of Reynolds number to reach flow conditions which are representative of fully developed turbulence ($Re_{\unicode[STIX]{x1D70F}}\approx 1000$). Numerical simulations are carried out over very long integration times to get adequate convergence of the flow statistics, and specifically to achieve high-fidelity representation of the secondary motions which arise. The intensity of the latter is found to be on the order of 1 %–2 % of the bulk velocity, and approximately unaffected by Reynolds number variation, at least in the range under scrutiny. The smallness of the mean convection terms in the streamwise vorticity equation points to a simple characterization of the secondary flows, which in the asymptotic high-$Re$ regime are approximated with good accuracy by eigenfunctions of the Laplace operator, in the core part of the duct. Despite their effect of redistributing the wall shear stress along the duct perimeter, we find that secondary motions do not have a large influence on the bulk flow properties, and the streamwise velocity field can be characterized with good accuracy as resulting from the superposition of four flat walls in isolation. As a consequence, we find that parametrizations based on the hydraulic diameter concept, and modifications thereof, are successful in predicting the duct friction coefficient.

The production and diffusion of vorticity in duct flow

1964 ◽  
Vol 19 (3) ◽  
pp. 375-394 ◽  
Author(s):  
E. Brundrett ◽  
W. D. Baines
Keyword(s):  
Reynolds Number ◽  
Turbulent Flow ◽  
Secondary Flow ◽  
Longitudinal Component ◽  
Secondary Flows ◽  
Duct Flow ◽  
Square Duct ◽  
Circular Pipes ◽  
Flow Circulation ◽  
Diffusion And Convection

Secondary flows in non-circular ducts are accompanied by a longitudinal component of vorticity. The equation of motion defining this component in a turbulent flow is composed of three terms giving the rates of production, diffusion and convection. Since the expression for production is the second derivative of Reynolds strees components, longitudinal vorticity cannot exist in laminar flow. For turbulent flow in a square duct the Reynolds stress tensor is examined in detail. Symmetry requirements alone provide relationships showing that the production is zero along all lines of symmetry. General characteristics of flow in circular pipes are sufficient to indicate where the production must be greatest. Experimental measurements verify this result and define the point density of production, diffusion and convection of vorticity. Data also indicate that the basic pattern of secondary flow is independent of Reynolds number, but that with increasing values of Reynolds number the flows penetrate the corners and approach the walls. A similar experimental investigation of a rectangular duct shows that the corner bisectors separate independent secondary flow circulation zones. Production of vorticity is again associated with the region near the bisector. However, there is some evidence that the secondary flow pattern is not so complex as inferred from the distortion of the main longitudinal flow.

Heat Transfer and Flowfield Measurements in the Leading Edge Region of a Stator Vane Endwall

1999 ◽  
Vol 121 (3) ◽  
pp. 558-568 ◽  
Author(s):  
M. B. Kang ◽  
A. Kohli ◽  
K. A. Thole
Keyword(s):  
Heat Transfer ◽  
Reynolds Number ◽  
Leading Edge ◽  
Reynolds Numbers ◽  
Secondary Flows ◽  
Edge Region ◽  
Leading Edge Vortex ◽  
Stator Vane ◽  
Edge Vortex ◽  
The Mean

The leading edge region of a first-stage stator vane experiences high heat transfer rates, especially near the endwall, making it very important to get a better understanding of the formation of the leading edge vortex. In order to improve numerical predictions of the complex endwall flow, benchmark quality experimental data are required. To this purpose, this study documents the endwall heat transfer and static pressure coefficient distribution of a modern stator vane for two different exit Reynolds numbers (Reex = 6 × 105 and 1.2 × 106). In addition, laser-Doppler velocimeter measurements of all three components of the mean and fluctuating velocities are presented for a plane in the leading edge region. Results indicate that the endwall heat transfer, pressure distribution, and flowfield characteristics change with Reynolds number. The endwall pressure distributions show that lower pressure coefficients occur at higher Reynolds numbers due to secondary flows. The stronger secondary flows cause enhanced heat transfer near the trailing edge of the vane at the higher Reynolds number. On the other hand, the mean velocity, turbulent kinetic energy, and vorticity results indicate that leading edge vortex is stronger and more turbulent at the lower Reynolds number. The Reynolds number also has an effect on the location of the separation point, which moves closer to the stator vane at lower Reynolds numbers.

Laminarization of Internal Flows Under the Combined Effect of Strong Curvature and Rotation

2008 ◽  
Vol 130 (9) ◽  
Author(s):  
K. M. Guleren ◽  
I. Afgan ◽  
A. Turan
Keyword(s):  
Secondary Flows ◽  
Combined Effect ◽  
Navier Stokes ◽  
Duct Flow ◽  
Internal Flows ◽  
Square Duct ◽  
Eddy Simulation ◽  
Numerical Predictions ◽  
Wide Range ◽  
Strong Curvature

The laminarization phenomenon for the flow under the combined effect of strong curvature and rotation is discussed based on numerical predictions of large-eddy simulation (LES). Initially, the laminarization process is presented for the fully developed flow inside a spanwise rotating straight square duct. LES predictions over a wide range of rotation numbers (Ro=0–5) show that the turbulent kinetic energy decreases monotonically apart from 0.2<Ro<0.5. Subsequently, a spanwise rotating U-duct flow is considered with Ro=±0.2. The interaction of curvature and Coriolis induced secondary flows enhances the turbulence for the negative rotating case, whereas this interaction ensues strong laminarization for the positive rotating case. Finally, the laminarization is presented in the impeller of a typical centrifugal compressor, rotating at a speed of Ω=1862rpm(Ro=0.6). The resulting LES predictions are observed to be better than those of Reynolds-averaged Navier-Stokes (RANS) in the regions where turbulence is significant. However, for the regions dominated by strong laminarization, RANS results are seen to approach those of LES and experiments.

Experimental Studies and Correlations of Radially Outward and Inward Air-Flow Heat Transfer in a Rotating Square Duct

1996 ◽  
Vol 118 (1) ◽  
pp. 23-30 ◽  
Author(s):  
C. R. Kuo ◽  
G. J. Hwang
Keyword(s):  
Heat Transfer ◽  
Reynolds Number ◽  
Experimental Studies ◽  
Leading Edge ◽  
Wall Heat Transfer ◽  
Square Duct ◽  
Reinforced Plastic ◽  
Flow Heat ◽  
The Mean ◽  
Heated Length

Experiments were conducted to investigate the convective heat transfer of radially outward and inward air flows in a uniformly heated rotating square duct. The interior duct surfaces, constructed by fiberglass-reinforced plastic, were plated with separated film heaters for distinguishing the local wall heat transfer rate. The duct hydraulic diameter, the actively heated length, and the mean rotation radius are 4, 120, and 180 mm, respectively. In the experiments, the parameters were the throughflow Reynolds number, Re = 1,000∼15,000; the rotation number, Ro = 0∼0.32; and the rotational buoyancy parameter, Ra* = 0∼0.5. For the outward flow the Coriolis-induced cross-stream secondary flow strongly enhanced the heat transfer on the leading edge. But for the radially inward flow the trend was reversed. When the throughflow Reynolds number was increased, the rotating-buoyancy decreased, then increased the heat transfer for the outward flow; however, the rotating-buoyancy always increased the heat transfer for the inward flow. The heat transfer data are correlated for the outward and inward flows for the ranges of parameters under study.

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.

scholarly journals Reynolds number dependence of mean flow structure in square duct turbulence

2010 ◽  
Vol 644 ◽  
pp. 107-122 ◽  
Author(s):  
ALFREDO PINELLI ◽  
MARKUS UHLMANN ◽  
ATSUSHI SEKIMOTO ◽  
GENTA KAWAHARA
Keyword(s):  
Reynolds Number ◽  
Secondary Flow ◽  
Buffer Layer ◽  
Turbulent Flows ◽  
Reynolds Numbers ◽  
The Other ◽  
Streamwise Vorticity ◽  
Square Duct ◽  
Layer Structures ◽  
The Mean

We have performed direct numerical simulations of turbulent flows in a square duct considering a range of Reynolds numbers spanning from a marginal state up to fully developed turbulent states at low Reynolds numbers. The main motivation stems from the relatively poor knowledge about the basic physical mechanisms that are responsible for one of the most outstanding features of this class of turbulent flows: Prandtl's secondary motion of the second kind. In particular, the focus is upon the role of flow structures in its generation and characterization when increasing the Reynolds number. We present a two-fold scenario. On the one hand, buffer layer structures determine the distribution of mean streamwise vorticity. On the other hand, the shape and the quantitative character of the mean secondary flow, defined through the mean cross-stream function, are influenced by motions taking place at larger scales. It is shown that high velocity streaks are preferentially located in the corner region (e.g. less than 50 wall units apart from a sidewall), flanked by low velocity ones. These locations are determined by the positioning of quasi-streamwise vortices with a preferential sign of rotation in agreement with the above described velocity streaks' positions. This preferential arrangement of the classical buffer layer structures determines the pattern of the mean streamwise vorticity that approaches the corners with increasing Reynolds number. On the other hand, the centre of the mean secondary flow, defined as the position of the extrema of the mean cross-stream function (computed using the mean streamwise vorticity), remains at a constant location departing from the mean streamwise vorticity field for larger Reynolds numbers, i.e. it scales in outer units. This paper also presents a detailed validation of the numerical technique including a comparison of the numerical results with data obtained from a companion experiment.

Dispersion of a ball settling through a quiescent neutrally buoyant suspension

1998 ◽  
Vol 361 ◽  
pp. 309-331 ◽  
Author(s):  
JAMES R. ABBOTT ◽  
ALAN L. GRAHAM ◽  
LISA A. MONDY ◽  
HOWARD BRENNER
Keyword(s):  
Settling Velocity ◽  
Volume Fraction ◽  
Functional Dependence ◽  
Vertical Component ◽  
Transversely Isotropic ◽  
Suspended Particles ◽  
Three Dimensions ◽  
Ball Diameter ◽  
The Mean ◽  
Neutrally Buoyant

Individual falling balls were allowed to settle through otherwise quiescent well-mixed suspensions of non-colloidal neutrally buoyant spheres dispersed in a Newtonian liquid. Balls were tracked in three dimensions to determine the variances in their positions about a mean uniform vertical settling path. The primary experimental parameters investigated were the size of the falling ball and the volume fraction and size of the suspended particles. Unlike the horizontal variances, the vertical variances were found to be affected by short-time deterministic behaviour relating to the instantaneous local configurational arrangement of the suspended particles. For sufficiently long intervals between successive observations, the trajectories of the balls were observed to disperse about their mean settling paths in a random manner. This points to the existence of a Gaussian hydrodynamic dispersivity that characterizes the linear temporal growth of the variance in the position of a falling ball. The functional dependence of these horizontal and vertical dispersivities upon the parameters investigated was established.The dispersivity dyadic was observed to be transversely isotropic with respect to the direction of gravity, with the vertical component at least 25 times larger than the horizontal component. The vertical dispersivity Dˆv (made dimensionless with the diameter of the suspended spheres and the mean settling velocity) was observed to decrease with increasing falling ball diameter, but to decrease less rapidly with concentration than theoretically predicted for very dilute suspensions; moreover, for falling balls equal in size to the suspended spheres, Dˆv increased linearly with increasing volume fraction ϕ of suspended solids.In addition to the above experiments performed on suspensions of spheres, previously published settling-velocity data on the fall of balls through neutrally buoyant suspensions of rods possessing an aspect ratio of 20 were re-analysed, and vertical dispersivities calculated therefrom. (These data, taken by several of the present investigators in conjunction with other researchers, had only been grossly analysed in prior publications to extract the mean settling velocity of the ball, no attempt having been made at the time to extract dispersivity data too.) The resulting vertical dispersivities, when rendered dimensionless with the rod length and mean settling velocity, showed no statistically significant dependence upon the falling-ball diameter; moreover, all other things being equal, these dispersivities were observed to increase with increasing rod concentration.

Sign in / Sign up

Export Citation Format

Share Document