Secondary cells and separation in developing laminar curved-pipe flows

Laminar flows through 180° curved bends of circular cross section are investigated numerically. For small curvature ratio, , defined as pipe radius over mean bend radius, the governing equations could be parabolized. The equations are solved for an range of from 0.04 to 0.143, a Dean number (De) range of from 277.5 to 1360, and for a uniform flow, a potential vortex, and a parabolic flow inlet condition. In all these studies a zero cross-stream flow at the inlet is assumed. A detailed study of the effects of , De, and inlet condition on the secondary flow pattern is carried out. Within the range of parameters investigated, up to three secondary cells are found in the cross-stream half-plane of a curved pipe. They are the Dean-type secondary cell, a secondary separation cell near the inner bend (closest to the center of curvature of the bend), and a third cell near the pipe center. The number of secondary cells in the cross-stream half-plane is greatly influenced by the inlet flow, and to a much lesser extent by and De. For example, only the Dean cell is found in a curved-pipe flow where and De are small and the inlet flow is either uniform or a potential vortex. When the inlet condition of the same case is changed to a parabolic flow, a three-cell structure results. Furthermore, as De increases to 1180, incipient axial flow separation begins at around 23° downstream of the curved-pipe entrance. The formation and extent of the separation and third cells are investigated together with their dependence on the parameters studied. This investigation further shows that, within the range of parameters examined, there is no secondary cell occurring near the outer bend, contrary to some earlier findings on fully developed curved-pipe flows.This work was supported by the Office of Naval Research under Grant No. N0014-81-K-0428 and by DTRC, Annapolis, Maryland, under Contract No. N00167-86-K-0075. Also, support in the form of an IPA awarded to RMCS during his sabbatical leave at DTRC, Annapolis, Maryland, in the spring of 1990 is gratefully acknowledged.     

Fluid flow in a helical pipe

Without simplifying the N-S equations of Germano's^[5], we study the flow in a helical circular pipe employing perturbation method. A third perturbation solution is fully presented. The first- second- and third-order effects of curvature κ and torsion τ on the secondary flow and axial velocity are discussed in detail. The first-order effect of curvature is to form two counter-rotating cells of the secondary flow and to push the maximum axial velocity to the outer bend. The two cells are pushed to the outer bend by the pure second-order effect of curvature. The combined higher-order (second-, third-) effects of curvature and torsion, are found to be an enlargement of the lower vortex of the secondary flow at expense of the upper one and a clockwise shift of the centers of the secondary vortices and the location of maximum axial velocity. When the axial pressure gradient is small enough or the torsion is sufficiently larger than the curvature, the location of the maximal axial velocity is near the inner bend. The equation of the volume flux is obtained from integrating the perturbation solutions of axial velocity. From the equation the validity range of the perturbation solutions in this paper can be obtained and the conclusion that the three terms of torsion have no effect on the volume flux can easily be drawn. When the axial pressure gradient is less than 22.67, the volume flux in a helical pipe is larger than that in a straight pipe.     

Experimental investigation on turbulent flow through a circular-sectioned 180° bend

The steady, developing turbulent flow in a circular-sectioned 180° bend has been investigated. The bend had a radius of 104 mm and a curvature radius ratio of 4.0 with long, straight upstream and downstream pipes. Measurements of the longitudinal, radial and circumferential components of mean velocity, and corresponding components of the Reynolds stress were obtained with a hot wire anemometer at a Reynolds number of 6×10⁴ and at various longitudinal stations. The velocity fields of the primary and secondary flows and the Reynolds stresses were illustrated in the form of contour map or vector diagram. Moreover, the mean quantities characterizing the bend flow, i.e., the deflection of the primary flow in the cross section, the intensity of the secondary flow and the turbulence energy, were shown in a graphic form against the longitudinal distances. In the section upstream from a bend angle of about 60°, both the flows through the 180° and the 90° bend are closely similar in their behavior. In the section from the bend angle of 90°, the high-velocity regions, however, occur near the upper and lower walls as a result of strong secondary flow and the turbulence with high level emerges in the central region of the bend. Just behind the bend exit, an additional pair of vortices appears in the outer part of the cross section owing to the transverse pressure difference. In the downstream tangent, the flow returns slowly to the proper flow in a straight pipe, but it needs a longer distance for recovery than in the 90° bend. Received: 23 April 1998/Accepted: 24 April 1999     

