Laminar-to-turbulent transition in pulsatile flow through a stenosis

Laminar-to-turbulent transition in pulsatile flow through a stenosis is studied by means of three-dimensional numerical simulations. The flow transition is associated with the occurrence of a flow instability initiating in the stenosis region. The instability is manifested by a three-dimensional symmetry-breaking and leads to asymmetric separation and intense swirling motion downstream of the stenosis. The above have profound effects on the wall shear stress (WSS). The simulations reveal that the asymmetric separation is extended several radii downstream of the stenosis with substantial WSS fluctuations, in both space and time, occurring in the poststenotic region.     

Microcontinuum model for pulsatile blood flow through a stenosed tube

The effects of polar nature of blood and pulsatility on flow through a stenosed tube have been analysed by assuming blood as a micropolar fluid. Linearized solutions of basic equations are obtained through consecutive applications of finite Hankel and Laplace transforms. The analytical expressions for axial and particle angular velocities, wall shear stress, resistance to flow and apparent viscosity have been obtained. The axial velocity profiles for Newtonian and micropolar fluids have been compared. The interesting observation of this analysis is velocity, in certain parts of cycle, for micropolar fluid is higher than Newtonain fluid. Variation of apparent viscosity eta a with tube radius shows both inverse Fahraeus-Lindqvist and Fahraeus-Lindqvist effects. Finally, the resistance to flow and wall shear stress for normal and diseased blood have been computed and compared.     

Phasic and spatial pressure measurements in a femoral artery branch model for pulsatile flow

Phasic and spatial time-averaged pressure distributions were measured in a 60-deg femoral artery branch model over a large range of branch flow ratios and at physiological Reynolds numbers of about 120 and 700. The results obtained with an in-vivo like flow wave form indicated spatial adverse time average pressure gradients in the branch vicinity which increased in magnitude with branch flow ratio, and the importance of the larger inertial effects at the higher Reynolds numbers. Pressure losses in the branch entrance region were relatively large, and corresponding flow resistances may limit branch flow, particularly at higher Reynolds numbers. The effect of branch flow was to reduce the pressure loss in the main lumen.     

Wall shear stress distribution in the human carotid siphon during pulsatile flow

Wall shear stress distribution in the carotid siphon, which is a multiple curved segment of the internal carotid artery, is investigated numerically under physiological flow conditions. The computer simulation of flow through the model segment is based on the time-dependent, three-dimensional Navier-Stokes equations, solved numerically with a finite element method. The study shows the behavior of the wall shear stress-vector field and identifies the zones of high and low wall shear stress values during the cardiac cycle.     

The flow of a viscous liquid in a converging tube

The problem of the viscous flow of an incompressible Newtonian liquid in a converging tapered tube has been solved in spherical polar coordinates. The method of the solution involves the Stokes' stream function and a technique introduced by Stokes in the study of a sphere oscillating in a fluid. The theory for the flow in a rigid tube includes: (1) the pulsatile flow with both radial and angular velocity components; (2) the steady state flow with both radial and angular velocity components and (3) the very slow steady state flow with only a radial velocity component present. For a tapered elastic tube, the velocity of the propagated pulse wave is determined. The solution given is in terms of the elastic constants of the system and the coordinates for this type of geometry. The pulse velocity is then related to the velocity in an elastic cylindrical tube with the necessary correction terms to account for the tapered tube. Supported in part by the American Heart Association (No. 62F4EG). This work was done during the tenure of an Established Investigatorship of the American Heart Association.     

Fluid particle motion and Lagrangian velocities for pulsatile flow through a femoral artery branch model

A flow visualization study using selective dye injection and frame by frame analysis of a movie provided qualitative and quantitative data on the motion of marked fluid particles in a 60 degree artery branch model for simulation of physiological femoral artery flow. Physical flow features observed included jetting of the branch flow into the main lumen during the brief reverse flow period, flow separation along the main lumen wall during the near zero flow phase of diastole when the core flow was in the downstream direction, and inference of flow separation conditions along the wall opposite the branch later in systole at higher branch flow ratios. There were many similarities between dye particle motions in pulsatile flow and the comparative steady flow observations.     

Numerical and experimental investigations of pulsatile blood flow pattern through a dysfunctional mechanical heart valve

