Numerical simulation of turbulent flow through series stenoses

The flow fields in the neighbourhoods of series vascular stenoses are studied numerically for the Reynolds numbers from 100 to 4000, diameter constriction ratios of 0.2–0.6 and spacing ratios of 1, 2, 3, 4 and ∞. In this study, it has been further verified that in the laminar flow region, the numerical predictions by k–ω turbulence model matched those by the laminar‐flow modelling very well. This suggests that the k–ω turbulence model is capable of the prediction of the laminar flow as well as the prediction of the turbulent stenotic flow with good accuracy. The extent of the spreading of the recirculation region from the first stenosis and its effects on the flow field downstream of the second stenosis depend on the stenosis spacing ratio, constriction ratio and the Reynolds number. For c₁ = 0.5 with c₂ ≤ c₁, the peak value of wall vorticity generated by the second stenosis is always less than that generated by the first stenosis. However, the maximum centreline velocity and turbulence intensity at the second stenosis are higher than those at the first stenosis. In contrast, for c₁ = 0.5 with c₂ = 0.6, the maximum values at the second stenosis are much higher than those at the first stenosis whether for centreline velocity and turbulence intensity or for wall vorticity. The peak values of the wall vorticity and the centreline disturbance intensity both grow up with the Reynolds number increasing. The present study shows that the more stenoses can result in a lower critical Reynolds number that means an earlier occurrence of turbulence for the stenotic flows. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

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.     

Numerical study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

Upstream roughness and Reynolds number effects on turbulent flow structure over forward facing step

This paper presents results of experiments conducted to investigate the effects of Reynolds number and upstream wall roughness on the turbulence structure in the recirculation and recovery regions of a smooth forward facing step. A reference smooth upstream wall and a rough upstream wall made from sand gains were studied. For the smooth upstream wall, experiments were conducted at Reynolds number based on the freestream velocity and step height (h), Re_h = 4940, 8400 and 8650. The rough wall experiments was performed at Re_h = 5100, 8200 and 8600 to closely match the corresponding Re_h experiment over the smooth wall. The reattachment lengths in the smooth wall experiments were L_r/h ≈ 2.2, but upstream roughness significantly reduced these values to L_r/h ≈ 1.3. The integral scales within the recirculation bubbles were independent of upstream roughness and Reynolds number; however, upstream roughness significantly increased the spatial coherence and integral scales outside the recirculation bubbles and in the recovery region. Irrespective of the upstream wall condition, the redeveloping boundary layer recovered at 25h from reattachment.     

An experimental investigation of the effect of a constriction on turbulent pipe flow

Measurements have been made of the distributions of the mean-velocity and the axial turbulence velocity component in a cross-section of a circular tube at various distances downstream from a number of different constrictions. Also spectral distributions of the turbulence velocity have been measured in the axis of the tube and in a point very close to the wall. The constrictions had a contraction ratio of 0.25 except one which had a ratio of 0.5. One of the constrictions was made of a thin rubber hose. When for this constriction the contraction ratio was reduced to a value smaller than 0.25, self-excited vibrations of the hose took place, producing an oscillating flow of the air in the tube. The Reynoldsnumber was kept at roughly 5,000. As could be expected, after 40 tube diameter distance downstream from the constrictions an almost complete recovery of the disturbed turbulent flow, as far as the distributions of the mean velocity and relative turbulence intensity are concerned, was obtained. Depending on the shape of the constriction even a shorter distance appeared to be sufficient. The flexible constriction then was in the non-vibrating condition. However, the spectral distributions showed in some cases still a difference with the undisturbed case, in particular in the low frequency range. If the flexible constriction was vibrating, the induced oscillations of the flow which showed up as discrete peaks in the spectral distributions, persisted over the entire length of the tube, again as expected.     

Near-wall investigation of backward-facing step flows

