Detailed visualization of pulsatile flow fields produced by modelled arterial stenoses

Transport physics and biorheology in the setting of hemostasis and thrombosis

The biophysics of blood flow can dictate the function of molecules and cells in the vasculature with consequent effects on hemostasis, thrombosis, embolism, and fibrinolysis. Flow and transport dynamics are distinct for (i) hemostasis vs. thrombosis and (ii) venous vs. arterial episodes. Intraclot transport changes dramatically the moment hemostasis is achieved or the moment a thrombus becomes fully occlusive. With platelet concentrations that are 50- to 200-fold greater than platelet-rich plasma, clots formed under flow have a different composition and structure compared with blood clotted statically in a tube. The platelet-rich, core/shell architecture is a prominent feature of self-limiting hemostatic clots formed under flow. Importantly, a critical threshold concentration of surface tissue factor is required for fibrin generation under flow. Once initiated by wall-derived tissue factor, thrombin generation and its spatial propagation within a clot can be modulated by γ'-fibrinogen incorporated into fibrin, engageability of activated factor (FIXa)/activated FVIIIa tenase within the clot, platelet-derived polyphosphate, transclot permeation, and reduction of porosity via platelet retraction. Fibrin imparts tremendous strength to a thrombus to resist embolism up to wall shear stresses of 2400 dyne cm(-2) . Extreme flows, as found in severe vessel stenosis or in mechanical assist devices, can cause von Willebrand factor self-association into massive fibers along with shear-induced platelet activation. Pathological von Willebrand factor fibers are A Disintegrin And Metalloprotease with ThromboSpondin-1 domain 13 resistant but are a substrate for fibrin generation due to FXIIa capture. Recently, microfluidic technologies have enhanced the ability to interrogate blood in the context of stenotic flows, acquired von Willebrand disease, hemophilia, traumatic bleeding, and drug action.

Effect of Pulsatile Flow Waveform and Womersley Number on the Flow in Stenosed Arterial Geometry

The salient hemodynamic flow features in a stenosed artery depend not only on the degree of stenosis, but also on its location in the circulatory system and the physiological condition of the body. The nature of pulsatile flow waveform and local Womersley number vary in different regions of the arterial system and at different physiological state, which affects the local hemodynamic wall parameters, for example, the wall shear stress (WSS) and oscillatory shear index (OSI). Herein, we have numerically investigated the effects of different waveforms and Womersley numbers on the flow pattern and hemodynamic parameters in an axisymmetric stenosed arterial geometry with 50% diametral occlusion. Temporal evolution of the streamlines and hemodynamic parameters are investigated, and the time-averaged hemodynamic wall parameters are compared. Presence of the stenosis is found to increase the OSI of the flow even at the far-downstream side of the artery. At larger Womersley numbers, the instantaneous flow field in the stenosed region is found to have a stronger influence on the flow profiles of the previous time levels. The study delineates how an approximation in the assumption of inlet pulsatility profile may lead to significantly different prediction of hemodynamic wall parameters.

Direct observation of von Willebrand factor elongation and fiber formation on collagen during acute whole blood exposure to pathological flow

OBJECTIVE

In severe stenosis, von Willebrand factor (vWF) experiences millisecond exposures to pathological wall shear rates (γ(w)). We sought to evaluate the deposition of vWF onto collagen surfaces under flow in these environments.

METHODS AND RESULTS

