Numerical studies of fluid flow through tubes with double constrictions

The flow fields in the neighbourhood of double constrictions in a circular cylindrical tube were studied numerically. The effects on the streamline, velocity and vorticity distributions as the flow passes through the constrictions in the tube were studied in the Reynolds number range 5–200. Double constrictions with dimensionless spacing ratios of 1, 2, 3 and ∞ were studied for a 50% constriction. It is noted that when the Reynolds number is below 10, no recirculation region is formed in the above constricted flow. For Reynolds numbers greater than 10, a recirculation region forms downstream of each of the constrictions. For constriction spacing ratios of 1, 2, and 3, when the Reynolds number is high, a recirculation region spreads between the valley of the constrictions. The recirculation region formed between the two constrictions has a diminishing effect on the generation of wall vorticity near the second constriction area. In general, the peak value of wall vorticity is found slightly upstream of each of the constrictions. When the Reynolds number is increased, the peak wall vorticity value increases and its location is moved upstream. Maximum wall vorticity generated by the first constriction is found to be always greater than the maximum wall vorticity generated by the second constriction. The extent of this spreading of the recirculation region from the first constriction and its effects on the second constriction depend on the constriction spacing ratio and the flow Reynolds number.     

NUMERICAL STUDY OF STEADY LAMINAR FLOW THROUGH TUBES WITH MULTIPLE CONSTRICTIONS USING CURVILINEAR CO-ORDINATES

The flow field through tubes with multiple axisymmetric constrictions in tubes was studied numerically. Two practical problem cases were considered and the numerical scheme was developed for both. In the first case there are one, two, three and four constrictions in the tube. The effects of the number of constrictions on wall shear stress, pressure drop, streamline, vorticity and velocity distributions as the flow passes through the tube were studied and the development of the periodicity characteristics was investigated. In the second case there were multiple constrictions in the tube equidistant from each other. For this case the governing equations were reformulated for a module at a sufficient distance downstream from the inlet where the entrance region effects could be ignored and flow field is assumed to repeat itself. The flow field solutions were obtained in this region. The governing equations were formulated in curvilinear co-ordinates and a finite volume discretization procedure was used to solve the problem. The computations were carried out over a range of Reynolds numbers between 50 to 250 for constrictions with 75 percent area reduction. The method is validated by comparing some of the solutions with experimental results.     

An experimental investigation of pulsatile flow through a smooth constriction

Measurements of flow disturbances in the downstream region of modeled stenoses in a rigid tube, with upstream pulsatile flow are reported. Experiments were conducted over physiologically relevant mean Reynolds numbers of 600; based on the tube diameter and the time-averaged value of upstream centerline velocity. Contoured constrictions with 25%, 50% and 75% area reductions were investigated and velocity data were obtained from ensemble averaging techniques (phase-locked waveform). Experimental data over extensive spatial regions of poststenotic fields were taken, employing a two-component laser Doppler velocimeter LDV. Constant time sampling techniques for performing data or frequency analyses were used to avoid velocity bias and to study the evolution of poststenotic flow disturbances. It is found that different types of flow disturbances exist downstream of the constriction. Data analysis methods with the aid of flow visualization allow accurate classification of the disturbances which are sensitive indicators of mild to moderate constrictions. Although the present study was motivated by a biological situation, sufficient data were reported in detail that they may also be used by investigators working in computational fluid dynamics.     

Laser Doppler velocimeter measurements to characterize turbulence in a constriction with sharp and rounded inlets

An important ingredient for successfully understanding and predicting turbulent flow behavior is an accurate characterization of the turbulence. While several models are available, it is not clear which turbulence model is appropriate for flow in a constriction geometry. To answer this question, experimental results were obtained using a laser Doppler velocimeter (LDV) and compared to predictions using three turbulence models in a computational fluid dynamics (CFD) code for two constrictions. Experimental results indicate that a small amount of rounding at the constriction inlet can have a significant effect on the velocity and turbulent kinetic energy distributions in the throat. Received: 11 September 2000 / Accepted: 29 October 2001     

Application of the PITM Method Using Inlet Synthetic Turbulence Generation for the Simulation of the Turbulent Flow in a Small Axisymmetric Contraction

