Lumped Parameter and Feedback Control Models of the Auto-regulatory Response in the Circle of Willis

The Circle of Willis (CoW) 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. A 1-dimensional model of the CoW has been created to simulate a series of possible clinical scenarios such as occlusions in afferent arteries, absent or string-like circulus vessels, or arterial infarctions. The model captures cerebral haemodynamic auto-regulation by using a proportional-integral-derivative (PID) controller to modify efferent resistances and maintain optimal efferent flowrates for a given circle geometry and afferent blood pressure. Results match limited clinical data and results obtained in prior studies to within 6%. In addition, a set of boundary conditions and geometry is presented for which the auto-regulated system cannot provide the necessary efferent flowrates and perfusion, representing a condition with increased risk of stroke and highlighting the importance of modelling the haemodynamics of the CoW. The system model created is computationally simple so it can be used to identify at-risk cerebral arterial geometries and conditions prior to surgery or other clinical procedures.     

Lumped parameter and feedback control models of the auto-regulatory response in the Circle of Willis

The Circle of Willis (CoW) 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. A 1-dimensional model of the CoW has been created to simulate a series of possible clinical scenarios such as occlusions in afferent arteries, absent or string-like circulus vessels, or arterial infarctions. The model captures cerebral haemodynamic auto-regulation by using a proportional-integral-derivative (PID) controller to modify efferent resistances and maintain optimal efferent flowrates for a given circle geometry and afferent blood pressure. Results match limited clinical data and results obtained in prior studies to within 6%. In addition, a set of boundary conditions and geometry is presented for which the auto-regulated system cannot provide the necessary efferent flowrates and perfusion, representing a condition with increased risk of stroke and highlighting the importance of modelling the haemodynamics of the CoW. The system model created is computationally simple so it can be used to identify at-risk cerebral arterial geometries and conditions prior to surgery or other clinical procedures.     

Metabolic model of autoregulation in the Circle of Willis

The Circle of Willis (CoW) is a ringlike structure of blood vessels found at the base of the brain. Its main function is to distribute oxygen-rich arterial blood to the cerebral mass. In a previous study, a one-dimensional (1D) model of the CoW was created to simulate a series of possible clinical scenarios such as occlusions in afferent arteries, absent or stringlike circulus vessels, or arterial infarctions (Moorhead et al., 2004, Comput. Methods Biomech. Biomed. Eng., 7(3), pp. 121-130). The model captured cerebral haemodynamic autoregulation by using a proportional-integral-derivative (PID) controller to modify efferent artery resistances. Although some good results and correlations were achieved, the model was too simple to capture all the transient dynamics of autoregulation. Hence a more physiologically accurate model has been created that additionally includes the oxygen dynamics that drive the autoregulatory response. Results very closely match accepted physiological response and limited clinical data. In addition, a set of boundary conditions and geometry is presented for which the autoregulated system cannot provide sufficient perfusion, representing a condition with increased risk of stroke and highlighting the importance of modeling the haemodynamics of the CoW. The system model created is computationally simple so it can be used to identify at-risk cerebral arterial geometries and conditions prior to surgery or other clinical procedures.     

The variability of the circle of Willis: univariate and bivariate analysis

The results are presented of a statistical analysis of the variability of the circle of Willis using univariate and bivariate methods. For this purpose 100 circles of Willis were available. From each circle 19 indexes of arterial size were determined, the basilar artery was measured in two places. Half the circumference was measured. This data yielded no evidence of differences between left- and right-sided vessels in the sample. An important source of variation is the general size of all vessels considered. When the data are cleared from this general size variation, correlation coefficients reveal interesting relations between the vessels. The posterior communicating arteries are strongly related to the ipsilateral carotid artery, whereas a strong inverse relationship exists with the basilar artery and the precommunicating part of the ipsilateral posterior cerebral artery. These relationships can be understood from the expected patterns of the blood flow in these vessels. Similar relationships can be found in the anterior part of the circle of Willis and in the vertebro-basilar junction. In a different manner, based on previous haemodynamic studies, the relation between blood flow and vessel size within the circle of Willis can be demonstrated by relating the ratios of the sizes of afferent and efferent arteries to the sizes of the posterior communicating arteries, an "intuitive" model. The supposed correlations of the outcome of this "intuitive" model with the size of the communicating arteries appeared to by highly significant. It is concluded that the variations of the circle of Willis are related to the individual variations of the blood flow in this arterial network.     

Numerical simulations of the blood flow in the patient-specific arterial cerebral circle region

