Wave propagation in a model of the arterial circulation

The propagation of the arterial pulse wave in the large systemic arteries has been calculated using a linearised method of characteristics analysis to follow the waves generated by the heart. The model includes anatomical and physiological data for the 55 largest arteries adjusted so that the bifurcating tree of arteries is well matched for forward travelling waves. The peripheral arteries in the model are terminated by resistance elements which are adjusted to produce a physiologically reasonable distribution of mean blood flow. In the model, the pressure and velocity wave generated by the contraction of the left ventricle propagates to the periphery where it is reflected. These reflected waves are re-reflected by each of the bifurcations that they encounter and a very complex pattern of waves is generated. The results of the calculations exhibit many of the features of the systemic arteries, including the increase of the pulse pressure with distance away from the heart as well as the initial decrease and then the large increase in the magnitude of back flow during late systole going from the ascending aorta to the abdominal aorta to the arteries of the leg. The model is then used to study the effects of the reflection or absorption of waves by the heart and the mechanisms leading to the incisura are investigated. Calculations are carried out with the total occlusion of different arterial segments in order to model experiments in which the effects of the occlusion of different arteries on pressure and flow in the ascending aorta were measured. Finally, the effects of changes in peripheral resistance on pressure and velocity waveforms are also studied. We conclude from these calculations that the complex pattern of wave propagation in the large arteries may be the most important determinant of arterial haemodynamics.     

Numerical simulation of pulsatile flow in a compliant curved tube model of a coronary artery

The endothelial cells (ECs) lining a blood vessel wall are exposed to both the wall shear stress (WSS) of blood flow and the circumferential strain (CS) of pulsing artery wall motion. These two forces and their interaction are believed to play a role in determining remodeling of the vessel wall and development of arterial disease (atherosclerosis). This study focused on the WSS and CS dynamic behavior in a compliant model of a coronary artery taking into account the curvature of the bending artery and physiological radial wall motion. A three-dimensional finite element model with transient flow and moving boundaries was set up to simulate pulsatile flow with physiological pressure and flow wave forms characteristic of the coronary arteries. The characteristic coronary artery curvature and flow conditions applied to the simulation were: aspect ratio (lambda) = 10, diameter variation (DV) = 6 percent, mean Reynolds number (Re) = 150, and unsteadiness parameter (alpha) = 3. The results show that mean WSS is about 50 percent lower on the inside wall than the outside wall while WSS oscillation is stronger on the inside wall. The stress phase angle (SPA) between CS and WSS, which characterizes the dynamics of the mechanical force pattern applied to the endothelial cell layer, shows that CS and WSS are more out of phase in the coronaries than in any other region of the circulation (-220 deg on the outside wall, -250 deg on the inside wall). This suggests that in addition to WSS, SPA may play a role in localization of coronary atherosclerosis.     

Influence of the distensibility of large arteries on the longitudinal impedance: application for the development of non-invasive techniques to the diagnosis of arterial diseases

Background

This study shows that the arterial longitudinal impedance constitutes a hemodynamic parameter of interest for performance characterization of large arteries in normal condition as well as in pathological situations. For this purpose, we solved the Navier?CStokes equations for an incompressible flow using the finite element analysis method and the Arbitrary Lagrangian Eulerian (ALE) formulation. The mathematical model assumes a two-dimensional flow and takes into account the nonlinear terms in the equations of fluid motion that express the convective acceleration, as well as the nonlinear deformation of the arterial wall. Several numerical simulations of the blood flow in large vessels have been performed to study the propagation along an arterial vessel of a pressure gradient pulse and a rate flow pulse. These simulations include various deformations of the wall artery leading to parietal displacements ranging from 0 (rigid wall) to 15% (very elastic wall) in order to consider physiological and pathological cases.

Results

The results show significant changes of the rate flow and the pressure gradient wave as a function of aosc, the relative variation in the radius of the artery over a cardiac cycle. These changes are notable beyond a critical value of aosc equal to 0.05. This critical value is also found in the evolution of the longitudinal impedance. So, above a variation of radius of 5%, the convective acceleration, created by the fluid-wall interactions, have an influence on the flow detectable on the longitudinal impedance.

Conclusions