We investigate the turbulence modeling of second moment closure used both in RANS and PITM methodologies from a fundamental point of view and its capacity to predict the flow in a low turbulence wind tunnel of small axisymmetric contraction designed by Uberoi and Wallis. This flow presents a complex phenomenon in physics of fluid turbulence. The anisotropy ratio of the turbulent stresses τ ₁₁/τ ₂₂ initially close to 1.4 returns to unity through the contraction, but surprisingly, this ratio gradually increases to its pre-contraction value in the uniform section downstream the contraction. This point constitutes the interesting paradox of the Uberoi and Wallis experiment. We perform numerical simulations of the turbulent flow in this wind tunnel using both a Reynolds stress model developed in RANS modeling and a subfilter scale stress model derived from the partially integrated transport modeling method. With the aim of reproducing the experimental grid turbulence resulting from the effects of the square-mesh biplane grid on the uniform wind tunnel stream, we develop a new analytical spectral method of generation of pseudo-random velocity fields in a cubic box. These velocity fields are then introduced in the channel using a matching numerical technique. Both RANS and PITM simulations are performed on several meshes to study the effects of the contraction on the mean velocity and turbulence. As a result, it is found that the RANS computation using the Reynolds stress model fails to reproduce the increase of anisotropy in the centerline of the channel after passing the contraction. In the contrary, the PITM simulation predicts fairly well this turbulent flow according to the experimental data, and especially, the “return to anisotropy” in the straight section of the channel downstream the contraction. This work shows that the PITM method used in conjunction with an analytical synthetic turbulence generation as inflow is well suited for simulating this flow, while allowing a drastic reduction of the computational resources.     

Static pressure effects on calibration of velocity transducers at nozzle exits

Calibrations of velocity transducers are generally performed at nozzle exits where the turbulence level is minimal. Such nozzles require a high contraction ratio and are usually short with a consequent effect on stream static pressure due to formation of a vena contracta. Calibrations performed with a Pitot tube and no static pressure correction can lead to a significant velocity error. Measurements of static pressure, mean velocity and turbulence intensity are presented for a typical high contraction ratio nozzle. These show two distinct flow regimes which can be described very simply by the Reynolds form of the Navier-Stokes equations.     

A study of floc breakup and formation in flowing concentrated fiber suspensions

An experimental study was carried out to examine the effects of fiber concentration and average flow velocity on the breakup of fiber flocs and the subsequent reflocculation process downstream of a converging/diverging flow section used as a turbulence generator. Flocculation measurements were made using fiber optic probes and are presented in terms of floc-size distributions. These were complemented with pressure fluctuation measurements used to indicate the state of fluidization of the suspension.     

An experimental study of nonrheometric flows of visco-elastic fluids: I. flow into a re-entrant tube

A laser anemometer has been used to study the developing flow both upstream and downstream from the entry plane in a re-entrant tube geometry. A 0.75% polyacrylamide/water solution was used and Reynolds numbers (based on wall conditions in the fully developed downstream flow) in the range 100–500 were obtained in 1.82-cm and 2.40-cm-diameter tubes.The shear stress-shear rate relationship for the fluid was measured using a cone and plate geometry in conjunction with a Weissenberg rheogoniometer. Theoretical fully developed velocity profiles were calculated numerically from these measurements. The measured fully developed velocity profiles were found to be in excellent agreement with those calculated.Velocity profiles measured at the tube entry plane showed the pronounced wall region distortion typically predicted by recent numerical solutions of the flow of purely-viscous fluids through an abrupt tube contraction.It was found that the major velocity rearrangements were achieved within only a few diameters (both upstream and downstream) of the entry plane. In particular, the velocity distribution near the tube wall varied negligibly over the relatively longer distance (many diameters) that it took for the centreline velocity to achieve its fully developed value. Entry lengths were found to be only about half those for purely-viscous fluids.Calculation of the time of flight along the central streamline confirmed that the major rearrangements of velocity suffered by the fluid occurred over a relatively short time period. This indicates that hereditary integral constitutive equations may have to be used in theoretical analyses of this type of flow situation.     

Unsteady laminar separated flow through constricted channel