The Cerebral Circle Region, also known as the Circle of Willis (CoW), is a loop of arteries that form arterial connections between supply arteries to distribute blood throughout the cerebral mass. Among the population, only 25% to 50% have a complete system of arteries forming the CoW. 3D time-varying simulations for three different patient-specific artery anatomies of CoW were performed in order to gain a better insight into the phenomena existing in the cerebral blood flow. The models reconstructed on the basis of computer tomography images start from the aorta and include the largest arteries that supply the CoW and the arteries of CoW. Velocity values measured during the ultrasound examination have been compared with the results of simulations. It is shown that the flow in the right anterior artery in some cases may be supplied from the left internal carotid artery via the anterior communicating artery. The investigations conducted show that the computational fluid dynamic tool, which provides high resolution in both time and space domains, can be used to support physicians in diagnosing patients of different ages and various anatomical arterial structures.     

Numerical models of auto-regulation and blood flow in the cerebral circulation

A two-dimensional time-dependent computational fluid dynamics model of the Circle of Willis has been developed. To simulate, not only the peripheral resistance of the cerebrovascular tree but also its auto-regulation function, a new "active" boundary condition has been defined and developed using control theory to provide a model of the feedback mechanism. The model was then used to simulate different common abnormalities of the Circle of Willis while a pressure drop, simulating a rapid compression of the right internal carotid artery, was imposed. Test results using a simple tube compared excellently with experiment. The total time-dependent flux for each efferent artery was tabulated and showed the important relationship between geometrical variations in the Circle of Willis and the auto-regulation of blood flow by vascular vaso-dilation and contraction. From this study, it was found that the worst case seemed to be that of a missing or dysfunctional right A1 segment of the anterior cerebral artery. The use of valid physiological models of the peripheral resistance allows for more realistic models of the blood flow in the Circle whilst allowing an easy extension to 3D patient specific simulations.     

Numerical Models of Auto-regulation and Blood Flow in the Cerebral Circulation

A two-dimensional time-dependent computational fluid dynamics model of the Circle of Willis has been developed. To simulate, not only the peripheral resistance of the cerebrovascular tree but also its auto-regulation function, a new "active" boundary condition has been defined and developed using control theory to provide a model of the feedback mechanism. The model was then used to simulate different common abnormalities of the Circle of Willis while a pressure drop, simulating a rapid compression of the right internal carotid artery, was imposed. Test results using a simple tube compared excellently with experiment. The total time-dependent flux for each efferent artery was tabulated and showed the important relationship between geometrical variations in the Circle of Willis and the auto-regulation of blood flow by vascular vaso-dilation and contraction. From this study, it was found that the worst case seemed to be that of a missing or dysfunctional right A1 segment of the anterior cerebral artery. The use of valid physiological models of the peripheral resistance allows for more realistic models of the blood flow in the Circle whilst allowing an easy extension to 3D patient specific simulations.     

The variability of the circulus arteriosus (Willisii): order or anarchy?

The variation of the circulus arteriosus is studied using multivariate methods. The data which form the basis of this study are 19 measurements of half the circumference of the arteries that form the circle of Willis and its afferent and efferent branches; 100 circles of Willis were measured for this purpose. Since the number of variables per individual is large, multivariate statistical techniques are the most appropriate method to gain insight in the relations of vessel sizes that exist within the circle of Willis. So a principal component analysis was performed on the data. The results clearly show a number of relations between vessel sizes. In general, inverse relationships were found of vessels that have (at least partially) an identical irrigation area: both internal carotid arteries and the ipsilateral posterior communicating artery show an intimate relationship and are together inversely related to the basilar artery and the precommunicating part of the posterior cerebral artery. Inverse relationships are also found for both vertebral arteries and both precommunicating parts of the anterior cerebral arteries. The homonymous efferent arteries appear to be closely related and show an independent variation. Together the first six principal components explain 69% of the variance. These results support a haemodynamical hypothesis on the explanation of the variability of the circle of Willis. Moreover, the differential growth in the head-neck region during the first two decades of life is postulated to be the origin of a part of the variation.     

Computational modelling of emboli travel trajectories in cerebral arteries: influence of microembolic particle size and density

