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 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.     

Nonlinear dynamics and conformational changes of linear entangled polymers using the Rolie-Poly model with an “effective” maximum contour length

The extended Rouse-CCR tube model for linear entangled polymers recently proposed by the authors (Kabanemi and Hétu, J Non Newtonian-Fluid Mech 160:113–121, 2009), designed to capture the progressive changes in the average internal structure (kinked state) of polymer chains, is here used to analyze, by means of a time-dependent three-dimensional finite element method, chain segment dynamics, pressure drop, and stability of flow through a 4:1:4 constriction in a tube. The model predicts an enhancement of the pressure drop in the stretch-dominated flow regime, which is also observed experimentally. This excess pressure drop was not associated with the onset of flow instability. The model also predicts kinked configurations within chain segments in the entry section to the constriction tube, at the inception of flow, and prior to the development of upstream vortices. It is also shown how these kinked configurations within chain segments influence pressure drop transients.     

动脉局部狭窄时脉动流的有限元分析

本文利用有限元方法研究动脉局部狭窄下的脉动流流场,重点考查在50％与80％面积狭窄下的速度分布、压力分布、壁面剪应力分布及流动分离情况。几何形状及边界条件均模拟相应的脉动流实验模型。采用测得的随时间变化的速度分布作为入口端条件,并利用罚函数和逆风格式等计算技巧得出了光滑的与实验基本相符的速度、压力波形。本文讨论了不同狭窄下速度、压力、壁面剪应力的分布形态,给出了脉动流中狭窄处局部流动分离的间歇性变化规律,并结合实验与临床应用进行了讨论。     

Pressure weighted upwinding for flow induced force predictions: application to iced surfaces

This paper addresses two main topics, namely the development of a pressure‐weighted upwinding method and its application to flow induced forces on iced cylinders. Although the near‐wall convective upwinding exhibits special applicability to iced surfaces, its capabilities extend more generally to other applications. By fully linking pressure and velocity at a sub‐element level near the wall, a higher order accuracy can be obtained. Also, a non‐physical de‐coupling between pressure and velocity can be prevented. The method is developed under the context of a control‐volume‐based finite element method for 2‐D, incompressible flows. Drag and lift coefficients are predicted, based on the pressure weighted upwinding near the wall. The numerical predictions are successfully compared against experimental data, including flow induced forces on iced cables. Copyright © 2004 John Wiley & Sons, Ltd.     

Adaptive Lagrangian–Eulerian computation of propagation and rupture of a liquid plug in a tube

Liquid plug propagation and rupture occurring in lung airways can have a detrimental effect on epithelial cells. In this study, a numerical simulation of a liquid plug in an infinite tube is conducted using an Eulerian–Lagrangian approach and the continuous interface method. A reconstruction scheme is developed to allow topological changes during plug rupture by altering the connectivity information about the interface mesh. Results prior to the rupture are in reasonable agreement with the study of Fujioka et al. in which a Lagrangian method is used. For unity non‐dimensional pressure drop and a Laplace number of 1000, rupture time is shown to be delayed as the initial precursor film thickness increases and rupture is not expected for thicknesses larger than 0.10 of tube radius. During the plug rupture process, a sudden increase of mechanical stresses on the tube wall is recorded, which can cause tissue damage. The peak values of those stresses increase as the initial precursor film thickness is reduced. After rupture, the peaks in mechanical stresses decrease in magnitude as the plug vanishes and the flow approaches a fully developed behavior. Increasing initial pressure drop is shown to linearly increase maximum variations in wall pressure and shear stress. Decreasing the pressure drop and increasing the Laplace number appear to delay rupture because it takes longer for a fluid with large inertial forces to respond to the small pressure drop. Copyright © 2010 John Wiley & Sons, Ltd.     

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 new positive‐definite regularization of incompressible Navier–Stokes equations discretized with Q1/P0 finite element