Turbulent pipe flow downstream a 90° pipe bend with and without superimposed swirl

In the present work, the turbulent flow downstream a 90° pipe bend is investigated by means of stereoscopic particle image velocimetry. In particular, the three dimensional flow field at the exit of the curved pipe is documented for non-swirling and swirling flow conditions, with the latter being generated through a unique axially rotating pipe flow facility. The non-swirling flow was examined through snapshot proper orthogonal decomposition (POD) with the aim to reveal the unsteady behaviour of the Dean vortices under turbulent flow conditions, the so-called “swirl-switching” phenomenon. In respect to the swirling turbulent pipe flow, covering a wide range of swirl strengths, POD has been employed to study the effect of varying strength of swirl on the Dean vortices as well as the interplay of swirling motion and Dean cells. Furthermore, the visualised large scale structures in turbulent swirling flows through the bend are found to incline and tear up with increasing swirl intensity. The present time-resolved, three component, experimental velocity field data will provide a unique and useful database for future studies; in particular for the CFD community.     

Investigation on Secondary Flow Characteristics in a Curved Annular Duct with Struts

Under the influence of duct curvature, cross-sectional area variation and internal struts, the internal flow field within a curved annular duct becomes rather complicated and contains strong secondary flow. In this paper, the secondary flow characteristics in an annular duct with struts are experimentally and numerically investigated. The results show that large pressure gradients exist on the bends of hub and shroud. Meanwhile, two counter-rotating vortex pairs appear both along the hub-side and shroud-side surfaces. The hub-side vortex pair of which the vortex cores travel downstream parallelly evolves from the horseshoe vortex which is induced by the leading edge of the upstream strut, whereas the shroud-side vortex pair originates from the strut trailing edge and the corresponding vortex cores develop in a divergent way. Additionally, the effects of the duct exit Mach number on the secondary flow characteristics are also studied. As the exit Mach number increases, the streamwise pressure gradients increase and lead to more intense vortices, higher total pressure loss and larger flow distortion.     

弯曲动脉的血流动力学数值分析

利用计算流体力学的理论和方法对弯曲动脉中的血流动力学进行数值分析，是研究心血管疾病流体动力学机理的一种行之有效的方法。本文将升主动脉、主动脉弓和降主动脉联系起来作为弯曲动脉几何模型，给出了血液流动的边界条件以及计算条件。根据生理脉动流条件，对狗的弯曲动脉几何模型内发展中的血液流动进行了有限元数值模拟，并利用可视化方法对血液流动的轴向速度、二次流、壁面切应力等计算结果进行了分析。研究结果表明，在弯管内侧壁处，同时存在主流方向和二次流方向的回流，此处容易形成涡流。弯管内侧壁比外侧壁的壁面切应力具有更强的脉动性。     

环形截面螺旋管道内二次流动特性的研究

从曲线柱坐标系下的N-S方程出发,以曲率和挠率为小参数,采用摄动法求解了环形截面螺旋管道内的黏性流动,给出了完全二阶摄动解,结果表明：当挠率为零时,二次流表现为上下对称的四个涡;当挠率不为零,涡的对称性遭到破坏,二次涡的强度和个数受De数和环形截面内外径之比δ的影响,轴向速度最大值在De数较小时靠近管道的内侧,随着De数的增加,其最大值向外侧移动。     

90°弯管内气固两相湍流流动的数值研究

利用两相湍流KET模型对90°弯管内气固两相湍流流动进行了数值模拟,得到了弯管内两相流动的一些规律,并提出用颗粒动理学压力来定性表征弯管内磨损严重部位,为管道抗磨损设计提供了一定的理论依据。     

Dean instability and secondary flow structure in curved rectangular ducts