Around 250,000 heart valve replacements are performed every year around the world. Due their higher durability, approximately 2/3 of these replacements use mechanical prosthetic heart valves (mainly bileaflet valves). Although very efficient, these valves can be subject to valve leaflet malfunctions. These malfunctions are usually the consequence of pannus ingrowth and/or thrombus formation and represent serious and potentially fatal complications. Hence, it is important to investigate the flow field downstream of a dysfunctional mechanical heart valve to better understand its impact on blood components (red blood cells, platelets and coagulation factors) and to improve the current diagnosis techniques. Therefore, the objective of this study will be to numerically and experimentally investigate the pulsatile turbulent flow downstream of a dysfunctional bileaflet mechanical heart valve in terms of velocity field, vortex formation and potential negative effect on blood components. The results show that the flow downstream of a dysfunctional valve was characterized by abnormally elevated velocities and shear stresses as well as large scale vortices. These characteristics can predispose to blood components damage. Furthermore, valve malfunction led to an underestimation of maximal transvalvular pressure gradient, using Doppler echocardiography, when compared to numerical results. This could be explained by the shifting of the maximal velocity towards the normally functioning leaflet. As a consequence, clinicians should try, when possible, to check the maximal velocity position not only at the central orifice but also through the lateral orifices. Finding the maximal velocity in the lateral orifice could be an indication of valve dysfunction.     

Vortex generation in pulsatile flow through arterial bifurcation models including the human carotid artery

Visualization experiments were performed to elucidate the complicated flow pattern in pulsatile flow through arterial bifurcations. Human common carotid arteries, which were made transparent, and glass-models simulating Y- and T-shaped bifurcations were used. Pulsatile flow with wave forms similar to those of arterial flow was generated with a piston pump, elastic tube, airchamber, and valves controlling the outflow resistance. Helically recirculating flow with a pattern similar to that of the horseshoe vortex produced around wall-based protuberances in circular tubes was observed in pulsatile flow through all the bifurcations used in the present study. This flow type, which we shall refer to as the horseshoe vortex, has also been demonstrated to occur at the human common carotid bifurcation in steady flow with Reynolds numbers above 100. Time-varying flows also produced the horseshoe vortex mostly during the decelerating phase. Fluid particles of dye solution approaching the bifurcation apex diverged, divided into two directions perpendicularly, and then showed helical motion representing the horseshoe vortex formation. While this helical flow was produced, the stagnation points appeared on the wall upstream of the apex. Their position was dependent upon the flow distribution ratio between the branches in the individual arteries. The region affected by the horseshoe vortex was smaller during pulsatile flow than during steady flow. Lowering the Reynolds number together with the Womersley number weakened the intensity of helical flow. A separation bubble, resulting from the divergence or wall roughness, was observed at the outer or inner wall of the branch vessels and made the flow more complicated.     

Developing pulsatile flow in a deployed coronary stent

A major consequence of stent implantation is restenosis that occurs due to neointimal formation. This patho-physiologic process of tissue growth may not be completely eliminated. Recent evidence suggests that there are several factors such as geometry and size of vessel, and stent design that alter hemodynamic parameters, including local wall shear stress distributions, all of which influence the restenosis process. The present three-dimensional analysis of developing pulsatile flow in a deployed coronary stent quantifies hemodynamic parameters and illustrates the changes in local wall shear stress distributions and their impact on restenosis. The present model evaluates the effect of entrance flow, where the stent is placed at the entrance region of a branched coronary artery. Stent geometry showed a complex three-dimensional variation of wall shear stress distributions within the stented region. Higher order of magnitude of wall shear stress of 530 dyn/cm2 is observed on the surface of cross-link intersections at the entrance of the stent. A low positive wall shear stress of 10 dyn/cm2 and a negative wall shear stress of -10 dyn/cm2 are seen at the immediate upstream and downstream regions of strut intersections, respectively. Modified oscillatory shear index is calculated which showed persistent recirculation at the downstream region of each strut intersection. The portions of the vessel where there is low and negative wall shear stress may represent locations of thrombus formation and platelet accumulation. The present results indicate that the immediate downstream regions of strut intersections are areas highly susceptible to restenosis, whereas a high shear stress at the strut intersection may cause platelet activation and free emboli formation.     

Developed pulsatile flow in a deployed coronary stent

The patho-physiologic process of restenosis and tissue growth may not be completely eliminated and is the primary concern of clinicians performing angioplasty and stent implantation procedures. Recent evidence suggests that the restenosis process is influenced by several factors: (1) geometry and size of vessel; (2) stent design; and (3) it's location that alter hemodynamic parameters, including local wall shear stress (WSS) distributions. The present three-dimensional (3D) analysis of pulsatile flow in a deployed coronary stent: (1) shows complex 3D variation of hemodynamic parameters; and (2) quantifies the changes in local WSS distributions for developed flow and compares with recently published WSS data for developing flow. Higher order of magnitude of WSS of 290 dyn/cm(2) is observed on the surface of cross-link intersections at the entrance of the stent for developed flow, which is about half of that for developing flow. Low WSS of 0.8 dyn/cm(2) and negative WSS of -8 dyn/cm(2) are seen at the immediate upstream and downstream regions of strut intersections. Persistent recirculation is observed at the downstream region of each strut cross-link and the regions of low and negative WSS may lead to patho-physiologic conditions near the stented region. The key finding of this study is that the location of stent in the coronary artery determines the developing or developed nature of the flow, which in turn, results in varied level of WSS.     

