Pressure drop and flow rate measurements in a human aortic bifurcation cast for steady and pulsatile flow

Pressure drop and flow rate measurements in a rigid cast of a human aortic bifurcation under both steady and physiological pulsatile flow conditions are reported. Integral momentum and mechanical energy balances are used to calculate impedance, spatially averaged wall shear stress and viscous dissipation rate from the data. In the daughter branches, steady flow impedance is within 30% of the Poiseuille flow prediction, while pulsatile flow impedance is within a factor of 2 of fully developed, oscillatory, straight tube flow theory (Womersley theory). Estimates of wall shear stress are in accord with measurements obtained from velocity profiles. Mean pressure drop and viscous dissipation rate are elevated in pulsatile flow relative to steady flow at the mean flow rate, and the exponents of their Reynolds number dependence are in accord with available theory.     

An LDA and flow visualization study of pulsatile flow in an aortic bifurcation model

The structure of pulsatile flow in a rigid aortic bifurcation model was studied by means of a flow visualization technique and three-dimensional laser-Doppler anemometry. The model was made of glass, having the same shape as that of the average human aortic bifurcation. It was installed into a mock circulatory loop that generated physiological pulsatile flow. Flow separation was observed during accelerated and decelerated flow periods. Double helical flow existed inside the flow separation in the early accelerated flow period. In the decelerated flow period, disturbed flow appeared behind the separation zone. Flow was strongly disturbed during the back flow period, and then was gradually stabilized in the forward flow period. The flow separation and the disturbances released from the flow separation zone greatly influenced near-wall velocities along the lateral wall. The wave form of the near-wall velocity in the flow separation zone was much different from that observed in the aortic portion and behind the separation zone; for example, the magnitude of the negative peak velocity in the direction of the tube axis was larger than that of the positive one, and mean velocity in a cycle was very low. This abnormal phasic change of the near-wall velocity may be associated with atherogenesis. The three-dimensional velocity measurement is very useful for the detailed analysis of near-wall velocity patterns.     

Three-dimensional visualization of velocity fields downstream of six mechanical aortic valves in a pulsatile flow model

Velocity fields downstream of 27 mm Bj?rk-Shiley Standard, Bj?rk-Shiley Convex-Concave, Bj?rk-Shiley Monostrut, Hall-Kaster (Medtronic-Hall), St. Jude Medical and Starr-Edwards Silastic Ball aortic valves were studied in a pulsatile mock circulation. Stroke volume was 70 cm3 and frequency 71 min-1 and 88 min-1. Fluid velocity was measured by a catheter mounted hot-film anemometer probe in a glycerol water mixture one and two diameters downstream of the aortic valve. Velocity fields were dynamically visualized by a three-dimensional technique and revealed qualitative independence of frequency. All profiles were flat in the acceleration phase of systole. From peak systole and throughout the systolic deceleration phase profiles characteristic of the individual valves appeared. The pivoting and tilting disc valves caused a skewed velocity profile with highest velocities downstream of the major orifice and lowest velocities downstream of the minor orifice. The differences between the three investigated Bj?rk-Shiley valves were remarkable. The St. Jude Medical valve generated velocity peaks downstream of the two major orifices and the central slit, and lower velocities in the hinge areas. A rather flat profile with central hollowing was seen downstream of the Starr-Edwards Ball valve. All velocity profiles were more or less dampened two diameters downstream.     

Steady pressure-flow relationship in a cast of the upper and central human airways

Pressure drops across the upper (larynx) and central airways of a human lung cast were measured at steady state inspiratory and expiratory flows. Air, HeO₂ and SF6-O₁ gas mixtures were used at tracheal Reynolds' numbers ranging from 145 to 30 000. The pressure-flow characteristics of the model were analysed using standard pressure-flow diagrams and Moody plots. We found that the asymmetry between inspiratory and expiratory resistances, observed in the central airways (larynx excluded), was markedly reduced in the presence of the larynx. However, static pressure differences were greater across the entire model of the upper and central airways than across the model of the five generations of the tracheo-bronchial tree (without larynx) at the same flow-rates. In addition, our results showed that the presence of the larynx tended to reduce the zone of fully developed laminar flow in the Moody diagram with the higher density gas, while extending the zone of turbulent flow even for the low density gas at low Reynold's numbers.     

Some flow visualization and laser-Doppler-velocity measurements in a true-to-scale elastic model of a human aortic arch--a new model technique.

Flow studies were done in an elastic true-to-scale silicone rubber model of an aortic arch to study further hemodynamic influences on atherosclerosis. The model was prepared from a cast of a young woman. A revised model technique was used. The model had a compliance similar to that of the human aortic arch. Velocity measurements were done in the model with a two component laser-Doppler-anemometer in steady and pulsatile flow using a calcium chloride solution with a viscosity of eta = 3.18 mPas and density of rho = 1.28 kg/m3 at 20 degrees C. The time average Reynolds numbers over a whole cycle in the ascending aorta was Re = 1350. The Womersley parameter for pulsatile flow was a = 20. The pulse wave velocity in the ascending aorta was about c = 5.4 m/sec. The secondary flow behavior was discussed for steady and pulsatile flow. Reverse flows were found, especially along the inner radius of the aortic arch in the descending aorta in steady and pulsatile flow and also in small areas of the ascending aorta and at the branches of the aortic arch. The formation of atherosclerotic plaques at preferred local flow regions is discussed.     