In the present paper unsteady Navier-Stokes equations have been solved numerically by finite-difference technique in staggered grid distribution for a flow through a channel with locally symmetric and asymmetric constrictions. A coordinate stretching has been made to map the infinite irregular geometry into a finite regular computational domain. Pressure and pressure-velocity corrections scheme have been developed. Convergence criteria (in terms of continuity equation) has been achieved after few time iterations. The critical Reynolds number for asymmetric flow through a symmetric constriction has been found. Critical values depend on the area reduction and the length of the constriction. The increment of Reynolds number grows the asymmetry of the flow. The root mean square (r.m.s.) centreline vertical velocity for asymmetric flow through a symmetric constriction has been drawn at different Reynolds numbers. For flow through symmetric constriction the centreline vertical velocity shows finite oscillation behind the constriction at high Reynolds number.     

Parametric vibrations of flexible hoses excited by a pulsating fluid flow,Part II: Experimental research

The paper presents the results of experimental studies of vibrations of an elastic hose which are induced by a pulsating fluid flow. It was found that there is a possibility of parametric resonances: principal and combination associated with certain modes of vibrations. The influence of frequency and the amplitude of pulsation, average flow velocity, pressure inside pipe, the length of the hose, and the temperature on the ranges of parametric vibrations were analysed. The character of vibrations in resonance ranges was determined by showing time histories and the results of spectral analyses. A flexible hose applied in high-pressure hydraulic systems was used as an object of research. The values of basic parameters which describe the hose׳s mechanical properties were identified experimentally. The results of the experiments were compared with the results of numerical simulations conducted on the basis of the methodology proposed in Part I of this paper.     

Decay of passive-scalar fluctuations in slightly stretched grid turbulence

Isotropic turbulence is closely approximated by stretching a grid flow through a short (1.36:1) secondary contraction. The flow is operated at small values of the Taylor microscale Reynolds number (about 25–55) and is slightly heated just downstream of the grid, so that the temperature serves as a passive scalar and the initial velocity/thermal length-scale ratio is about 1. For the same grid, the contraction reduces the skewness and kurtosis of the thermal fluctuations and their derivative. The thermal fluctuations and their mean dissipation rates follow a power-law rate of decay that depends on the geometry of the grid. Comparison with velocity measurements shows that, for three different grids, the ratio between the temperature and velocity power-law exponents closely matches the velocity/thermal timescale ratio. For the present measurements, the timescale ratio is slightly larger than 1 but does not exceed 1.2, in accordance with the proposal by Corrsin (J Aeronaut Sci 18(6):417–423, 1951b).     

Free-stream turbulence effects on the instantaneous pressure and forces on cylinders of rectangular cross section

Simultaneous measurements of instantaneous pressure distributions on rectangular cylinders of length to height ratio(B/D) of 1.0, 2.5 and 3.0 in smooth nonturbulent and homogeneous turbulent flows were made and the data were analyzed by phase averaging and spectral analysis in addition to more conventional methods. The turbulence in the inflow stream is nearly homogeneous and isotropic with the intensity and the scale of 5% and 1.2-1.5 times the cylinder height, respectively. The main effects of the turbulence in the inflow free stream of this scale and intensity are to laterally move the separated shear flow off the upstream corners and cause intermittent reattachment on the side surfaces of cylinders of B/D of 2.5 and larger. For the cylinder with smaller B/D, the flow does not reattach with or without turbulence in the free stream, and the instantaneous surface pressure distributions fluctuate quite periodically at a frequency corresponding to the Strouhal frequency of the vortex shedding. The effects of the free-stream turbulence appear in the increased fluctuation on the front surface as buffeting due to the impinging turbulence. When the separated shear layers reattach due to the influence of the free-stream turbulence, the reattachment point moves intermittently, the pressure distributions downstream of the reattachment fluctuate periodically, and a mild peak is formed in the spectra at a frequency much larger than the Strouhal frequency.     

Investigation of turbulence effects on forced convection in a composite porous/fluid duct: Constant wall flux and constant wall temperature cases