Ischaemic stroke is responsible for up to 80 % of stroke cases. Prevention of the reoccurrence of ischaemic attack or stroke for patients who survived the first symptoms is the major treatment target. Accurate diagnosis of the emboli source for a specific infarction lesion is very important for a better treatment for the patient. However, due to the complex blood flow patterns in the cerebral arterial network, little is known so far of the embolic particle flow trajectory and its behaviour in such a complex flow field. The present study aims to study the trajectories of embolic particles released from carotid arteries and basilar artery in a cerebral arterial network and the influence of particle size, mass and release location to the particle distributions, by computational modelling. The cerebral arterial network model, which includes major arteries in the circle of Willis and several generations of branches from them, was generated from MRI images. Particles with diameters of 200, 500 and 800  $\upmu \hbox {m}$ and densities of 800, 1,030 and 1,300 $\hbox {kg/m}^{3}$ were released in the vessel’s central and near-wall regions. A fully coupled scheme of particle and blood flow in a computational fluid dynamics software ANASYS CFX 13 was used in the simulations. The results show that heavy particles (density large than blood or a diameter larger than 500  $\upmu \hbox {m}$ ) normally have small travel speeds in arteries; larger or lighter embolic particles are more likely to travel to large branches in cerebral arteries. In certain cases, all large particles go to the middle cerebral arteries; large particles with higher travel speeds in large arteries are likely to travel at more complex and tortuous trajectories; emboli raised from the basilar artery will only exit the model from branches of basilar artery and posterior cerebral arteries. A modified Circle of Willis configuration can have significant influence on particle distributions. The local branch patterns of internal carotid artery to middle cerebral artery and anterior communicating artery can have large impact on such distributions.     

1D simulation of blood flow characteristics in the circle of Willis using THINkS

One-dimensional (1D) simulation of the complete vascular network, so called THINkS (Total Human Intravascular Network Simulation) is developed to investigate changes of blood flow characteristics caused by the variation of CoW. THINkS contains 158 major veins, 85 major arteries, and 77 venous and 43 arterial junctions. THINkS is validated with available in vivo blood flow waveform data. The overall trends of flow rates in variations of the CoW, such as the missing anterior cerebral artery (missing-A1) or missing posterior cerebral artery (missing-P1), are confirmed by in vivo experimental data. It is demonstrated that the CoW has the ability to shunt blood flow to different areas in the brain. Flow rates in efferent arteries remain unaffected under the variation of CoW, while the flow rates in afferent vessels can be subject to substantial changes. The redistribution of blood flow can cause particular vessels to undergo extra flow rate and hemodynamic stresses.     

A mathematical model of the flow in the circle of Willis

A mathematical model of the flow in the circle of Willis has been designed and the effects of (a) the large anatomical variation of the communicating arteries and (b) physiological changes of the resistances of the vertebral arteries have been studied. The influence of the posterior perforating arteries on the flow in the posterior communicating arteries has been investigated as well, with special attention being paid to the possible occurrence of a 'dead point'. In the model, the influence of diameters of the communicating arteries on the flow in the afferent vessels and the segments of the circle turns out to be considerable, especially in the range of the anatomical variation of the diameters. Within this range flow reductions due to an increased resistance of the vertebral artery will be compensated for by the system. Assuming that the values and ratios of the peripheral resistances are within the physiological range, a dead point is not to be expected in the flow in the posterior communicating arteries.     

Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows

Blood flow in the circle of Willis (CoW) is modelled using the 1-D equations of pressure and flow wave propagation in compliant vessels. The model starts at the left ventricle and includes the largest arteries that supply the CoW. Based on published physiological data, it is able to capture the main features of pulse wave propagation along the aorta, at the brachiocephalic bifurcation and throughout the cerebral arteries. The collateral ability of the complete CoW and its most frequent anatomical variations is studied in normal conditions and after occlusion of a carotid or vertebral artery (VA). Our results suggest that the system does not require collateral pathways through the communicating arteries to adequately perfuse the brain of normal subjects. The communicating arteries become important in cases of missing or occluded vessels, the anterior communicating artery (ACoA) being a more critical collateral pathway than the posterior communicating arteries (PCoAs) if an internal carotid artery (ICA) is occluded. Occlusions of the VAs proved to be far less critical than occlusions of the ICAs. The worst scenario in terms of reduction in the mean cerebral outflows is a CoW without the first segment of an anterior cerebral artery combined with an occlusion of the contralateral ICA. Furthermore, in patients without any severe occlusion of a carotid or VA, the direction of flow measured at the communicating arteries corresponds to the side of the CoW with an absent or occluded artery. Finally, we study the effect of partial occlusions of the communicating arteries on the cerebral flows, which again confirms that the ACoA is a more important collateral pathway than the PCoAs if an ICA is occluded.     

Blood vessel size of circulus arteriosus cerebri (circle of Willis): a statistical research on 100 human subjects

