Computation of steady flow in a two-dimensional arterial model

The Navier-Stokes equations are solved numerically for steady flow through a double-branched two-dimensional section of a three-dimensional model of a canine aorta for which experimental data is available. The numerical scheme involves transforming the physical coordinates to a curvilinear boundary-fitted coordinate system and performing finite-difference computations in the transformed system. Shear stress at the wall is calculated for a Reynolds number of a 1,000 with branch-to-main aortic flow rate ratios as a parameter. The results are compared with the aforementioned experimental data and show reasonable qualitative agreement.     

Disorder distal to modeled stenoses in steady and pulsatile flow

Measurements of the velocity and energy spectra were made in the distal region of modeled stenoses in a rigid tube with both steady and pulsatile water flows. Reynolds numbers of 318–2540 and a pulsatile flow frequency parameter of 15 were employed. The effects of the degree of stenosis, the stenosis geometry and the presence or absence of the downstream confining wall on the development of flow disturbances were investigated. Visualization of the distal flow patterns in stenotic and free jets illustrated the existence of complex fields which included vortex shedding, highly turbulent regions, and recirculation zones. Significant flow disorder was created by a mild stenosis in pulsatile, but not in steady, flow. Nondimensionalization employing the stenosis diameter and flow velocity in the throat of the constriction correlates the vortex shedding frequency and energy spectra within a limited postestenotic region.     

Comparison of steady and pulsatile flow near the ventral and dorsal walls of casts of human aortic bifurcations

Steady and pulsatile flows were passed through casts of human aortic bifurcations and, by means of a laser Doppler anemometer, fluid velocities were measured at selected sites near the ventral and dorsal walls. At these sites, in the vicinity of the bifurcation, the influence of secondary flow is significant and therefore an appreciation of the phasic variation of secondary flow patterns is important. Results are presented comparing the flow direction in both steady and pulsatile flow at sites in three casts. The common features of the flow at these sites were the persistence of the flow direction during the accelerating and decelerating phases of the pulsatile cycle, and the consistently smaller angle (measured from the inlet centerline) of the pulsatile flow direction as compared to the angle of the flow direction in steady flow.     

A turbulence model for pulsatile arterial flows

Difficulties in predicting the behavior of some high Reynolds number flows in the circulatory system stem in part from the severe requirements placed on the turbulence model chosen to close the time-averaged equations of fluid motion. In particular, the successful turbulence model is required to (a) correctly capture the "nonequilibrium" effects wrought by the interactions of the organized mean-flow unsteadiness with the random turbulence, (b) correctly reproduce the effects of the laminar-turbulent transitional behavior that occurs at various phases of the cardiac cycle, and (c) yield good predictions of the near-wall flow behavior in conditions where the universal logarithmic law of the wall is known to be not valid. These requirements are not immediately met by standard models of turbulence that have been developed largely with reference to data from steady, fully turbulent flows in approximate local equilibrium. The purpose of this paper is to report on the development of a turbulence model suited for use in arterial flows. The model is of the two-equation eddy-viscosity variety with dependent variables that are zero-valued at a solid wall and vary linearly with distance from it. The effects of transition are introduced by coupling this model to the local value of the intermittency and obtaining the latter from the solution of a modeled transport equation. Comparisons with measurements obtained in oscillatory transitional flows in circular tubes show that the model produces substantial improvements over existing closures. Further pulsatile-flow predictions, driven by a mean-flow wave form obtained in a diseased human carotid artery, indicate that the intermittency-modified model yields much reduced levels of wall shear stress compared to the original, unmodified model. This result, which is attributed to the rapid growth in the thickness of the viscous sublayer arising from the severe acceleration of systole, argues in favor of the use of the model for the prediction of arterial flows.     

Gas transport in serpentine microporous tubes under steady and pulsatile blood flow conditions

A serpentine gas exchange unit was built with cylindrical tubular microporous membranes featuring periodic arcs with a fixed curvature ratio (ratio of tube radius to radius of curvature) of 1/14 and circular angles between 30 and 360 deg. Oxygen transfer was measured under steady and pulsatile blood flow conditions in vitro and ex vivo to assess the design features which most effectively augment gas transfer. Under steady blood flow conditions, oxygen transfer increased with circular angles beyond 70 deg. Under pulsatile conditions, a wide range of geometrical and fluid mechanical parameters could be combined to enhance gas transfer performance, which eventually depended upon the secondary Reynolds number and the Womersley parameter.     

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.     