Distinct from viscometry experiments that last many seconds, we deployed microfluidic devices for single-pass perfusion of whole blood or platelet-free plasma over fibrillar type 1 collagen (<50 ms transit time) at pathological γ(w) or spatial wall shear rate gradients (grad γ(w)). Using fluorescent anti-vWF, long thick vWF fibers (>20 μm) bound to collagen were visualized at constant γ(w)>30000 s(-1) during perfusion of platelet-free plasma, a process enhanced by EDTA. Rapid acceleration or deceleration of EDTA platelet-free plasma at grad γ(w)=±1.1×10(5) to ±4.3×10(7) s(-1)/cm did not promote vWF deposition. At 19400 s(-1), EDTA blood perfusion resulted in rolling vWF-platelet nets, although blood perfusion (normal Ca(2+)) generated large vWF/platelet deposits that repeatedly embolized and were blocked by anti-glycoprotein Ib or the α(IIb)β(3) inhibitor GR144053 and did not require grad γ(w). Blood perfusion at venous shear rate (200 s(-1)) produced a stable platelet deposit that was a substrate for massive but unstable vWF-platelet aggregates when flow was increased to 7800 s(-1).

CONCLUSIONS

Triggered by collagen and enhanced by platelet glycoprotein Ib and α(IIb)β(3), vWF fiber formation occurred during acute exposures to pathological γ(w) and did not require gradients in wall shear rate.

Neutrophil adhesion on endothelial cells in a novel asymmetric stenosis model: effect of wall shear stress gradients

Leukocytes play a pivotal role in the progression of atherosclerosis. A novel three-dimensional in vitro asymmetric stenosis model was used to better investigate the role of local hemodynamics in the adhesion of leukocytes to an established plaque. The adhesion of a human promyelocytic cell line (NB4) on a human abdominal aortic endothelial cell (EC) monolayer was quantified. NB4 cells were circulated over TNF-alpha stimulated and nonstimulated ECs for 1 or 6 h at 1.25 or 6.25 dynes/cm(2) and compared to static conditions. Cytokine stimulation increased significantly EC expression of intercellular adhesion molecule and vascular cell adhesion molecule. Under static conditions, neutrophils adhered overall more than under flow, with decreased adhesion with increasing shear. Adhesion was significantly higher in the recirculation region distal to the stenosis than in the inlet. Preshearing the ECs decreased the expression of cell adhesion molecules in inflamed endothelium and significantly decreased adhesion. However, the ratio of adhesion between the recirculation zone and the inlet increased, hence exhibiting an increased regional difference. This work suggests an important role for neutrophil-EC interactions in the atherosclerotic process, especially in wall shear stress gradient regions. This is important clinically, potentially helping to explain plaque stability.

Haemodynamics and blood flow measured using ultrasound imaging

Visualization of, and measurements related to, haemodynamic phenomena in arteries may be made using ultrasound systems. Most ultrasound technology relies on simple measurements of blood velocity taken from a single site, such as the peak systolic velocity for assessment of the degree of lumen reduction caused by an arterial stenosis. Real-time two-dimensional (2D) flow field visualization is possible using several methods, such as colour flow, blood flow imaging, and echo particle image velocimetry; these have applications in the examination of the flow field in diseased arteries and in heart chambers. Three-dimensional (3D) and four-dimensional ultrasound systems have been described. These have been used to provide 2D velocity profile data for the estimation of volumetric flow. However, they are limited for haemodynamic evaluation in that they provide only one component of the velocity. The provision of all seven components (three space, three velocity, and one time) is possible using image-guided modelling, in which 3D ultrasound is combined with computational fluid dynamics. This method also allows estimation of turbulence data and of relevant quantities such as the wall shear stress.

Simulation and validation of arterial ultrasound imaging and blood flow

We reviewed the simulation and validation of arterial ultrasound imaging and blood flow assessment. The physical process of ultrasound imaging and measurement is complex, especially in disease. Simulation of physiological flow in a phantom with tissue equivalence of soft tissue, vessel wall and blood is now achievable. Outstanding issues are concerned with production of anatomical models, simulation of arterial disease, refinement of blood mimics to account for non-Newtonian behavior and validation of velocity measurements against an independent technique such as particle image velocimetry. String and belt phantoms offer simplicity of design, especially for evaluation of velocity estimators, and have a role as portable test objects. Electronic injection and vibrating test objects produce nonphysiologic Doppler signals, and their role is limited. Computational models of the ultrasound imaging and measurement process offer considerable flexibility in their ability to alter multiple parameters of both the propagation medium and ultrasound instrument. For these models, outstanding issues are concerned with the inclusion of different tissue types, multilayer arteries, inhomogeneous tissues and diseased tissues.