Incorporating the microbial loop in a simple plankton model

The microbial loop, which recycles nutrients in the upper layers of the ocean, is an integral part of plankton dynamics. The usual method for modelling the complex patterns involved has been to consider the ''Z'' in N/P/Z (nutrient/phytoplankton/zooplankton) models as containing all possible grazers on P and, implicitly, relegate the carnivorous metazoans to the loss term on Z. I propose the opposite approach (to define Z explicitly as the metazoans responsible for export fluxes) and to simulate the effects of the microbial loop implicitly in terms of grazing and excretion rates. The reasons for taking this alternative route are (i) the importance of copepods in the carbon/nitrogen flux from the euphotic zone to deeper water compared with (ii) the predominantly internal role of the microzooplankton in recycling nutrients; and (iii) the problems of sampling the microbial component, compared with sampling larger metazoans. Finally, there is the need to keep plankton models as simple as possible for later use in coupled physical/biological systems.     

An experimental study of pressure losses in pulsatile flows through rigid and pulsating stenosis

The time-dependent pressure curves of a pulsatile flow across rigid and pulsating stenoses were investigated experimentally in a laboratory simulator of the outflow tract of the heart right ventricle. The experiments were performed within the range of physiological conditions of frequency and flow rate. The experimental setup consisted of a closed flow system which was operated by a pulsatile pump, and a test chamber which enabled checking different modes of stenosis. Rigid constrictions were simulated by means of axisymmetric blunt-ended annular plugs with moderate-to-severe area reductions. The pulsating stenosis consisted of a short starling resistor device operated by a pulsating external pressure which was synchronized by the pulsatile flow. It was found that the shape of the time-dependent pressure curve upstream of the stenosis was different in the case of rigid stenosis than in the pulsating one. Potential clinical applications of the work may relate to diagnosis of the type of stenosis in the congenital heart disease known as Tetralogy of Fallot.     

Branch junctions and the flow of water through xylem in Douglas-fir and ponderosa pine stems

Water flowing through the xylem from the roots to the leaves of most plants must pass through junctions where branches have developed from the main stem. These junctions have been studied as both flow constrictions and components of a hydraulic segmentation mechanism to protect the main axes of the plant. The hydraulic nature of the branch junction also affects the degree to which branches interact and can respond to changes in flow to other branches. The junctions from shoots of two conifer species were studied, with particular emphasis on the coupling between the downstream branches. Flow was observed qualitatively by forcing stain through the junctions and the resulting patterns showed that flow into a branch was confined to just part of the subtending xylem until a considerable distance below the junction. Junctions were studied quantitatively by measuring flow rates in a branch before and after flow was stopped in an adjacent branch and by measuring the hydraulic resistance of the components of the junction. Following flow stoppage in the adjacent branch, flow into the remaining branch increased, but considerably less than predicted based on a simple resistance analogue for the branch junction that assumes the two branches are fully coupled. The branches downstream from a junction, therefore, appear to be limited in their interconnectedness and hence in their ability to interact.     

Interaction of peristaltic flow with pulsatile flow in a circular cylindrical tube

The effect of pulsatile flow on peristaltic transport in a circular cylindrical tube is analysed. The flow of a Newtonian viscous incompressible fluid in a flexible circular cylindrical tube on which an axisymmetric travelling sinusoidal wave is imposed, is considered. The initial flow in the tube is induced by an arbitrary periodic pressure gradient. A perturbation solution with amplitude ratio (wave amplitude/tube radius) as a parameter is obtained when the frequency of the travelling wave and that of the imposed pressure gradient are equal. The interaction effects of periodic wall induced flow and periodic pressure imposed flow are visualized through the presence of substantially different components of steady and higher harmonic oscillating flow in the first order flow solution. Numerical results show a strong variation of steady state velocity profiles with boundary wave number and Reynolds number and a strong phase shift behaviour of the flow in the radial direction.     

Waveform dependence of pulsatile flow in a stenosed channel