The interpretation of the evolution of the longitudinal impedance shows that it could be a mean to test the performance of large arteries and can contribute to the diagnosis of parietal lesions of large arteries. For a blood vessel with a wall displacement higher than 5% similar to those of large arteries like the aorta, the longitudinal impedance is substantially greater than that obtained in the absence of wall displacement. This study also explains the effects of convective acceleration, on the shape of the decline of the pressure gradient wave and shows that they should not be neglected when the variation in radius is greater than 5%.     

The role of the surrounding tissue in the propagation of waves through the arterial system.

A theoretical analysis of the flow in arteries is presented, taking into consideration the role played by the surrounding tissues in determining the speed of propagatoion and the damping of the blood pressure pulse. This study was undertaken (a) to exhibit a method of computing the flow properties with a more nearly accurate model, (b) to see if the displacement on the skin would be related to the arterial wall displacement, and hence to pressure, velocity and flow rate of blood in the artery, and if it is likely to be measurable. It was found that the pressure of the 'viscous' part in the surrounding tissue increases the pulse velocity and the damping of the wave over the values found by other models which considered only thick-walled elastic tubes with no surrounding tissue. This study also shows that measurements on the skin can provide information about changes in arterial circulation due to diseases such as: edema, arteriosclerosis and others where the Young's modulus for either the arterial wall or the surrounding tissues is altered.     

Validation of a fluid–structure interaction numerical model for predicting flow transients in arteries

Fluid–structure interaction (FSI) numerical models are now widely used in predicting blood flow transients. This is because of the importance of the interaction between the flowing blood and the deforming arterial wall to blood flow behaviour. Unfortunately, most of these FSI models lack rigorous validation and, thus, cannot guarantee the accuracy of their predictions. This paper presents the comprehensive validation of a two-way coupled FSI numerical model, developed to predict flow transients in compliant conduits such as arteries. The model is validated using analytical solutions and experiments conducted on polyurethane mock artery. Flow parameters such as pressure and axial stress (and precursor) wave speeds, wall deformations and oscillating frequency, fluid velocity and Poisson coupling effects, were used as the basis of this validation. Results show very good comparison between numerical predictions, analytical solutions and experimental data. The agreement between the three approaches is generally over 95%. The model also shows accurate prediction of Poisson coupling effects in unsteady flows through flexible pipes, which up to this stage have only being predicted analytically. Therefore, this numerical model can accurately predict flow transients in compliant vessels such as arteries.     

Numerical simulation of local blood flow in the carotid and cerebral arteries under altered gravity

A computational fluid dynamics (CFD) approach was presented to model the blood flows in the carotid bifurcation and the brain arteries under altered gravity. Physical models required for CFD simulation were introduced including a model for arterial wall motion due to fluid-wall interactions, a shear thinning fluid model of blood, a vascular bed model for outflow boundary conditions, and a model for autoregulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models were solved iteratively using the pseudocompressibility method and dual time stepping. Gravity source terms were added to the Navier-Stokes equations to take the effect of gravity into account. For the treatment of complex geometry, a chimera overset grid technique was adopted to obtain connectivity between arterial branches. For code validation, computed results were compared with experimental data for both steady-state and time-dependent flows. This computational approach was then applied to blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other using an anatomical data set. A three-dimensional Circle of Willis configuration was reconstructed from subject-specific magnetic resonance images using an image segmentation method. Through the numerical simulation of blood flow in two model problems, namely, the carotid bifurcation and the brain arteries, it was observed that the altered gravity has considerable effects on arterial contraction/dilatation and consequent changes in flow conditions.     

Surgical Technique for the Implantation of Tissue Engineered Vascular Grafts and Subsequent In Vivo Monitoring

The development of Tissue Engineered Vessels (TEVs) is advanced by the ability to routinely and effectively implant TEVs (4-5 mm in diameter) into a large animal model. A step by-step protocol for inter-positional placement of the TEV and real-time digital assessment of the TEV and native carotid arteries is described here. In vivo monitoring is made possible by the implantation of flow probes, catheters and ultrasonic crystals (capable of recording dynamic diameter changes of implanted TEVs and native carotid arteries) at the time of surgery. Once implanted, researchers can calculate arterial blood flow patterns, invasive blood pressure and artery diameter yielding parameters such as pulse wave velocity, augmentation index, pulse pressures and compliance. Data acquisition is accomplished using a single computer program for analysis throughout the duration of the experiment. Such invaluable data provides insight into TEV matrix remodeling, its resemblance to native/sham controls and overall TEV performance in vivo.     