Structure of pulsatile flow in a model of elastic cerebral aneurysm

This study examines the effect of aneurysmal wall elasticity on the structure of flow within an elastic aneurysm during pulsatile flow. We visualized flow structure in a model of an elastic saccular aneurysm located at the bifurcation of the anterior cerebral artery and extending to the anterior communicating artery, and measured changes in the diameter of the aneurysm wall during pulsatile flow using particle imaging velocimetry (PIV). We similarly measured these features during steady flow by PIV and found that dilation of the aneurysmal wall absorbed the dynamic energy within the aneurysm. Accordingly, aneurysm wall elasticity functions as a biocompatible reaction that relieves wall shear stress acting on the vascular wall during pulsatile flow, and should thus inhibit the development and rupture of an aneurysm.     

Comparison of LES of steady transitional flow in an idealized stenosed axisymmetric artery model with a RANS transitional model

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y(+) less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y(+) less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.     

A model of pulsatile flow in a uniform deformable vessel.

Simulations of blood flow in natural and artificial conduits usually require large computers for numerical solution of the Navier-Stokes equations. Often, physical insight into the fluid dynamics is lost when the solution is purely numerical. An alternative to solving the most general form of the Navier-Stokes equations is described here, wherein a functional form of the solution is assumed in order to simplify the required computations. The assumed forms for the axial pressure gradient and velocity profile are chosen such that conservation of mass is satisfied for fully established pulsatile flow in a straight, deformable vessel. The resulting equations are cast in finite-difference form and solved explicitly. Results for the limiting cases of rigid wall and zero applied pressure are found to be in good agreement with analytical solutions. Comparison with the experimental results of Klanchar et al. [Circ. Res. 66, 1624-1635 (1990]) also shows good agreement. Application of the model to realistic physiological parameter values provides insight as to the influence of the pulsatile nature of the flow field on wall shear development in the presence of a moving wall boundary. Specifically, the model illustrates the dependence of flow rate and shear rate on the amplitude of the vessel wall motion and the phase difference between the applied pressure difference and the oscillations of the vessel radius. The present model can serve as a useful tool for experimentalists interested in quantifying the magnitude and character of velocity profiles and shearing forces in natural and artificial biologic conduits.     

Differential effects of pulsatile versus steady flow on coronary endothelial membrane potential

Steady shear stress stimulates transient hyperpolarization coupled to calcium-sensitive potassium (KCa) channels and sustained depolarization linked to chloride-selective channels. Physiological flow is pulsatile not static, and whereas in vivo data suggest phasic shear stress may preferentially activate KCa channels, its differential effects on both currents remain largely unknown. To determine this interaction, coronary endothelial cells were cultured in glass capillary flow tubes, loaded with the voltage-sensitive dye bis-(1,3-dibutylbarbituric acid)trimethine oxonol, and exposed to constant or pulsatile shear stress. The latter was generated by a custom servoperfusion system employing physiological pressure and flow waveforms. Steady shear induced a sustained depolarization inhibited by the Cl- channel blocker DIDS. Even after exposure to steady flow, subsequent transition to pulsatile shear stress further stimulated DIDS-sensitive depolarization. DIDS pretreatment "unmasked" a pulsatile flow-induced hyperpolarization of which magnitude was further enhanced by nifedipine, which augments epoxygenase synthesis. Pulse-shear hyperpolarization was fully blocked by KCa channel inhibition (charybdotoxin + apamin), although these agents had no influence on membrane potential altered by steady flow. Thus KCa-dependent hyperpolarization is preferentially stimulated by pulsatile over steady flow, whereas both can stimulate Cl--dependent depolarization. This supports studies showing greater potency of pulsatile flow for triggering KCa-dependent vasorelaxation.     

Low Reynolds number steady state flow through a branching network of rigid vessels II. A finite element mixture model

This research aims at formulating and verifying a finite element mixture formulation for blood perfusion. The equations derived in a companion paper [3] are discretized according to the Galerkin method. A flow experiment in a rigid model of a vascular tree of about 500 vessels is performed in order to verify the finite element mixture formulation. Although the comparison of numerical results and experimental measurements is not conclusive as far as the validity of the theory is concerned, the results do suggest that the finite element model has predictive power in the case of low Reynolds number steady-state flow of a Newtonian fluid in a rigid vascular tree.     

Finite element simulation of pulsatile flow through arterial stenosis.