The electrodiffusion technique has been used to investigate reattaching and recirculating flows behind a backward-facing step. The instantaneous wall shear rate vectors were determined using the current signals provided by a three-segment electrodiffusion probe. The near-wall extents of two counter-rotating recirculation zones located behind the step were determined under turbulent flow conditions in a water channel. The near-wall flow inside these recirculation zones was found to be very unsteady, with strong low-frequency fluctuations. The streamwise profiles of the wall shear stress were measured at several values of the Reynolds number and a high level of skin friction was obtained in the reverse-flow region. The strong dependence of the peak value of skin friction on the Reynolds number confirms the viscous-dominated character of the reverse flow appearing inside the recirculation zone. Received: 22 May 2000/Accepted: 28 March 2001     

NUMERICAL STUDIES OF TRANSITIONAL TURBULENT PULSATILE FLOW IN PIPES WITH RING–TYPE CONSTRICTIONS

Pulsatile flows in the vicinity of mechanical ring-type constrictions in pipes were studied for transitional turbulent flow with a Reynolds number (Re) of the order of 10⁴. The Womersley number (Nw) is in the range 30–50, with a corresponding Strouhal number (St) range of 0·0143–0·0398. The pulsatile flows considered are a pure sinusoidal flow, a physiological flow and an experimental pulsatile flow profile for mechanical aortic valve flow simulations. Transitional laminar and turbulent flow characteristics in an alternating manner within the pulsatile flow fields were studied numerically. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity and turbulent shear stress are smaller during the acceleration phase than during the deceleration period. Various parametric equations have been formulated through numerical experimentation to better describe the relationships between the instantaneous flow rate (Q), the pressure loss (ΔP), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max) for turbulent pulsatile flow in the vicinity of constrictions in the vascular tube. An elliptic relationship has been found to exist between the instantaneous flow rate and the instantaneous pressure gradient. Other linear and quadratic relations between various flow parameters were also obtained.     

Numerical simulation of two‐dimensional laminar incompressible offset jet flows

Two‐dimensional transient laminar incompressible offset jet is simulated numerically to gain insight into convective recirculation and flow processes induced by an offset jet. The behaviour of the jet with respect to offset ratio (OR) and Reynolds number (Re) are described in detail. The transient development of the velocity is simulated for various regions: recirculation, impingement and wall jet development. It is found that the reattachment length is dependent on both Re and OR for the range considered. Simulations are made to show the effect of entrainment on recirculation eddy. A detailed study of u velocity decay is reported. The decay rate of horizontal velocity component (u) is linear in impingement region. It is found that at high OR, velocity decay depends on Re only. Velocity profile in the wall jet region shows good agreement with experimental as well as similarity solutions. It is found that the effect of Re and OR are significant to bottom wall vorticity up to impingement region. Far downstream bottom wall vorticity is independent of OR. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical study of physiological turbulent flows through series arterial stenoses

A numerical investigation on the characteristics of transitional turbulent flow over series bell‐shape stenoses for a physiological pulsatile flow is presented in the present study. The flow behaviours for the physiological pulsatile flow are studied by considering the effects of the Reynolds number, Womersley number, constriction ratio and spacing ratio of the stenoses on the pulsatile turbulent flow fields. Especially, the mutual influences between the double stenoses under different flow conditions are considered. The numerical results show that the variation of these flow parameters puts significant impacts on the flow developments in the arteries with series stenoses. The double stenoses lead to the higher peak turbulence disturbance and the greater area with comparatively high turbulence intensity distal to the stenoses in comparison with the single stenosis. The analysis shows that for the physiological pulsatile flow, the downstream stenosis usually does not have perceptible influences on the upstream flow fields. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of Reynolds number on physiological‐type pulsatile flows in a pipe with ring‐type constrictions