Can the modified Allen's test always detect sufficient collateral flow in the hand? A computational study

Blood flow in the largest arteries of the arm up to the digital arteries is numerically modelled using the one-dimensional equations of pressure and flow wave propagation in compliant vessels. The model can be applied to different anatomies of arterial networks and can simulate compression of arteries, these allowing us to simulate the modified Allen's test (MAT) and to assess its suitability for the detection of sufficient collateral flow in the hand if radial blood supply is interrupted. The test measures blood flow in the superficial palmar arch before and during compression of the radial artery. The absence of reversal flow in the palmar arch with the compression indicates insufficient collateral flow and is referred to as a positive MAT. This study shows that small calibres of the superficial palmar arch and insufficient compression of the radial artery can lead to false-positive results. Measurement of the drop in digital systolic pressures with compression of the radial artery has proved to be a more sensitive test to predict the presence of sufficient ulnar collateral flow in networks with small calibres of the superficial palmar arch. However, this study also shows that digital pressure measurements can fail in detecting enough collateral flow if the radial artery is insufficiently compressed.     

Can the modified Allen's test always detect sufficient collateral flow in the hand? A computational study

Blood flow in the largest arteries of the arm up to the digital arteries is numerically modelled using the one-dimensional equations of pressure and flow wave propagation in compliant vessels. The model can be applied to different anatomies of arterial networks and can simulate compression of arteries, these allowing us to simulate the modified Allen's test (MAT) and to assess its suitability for the detection of sufficient collateral flow in the hand if radial blood supply is interrupted. The test measures blood flow in the superficial palmar arch before and during compression of the radial artery. The absence of reversal flow in the palmar arch with the compression indicates insufficient collateral flow and is referred to as a positive MAT. This study shows that small calibres of the superficial palmar arch and insufficient compression of the radial artery can lead to false-positive results. Measurement of the drop in digital systolic pressures with compression of the radial artery has proved to be a more sensitive test to predict the presence of sufficient ulnar collateral flow in networks with small calibres of the superficial palmar arch. However, this study also shows that digital pressure measurements can fail in detecting enough collateral flow if the radial artery is insufficiently compressed.     

Comparison of left anterior descending coronary artery hemodynamics before and after angioplasty

Coronary artery disease (CAD) is characterized by the progression of atherosclerosis, a complex pathological process involving the initiation, deposition, development, and breakdown of the plaque. The blood flow mechanics in arteries play a critical role in the targeted locations and progression of atherosclerotic plaque. In coronary arteries with motion during the cardiac contraction and relaxation, the hemodynamic flow field is substantially different from the other arterial sites with predilection of atherosclerosis. In this study, our efforts focused on the effects of arterial motion and local geometry on the hemodynamics of a left anterior descending (LAD) coronary artery before and after clinical intervention to treat the disease. Three-dimensional (3D) arterial segments were reconstructed at 10 phases of the cardiac cycle for both pre- and postintervention based on the fusion of intravascular ultrasound (IVUS) and biplane angiographic images. An arbitrary Lagrangian-Eulerian formulation was used for the computational fluid dynamic analysis. The measured arterial translation was observed to be larger during systole after intervention and more out-of-plane motion was observed before intervention, indicating substantial alterations in the cardiac contraction after angioplasty. The time averaged axial wall shear stress ranged from -0.2 to 9.5 Pa before intervention compared to -0.02 to 3.53 Pa after intervention. Substantial oscillatory shear stress was present in the preintervention flow dynamics compared to that in the postintervention case.     

A Mathematical Study on Three Layered Oscillatory Blood Flow Through Stenosed Arteries

A mathematical model is constructed to examine the characteristics of three layered blood flow through the oscillatory cylindrical tube (stenosed arteries).The proposed model basically consists three layers of blood (viscous fluids with different viscosities) named as core layer (red blood cells),intermediate layer (platelets/white blood cells) and peripheral layer (plasma).The analysis was restricted to propagation of small-amplitude harmonic waves,generated due to blood flow whose wave length is larger compared to the radius of the arterial segment.The impacts of viscosity of fluid in peripheral layer and intermediate layer on the interfaces,average flow rate,mechanical efficiency,trapping and reflux are discussed with the help of numerical and computational results.This model is the generalized form of the preceding models.On the basis of present discussion,it is found that the size of intermediate and peripheral layers reduces in expanded region and enhances in contracted region with the increasing viscosity of fluid in peripheral layer,whereas,opposite effect is observed for viscosity of fluid in intermediate layer.Final conclusion is that the average flow rate and mechanical efficiency increase with the increasing viscosity of fluid in both layers,however,the effects of the viscosity of fluid in both layers on trapping and reflux are opposite to each other.     