The brain weight of 100 fresh cadavers of Italian subjects (60 males and 40 females), aged between 17 and 84 years, was obtained and the corrected circumference of the following blood vessels was measured: basilar artery, internal carotid arteries, anterior and posterior cerebral arteries, and anterior and posterior communicating arteries. The cerebral 'potential flow' was expressed in each case by adding the circumference of the basilar artery to that of the internal carotid arteries. Moreover, the sides and the perimeter of the circle of Willis as well as the length of the basilar artery were calculated. The statistical analysis of the obtained data yielded the following main results: (1) the brain weight decreases with aging, is lower in females than in males and and is statistically correlated neither with the circumferences of the considered arteries and the 'potential flow' nor with the perimeter of the arterial polygon (circle of Willis); (2) the arteries of the left side appear to be larger than those of the right one; (3) no significant difference exists in the circumference and length of the arteries between males and females; (4) the increase of the perimeter of the arterial polygon is achieved by means of a harmonious increase of all its sides; (5) the anterior and posterior communicating arteries have an anarchic pattern, because of the relatively frequent anomalies and the lack of a correlation between their circumference and that of the vessel of origin or of outlet.     

Study of the collateral capacity of the circle of Willis of patients with severe carotid artery stenosis by 3D computational modeling

This numerical study aims to investigate the capacity of the circle of Wills (CoW) to provide collateral blood supply for patients with unilateral carotid arterial stenosis. The basic 3D geometry of the CoW was reconstructed based on a magnetic resonance angiogram of a normal human subject. A total of 52 computational fluid dynamics simulations were performed for four geometry configurations of the CoW with an artificially inserted axisymmetric stenosis of different luminal area reductions in an internal carotid artery (ICA) under a variety of boundary conditions. The CoW geometric configurations included (a) a normal CoW with all communicating arteries; (b) as model (a) but with enlarged communicating arterial diameters; (c) as (a) but with the ipsilateral posterior communicating artery missing, and (d) as (c) but with enlarged communicating arteries. It is found that the blood perfusion pressure drop between the ipsilateral ICA and the middle cerebral artery (MCA) only becomes significant when the degree of stenosis is greater than 86%. The cerebral autoregulation range varied significantly between the different CoW configurations for the severe stenosis cases. Without causing the flow rates to decrease at the efferent arterial ends, the mean perfusion pressure in the ipsilateral ICA can drop from 100 to 73, 67, 92 and 84mmHg for the CoW models (a)-(d) with 96% luminal area reduction stenosis, respectively. The additional pathways are able to raise the ipsilateral MCA pressure significantly without reducing the total flow perfusion. Cerebral autoregulation effects were not directly included in the study. Therefore, the findings in the study should be interpreted with cautions when comes to the biological and clinical significance.     

Vascular anatomy of the gills of the stingarees Urolophus mucosus and U. paucimaculatus (Urolophidae,Elasmobranchii)

The branchial vascular anatomy of Urolophus mucosus and U. paucimaculatus was studied by scanning electron microscopical examination of critical-point-dried tissue or of vascular corrosion casts. The vasculature could be divided into arterioarterial and arteriovenous pathways, which channel the flow of blood through the gills. The arterioarterial pathway consists of an afferent branchial artery which gives rise to afferent distributing arteries that run through the tissues of the interbranchial septum and supply the afferent filament arteries of several filaments. Afferent filament arteries open regularly into a corpus cavernosum in the core of the filament; unlike other elasmobranchs no septal corpora cavernosa are found. At the tip of the filament, channels of the corpus cavernosum connect to a channel which passes across the distal end of the filament from afferent to efferent side. This channel always connects to the afferent filament artery, and in many filaments it connects to the efferent filament artery as well. In addition, a vascular arcade connects all the afferent filament arteries along the entire length of each hemibranch. The filament corpus cavernosum supplies the secondary lamellae. The lamellae drain into efferent lamellar arterioles which in turn drain into the efferent filament artery and the efferent branchial artery. The vascular anatomy of the arteriovenous pathway is similar to that described in other elasmobranchs and consists of arteriovenous anastomoses, found only arising from efferent arterial circulation, and the venolymphatic system, which is composed of the central venous sinus and the companion vessels.     

Hemodynamic role of the circle of Willis in stenoses of internal carotid arteries. An analytical solution of a linear model