A numerical and experimental study of periodic flow in a model of a corrugated vessel with application to stented arteries

The investigation presented in this paper has attempted to study, from a fluid dynamics point of view, the consequences of the presence of a stent on the flow of blood. The method adopted is mainly by numerical simulation using a finite element technique, but visualization experiments and an analytical study have also been carried out. The flow of blood is treated as being transient, laminar and Newtonian within a rigid section of the stented vessel. The flow is driven by an imposed pressure gradient in the form of a physiological waveform. A periodic boundary condition has been applied. Particle trajectories, fluid shear rate contours and wall shear stress maps, together with the use of a simplified form of particle image velocimetry in the experimental work, have been used to interpret the results. It was found that the vessel wall experienced oscillating levels of wall shear stress, in particular extensive exposure to low shear stress. Regions of flow recirculation, points of flow separation and reattachment constantly move in regions of low shear stress. Vortices form and then are destroyed rapidly by reversing flow. Their potential physiological significance to stented arteries is discussed.

Relation of the flow field distal to a moderate stenosis to the Doppler power

An experimental investigation was undertaken to establish how different flow regimes affect the Doppler signal. A rigid tube model consisting of a 70% asymmetric area stenosis was used with steady and pulsatile flow conditions. The characteristics of the flow field at various sites was determined using a photochromic flow visualization method. Continuous-wave Doppler measurements were made using a 41% suspension of human red blood cells (RBCs) in saline as well as a dilute suspension of 4% fixed RBCs. For steady flow, the photochromic results indicated that for Reynolds numbers (Re) of 545 and 1410, turbulence was generated and the length of the turbulent region was found to increase with increasing Re. Under pulsatile flow conditions, turbulence was triggered around peak systole and began to dissipate in late deceleration, and by the end of diastole the flow field almost relaminarized. During the turbulent phase of the flow cycle, the poststenotic flow field was seen to consist of four distinct flow regimes similar to those observed for steady flow. For higher Womersley parameters and Reynolds numbers the turbulent zone was found to be larger and to occupy a greater fraction of the flow cycle. These flow visualization results were compared with the Doppler power measurements made at the same locations and under similar flow conditions. At physiological hematocrits (41%) the onset of turbulence for both steady and pulsatile flow increased the backscattered Doppler power. The location of the peak Doppler power coincided with the region of maximum turbulence observed using the photochromic technique.

A transputer-based physiological signal processing system. Part 2--System testing and investigation of flow through models of very small arterial stenoses

This paper describes the performance testing of a novel transputer-based physiological digital signal processing (DSP) unit and its application in the interpretation of pulsed Doppler ultrasound signals, obtained from models of arterial stenoses. The first test used the DSP unit as a stand-alone spectrum analyser using (1) sinusoidal frequencies (50 Hz to 10 kHz) and (2) filtered white noise (centre frequency 3 kHz, bandwidth 2.5 kHz). For the second test, the DSP unit was attached to a 30-channel multi-gate Doppler ultrasound scanner (transmitting a 4.8 MHz pulse with a repetition frequency of 4.8 kHz) and a vessel tracking unit. The Doppler ultrasound signals obtained from steady flow (100-600 ml/min) in a rigid acrylic tube (internal diameter 6 mm) were then analyzed by the DSP unit and a commercially available system. Lastly, an in vitro investigation into the flow disturbances around very small stenoses (2-25% cross-sectional area reduction), using steady flow (100-600 ml/min), was undertaken. The results indicated that the system was capable of detecting stenoses as small as 5% cross-sectional area reduction.

A transputer-based physiological signal processing system. Part 1--System design