A new regularization method is proposed for the Galerkin approximation of the incompressible Navier–Stokes equations with Q1/P0 element, by newly introducing a square‐type linear form into the variational divergence‐free constraint regularized with the global pressure jump (GPJ) method. The addition of the square‐type linear form is intended to eliminate the hydrostatic pressure mode appearing in confined flows, and to make the discretized matrix positive definite and then non‐singular without the pressure pegging trick. Effects of the free parameters for the regularization on the solutions are numerically examined with a 2‐D driven cavity flow problem. Furthermore, the convergences in the conjugate gradient iteration for the solution of the pressure Poisson equation are compared among the mixed method, the GPJ method and the present method for both leaky and non‐leaky 3‐D driven cavity flows. Finally, the non‐leaky 3‐D cavity flows at different Re numbers are solved to compare with the literature data and to demonstrate the accuracy of the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

A Parallel 3D Unsteady Flow Analysis in an Asymmetrically Constricted Vessel

In this study, parallel computation of unsteady incompressible flow in an asymmetrically constricted 3D vessel has been presented. A time accurate cell centered finite volume method (FVM) in conjunction with pseudo-compressibility technique and Roe's flux difference splitting of nonlinear terms has been employed for solving the Navier-Stokes (NS) equations on the multiple instruction multiple data (MIMD) machine VPP700. The influence of Reynolds' number ( Re ) and the Strouhal number ( St ) on flow dynamic factors like wall pressure (WP), wall shear stress (WSS), central axis velocity (CAV), etc., have been analyzed. Three-dimensional (3D) features in the formation and detachment of separation zones, which are sensitive to both Re and St have been noticed on the diverging wall of the constriction.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of laminar flow and heat transfer over banks of staggered cylinders

A study is conducted to investigate forced convective flow and heat transfer over a bank of staggered cylinders. Using a novel numerical formulation based on a non‐orthogonal collocated grid in a physical plane, the effects of Reynolds number and cylinder spacing on the flow and heat transfer behaviour are systematically studied. It is observed that both the Reynolds number and cylinder spacing influence the recirculatory vortex formation and growth in the region between the cylinders; in turn, the rates of heat transfer between the fluid and the staggered cylinders are affected. As the cylinder spacing decreases, the size and length of eddies reduce. For sufficiently small spacings, eddy formation is completely suppressed even at high Reynolds number. Pressure drop and Nusselt number predictions based on numerical study are in excellent agreement with available correlations. The study provides useful insight on the detailed flow and heat transfer phenomena for the case of a bank of staggered cylinders. Copyright © 2002 John Wiley & Sons, Ltd.     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

Topology optimization of microfluidic mixers

This paper demonstrates the application of the topology optimization method as a general and systematic approach for microfluidic mixer design. The mixing process is modeled as convection dominated transport in low Reynolds number incompressible flow. The mixer performance is maximized by altering the layout of flow/non‐flow regions subject to a constraint on the pressure drop between inlet and outlet. For a square cross‐sectioned pipe the mixing is increased by 70% compared with a straight pipe at the cost of a 2.5 fold increase in pressure drop. Another example where only the bottom profile of the channel is a design domain results in intricate herring bone patterns that confirm findings from the literature. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Blockage of constrictions by particles in fluid–solid transport

Blockage is an important phenomenon in particulate flow. Work was undertaken to provide a better understanding of key hydrodynamic multiphase flow factors which cause, or contribute to, stalling and blockage in particulate feeding systems such as those used for feeding biomass into reactors. Rubber and plastic particles were hydraulically conveyed along a horizontal rectangular duct leading to constrictions of different geometries. Experimental results showed that large size, irregular shape, high volumetric concentrations of particles, small constriction dimensions and particle compressibility all increased the likelihood of blockage. Reynolds number also had a significant effect on particle behaviour and blockage propensity. The pressure drop needed to break a blockage is also considered, based on a simple horizontal packed bed model.     

Linear and non‐linear simulations of feedback control in plane Poiseuille flow