Forced convection heat transfer in composite porous/fluid domains is of great practical and theoretical significance. However, research in this area traditionally addressed only the laminar flow case in both homogeneous fluid and porous regions of the domain. This paper investigates the interaction between turbulent flow in the center of a circular tube filled with a homogeneous fluid and laminar flow in the porous layer adjacent to the tube wall. A two-layer algebraic turbulence model suggested by Cebeci and Smith is utilized for the flow in the central region of the tube. The effects of turbulence in the central region on velocity and temperature distributions as well as on the Nusselt number are analyzed.     

From single-scale turbulence models to multiple-scale and subgrid-scale models by Fourier transform

A theoretical method based on mathematical physics formalism that allows transposition of turbulence modeling methods from URANS (unsteady Reynolds averaged Navier–Stokes) models, to multiple-scale models and large eddy simulations (LES) is presented. The method is based on the spectral Fourier transform of the dynamic equation of the two-point fluctuating velocity correlations with an extension to the case of non-homogenous turbulence. The resulting equation describes the evolution of the spectral velocity correlation tensor in wave vector space. Then, we show that the full wave number integration of the spectral equation allows one to recover usual one-point statistical closure whereas the partial integration based on spectrum splitting gives rise to partial integrated transport models (PITM). This latter approach, depending on the type of spectral partitioning used, can yield either a statistical multiple-scale model or subfilter transport models used in LES or hybrid methods, providing some appropriate approximations are made. Closure hypotheses underlying these models are then discussed by reference to physical considerations with emphasis on identification of tensorial fluxes that represent turbulent energy transfer or dissipation. Some experiments such as the homogeneous axisymmetric contraction, the decay of isotropic turbulence, the pulsed turbulent channel flow and a wall injection induced flow are then considered as typical possible applications for illustrating the potentials of these models.      

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

Measurement of bubble size, velocity and concentration in flashing flow behind a sudden constriction

In the present paper the results of investigations in flashing flow behind a sudden constriction in vertical upflow are described. Flow visualization, laser-Doppler and phase-Doppler anemometry have been used to measure local bubble and fluid velocities, local bubble sizes and void fractions. The measurements were performed in the midplane of a two-dimensional channel with a 2:1 stepwise constriction.It was found that bubble nucleation takes place in the recirculation zone immediately behind the constriction, which is the location of the lowest static pressure. These bubbles are transported downstream by the mean flow field, while undergoing further growth. No additional nucleation was observed downstream of the recirculation zone. A periodic, cloudwise behaviour of the bubble formation was found which could be explained by the interaction between the bubble growth and the mean flow field. This interaction results in strong disturbances of the mean flow field, which show up as an increase of the fluctuating bubble velocity by a factor of 3 compared to single-phase measurements in a region of 10 step heights behind the constriction. However, these fluctuations appear more like a periodic change in the mean velocity rather than a higher turbulence level. The measured arithmetic mean bubble diameters rise from approx. 50 μm in the recirculation region to about 70–80 μm 50 step heights downstream. Maximum local bubble number density and void fraction were found to be 160001/cm³ and 0.8%, respectively.     

Fluid-elastic instability in tube arrays subjected to air-water and steam-water cross-flow

Flow induced vibrations in heat exchanger tubes have led to numerous accidents and economic losses in the past. Efforts have been made to systematically study the cause of these vibrations and develop remedial design criteria for their avoidance. In this research, experiments were systematically carried out with air-water and steam-water cross-flow over horizontal tubes. A normal square tube array of pitch-to-diameter ratio of 1.4 was used in the experiments. The tubes were suspended from piano wires and strain gauges were used to measure the vibrations. Tubes made of aluminum; stainless steel and brass were systematically tested by maintaining approximately the same stiffness in the tube-wire systems. Instability was clearly seen in single phase and two-phase flow and the critical flow velocity was found to be proportional to tube mass. The present study shows that fully flexible arrays become unstable at a lower flow velocity when compared to a single flexible tube surrounded by rigid tubes. It is also found that tubes are more stable in steam-water flow as compared to air-water flow. Nucleate boiling on the tube surface is also found to have a stabilizing effect on fluid-elastic instability.     

An experimental study of nonrheometric flows of viscoelastic fluids: II. Flow in a liquid jet exiting from a vertical tube