Flow-tissue interaction with compliance mismatch in a model stented artery

The insertion of an endovascular prosthesis is known to have a thrombogenic effect that is also a consequence of the interaction between the flowing blood and the stented arterial segment; in fact the prosthesis induces a compliance mismatch and a possible small expansion along the vessel that eventually gives rise to an anomalous distribution of wall shear stresses. The fluid dynamics inside a rectilinear elastic vessel with compliance and section variation is studied here numerically. A recently introduced perturbative approach is employed to model the interaction between the fluid and the elastic tissue; this approximate technique is first validated by comparison with a complete solution within a simple one-dimensional model of the same system. Then it is applied to an axisymmetric model in order to evaluate the flow dynamics and the distribution of wall shear stress in the stented vessel. Compliance mismatch is shown to produce more intense negative wall shear stresses in the stented segment while rapid variations of wall shear stress are found at the stent ends. These effects are enhanced when the prosthesis is accompanied by a small increase of the vessel lumen.     

A data-driven model to study utero-ovarian blood flow physiology during pregnancy

In this paper, we describe a mathematical model of the cardiovascular system in human pregnancy. An automated, closed-loop 1D–0D modelling framework was developed, and we demonstrate its efficacy in (1) reproducing measured multi-variate cardiovascular variables (pulse pressure, total peripheral resistance and cardiac output) and (2) providing automated estimates of variables that have not been measured (uterine arterial and venous blood flow, pulse wave velocity, pulsatility index). This is the first model capable of estimating volumetric blood flow to the uterus via the utero-ovarian communicating arteries. It is also the first model capable of capturing wave propagation phenomena in the utero-ovarian circulation, which are important for the accurate estimation of arterial stiffness in contemporary obstetric practice. The model will provide a basis for future studies aiming to elucidate the physiological mechanisms underlying the dynamic properties (changing shapes) of vascular flow waveforms that are observed with advancing gestation. This in turn will facilitate the development of methods for the earlier detection of pathologies that have an influence on vascular structure and behaviour.

     

Analysis of blood flow through a model of the human arterial system under periodic body acceleration

The human system may be subjected to a body acceleration deliberated for example by making subjects lie down on vibrating tables or more frequently unintentionally, for example during travel in water and land or in air and space. The present study is concerned with the effects of externally imposed body accelerations on blood flow in a branched system of arteries. A finite-element model of flow in the arterial system subject to periodic body accelerations is presented. Computational results on the flow rates through selected arteries and the corresponding inlet and outlet pressures under different conditions (magnitude, frequency and direction) of applied acceleration are presented.     

One-dimensional and three-dimensional models of cerebrovascular flow

The Circle of Willis is a ring-like structure of blood vessels found beneath the hypothalamus at the base of the brain. Its main function is to distribute oxygen-rich arterial blood to the cerebral mass. One-dimensional (1D) and three-dimensional (3D) computational fluid dynamics (CFD) models of the Circle of Willis have been created to provide a simulation tool which can potentially be used to identify at-risk cerebral arterial geometries and conditions and replicate clinical scenarios, such as occlusions in afferent arteries and absent circulus vessels. Both models capture cerebral haemodynamic autoregulation using a proportional-integral (PI) controller to modify efferent artery resistances to maintain optimal efferent flow rates for a given circle geometry and afferent blood pressure. The models can be used to identify at-risk cerebral arterial geometries and conditions prior to surgery or other clinical procedures. The 1D model is particularly relevant in this instance, with its fast solution time suitable for real-time clinical decisions. Results show the excellent correlation between models for the transient efferent flux profile. The assumption of strictly Poiseuille flow in the 1D model allows more flow through the geometrically extreme communicating arteries than the 3D model. This discrepancy was overcome by increasing the resistance to flow in the anterior communicating artery in the 1D model to better match the resistance seen in the 3D results.     

Small disturbance in arterial blood flow