The effects of Reynolds number on the physiological‐type of laminar pulsatile flow fields within the vicinity of mechanical ring‐type constriction in small pipes were studied numerically. The parameters considered are: the Reynolds number (Re) in the range of 50–1500; Strouhal number (St) in the range of 0.00156–3.98; Womersley number (Nw) from 0.0 to 50.0. The pulsatile flows considered were physiological‐type of simulated flows. Within a pulsating cycle, detailed flow characteristics were studied through the pulsating contours of streamline (ψ), vorticity (Ω), shear stress (τ) and isobar. The relations between the instantaneous flow rate (Q) and instantaneous pressure gradients (dp/dz) are observed to be elliptic. The relations between the instantaneous flow rate (Q) and pressure loss (P_loss) are quadratic. Linear relations were observed between the instantaneous flow rate (Q) and the maximum velocity, maximum vorticity and maximum shear stress. The Reynolds number of the flow in a pulsating cycle was found to have significant effects on the recirculation length and the pressure gradient within the pulsatile flow regime. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Experimental study and large eddy simulation for the turbulent flow around four cylinders in an in-line square configuration

The turbulent flows around four cylinders in an in-line square configuration with different spacing ratios of 1.5, 2.5, 3.5 and 5.0 have been investigated experimentally at subcritical Reynolds numbers from 11,000 to 20,000. The mean and fluctuating velocity distributions were obtained using the laser Doppler anemometry (LDA) measurement. The digital particle image velocimetry (DPIV) was employed to characterize the full field vorticity and velocity distributions as well as other turbulent quantities. The experimental study indicated that several distinct flow patterns exist depending on the spacing ratio and subcritical Reynolds number for turbulent flow. The three-dimensional numerical simulations were also carried out using the large eddy simulation (LES) at Reynolds number of 15,000 with the spacing ratio of 1.5 and 3.5. The results show that the LES numerical predictions are in good agreement with the experimental measurements. Therefore, the three-dimensional vortex structures and the full field instantaneous and mean quantities of the flow field such as velocity field, vorticity field, etc., which are very difficult to obtain experimentally, can be extracted from the simulation results for the deepening of our understanding on the complex flow phenomena around four cylinders in in-line configuration.     

An experimental investigation of heat transfer to pulsating pipe air flow with different amplitudes

Heat transfer characteristics to both laminar and turbulent pulsating pipe flows under different conditions of Reynolds number, pulsation frequency, pulsator location and tube diameter were experimentally investigated. The tube wall of uniform heat flux condition was considered for both cases. Reynolds number varied from 750 to 12,320 while the frequency of pulsation ranged from 1 to 10 Hz. With locating the pulsator upstream of the inlet of the test section tube, results showed an increase in heat transfer rate due to pulsation by as much as 30% with flow Reynolds number of 1,643 and pulsation frequency of 1 Hz, depending on the upstream location of the pulsator valve. Closer the valve to the tested section inlet, the better improvement in the heat transfer coefficient is achieved. Upon comparing the heat transfer results of the upstream and the downstream pulsation, at Reynolds number of 1,366 and 1,643, low values of the relative mean Nusselt number were obtained with the upstream pulsation. Comparing the heat transfer results of the two studied test sections tubes for Reynolds number range from 8,000 to 12,000 and pulsation frequency range from 1.0 to 10 Hz showed that more improvement in heat transfer rate was observed with a larger tube diameter. For Reynolds number ranging from 8,000 to 12,000 and pulsation frequency of 10 Hz, an improvement in the relative mean Nusselt number of about 50% was obtained at Reynolds number of 8,000 for the large test section diameter of 50 mm. While, for the small test section diameter of 15 mm, at same conditions of Reynolds number and frequency, a reduction in the relative mean Nusselt number of up to 10% was obtained.     

Skin-friction drag reduction in turbulent channel flow based on streamwise shear control