The problem of blood flow through a stenosis is solved using the incompressible Navier-Stokes equations in a rigid circular tube presenting a partial occlusion. Calculations are based on a Galerkin finite element method. The time marching scheme employs a predictor-corrector technique using a variable time step. Results are obtained for steady and physiological pulsatile flows. Computational experiments analyse the effect of varying the degree of stenosis, the stricture length, the Reynolds number and Womersley number. The method gives results which agree well with previous computations for steady flows and experimental findings for steady and pulsatile flows.     

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.     

Application of large-eddy simulation to the study of pulsatile flow in a modeled arterial stenosis

The technique of large-eddy simulation (LES) has been applied to the study of pulsatile flow through a modeled arterial stenosis. A simple stenosis model has been used that consists of a one-sided 50 percent semicircular constriction in a planar channel. The inlet volume flux is varied sinusoidally in time in a manner similar to the laminar flow simulations of Tutty (1992). LES is used to compute flow at a peak Reynolds number of 2000 and a Strouhal number of 0.024. At this Reynolds number, the flow downstream of the stenosis transitions to turbulence and exhibits all the classic features of post-stenotic flow as described by Khalifa and Giddens (1981) and Lieber and Giddens (1990). These include the periodic shedding of shear layer vortices and transition to turbulence downstream of the stenosis. Computed frequency spectra indicate that the vortex shedding occurs at a distinct high frequency, and the potential implication of this for noninvasive diagnosis of arterial stenoses is discussed. A variety of statistics have been also extracted and a number of other physical features of the flow are described in order to demonstrate the usefulness of LES for the study of post-stenotic flows.     

Numerical simulations of pulsatile blood flow using a new constitutive model

In the present paper we use a new constitutive equation for whole human blood [R.G. Owens, A new microstructure-based constitutive model for human blood, J. Non-Newtonian Fluid Mech. (2006), to appear] to investigate the steady, oscillatory and pulsatile flow of blood in a straight, rigid walled tube at modest Womersley numbers. Comparisons are made with the experimental results of Thurston [Elastic effects in pulsatile blood flow, Microvasc. Res. 9 (1975), 145-157] for the pressure drop per unit length against volume flow rate and oscillatory flow rate amplitude. Agreement in all cases is very good. In the presentation of the numerical and experimental results we discuss the microstructural changes in the blood that account for its rheological behaviour in this simple class of flows. In this context, the concept of an apparent complex viscosity proves to be useful.     

Viscous energy dissipation for steady flow in models of arterial bifurcations

In this paper the viscous energy dissipation in a series of glass model symmetric bifurcations—typical of human vascular branching—was studied. The bifurcations studied have included angles of 75°, 100° and 125° and total output/input area ratios of 0.73, 1.07 and 1.33. The flowrate range studied corresponded to parent tube Reynolds numbers in the range 100–1000.Pressure and flow measurements were made using a highly sensitive variable reluctance pressure transducer and electromagnetic flowmeter. The measurements were made in such a way as to indicate the net effect of the bifurcation.It was found that a dimensionless form of the viscous dissipation had a constant value for all the geometries investigated up to a Reynolds number of 800. Above this, for an angle of 125°, the measured energy dissipation increased. An analytical model based on entry flow principles showed good agreement with measured values except at an area ratio of 0.73. The reasons for this are discussed with the physiological implications of the results.     

Platelet near-wall excess in porcine whole blood in artery-sized tubes under steady and pulsatile flow conditions

Platelet margination (enhanced platelet concentration in the near wall region of a blood vessel) has been well documented in small vessels. In artery-sized vessels margination has only been demonstrated in one study, using ghost cell suspensions and under relatively non-physiologic conditions of steady flow and 50 cm development length. Local sampling experiments were performed to confirm platelet margination in artery-sized stainless steel tubes, for a typical anatomical length and under pulsatile flow, using fresh EDTA-anticoagulated porcine whole blood (N=21). Experiments were designed using three-dimensional Computational Fluid Dynamics (CFD) to model the sample region with greater fidelity. Steady flow experiments in 50 cm long tubes verify published laser Doppler measurements of platelet margination in 3 mm ID tubes at normal arterial shear rate (500 s(-1). Margination persists under pulsatile flow conditions (63.8 pulses/min), but in steady flow at length of 10 cm, margination is reduced. Platelet margination ratio (the ratio of the platelet concentration near the wall to bulk average platelet count) ranges from 1.21 to 2.48. No significant effects of calculated sampling thickness (20 microm and 50 microm) or pulsatility were detected. Hematocrit margination ratio is 0.68 to 0.90. Two model platelet concentration profiles are fit to the experimental results.