The pressure driven, fully developed turbulent flow of incompressible viscous fluid (water) in 120° curved ducts of rectangular cross-section is investigated experimentally and numerically. Three different types of curved duct (A-CL, B-SL and C-IL) with continuously varying curvature conform to blade profile as the inner and outer curvature walls to simplify and guide the impeller design of pumps. After validating the numerical method against Particle Image Velocimetry (PIV) measurements, the flow development in the ducts is analyzed in detail by Computational Fluid Dynamics (CFD) for a wide range of Reynolds numbers (Re = 2.4 × 10⁴–1.4 × 10⁵) and aspect ratios (Ar > 1.0, =1.0 and <1.0). The results clearly depict the existence of multiple Dean vortices along the duct: while the axial velocity profile is more related to an inner Dean vortex (called split base vortex), the wall pressure is more influenced by the Dean vortex attached to the inner curvature wall (called ICW Dean vortex). The induced multiple Dean vortices and the secondary flow patterns in the duct cannot be faithfully predicted by using traditional techniques. Therefore, a new criterion based on the vortex core velocities is proposed. With this approach, the effects of Re, Cr and Ar on the Dean instabilities in curved ducts are carefully studied. Decreasing Re promotes the generation of Dean vortices closer to the duct inlet, a trend that is as opposed to laminar flow. In addition, a new pair of vortices called entrainment Dean vortex occurs near the outlet of the curved duct with Ar = 1.0, which has not been previously reported in the literature.     

Low Reynolds number axisymmetric turbulent boundary layer on a cylinder

In the present study, an axisymmetric turbulent boundary layer growing on a cylinder is investigated experimentally using hot wire anemometry. The combined effects of transverse curvature as well as low Reynolds number on the mean and turbulent flow quantities are studied. The measurements include the mean velocity, turbulence intensity, skewness and flatness factors in addition to wall shear stress. The results are presented separately for the near wall region and the outer region using dimensionless parameters suitable for each case. They are also compared with the results available in the open literature.The present investigation revealed that the mean velocity in near wall region is similar to other simple turbulent flows (flat plate boundary layer, pipe and channel flows); but it differs in the logarithmic and outer regions. Further, for dimensionless moments of higher orders, such as skewness and flatness factors, the main effects of the low Reynolds number and the transverse curvature are present in the near wall region as well as the outer region.     

Vortical patterns in turbulent flow downstream a 90° curved pipe at high Womersley numbers

The present experimental work focuses on highly pulsatile, i.e. inertia dominated, turbulent flow downstream a curved pipe and aims at investigating the vortical characteristics of such a flow. The flow parameters (Dean and Womersley number) investigated are of the same order as those met in the internal combustion engine environment. The technique employed is time-resolved stereoscopic particle image velocimetry at different cross-sections downstream the pipe bend. These measurements allow the large-scale structures that are formed to be analyzed by means of proper orthogonal decomposition. The flow field changes drastically during a pulsatile cycle, varying from a uniform flow direction across the pipe section from the inside to the outside of the bend to vortical patterns consisting of two counter-rotating cells. This study characterizes and describes pulsatile curved pipe flow at Womersley numbers much higher than previously reported in the literature. Furthermore, the oscillatory behaviour of the Dean cells for the steady flow – the so-called ‘swirl switching’ – is investigated for different downstream stations from the bend exit and it is shown that this motion does not appear in the immediate vicinity of the bend, but only further downstream.     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

A criterion for detection of the onset of Dean instability in Newtonian fluids

Dean instability for Newtonian fluids in laminar secondary flow in 180° curved channels was studied experimentally and numerically. The numerical study used Fluent CFD code to solve the Navier–Stokes equations, focusing on flow development conditions and the parameters influencing Dean instability. An accurate criterion based on the radial gradient of the axial velocity was defined that allows detection of the instability threshold, and this criterion is used to optimize the grid geometry. The effects on Dean instability of the curvature ratio (from 5.5 to 20) and aspect ratio (from 0.5 to 12) are studied. In particular, we show that the critical value of the Dean number decreases with the increasing duct curvature ratio. The variation of the critical Dean number with duct aspect ratio is less regular.In the experimental study, flows were visualized in several tangential positions of a 180° curved channel with aspect ratio 8 and curvature ratio 10. The flow is hydrodynamically developed at the entrance to the curved channel. The critical Dean number is detected and the development of secondary flow vortices by additional counter-rotating vortex pairs is observed. A diagram of different critical Dean numbers is established.     

Direct numerical simulation of low Reynolds number flows in an open‐channel with sidewalls