A laser anemometer has been used to study the region of accelerating shear flow near the exit of a vertical tube. It is in this region that the transition between steady laminar shear flow in the upstream tube and elongational flow in the downstream liquid jet takes place.Downstream velocity profiles were measured for solutions of 0.9% polyacrylamide in 85% glycerol/water and 0.9% polyacrylamide in water. Reynolds numbers (based on wall conditions in the fully developed upstream flow) ranged from 45 to 310 and Froude numbers from 0.294 to 4.11. Tubes, having sharpedged and rounded exit corners, with diameters of 1.25 cm and 1.90 cm were usedUpstream velocity profiles were measured for a solution of 0.9% polyacrylamide in water. Reynolds numbers ranged from 16 to 670. Only tubes having sharp-edged exit corners were used.It was found that the transition region did not extend upstream into the tube but was confined to the downstream jet. The transition took place over a distance of about 3–5 tube diameters depending upon the value of the Froude number. The axial distance downstream from the tube exit plane at which the velocity profile first became flat increased with increasing Froude number. The magnitude of the jet velocity at this point decreased with increasing Froude number.The condition of the tube exit corner was found to influence the flow in the transition region. Downstream velocity profiles obtained using tubes having rounded exit corners initially develop more slowly than, but soon catch up with and eventually overtake, the corresponding profiles obtained using tubes with sharp-edged exit corners.Downstream velocity profiles obtained for the 0.9% polyacrylamide in 85% glycerol/water solution were found to develop smoothly. The transition from steady shear flow in the tube to elongational flow in the jet took place through the combined processes of acceleration of the outer layers of the jet due to radial transfer of momentum with adjacent inner layers, the process spreading steadily inwards with increasing axial distance from the tube exit plane, and acceleration of the whole due to gravity. However, the velocity profiles obtained for the 0.9% polyacrylamide in water solution did not always develop so smoothly. At a Reynolds number of 310 and Froude number of 2.06 the radial momentum transfer process was restricted to a narrow outer region of the jet until a downstream axial distance of about 2 tube diameters was reached. Thereafter, the transition to a flat profile took place smoothly.     

The development of an axisymmetric curved turbulent wall jet

An experimental study has been carried out of the low speed Coanda wall jet with both streamwise and axisymmetric curvature. A single component laser Doppler technique was used, and by taking several orientations at a given point, values of the three mean velocities and five of the six Reynolds stresses were obtained. The lateral divergence and convex streamwise curvature both enhanced the turbulence in the outer part of the jet compared with a plane two-dimensional wall jet. The inner layer exhibited a large separation of the positions of maximum velocity and zero shear stress. It was found that the streamwise mean velocity profile became established very rapidly downstream of the slot exit. The profile appeared fairly similar at later downstream positions, but the mean radial velocity and turbulence parameters showed the expected nonself preservation of the flow. Removal of the streamwise curvature resulted in a general return of the jet conditions toward those expected of a plane wall jet. The range and accuracy of the data may be used for developing turbulence models and computational techniques for this type of flow.     

Near Field Round Jet Flow Downstream from an Abrupt Contraction Nozzle with Tube Extension

Round air jet development downstream from an abrupt contraction coupled to a uniform circular tube extension with length to diameter ratio L/D?=?1.2 and L/D?=?53.2 is studied experimentally. Smoke visualisation and systematic hot film velocity measurements are performed for low to moderate Reynolds numbers 1130?<?Re _b ?<?11320. Mean and turbulent velocity profiles are quantified at the tube exit and along the centerline from the tube exit down to 20 times the diameter D. Flow development is seen to be determined by the underlying jet structure at the tube exit which depends on Reynolds number, initial velocity statistics at the tube exit and the presence/absence of coherent structures. It is shown that the tube extension ratio L/D as well as the sharp edged abrupt contraction influence the initial jet structure at the tube exit. For both L/D ratios, the presence of the abrupt contraction results in transitional jet flow in the range 2000?<?Re _b ?<?4000 and in flow features associated with forced jets and high Reynolds numbers Re _b ?>?10⁴. The tube extension ratio L/D downstream from the abrupt contraction determines the shear layer roll up so that for L/D?=?1.2 flow visualisation suggests the occurrence of toroidal vortices for Re _b ?<?4000 whereas helical vortices are associated with the transitional regime for L/D?=?53.2. Found flow features are compared to features reported in literature for smooth contraction nozzles and long pipe flow.