The behaviour of a small disturbance in an arterial blood flow has been studied analytically. The growth equation governing growth or decay of a disturbance has been obtained and solved. The behaviour of wave amplitude has been investigated as the wave propagates in time. The application of results to the human arterial system shows that the shock waves are not expected under normal physiological conditions. In the case of a pathologically increased pressure rise at the root of aorta, shock-like transitions may develop in the periphery. It is observed that the friction effects are to resist the tendency of shock formation in arteries.     

Dynamic response of wall shear stress on the stenosed artery

The present study deals with an appropriate mathematical model of an artery in the presence of constriction in which the generated wall shear stress due to blood flow is analysed. The geometry of the stenosed arterial segment in the diseased state, causing malfunction of the cardiovascular system, is formed mathematically. The flowing blood contained in the stenosed artery is treated as non-Newtonian and the flow is considered to be two-dimensional. The motion of the arterial wall and its effect on local fluid mechanics is not ruled out from the present pursuit. The flow analysis applies the time-dependent, two-dimensional incompressible nonlinear Navier–Stokes equations for non-Newtonian fluids. The flow-field can be obtained primarily following the radial coordinate transformation, using the appropriate boundary conditions and finally adopting a suitable finite difference scheme numerically. The influences of flow unsteadiness, the arterial wall distensibility and the presence of stenosis on the flow-field and the wall shear stresses are quantified in order to indicate the susceptibility to atherosclerotic lesions and thereby to validate the applicability of the present theoretical model.     

Mechanics of blood flow

The historical development of the mechanics of blood flow can be traced from ancient times, to Leonardo da Vinci and Leonhard Euler and up to the present times with increasing biological knowledge and mathematical analysis. In the last two decades, quantitative and numerical methods have steadily given more complete and precise understanding. In the arterial system wave propagation computations based on nonlinear one-dimensional modeling have given the best representation of pulse wave propagation. In the veins, the theory of unsteady flow in collapsible tubes has recently been extensively developed. In the last decade, progress has been made in describing the blood flow at junctions, through stenoses, in bends and in capillary blood vessels. The rheological behavior of individual red blood cells has been explored. A working model consists of an elastic membrane filled with viscous fluid. This model forms a basis for understanding the viscous and viscoelastic behavior of blood.     

Ultrasound techniques for measurement of blood flow and tissue motion

This article will review the ability of ultrasound techniques to provide 3D information on arterial geometry, blood flow and tissue motion.3D systems. 3D datasets can be obtained by sequential acquisition of 2D slices. Ideally a transducer is required in which there is full control of beam steering within a 3D volume. This requires a 2D array consisting of several thousand elements. Prototype 2D arrays have been built which provide several 3D datasets per second.Blood velocity measurement. Current Doppler systems estimate only the component of velocity in the direction of the Doppler beam. Lack of knowledge of the direction of blood motion and also other effects associated with 'spectral broadening' limit the accuracy of velocity measurement. Improved accuracy can be obtained using vector Doppler systems using 2 or 3 beam directions; this approach is referred to as 'vector Doppler'.Tissue motion. Doppler techniques can also be used to detect tissue motion (Tissue Doppler Imaging or TDI). Motion of the artery wall can be calculated from the TDI images. It is possible to estimate simultaneously motion for adjacent diameters within the longitudinal plane, and to visualise the relative motion at different parts of the wall.     

Effect of initial stresses on the wave propagation in arteries

A theoretical analysis for the problem of wave propagation in arteries is presented. Blood is treated as a Newtonian, viscous incompressible fluid. On the basis of information derived from experimental investigations on the mechanical properties of wall tissues, the arterial wall is considered to be nonlinearly viscoelastic and orthotropic. The analysis is carried out for a cylindrical artery, under the purview of membrane theory, by taking account the effect of initial stresses. The motion of the wall and that of the fluid are assumed to be axisymmetric. Particular emphasis has been paid to the propagation of small amplitude harmonic waves whose wavelength is large compared to the radius of the vessel. By employing the equations of motion of the fluid and those for the wall, together with the equations of continuity, a frequency equation is derived by exploiting the conditions of continuity of the velocity of the arterial wall and that of blood on the endosteal surface of the wall. In order to illustrate the validity of the derived analytical expressions a quantitative analysis is made for the variations of the phase velocities as well as the transmission coefficient with frequency corresponding to different transmural pressures which relate to different initial stresses. Computational results indicate that phase velocities increase with the increase of transmural pressures.