This paper examines the performance of optimal linear quadratic state and output feedback controllers in stabilizing two‐dimensional perturbations in a plane Poiseuille flow. The synthesis of the controllers is based on a linearized model of the flow using a new set of interpolating polynomials in the wall‐normal direction, which automatically satisfy the homogeneous Dirichlet and Neumann boundary conditions at the walls and eliminate spurious eigenvalues. The controllers are implemented into a non‐linear Navier–Stokes solver, which is modified to compute the evolution of the flow perturbations. Two cases are examined, one with small initial disturbances that do not violate the linearity assumptions and the other with much larger disturbances that trigger the non‐linear convection terms. For the smallest disturbances, the solver accurately reproduced the results of the linear simulations of open‐ and closed‐loop systems. The simulations for the larger disturbances without control showed a rapid initial growth but the flow soon reached a saturated state in agreement with previous findings in the literature. The large initial growth is a consequence of the non‐normal nature of the system dynamics. The state feedback and output feedback controllers were able to reduce significantly the perturbation energy. For the larger disturbances, the energy calculated from the state variables is well below the energy evaluated by direct integration of the velocity field. This is probably due to the non‐linear terms transferring energy to harmonics of the considered wavenumber, which are not sensed by the linear controller. Copyright © 2008 John Wiley & Sons, Ltd.     

Multiscale computations of thin films in multiphase flows

In multiphase flows thin films are often encountered when fluid masses collide. These films can become very thin and in direct numerical simulations (DNS) it is often impractical to resolve their thickness fully, even with adaptive grid refining. Here we examine the collision of a fluid drop with a wall and develop a multiscale approach to compute the flow in the film between the drop and the wall. By using a semi-analytical model for the flow in the film we capture the evolution of films thinner than the grid spacing reasonably well.     

Unsteady Two-Dimensional Blood Flow in Porous Artery with Multi-Irregular Stenoses

The flow characteristics of an unsteady axisymmetric two-dimensional (2D) blood flow in a diseased porous arterial segment with flexible walls are investigated. The arterial walls mimic the irregular constrictions whereas the lumen containing the thrombus, cholesterol, and fatty plaques represents the porous medium. The governing equations with appropriate initial and boundary conditions are solved numerically using MAC method. The discretization is done on staggered grid with non-uniform grid size and pressure-poisson equation is solved following SOR method. The pressure and velocity corrections are made cyclically until the steady state is achieved. It is observed that for decreasing permeability, flow is highly decelerated while pressure drop and wall shear stress increases. The separation zones and re-circulation regions are found for severe stenoses. Flow separation and re-circulation diminishes for decreasing permeability of the porous medium. Comparisons are provided with published experimental and numerical results.     

Two‐phase pressure drop caused by sudden flow area contraction/expansion in small circular pipes

Theoretical studies have been made to determine the pressure drops caused by abrupt flow area expansion/contraction in small circular pipes for two‐phase flow of air and water mixtures at room temperature and near atmospheric pressure. Two‐phase computational fluid dynamics (CFD) calculations, using Eulerian–Eulerian model (with the air phase being compressible for pipe contraction case) are employed to calculate the pressure drop across sudden expansion and contraction. The pressure drop is determined by extrapolating the computed pressure profiles upstream and downstream of the expansion/contraction. The larger and smaller tube diameters are 1.6 and 0.84 mm, respectively. Computations have been performed with single‐phase water and air, and two‐phase mixtures in a range of Reynolds number (considering all‐liquid flow) from 1000 to 12 000 and flow quality from 1.2 × 10^?3 to 1.6 × 10^?2. The numerical results are validated against experimental data from the literature and are found to be in good agreement. The expansion and contraction loss coefficients are found to be different for single‐phase flow of air and water, and they agreed reasonably well with the commonly used theoretical predictions. Based on the numerical results as well as experimental data, correlations are developed for two‐phase flow pressure drops caused by the flow area contraction as well as expansion. Copyright © 2010 John Wiley & Sons, Ltd.