This paper, the first of two, details the design and in-vitro testing of a transputer-based physiological signal processing system. The heart of the system is a transputer-based digital signal processing (DSP) board which can act as a stand-alone spectrum analyser, designed to operate in the audio-frequency band up to 25 kHz. The board comprises a T800 processor, two A100 transversal filters, 12 bit A-D circuitry capable of sampling up to 48 kHz, memory and address mapper. The initial application of the system is for the detection of early arterial disease. For this the DSP board is harnessed to the front end of a multigate pulsed Doppler ultrasound scanner operating at 4.8 MHz insonation frequency and incorporating a vessel wall tracking unit. The complete system performs a Fourier transform on the backscattered signals, providing spectral information on discrete areas of flow (0.6 mm3) across the vessel lumen in real time. This first paper describes the hardware, and the second describes the performance testing of the system on the bench and an assessment of its ability to detect low grade stenoses during steady flow.

Detection of early atherosclerosis by analysis of ultrasonic Doppler signals produced by mural flow disturbances

The paper describes an in vitro study using a multi-gate Doppler ultrasound system to investigate flow disturbances in a blood analogue caused by small stenoses (2-25% cross-sectional area reduction), using steady flow (100-600 ml min-1) in a 6 mm diameter rigid artery model. The results indicate that stenoses greater than 5% were detectable.

Simulation of three-dimensional pulsatile flow through an asymmetric stenosis

The main objective of this work was to use desktop workstations to evaluate the computer code HEMO as a tool for predicting coronary blood flows. The flows are usually characterised by complex vortical structures and transitional effects, and as such present challenging computational problems. As the results of the computations shown in the paper demonstrate, we can predict realistic pulsatile flows in constricted tubes using the Sun Sparcstation 1+ in a matter of hours. The results shown in the paper have also demonstrated that the computer simulations can be very useful as a complementary tool for experimental investigations.

Spatial and temporal variations of wall shear stress within an end-to-side arterial anastomosis model

Wall shear stress has been strongly implicated in the initiation of fibrous intimal hyperplasia that leads to arterial bypass graft failure. In this study, the photochromic tracer technique was used to determine the instantaneous value of the wall shear stress in order to provide a detailed description of its spatial and temporal variations within a 45 degrees end-to-side anastomosis model. At the heel and on the bed across from the heel, the shear stress was close to zero throughout the sinusoidal flow cycle. Flow separation was produced just beyond the toe where the peak phase-averaged shear stress was 6 dyn cm-2, whilst on the bed across from the toe it was 27 dyn cm-2. The stagnation point was seen to fluctuate quite sharply on the bed over a distance of about one tube radius upstream from the site across from the toe. This led to large spatial gradients of the wall shear stress with a peak value of 580 dyn cm-3. In addition, the sudden motion of the stagnation point around peak flow produced sharp temporal gradients of the wall shear stress with a peak absolute value of 3400 dyn cm-2 s-1. When compared to the sites where intimal hyperplasia tends to occur, a strong correlation is seen with low wall shear stress at the heel and toe, and with the sharp temporal variations of the magnitude and spatial gradient of the wall shear stress on the bed across from the junction.

On the Doppler signal from a steady flow asymmetrical stenosis model: effects of turbulence

A steady flow model with a 70% (by area) asymmetrical stenosis was used to examine how changing flow regimes (laminar to turbulent) affect the Doppler signal. Human red blood cells (RBCs) (Hct = 42%) in saline were employed at a flow rate corresponding to a Reynold's number of approximately 545. A dilute suspension of 4% fixed RBCs was also used for the purpose of backscattered power comparison. Measurements of the Doppler signal enabled the backscattered power, time domain statistics, frequency spectra, frequency domain statistics, various spectral indices, autocorrelation function and decorrelation time to be calculated as a function of distance from the stenosis. It is shown that the characteristics of the Doppler signal measured at each site provide information on the nature of the insonated flow field and these correlate well with those expected. The results demonstrate that the onset of turbulence not only affects the Doppler spectrum but also has a profound effect on the signal power, the decorrelation time and the signal statistics.