Bloodflow in arteries often shows a rich variety of vortical flows, which are dominated by the complex geometry of blood vessels, the dynamic pulsation of blood flow, and the complicated boundary conditions. With a two-dimensional model of unsteady flow in a stenosed channel, the pulsatile influence on such vortical fluid dynamics has been numerically studied in terms of waveform dependence on physiological pulsation. Results are presented for unsteady flows downstream of the stenosed portion with variation in the wavefiorms of systole and diastole. Overall, a train of propagating vortex waves is observed for all the cases, but it shows great sensitivity to the waveforms. The generation and development of the vortex waves may be linked to the presence of an adverse pressure gradient within a specific interval between two points of inflection of the systolic waveform. The adverse pressure gradient consists of a global pressure gradient that is found to be closely related to the dynamnics of' the pulsation, and a local pressure gradient, which is obsented to be dominated by the nonlinear vortex dynamics.     

Simple pulsatile flow in an artery with a constriction

A smooth isolated, axisymmetric occlusion in a straight vascular tube is a tractable problem for pulsatile flow calculations via finite-difference approximations to the Navier-Stokes equation. Steady flow depends on the Reynolds number and two geometric parameters which describe the stenosis. The mere addition of a simple harmonic to the mean flow adds two more parameters. One is the reduced frequency, or Strokes number, and the other epsilon, the ratio of unsteady to steady flux. After describing steady stenosis flow examples, the dynamic patterns of pulsatile flow are illustrated indicating the inadequacy of basing hypotheses of atherosclerosis on mean (steady) flow.     

Local overall volumetric gas-liquid mass transfer coefficients in gas-liquid-solid reversed flow jet loop bioreactor with a non-Newtonian fluid

The local overall volumetric gas-liquid mass transfer coefficients at the specified point in a gas-liquid-solid three-phase reversed flow jet loop bioreactor (JLB) with a non-Newtonian fluid was experimentally investigated by a transient gassing-in method. The effects of liquid jet flow rate, gas jet flow rate, particle density, particle diameter, solids loading, nozzle diameter and CMC concentration on the local overall volumetric gas-liquid mass transfer coefficient (K(L)a) profiles were discussed. It was observed that local overall K(L)a profiles in the three-phase reversed flow JLB with non-Newtonian fluid increased with the increase of gas jet flow rate, liquid jet flow rate, particle density and particle diameter, but decreased with the increase of the nozzle diameter and CMC concentration. The presence of solids at a low concentration increased the local overall K(L)a profiles, and the optimum of solids loading for a maximum profile of the local overall K(L)a was found to be 0.18x10(-3)m(3) corresponding to a solids volume fraction, varepsilon(S)=2.8%.     

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery.

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.     

In-vitro pulsatile flow visualization studies in a pulmonary artery model

In-vitro pulsatile flow visualization studies were conducted in an adult-sized pulmonary artery model to observe the effects of valvular pulmonic stenosis on the flow fields of the main, left and right pulmonary arteries. The flow patterns revealed that as the degree of stenosis increased, the jet-type flow created by the valve became narrower, and it impinged on the far (distal) wall of the left pulmonary artery further downstream from the junction of the bifurcation. This in turn led to larger regions of disturbed turbulent flow, as well as helical-type secondary flow motions in the left pulmonary artery, compared to the right pulmonary artery. The flow field in the main pulmonary artery also became more disturbed and turbulent, especially during peak systole and the deceleration phase. The flow visualization observations have been valuable in helping to conduct further quantitative studies such as pressure and velocity field mapping. Such studies are important to understanding the fluid mechanics characteristics of the main pulmonary artery and its two major branches.     

Two-equation turbulence modeling of pulsatile flow in a stenosed tube

The study of pulsatile flow in stenosed vessels is of particular importance because of its significance in relation to blood flow in human pathophysiology. To date, however, there have been few comprehensive publications detailing systematic numerical simulations of turbulent pulsatile flow through stenotic tubes evaluated against comparable experiments. In this paper, two-equation turbulence modeling has been explored for sinusoidally pulsatile flow in 75% and 90% area reduction stenosed vessels, which undergoes a transition from laminar to turbulent flow as well as relaminarization. Wilcox's standard k-omega model and a transitional variant of the same model are employed for the numerical simulations. Steady flow through the stenosed tubes was considered first to establish the grid resolution and the correct inlet conditions on the basis of comprehensive comparisons of the detailed velocity and turbulence fields to experimental data. Inlet conditions based on Womersley flow were imposed at the inlet for all pulsatile cases and the results were compared to experimental data from the literature. In general, the transitional version of the k-omega model is shown to give a better overall representation of both steady and pulsatile flow. The standard model consistently over predicts turbulence at and downstream of the stenosis, which leads to premature recovery of the flow. While the transitional model often under-predicts the magnitude of the turbulence, the trends are well-described and the velocity field is superior to that predicted using the standard model. On the basis of this study, there appears to be some promise for simulating physiological pulsatile flows using a relatively simple two-equation turbulence model.