A direct numerical simulation of low Reynolds number turbulent flows in an open‐channel with sidewalls is presented. Mean flow and turbulence structures are described and compared with both simulated and measured data available from the literature. The simulation results show that secondary flows are generated near the walls and free surface. In particular, at the upper corner of the channel, a small vortex called inner secondary flows is simulated. The results show that the inner secondary flows, counter‐rotating to outer secondary flows away from the sidewall, increase the shear velocity near the free surface. The secondary flows observed in turbulent open‐channel flows are related to the production of Reynolds shear stress. A quadrant analysis shows that sweeps and ejections are dominant in the regions where secondary flows rush in toward the wall and eject from the wall, respectively. A conditional quadrant analysis also reveals that the production of Reynolds shear stress and the secondary flow patterns are determined by the directional tendency of the dominant coherent structures. Copyright © 2009 John Wiley & Sons, Ltd.     

Laser-Doppler measurements of laminar and turbulent flow in a pipe bend

Laser-Doppler measurements are reported for laminar and turbulent flow through a 90° bend of circular cross-section with mean radius of curvature equal to 2.8 times the diameter. The measurements were made in cross-stream planes 0.58 diameters upstream of the bend inlet plane, in 30, 60 and 75° planes in the bend and in planes one and six diameters downstream of the exit plane. Three sets of data were obtained: for laminar flow at Reynolds numbers of 500 and 1093 and for turbulent flow at the maximum obtainable Reynolds number of 43 000. The results show the development of strong pressure-driven secondary flows in the form of a pair of counter-rotating vortices in the streamwise direction. The strength and character of the secondary flows were found to depend on the thickness and nature of the inlet boundary layers, inlet conditions which could not be varied independently of Reynolds number. The quantitative anemometer measurements are supported by flow visualization studies. Refractive index matching at the fluid-wall interface was not used; the measurements consist, therefore, of streamwise components of mean and fluctuating velocities only, supplemented by wall pressure measurements for the turbulent flow. The displacement of the laser measurement volume due to refraction is allowed for in simple geometrical calculations. The results are intenden for use as benchmark data for calibrating flow calculation methods.     

Evolution of turbulence characteristics from straight to curved pipes

Fully developed, statistically steady turbulent flow in straight and curved pipes at moderate Reynolds numbers is studied in detail using direct numerical simulations (DNS) based on a spectral element discretisation. After the validation of data and setup against existing DNS results, a comparative study of turbulent characteristics at different bulk Reynolds numbers Re_b = 5300 and 11,700, and various curvature parameters κ = 0, 0.01, 0.1 is presented. In particular, complete Reynolds-stress budgets are reported for the first time. Instantaneous visualisations reveal partial relaminarisation along the inner surface of the curved pipe at the highest curvature, whereas developed turbulence is always maintained at the outer side. The mean flow shows asymmetry in the axial velocity profile and distinct Dean vortices as secondary motions. For strong curvature a distinct bulge appears close to the pipe centre, which has previously been observed in laminar and transitional curved pipes at lower Re_b only. On the other hand, mild curvature allows the interesting observation of a friction factor which is lower than in a straight pipe for the same flow rate.All statistical data, including mean profile, fluctuations and the Reynolds-stress budgets, is available for development and validation of turbulence models in curved geometries.     

A pressure probe for the detection of the flow direction close to walls. Case study: flow through a bend

We report on the development of a new pressure probe that detects the flow direction in wall-bound flow. Two pressure differences are measured and combined in a pressure coefficient which is proportional to the flow direction in a ±20° range. The probe is applied to the secondary flow vortex pair generated in a 90° pipe bend, with excellent results. The vortex pair, and its downstream decay, are identified. Furthermore, the stability of the vortex pair is found to depend sensitively on the upstream conditions. When these are fixed, the vortices stay put; in other words they are spatially stabilized. The consequences for installation effects on flow metering are discussed.     

Splitting,merging and spanwise wavenumber selection of Dean vortex pairs