Quantitative flow visualization: toward a comprehensive flow diagnostic tool

Quantitative flow visualization has many roots and has takenseveral approaches. The advent of digital image processing hasmade it possible to practically extract useful information fromevery kind of flow image. In a direct approach, the image intensityor color (wavelength or frequency) can be used as an indicationof concentration, density and temperature fields or gradientsof these scalar fields in the flow (Merzkirch, 1987). For whole-fieldvelocity measurement, the method of choice by experimental fluidmechanicians has been the technique of Particle Image Velocimetry(DPIV). This paper presents a novel approach to extend the DPIVtechnique from a planar method to a full three-dimensional volumemapping technique useful in both engineering and biologicalapplications.     

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.     

Steady flow through a constricted cylinder by multiparticle collision dynamics

The flow characterization of blood through healthy and diseased flow geometries is of interest to researchers and clinicians alike, as it may allow for early detection, and monitoring, of cardiovascular disease. In this paper, we use a numerically efficient particle-based flow model called multiparticle collision dynamics (MPC for short) to study the effect of compressibility and slip of flow of a Newtonian fluid through a cylinder with a local constriction. We use a cumulative averaging method to compare our MPC results to the finite-element solution of the incompressible no-slip Navier-Stokes equations in the same geometry. We concentrate on low Reynolds number flows [ $Re \in (4,30)$ ] and quantify important differences observed between the MPC results and the Navier-Stokes solution in constricted geometries. In particular, our results show that upstream recirculating zones can form with the inclusion of slip and compressibility, which are not observed in the flow of an incompressible Newtonian fluid using the no-slip assumption, but have been observed experimentally for blood. Important flow features are also presented that could be used as indicators to observe compressibility and slip in experimental data where near-wall data may be difficult to obtain. Finally, we found that the cumulative averaging method used is ideal for steady particle-based flow methods, as macroscopic no-slip is readily obtained using the MPC bounce-back rule. Generally, a small spurious slip is seen using other averaging methods such as weighted spatial averages or averages over several runs, and the bounce-back rule has to be modified so as to achieve macroscopic no-slip. No modifications of the bounce-back rule were required for our simulations.     

Pulsatile blood flow in a rigid pulmonary alveolar sheet with porous walls

This paper deals with the pulsatile blood flow in the lung alveolar sheets by idealizing each of them as a channel covered by porous media. As the blood flow in the lung is of low Reynolds number, a creeping flow is assumed in the channel. The analytical and numerical results for the velocity and pressure distribution in the porous medium are presented. The effect of an imposed slip condition is also studied. Comparisons with the corresponding results for the steady-state case are made at the end.     

Pulsed magnetohydrodynamic blood flow in a rigid vessel under physiological pressure gradient

Blood flow in a steady magnetic field has been of great interest over recent years. Many researchers have examined the effects of magnetic fields on velocity profiles and arterial pressure, and major studies have focused on steady or sinusoidal flows. In this paper, we present a solution for pulsed magnetohydrodynamic blood flow with a somewhat realistic physiological pressure wave obtained using a Windkessel lumped model. A pressure gradient is derived along a rigid vessel placed at the output of a compliant module which receives the ventricle outflow. Then, velocity profile and flow rate expressions are derived in the rigid vessel in the presence of a steady transverse magnetic field. As expected, results showed flow retardation and flattening. The adaptability of our solution approach allowed a comparison with previously addressed flow cases and calculations presented a good coherence with those well established solutions.     

A simulated dye method for flow visualization with a computational model for blood flow

A numerical dye method for the visualization of unsteady three-dimensional flow calculations is introduced by coupling the unsteady convection-diffusion equation to the Navier-Stokes equation for mass and momentum. This system of equations is descretized using a finite volume projection-like algorithm with generalized coordinates and overset grids. A powerful pressure prediction method is used to accelerate the convergence of the Pressure Poisson equation. To demonstrate the visualization technique, blood flow through the aortic arch region and the three main arterial branches is computed using various Womersley numbers. In this technique, parcels of fluid are followed in time as a function of the cardiac cycle without having to track individual particles, which in turn aids us to better understand some important aspects of the three-dimensionality of the developing unsteady flow. Using this numerical dye method we analyze the strength of the cross flow during the cardiac cycle, the relationship between the penetration of blood into the aortic branches from its relative position in the ascending aortic region and the effects of the Womersley parameter. This technique can be very useful in the design and development of stents where the topology of the device would require understanding where the blood emanating from the heart ends up at the end of the cardiac cycle. Moreover, this method could be useful in investigating the influence of flow and geometry on the local introduction of medication.     

Steady flow through a double converging-diverging tube model for mild coronary stenoses

Velocity profiles and the pressure drop across two mild (62 percent) coronary stenoses in series have been investigated numerically and experimentally in a perspex-tube model. The mean flow rate was varied to correspond to a Reynolds number range of 50-400. The pressure drop across two identical (62 percent) stenoses show that for low Reynolds numbers the total effect of two stenoses equals that of two single stenoses. A reduction of 10 percent is found for the higher Reynolds numbers investigated. Numerical and experimental results obtained for the velocity profiles agree very well. The effect of varying the converging angle of a single mild (62 percent) coronary stenosis on the fluid flow has been determined numerically using a finite element method. Pressure-flow relation, especially with respect to relative short stenoses, is discussed.