A mathematical model of blood flow through the circle of Willis was developed, within a linear framework. Comprehensive analytical solutions, including a remarkably small number of parameters, were derived in the cases of obstructive lesions of extracranial carotid arteries. The influence of these lesions and the role of anterior and posterior communicating arteries on the blood pressure at the entry of the cerebral territories were quantified and analyzed emphasizing that the responses of the system of Willis to obstructive carotid lesions are extremely varied, depending on the communicating artery anatomy. Comparison with numerical results obtained by using a non-linear model showed no physiologically significant differences. Such a model might be an essential tool for an accurate assessment of the cerebral hemodynamics in carotid diseases.     

Cross-flow at the anterior communicating artery and its implication in cerebral aneurysm formation

The anterior communicating artery (ACoA) is an important element of the circle of Willis. While the artery itself is short and small, a large number of intracranial aneurysms can be found at the ACoA. Four subject-specific ACoA models are constructed from 3D rotational angiographic images. The ACoA of these models ranged from 1.7 to 2.7 mm in diameter and 1.5 to 5.7 mm in length. Pulsatile flows through these four ACoA models are studied numerically. Blood is found to move in two opposite directions simultaneously within the ACoA, giving a much higher wall shear at the ACoA. These two opposite flow streams produce a cross-flow that is dependent on the flow rates at the anterior cerebral arteries and internal carotid arteries (ICAs). A larger and shorter ACoA allows flow through the ACoA easily, leading to a greater cross-flow and higher hemodynamic forces on the artery. This cross-flow may disappear when there is a sufficient net flow for a smaller and longer ACoA. Wall shear stress can be as high as 185 Pa at smaller ACoAs, but it can be lowered by asymmetric waveforms at the ICAs. A functional circle of Willis also promotes cross-flow at both the ACoA and posterior communicating 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.     

The development of gill arches and gill blood vessels of the rainbow trout,Salmo gairdneri

Gill development begins on the sixth day of incubation at 10^°C and is complete by 31 days (hatching). Gill arches are formed by fusion and perforation of ectoderm and endoderm across the pharyngeal wall. A primary branchial artery forms within each arch and a second branchial artery forms as a branch from its ventral end. A series of filament loop vessels forms connecting the two arteries and when several are patent a unidirectional blood flow is established via afferent (second) branchial artery, filament loop vessels to efferent (primary) branchial artery. Part of the efferent branchial artery just above its junction with the afferent branchial artery constricts and occludes. It is suggested that this change in the pattern of blood flow is dependent on differences in resistance of the two branchial arteries. A later extension of the gill ventrally is thought not to be homologous with similar regions in elasmobranchs and Acipenser.     

Experimental study of hemodynamics in the circle of willis

Background

The Circle of Willis (CoW) is an important collateral pathway of the cerebral blood flow. An experimental study of the cerebral blood flow (CBF) distribution in different anatomical variations may help to a better understanding of the collateral mechanism of the CoW.

Methods

An in-vitro test rig was developed to simulate the physiological cerebral blood flow in the CoW. Ten anatomical variations were considered in this study, include a set of different degrees of stenosis in L-ICA and L-ICA occlusion coexist with common anatomical variations. Volume flow rates of efferent arteries and pressure signals at the end of communicating arteries of each case were recorded. Physiological pressure waveforms were applied as inlet boundary condition.

Results

In the development of L-ICA stenosis, the total CBF decreases with the increase of stenosis degree. The blood supply of ipsilateral middle cerebral artery (MCA) was affected most by the stenosis of L-ICA. Anterior communicating artery (ACoA) and ipsilateral posterior communicating artery (PCoA) function as important collateral pathways of cerebral collateral circulation when unilateral stenosis occurred. The blood supply of anterior cerebral circulation was compensated by the posterior cerebral circulation through ipsilateral PCoA when L-ICA stenosis degree is greater than 40% and the affected side was compensated immediately by the unaffected side through ACoA. Blood flow of the anterior circulation and the total CBF reached the minimum among all cases studied when L-ICA occlusion coexist with the absence of PCoA.

Conclusion

The results demonstrated the flow distribution patterns of the CoW under anatomical variations and clarified the collateral mechanism of the CoW. The flow ACoA is the most sensitive indexes to the morphology change of ipsilateral ICA. The relative independence of the circulation in anterior and posterior sections of the CoW is not broken and the function of ipsilateral PCoA is not activated until a severe stenosis of unilateral ICA occurs. PCoA is the most important collateral pathway of the collateral circulation and the missing of PCoA has the highest risk of stroke when the ipsilateral ICA has severe stenosis. These findings may provide the basis for future therapeutic and diagnosis applications.