Splitting, merging and spanwise wavenumber selection are studied during the initial development of Dean vortex pairs in a channel with mild curvature, an aspect ratio of 40, and an inner to outer radius ratio of 0.979. Two types of splitting events, and four types of merging events are evident from flow visualizations at Dean numbers from 75 to 220. These events are described in detail along with observations that occurrances of these different events are tied to spatial and temporal variations of spanwise wavenumbers of vortex pair spacing. Also discussed are frequency spectra of different events, recurrent phenomena, and the part played by splitting and merging in laminar to turbulent transition in curved channels.This study was sponsored by the Propulsion Directorate, U.S. Army Aviation Research and Technology Activity-AVSCOM, NASA-Defense Purchase Request 030030-P. Professor C. S. Subramanian provided assistance in setting up the traversing system and data acquisition systems used for this study. Mr. R. E. Hughes assisted in obtaining some of the results presented.     

Electrokinetic flow and energy conversion in a curved microtube

Curved channels are ubiquitous in microfluidic systems. The pressuredriven electrokinetic flow and energy conversion in a curved microtube are investigated analytically by using a perturbation analysis method under the assumptions of the small curvature ratio and the Reynolds number. The results indicate that the curvature of the microtube leads to a skewed pattern in the distribution of the electrical double layer (EDL) potential. The EDL potential at the outer side of the bend is larger than that at the inner side of the bend. The curvature shows an inhibitory effect on the magnitude of the streaming potential field induced by the pressure-driven flow. Since the spanwise pressure gradient is dominant over the inertial force, the resulting axial velocity profile is skewed into the inner region of the curved channel. Furthermore, the flow rate in a curved microtube could be larger than that in a straight one with the same pressure gradient and shape of cross section. The asymptotic solutions of the axial velocity and flow rate in the absence of the electrokinetic effect are in agreement with the classical results for low Reynolds number flows. Remarkably, the curved geometry could be beneficial to improving the electrokinetic energy conversion (EKEC) efficiency.     

Swirling turbulent flow through a curved pipe

An experimental study of a swirling turbulent flow through a curved pipe with a pipe-to-mean-bend radius ratio of 0.077 and a flow Reynolds number based on pipe diameter and mean bulk velocity of 50,000 has been carried out. A rotating section, six pipe diameters long, is set up at six diameters upstream of the curved bend entrance. The rotating section is designed to provide a solid-body rotation to the flow. At the entrance of the rotating section, a fully-developed turbulent pipe flow is established. This study reports on the flow characteristics for the case where the swirl number, defined as the ratio of the pipe circumferential velocity to mean bulk velocity, is one. Wall static pressures, mean velocities, Reynolds stresses and wall shear distribution around the pipe are measured using pressure transducers, rotating-wires and surface hot-film gauges. The measurements are used to analyze the competing effects of swirl and bend curvature on curved-pipe flows, particularly their influence on the secondary flow pattern in the crossstream plane of the curved pipe. At this swirl number, all measured data indicate that, besides the decaying combined free and forced vortex, there are no secondary cells present in the cross-stream plane of the curved pipe. Consequently, the flow displays characteristics of axial symmetry and the turbulent normal stress distributions are more uniform across the pipe compared to fully-developed pipe flows.List of symbols B calibration constant - e bridge voltage - e ₀ bridge voltage at zero flow - C _f total skin friction coefficient, = 2 _w/ W ₀ ² - D pipe diameter, = 7.62 cm - De Dean number, = ^1/2 Re - M angular momentum - n calibration constant - N _s swirl number, = D/2 W ₀ - r radial coordinate - R mean bend radius of curvature, = 49.5 cm - Re pipe Reynolds number, = DW ₀/ - S axial coordinate along the upstream (measured negative) and downstream (measured positive) tangent - U, V, W mean velocities along the radial, tangential and axial directions, respectively - u, v, w mean fluctuating velocities along the radial, tangential and axial directions, respectively - u, v, w root mean square normal stress along the radial, tangential and axial directions, respectively - v ^{ov2}, u^{ov2} normal stress along the tangential and radial direction, respectively - W ₀ mean bulk velocity, 10 m/s - W _c W measured at pipe axis - W total wall friction velocity, - total wall friction velocity measured at S/D = -18 - ,v vw, w7#x016B; turbulent shear stresses - pipe-to-mean-bend radius ratio, = D/2 R = 0.077 - axial coordinate measured from bend entrance - fluid kinematic viscosity - fluid density - _w mean total wall shear stress - instantaneous total wall shear - azimuthal coordinate measured zero from pipe hori zontal diameter near outer bend - angular speed of the rotating section