It is known that stretching and intensification of a hairpin vortex by mean shear play an important role to create a hairpin vortex packet, which generates the large Reynolds shear stress associated with skin-friction drag in wall-bounded turbulent flows. In order to suppress the mean shear at the wall for high efficient drag reduction (DR), in the present study, we explore an active flow control concept using streamwise shear control (SSC) at the wall. The longitudinal control surface is periodically spanwise-arranged with no-control surface while varying the structural spacing, and an amplitude parameter for imposing the strength of the actuating streamwise velocity at the wall is introduced to further enhance the skin-friction DR. Significant DR is observed with an increase in the two parameters with an accompanying reduction of the Reynolds stresses and vorticity fluctuations, although a further increase in the parameters amplifies the turbulence activity in the near-wall region. In order to study the direct relationship between turbulent vortical structures and DR under the SSC, temporal evolution with initial eddies extracted by conditional averages for Reynolds-stress-maximizing Q2 events are examined. It is shown that the generation of new vortices is dramatically inhibited with an increase in the parameters throughout the flow, causing fewer vortices to be generated under the control. However, when the structural spacing is sufficiently large, the generation of new vortex is not suppressed over the no-control surface in the near-wall region, resulting in an increase of the second- and fourth-quadrant Reynolds shear stresses. Although strong actuating velocity intensifies the near-wall turbulence, the increase in the turbulence activity is attributed to the generation of counter-clockwise near-wall vortices by the increased vortex transport.     

Pulsatile flows and wall-shear stresses in models simulating normal and stenosed aortic arches

Pulsatile aqueous glycerol solution flows in the models simulating normal and stenosed human aortic arches are measured by means of particle image velocimetry. Three transparent models were used: normal, 25% stenosed, and 50% stenosed aortic arches. The Womersley parameter, Dean number, and time-averaged Reynolds number are 17.31, 725, and 1,081, respectively. The Reynolds numbers based on the peak velocities of the normal, 25% stenosed, and 50% stenosed aortic arches are 2,484, 3,456, and 3,931, respectively. The study presents the temporal/spatial evolution processes of the flow pattern, velocity distribution, and wall-shear stress during the systolic and diastolic phases. It is found that the flow pattern evolving in the central plane of normal and stenosed aortic arches exhibits (1) a separation bubble around the inner arch, (2) a recirculation vortex around the outer arch wall upstream of the junction of the brachiocephalic artery, (3) an accelerated main stream around the outer arch wall near the junctions of the left carotid and the left subclavian arteries, and (4) the vortices around the entrances of the three main branches. The study identifies and discusses the reasons for the flow physics’ contribution to the formation of these features. The oscillating wall-shear stress distributions are closely related to the featured flow structures. On the outer wall of normal and slightly stenosed aortas, large wall-shear stresses appear in the regions upstream of the junction of the brachiocephalic artery as well as the corner near the junctions of the left carotid artery and the left subclavian artery. On the inner wall, the largest wall-shear stress appears in the region where the boundary layer separates.     

Finite difference and finite element methods for mhd channel flows

A Galerkin finite element method and two finite difference techniques of the control volume variety have been used to study magnetohydrodynamic channel flows as a function of the Reynolds number, interaction parameter, electrode length and wall conductivity. The finite element and finite difference formulations use unequally spaced grids to accurately resolve the flow field near the channel wall and electrode edges where steep flow gradients are expected. It is shown that the axial velocity profiles are distorted into M-shapes by the applied electromagnetic field and that the distortion increases as the Reynolds number, interaction parameter and electrode length are increased. It is also shown that the finite element method predicts larger electromagnetic pinch effects at the electrode entrance and exit and larger pressure rises along the electrodes than the primitive-variable and streamfunction–vorticity finite difference formulations. However, the primitive-variable formulation predicts steeper axial velocity gradients at the channel walls and lower axial velocities at the channel centreline than the streamfunction–vorticity finite difference and the finite element methods. The differences between the results of the finite difference and finite element methods are attributed to the different grids used in the calculations and to the methods used to evaluate the pressure field. In particular, the computation of the velocity field from the streamfunction–vorticity formulation introduces computational noise, which is somewhat smoothed out when the pressure field is calculated by integrating the Navier–Stokes equations. It is also shown that the wall electric potential increases as the wall conductivity increases and that, at sufficiently high interaction parameters, recirculation zones may be created at the channel centreline, whereas the flow near the wall may show jet-like characteristics.     

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.