Skip to main content

A computational model study of the influence of the anatomy of the circle of willis on cerebral hyperperfusion following carotid artery surgery

Abstract

Background

Cerebral hyperperfusion syndrome develops in a small subset of patients following carotid artery surgery (CAS) performed to treat severe carotid artery stenosis. This syndrome has been found to have a close correlation with cerebral hyperperfusion occurring after CAS. The purpose of this study is to investigate whether and how the anatomy of the Circle of Willis (CoW) of the cerebral circulation influences post-CAS cerebral hyperperfusion.

Methods

A computational model of the cerebral circulation coupled with the global cardiovascular system has been developed to investigate hemodynamic events associated with CAS. Nine topological structures of the CoW were investigated in combination with various distribution patterns of stenosis in the feeding arteries of the cerebral circulation.

Results

The occurrence of post-CAS cerebral hyperperfusion was predicted for the CoW structures that have poor collateral pathways between the stenosed cerebral feeding arteries and the remaining normal feeding arteries. The risk and the localization of post-CAS hyperperfusion were determined jointly by the anatomy of the CoW and the distribution pattern of stenosis in the cerebral feeding arteries. The presence of basilar artery stenosis or contralateral ICA stenosis increased the risk of post-CAS hyperperfusion and enlarged the cerebral region affected by hyperperfusion. For a certain CoW structure, the diameters of the cerebral communicating arteries and the severity of carotid artery stenosis both had a significant influence on the computed post-CAS cerebral hyperperfusion rates. Moreover, post-CAS cerebral hyperperfusion was predicted to be accompanied with an excessively high capillary transmural pressure.

Conclusions

This study demonstrated the importance of considering the anatomy of the CoW in assessing the risk of post-CAS cerebral hyperperfusion. Particularly, since the anatomy of the CoW and the distribution pattern of stenosis in the cerebral feeding arteries jointly determine the risk and localization of post-CAS cerebral hyperperfusion, a patient-specific hemodynamic analysis aimed to help physicians identify patients at high risk of cerebral hyperperfusion should account for the combined effect of the anatomy of cerebral arteries and cerebral feeding artery stenoses on cerebral hemodynamics.

Background

Extracranial internal carotid artery (ICA) stenosis accounts for 15-20% of ischemic strokes and is usually treated by carotid artery surgery (CAS) such as carotid endarterectomy or stenting [1, 2]. A potential risky problem with CAS is that cerebral hyperperfusion syndrome (CHS) (characterized by ipsilateral headache, seizure or intracranial hemorrhage (ICH)) develops in a small subset (0.75-3%) of patients following successful CAS [2]. Although rare, CHS can lead to significant morbidity and mortality if not correctly recognized and treated [2, 3].

The most pronounced hemodynamic event associated with CAS is a sudden increase in cerebral blood flow (CBF). Generally, an over 100% increase in CBF after CAS compared to the pre-CAS value is considered as hyperperfusion [4]. Post-CAS hyperperfusion has been observed in 9-14% of patients in clinical studies [2, 4] and suggested to be an important hemodynamic factor underlying CHS, for instance, the risk of developing CHS is 10 times higher in patients with hyperperfusion than those without [2], and ICH develops in 3.3% of patients with hyperperfusion vs. only 0.24% of those without [4]. In fact, there is evidence that identifying patients at high risk of hyperperfusion and treating them early help to reduce the incidence of ICH and lead to better prognosis [5, 6].

Clinical studies [2, 7] have identified some risk factors for CHS or cerebral hyperperfusion, such as hypertension, high-grade ICA stenosis, decreased cerebral vasoreactivity and contralateral ICA stenosis. Most of these factors are associated closely with hemodynamics in the cerebral circulation. In fact, emerging evidence supports that a pre-operative evaluation of cerebral hemodynamic status may help to identify patients at high risk of post-operative hyperperfusion [5, 6, 8]. The cerebral circulation possesses many collateral vessels which play an important role in maintaining cerebral perfusion in case occlusive disease develops in the feeding arteries of the cerebral circulation [9]. Cerebral collateral vessels are commonly divided into primary and secondary collateral pathways, with the former constituted mainly by the Circle of Willis (CoW), while the latter by the ophthalmic artery and leptomeningeal vessels [9]. Many studies have demonstrated that the status of primary collateral flows is a determinant factor for clinical symptoms and outcomes of intervention in patients with severe ICA stenosis [1012]. Although the secondary collaterals may also play some roles in compensating for severe ischemia [13, 14], their compensatory capability seems to be limited [15].

Theoretically, a complete CoW is able to maintain sufficient cerebral perfusion when any single cerebral feeding artery is occluded. However, this ability can be impaired by a topological variation in CoW and coexistence of stenoses in multiple feeding arteries. In fact, a complete CoW structure exists in only about 50% of the population, with various incomplete CoW structures existing in the remaining population [16, 17]. The role of the anatomy of the CoW in regulating cerebral blood flows has been well described [17]; whereas, it remains unclear how the anatomy of the CoW influences post-CAS hyperperfusion, particularly when occlusive disease is present in multiple cerebral feeding arteries. To answer this question, we have developed a novel computational model of the cerebral circulation which is capable of describing cerebral hemodynamics under various physiological/pathological conditions.

Methods

Cardiovascular model

A cardiovascular model developed in our previous studies [18] has been extended to include the cerebral circulation (see Figure 1(A, B)). The cardiovascular model describes pulse wave propagation in the largest 83 arteries with a one-dimensional (1-D) sub-model while describing hemodynamics in the remaining cardiovascular system (including the microcirculation, the pulmonary circulation and the heart) with a lumped parameter (0-D) sub-model. The intracranial cerebral artery network is comprised by the 18 largest arteries (see Figure 1(B)) and the downstream vascular system corresponding to each cerebral efferent artery is divided into arteriolar, capillary, venular and venous compartments (see Figure 2). To reduce the degree of the complexity of the cerebral circulation model, we have neglected the smaller cerebral arteries and the secondary collateral vessels. Cerebral venous flows in different cerebral regions were assumed to converge to the neck veins. An extravascular pressure was applied to the intracranial vessels to account for the critical closing pressure (p e in Figure 2) at which blood flow to the brain stops. In addition, the cerebral arteriolar resistances and venous resistances were set to be changeable with cerebral hemodynamic conditions. These characteristics enabled the model to simulate cerebral hemodynamics under various physiological/pathological conditions.

Figure 1
figure 1

Schematic description of a computational model of the cardiovascular system (Panel A). The cerebral circulation is modeled as a portion of the cardiovascular system (Panel B). The cerebral feeding arteries are No. 40/47 - left/right ICA and No. 56 - Basilar artery (BA), communicating arteries are No. 59/69 - left/right posterior communicating artery (PCoA) and No. 64 - anterior communicating artery (ACoA), and cerebral efferent arteries are No. 58/70 - the right/left posterior cerebral artery II (PCA II), No. 61/67 - right/left middle cerebral artery (MCA) and No. 63/65 - right/left anterior cerebral artery II (ACA II). Panel C illustrates the nine typical CoW structures. The value given under each structure indicates the appearance frequency of the structure in the population. The arrows denote the locations of missing or fetal-type arteries.

Figure 2
figure 2

Lumped parameter modeling of the distal vascular systems of the intracranial cerebral circulation. The vascular systems distal to each cerebral efferent artery are divided into four compartments (namely, arteriolar, capillary, venular and venous compartments) with each compartment represented by a certain assembly of lumped parameters. The flows through all the distal vasculatures are assumed to converge into the neck veins. Please see the text for more details.

1-D governing equations for pulse wave propagation in the arteries

Blood flows in an artery were described by the 1-D mass and momentum conservation equations derived by integrating the three-dimensional mass conservation and Navier-Stokes equations over the cross-section of the artery [19].

A t + A U z = 0 ,
(1)
U t + z U 2 2 + P ρ = K υ U A
(2)

The system of Eqs. 1 and 2 was completed by a pressure-area relationship that has been previously used to describe the deforming response of arterial wall to changes in transmural pressure [1723].

P = P 0 + P e + β A - A 0 , with β = π E h 0 r 0 ( 1 - σ 2 ) .
(3)

Here, t is the time, z the axial coordinate along the artery; and ρ the blood density (ρ ≈ 1.06 g/cm3); A, U and P represent the lumen area, mean flow velocity and intravascular blood pressure, respectively; and K υ is the coefficient of the viscous term; P 0 is the reference pressure at A = A 0 and was set to be 85 mmHg; P e is the extravascular pressure; E is the Young's modulus; h 0 the wall thickness; r 0 the radius of the artery at the reference pressure; and σ the Poisson's ratio, here taken to be 0.5 by assuming arterial wall to be incompressible.

It is noted that, the cross-sectional velocity profile of blood flow changes transiently over a cardiac cycle and varies along the arterial system; this raises an issue as to how to correctly model the convective and viscous terms when reducing a three-dimensional blood flow model into a 1-D model. At this point, many modeling methods have been proposed based on various assumptions [17, 18, 20, 2427]. In this study, we employed a relatively simple modeling method in which the coefficient of the viscous term in Eq.2 is taken to be -8πυ (υ being the kinematic viscosity of blood ≈ 0.045 cm2/s) based on a Poiseuille flow assumption [20, 21], while the 1-D convective term is derived by assuming a flat velocity profile [17, 1922]. These assumptions led to several simplifications in the numerical treatment of flow conditions at the boundaries [19, 21]. Meanwhile, the error induced by the assumptions in the prediction of blood flow distribution in the cardiovascular system should be negligible since normal large arteries generate fairly less blood pressure loss in comparison with the downstream resistant micro-vasculatures [18].

From Eq. 3, arterial transmural pressure is related linearly to the change in arterial radius relative to its reference value. According to the data reported in previous experimental studies [28, 29], the linear relation is acceptable when arterial transmural pressure varies within the physiological range (e.g., from diastolic to systolic pressure). Previous computational studies [1723] have indeed demonstrated that employing the relation does not prevent a reasonable prediction of pulse wave propagation in large arteries. However, it should be noted that a linear pressure-radius relation fails to be proper when arterial transmural pressure varies beyond the general physiological range. Experimental studies [30] have demonstrated that when transmural pressure is reduced progressively from an over systolic to minus value, the pattern of arterial wall deformation changes from stretching to buckling and collapsing, exhibiting a highly non-linear pressure-radius relation.

Flows in different arteries were linked by imposing the conservation of mass and continuity of total pressure at the bifurcations [1723].

Stenosis model

Eq. 2 cannot fully account for the pressure drop induced by an arterial stenosis. To compensate for this limitation, an experiment-based empirical stenosis model [31] has been incorporated to relate the stenosis-induced pressure drop to the geometry of stenosis:

Δ P = K v μ A 0 D 0 Q + K t ρ 2 A 0 2 A 0 A s - 1 2 Q Q + K u ρ L s A 0 Q ° ,
(4)

where ΔP and Q denote pressure drop and flow rate through the stenosis, respectively; Q is the time derivative of Q, A 0 and A s refer to the cross-sectional areas of the normal and stenotic segments, respectively, L s represents the stenosis length, and μ is the blood viscosity. Further, K v, K t and K u are empirical coefficients, with K v = 32(0.83L s + 1.64D s)×(A 0/A s)2/D 0, K t = 1.52, and K u = 1.2, where D 0 and D s are the diameters corresponding to A 0 and A s. The degree (severity) of stenosis is defined as the percentage reduction in arterial diameter (= (1-D s/D 0) ×100%).

Governing equations for the 0-D sub-model of the cerebral circulation

Following the general 0-D modeling method [19, 32, 33], the viscous resistance, blood inertia and compliance of each vascular segment were mimicked respectively by three electric components (resistor (R), inductor (L) and capacitor (C)). In analogy to the principles of electric circuit, the governing equations were formulated by imposing mass and momentum conservation along the flow pathway (from arterioles to veins) (see Figure 2).

At a 'capacitor' component, mass conservation reads

d V i , j d t = Q i , j - Q i + 1 , j ,
(5)

and at an 'inductor' component, momentum conservation reads

d Q i , j d t = P i - 1 , j - Q i , j R i , j - P i , j L i , j ,
(6)

where V, Q and P represent blood volume, flow rate and blood pressure, respectively; P is related to V by P = V/C + P e, with P e being the critical closing pressure, here taken to be 18.1 mmHg in accord with the cerebral autoregulation curve (see Figure 3). The subscript 'i' is a sequence number that is increased from the arterial side toward the venous side (i I = [1, 4]); whereas 'j' denotes the labeled number of the cerebral efferent arteries in the 1-D sub-model (j J = {58, 61, 63, 65, 67, 70}) (see Figure 2).

Figure 3
figure 3

Relationship between cerebral perfusion pressure and perfusion rate. The continuous line represents the fitted fourth-degree polynomial function based on the experimental data.

Specially, cerebral vein may collapse due to the effect of extravascular pressure (P e), resulting in a varying venous resistance. This has been accounted for by modeling cerebral venous resistance (R 5, j in Figure 2) as a function of cerebral venous blood pressure (P 4, j ), downstream neck venous pressure (P nv) and extravascular pressure (P e) [34, 35].

R 5, j = R v0, j P 4, j - P nv P 4, j - P e ,
(7)

where R v0, j is a constant venous resistance component.

Cerebral venous flow (Q 5, j in Figure 2) can be calculated as

Q 5, j = P 4, j - P nv R 5, j .
(8)

Substituting Eq. 7 into Eq. 8, one gets

Q 5, j = P 4, j - P e R v0, j .
(9)

From Eq. 9, cerebral venous flow is independent of the extracranial venous pressure (P nv) as far as P e is larger than P nv [36].

0-D modeling of other portions of the cardiovascular system, such as the pulmonary circulation, the heart, has been described in detail in [18, 19].

Numerical methods

The equations system of the cardiovascular model consists of a 1D partial differential-algebraic sub-system coupled with a 0D ordinary differential-algebraic sub-system. The two sub-systems were solved numerically using the two-step Lax-Wendroff method and a fourth-order Runge-Kutta method, respectively. The solutions of the sub-systems were then linked at the 0-1D interfaces where mass and momentum conservation is imposed. More details on the numerical methods employed to treat flow conditions at the bifurcations and the 0-1D interfaces have been given elsewhere [19].

Cerebral autoregulation

Generally, cerebral autoregulation is a dynamic process, for instance, it takes several to several tens of seconds to restore cerebral perfusion upon an abrupt change in perfusion pressure [3739]. For the present problem, since a chronic ICA stenotic disease imposes a long-term influence on cerebral perfusion, a static cerebral pressure-perfusion rate relationship (herein termed cerebral autoregulation curve) would be sufficient, and which has herein been constructed by fitting a fourth-degree polynomial function to available experimental data [4042] (see Figure 3). Note that the data points that deviate apparently from the neighboring data points have been considered as noise data and removed from the input data of function fitting. The perfusion rate represents the ratio of real flow rate to normal flow rate. The autoregulation curve was incorporated into the model by adjusting the cerebral arteriolar resistances via a negative feedback process.

R n + 1 = R n 1 - α Q ̄ T - Q ̄ n Q ̄ T
(10)

Here R represents the arteriolar resistance corresponding to each cerebral efferent artery and Q the mean flow rate averaged over a cardiac cycle. Since resistance adjustment was performed at intervals of a cardiac cycle, the upper subscript 'n' denotes current cardiac cycle, whereas 'n+1' indicates the next cardiac cycle. QT denotes the target flow rate calculated from the cerebral autoregulation curve. α is the under-relaxation factor used to stabilize the numerical simulation (here taken to be 0.9).

Anatomical variations in CoW

Based on the data collected from the literature [16, 17, 43, 44], we categorized the frequently observed anatomical variations in CoW into nine types (see Figure 1(C)). Each type has a specific frequency of appearance in the population, with type 1 being the most prevalent structure.

Physiological data

The geometrical parameters of the cerebral arteries were assigned based on the data reported in the literature [17, 24, 27] and the in vivo data available in our lab. The elastic parameters of the arteries have been estimated according to an elastic modulus-artery radius relationship constructed based on experimental data [26]. The reference flow rate through each cerebral efferent artery at normal perfusion pressure was assigned based on previously reported data [24, 27, 43, 45], giving a flow division among PCA II, MCA and ACA II of 0.96:2:1 and a total cerebral flow rate of about 12 ml/s. The assigned data for the cerebral circulation model are summarized in Table 1. The data for other portions of the cardiovascular system have been given in our previous studies [18, 19].

Table 1 Physiological data of the cerebral circulation

Computation conditions

Each set of computation comprised three steps, with the computation for each step being continuously run for 30 cardiac cycles to guarantee the convergence of computation (inter-cardiac cycle error for mean flow rate within 0.1%): during the first 30 cardiac cycles (Step I), the reference cerebral arteriolar resistances are estimated under normal perfusion conditions (in the absence of artery stenosis); during the second 30 cardiac cycles (Step II), stenoses are introduced in certain cerebral feeding arteries (including at least one or both of the ICA) and the cerebral arteriolar resistances are further modified to match the cerebral autoregulation curve; and at the beginning of the last 30 cardiac cycles (Step III), an ICA stenosis is suddenly removed to simulate CAS. The post-CAS hyperperfusion rate in each cerebral efferent artery is calculated at the end of step III. It is noted that we have herein assumed that cerebral arteriolar resistances do not change immediately after CAS in order to simulate the largest post-CAS hyperperfusion rate.

The nine types of CoW structure illustrated in Figure 1(C) were studied. Each type of CoW structure was further investigated in combination with three distribution patterns of stenosis in the cerebral feeding arteries: (1) unilateral ICA stenosis, (2) bilateral ICA stenosis, and (3) coexisting unilateral ICA stenosis and basilar artery (BA) stenosis. The degree of each stenosis was set uniformly to be 75% to represent a severe stenotic condition. Heart rate has been fixed at 60BPM in all the computations.

Definition of hyperperfusion rate

Hyperperfusion rate (C H) was defined as the percentage change of post-CAS flow rate relative to the pre-CAS value [4].

C H = Q ̄ a Q ̄ b - 1 × 1 0 0 %
(11)

Here, Q b and Q a refer respectively to the mean flow rates before and after CAS.

Results

Computed flow rates through the left/right ICA and the BA for different CoW structures under normal conditions

The computed mean flow rates through the left/right ICA and the BA for three CoW structures (Types 1 to 3) under normal conditions (in the absence of feeding artery stenosis) are given in Table 2 in comparison with the corresponding in vivo data [43]. It was observed that the pattern of flow division among the three cerebral feeding arteries depends strongly on the CoW structure. And the computations reasonably captured the flow division patterns described by the in vivo data.

Table 2 Mean flow rates through the cerebral feeding arteries computed for three CoW structures under normal conditions (compared with measured data [43])

Hemodynamics before and after CAS

The type 6 CoW structure was taken as an example to illustrate the hemodynamic changes corresponding to onset of ICA stenosis and CAS. Herein, a stenosis was generated in the left ICA at the beginning of Step II and removed at the beginning of Step III. Before CAS, the cerebral distal resistance corresponding to each cerebral efferent artery was gradually modified from its initial value to regulate the flow rate toward the target value (see Figure 4). After CAS, the flow rates through the left MCA and ACA II were remarkably increased; whereas the flow rates through other cerebral efferent arteries were less changed. The computed hyperperfusion rates for the left MCA and ACA II were both larger than 100% (being 171% and 173%, respectively), indicating the occurrence of post-CAS hyperperfusion in the left MCA and ACA II territories.

Figure 4
figure 4

Computed changes in distal resistances of (Panel A) and mean flow rates through (Panel B) the efferent arteries of the type 6 CoW structure before and after CAS (unilateral left ICA stenosis). In the absence of ICA stenosis, the distal resistance corresponding to each cerebral efferent artery is gradually changed from its initial value (roughly assigned at the beginning of computation) to a reference value (step I). The distal resistances of the left MCA and ACA II are further reduced by about 50% following the onset of the left ICA stenosis (step II). After CAS is implemented at the beginning of step III, marked flow overshoots through the left MCA and ACA II appear.

The pre- and post-CAS blood flow/pressure waves in the six cerebral efferent arteries are plotted together with the results for normal conditions (in the absence of ICA stenosis) in Figure 5. As expected, the pressure/flow waves in the hyperperfusion-affected arteries (left MCA, ACA II) changed significantly with the onset and removal of the ICA stenosis. In contrast, the pressure/flow waves in the remaining cerebral efferent arteries showed little change.

Figure 5
figure 5

Computed changes in pressure/flow waves in the efferent arteries of the type 6 CoW structure before and after CAS: continuous line (normal), dashed line (in the presence of left ICA stenosis), dash-dotted line (after removal of the stenosis).

Figure 6 shows the transmural pressures distal to the left MCA (A) and right MCA (B). The transmural pressures distal to the left MCA were reduced as the ICA stenosis is present and increased after the stenosis is removed via CAS; whereas, those distal to the right MCA showed little change during the entire process.

Figure 6
figure 6

Computed changes in transmural pressures distal to the left MCA (Panel A) and right MCA (Panel B) of the type 6 CoW structure before and after CAS (unilateral left ICA stenosis). The onset and removal of the left ICA stenosis induce pronounced changes in transmural pressures distal to the left MCA in which post-CAS hyperperfusion occurs, but have little influence on transmural pressures distal to the right MCA.

Hyperperfusion rates in the case of unilateral ICA stenosis

Figure 7 shows the computed post-CAS hyperperfusion rates for the nine types of CoW structure in the cases of unilateral left ICA stenosis (A) and right ICA stenosis (B). The type 6 and type 4 CoW structures were found to induce hyperperfusion in the left ICA stenosis case and the right ICA stenosis case, respectively.

Figure 7
figure 7

Polar plots of the computed post-CAS hyperperfusion rates in the cerebral efferent arteries in the cases of unilateral left ICA stenosis (Panel A) and right ICA stenosis (Panel B).

Hyperperfusion rates in the case of bilateral ICA stenosis

In the case of bilateral ICA stenosis, we selectively removed the left or right ICA stenosis at each time. In this case, the types 4, 6, 8 and 9 CoW structures were found to induce post-CAS cerebral hyperperfusion (see Figure 8). Unlike the cases of the types 4, 6 structures, cerebral hyperperfusion induced by the types 8, 9 structures was not sensitive to the location side of the removed ICA stenosis and present in both hemispheres.

Figure 8
figure 8

Polar plots of the computed post-CAS hyperperfusion rates in the cerebral efferent arteries in the case of bilateral ICA stenosis: removal of the left ICA stenosis (Panel A), removal of the right ICA stenosis (Panel B).

Hyperperfusion rates in the case of coexisting BA stenosis and unilateral ICA stenosis

When BA stenosis existed concurrently with unilateral ICA stenosis, the types 4, 6 CoW structures that have been found to induce post-CAS hyperperfusion either in the case of left ICA stenosis or in the case of right ICA stenosis induced hyperperfusion in both cases (see Figure 9). Moreover, there was a significant enlargement of the cerebral region affected by post-CAS hyperperfusion (compare Figure 7 with Figure 9).

Figure 9
figure 9

Polar plots of the computed post-CAS hyperperfusion rates in the cerebral efferent arteries in the case of coexisting BA stenosis and unilateral ICA stenosis: left ICA stenosis (Panel A), right ICA stenosis (Panel B).

Discussion

The pre-CAS status of cerebral hemodynamics has been found to be an important factor for assessing post-CAS cerebral hyperperfusion in patients with severe ICA stenosis [5, 6, 8]. Cerebral hemodynamics may be evaluated directly by measuring intra-arterial blood flow using magnetic resonance angiography or transcranial Doppler sonography [10, 12, 46, 47] or indirectly via cerebral vasoreativity test [3], measurement of brain temperature [48] or brain oxygenation [49]. Despite the existence of these methods, an accurate measurement of blood flow rates in all the major cerebral arteries is yet difficult in clinical settings, which considerably hampers a full understanding of the collateral function of the cerebral artery network in pathological conditions. In contrast, the geometry of large cerebral arteries can nowadays be measured with satisfactory accuracy in clinical settings [11]. Computational hemodynamic modeling offers an alternative way to assess cerebral hemodynamics based on available geometrical data of cerebral arteries. A significant advantage of a computational model is that it allows us not only to quantify the blood flow rate in any cerebral artery of interest but also to evaluate the role of the entire cerebral artery network in regulating cerebral blood flows under various physiological/pathological conditions.

In this context, we have developed a 0-1D multi-scale model of the cerebral circulation (coupled with the global cardiovascular system) and applied it to investigate the influence of the anatomy of the CoW on post-CAS cerebral hyperperfusion. We should stress that although fully three-dimensional (3-D) modeling of the cerebral arteries can provide a more accurate and detailed description of blood flows compared to 0-D or 1-D modeling [45, 50]; it is not practical for the present study due to its high computational cost. In this study, each set of computation has to be run continually for tens of cardiac cycles, for which reason a modeling method that incurs lower computational cost would be more favorable. Reducing 3-D modeling into 1-D modeling significantly reduces the required computational effort, but at the expense of the loss of some geometric information, such as local artery surface shape, curvature and bifurcation structure. There is evidence that computation results obtained with 1-D and 3-D models of the cerebral arterial network are in good agreement in terms of mass-flow distribution and pressure drop along arteries [50]. The purpose of the present study determines that we are interested in mass-flow distribution rather than in local flow patterns; therefore, 1-D modeling should be a choice with a good balance between computational demand and physical detail for the description of the cerebral arterial network. For the peripheral portion of the cerebral circulation, its extreme complexity determines that 0-D modeling is the only practical way. Particularly, 0-D modeling allows us to readily account for certain physiological or pathological conditions by modifying model parameters.

The computed results presented in Figure 7, 8, 9 indicate that (1) the anatomy of the CoW and the distribution pattern of stenosis in the cerebral feeding arteries jointly determine the risk and the localization of post-CAS cerebral hyperperfusion; (2) the existence of BA stenosis or contralateral ICA stenosis tends to increase the risk of post-CAS hyperperfusion and enlarge the cerebral region affected by hyperperfusion; and (3) some CoW structures may induce post-CAS hyperperfusion in both hemispheres under certain conditions, such as the types 8, 9 CoW structures combined with bilateral ICA stenosis.

CoW structures susceptible to post-CAS hyperperfusion

As discussed above, the risk of a CoW structure for inducing post-CAS hyperperfusion should be always assessed in conjunction with the location of stenosis in the cerebral feeding arteries. According to the computed results, high-risk CoW structures are those (types 4 and 6) lacking collateral pathways from the contralateral ICA in the case of unilateral ICA stenosis and those lacking collateral pathways either from the contralateral ICA (types 4 and 6) or from the BA (types 8 and 9) in the case of bilateral ICA stenosis.

Influences of the diameters of the cerebral communicating arteries and the severity of ICA stenosis on post-CAS hyperperfusion

Although the major objective of this study is to investigate the influence of the anatomy of the CoW on post-CAS hyperperfusion, we should point out that the diameters of the cerebral communicating arteries and the severity of cerebral feeding artery stenosis are important factors for assessing post-CAS hyperperfusion as well. To confirm this viewpoint, we additionally carried out two sets of computation. In both sets of computations, the stenosis was located in the left ICA. In the first set of computation, the diameters of the ACoA and PCoAs of the type 1 CoW structure were reduced step by step to 10% of their reference values; whereas in the second set of computation, the severity of the ICA stenosis was increased gradually from 60% to 80% (herein the type 6 CoW structure was studied). The computed post-CAS hyperperfusion rates in the left MCA for the two cases are plotted in panels A and B of Figure 10, respectively. It was observed that even the complete CoW structure may induce hyperperfusion if the communicating arteries are small, and the value of hyperperfusion rate is highly sensitive to the severity of the ICA stenosis, especially when the severity is greater than 70%.

Figure 10
figure 10

Computed changes in hyperperfusion rate in the L. MCA with changes in: (A) the diameters of the cerebral communicating arteries (type 1 CoW structure), and (B) the degree of the ICA stenosis (type 6 CoW structure). In both cases, the stenosis is located in the left ICA. The diameter values shown in the horizontal axis of Figure 10(A) are normalized by the default values of the communicating arteries.

Sensitivity of cerebral arterial territory to post-CAS hyperperfusion

In the case of unilateral ICA stenosis, there is no apparent difference in hyperperfusion sensitivity between the ACA territory and the PCA territory; whereas, in the case of bilateral ICA stenosis, the ACA territory is more frequently affected by hyperperfusion in comparison with the PCA territory, with the appearance frequencies of hyperperfusion in these territories being 9 times vs. 3 times. The results support the clinical finding that a pronounced increase of blood flow velocity is often observed in the anterior part of the CoW immediately after CAS in the case of bilateral ICA stenosis [51]. This phenomenon can be explained from the fact that the PCA territory often possesses richer collateral pathways from the BA than the ACA territory, and the influence of this difference on cerebral perfusion is enhanced by the presence of bilateral ICA stenosis which reduces blood flows through both ICA, making cerebral perfusion rely more strongly on the blood flow supplied by the BA.

Change in transmural pressure after CAS

The most pronounced changes in transmural pressure after CAS were observed in the microcirculations of the cerebral territories subjected to post-CAS hyperperfusion. For instance, the proximal (arteriolar limb) capillary pressure in the left MCA territory (type 6 CoW structure with L. ICA stenosis) increased from a pre-CAS value of 18.6 mmHg to a post-CAS value of 46.9 mmHg (see Figure 6(A)). This phenomenon is attributable to a reduction in arteriolar resistance under pre-CAS ischemic conditions, which leads to a shift of pressure distribution from the arterioles toward the capillaries, ultimately resulting in a high capillary pressure when the proximal arterial perfusion pressure is recovered after CAS. The prediction is consistent with the results of the ischemia-reperfusion experiments on the isolated dog hind limb [52]; whereas, whether the similar phenomenon occurs in vivo in the human cerebral circulation remains not well known. If it were true, it might augment capillary leakage, increase the risk of edema, and hence be another causative factor for post-CAS CHS in addition to increased blood flow rate [53].

Limitations

A major limitation of this study is the absence of a sufficient comparison between model predictions and in vivo measurements. Actually, so far, we are not aware of any in vivo studies that systemically investigate the relationship between post-CAS hyperperfusion and the anatomy of the cerebral artery network with account of the distribution pattern of stenosis in the cerebral feeding arteries. At this point, further in vivo studies would be required to confirm the findings of the present study. Another limitation of this study may arise from the exclusion of the secondary cerebral collateral vessels from the present model, which potentially makes the model overestimate post-CAS hyperperfusion rate. Moreover, the CoW structures investigated in this study are limited to those illustrated in Figure 1(C), other CoW structures, such as those described elsewhere [14, 18, 44], would deserve further studies. Finally, since we did not take into account cerebral autoregulation in post-CAS computation by assuming cerebral distal resistances to remain constant after CAS, the predicted results may represent the largest values of post-CAS hyperperfusion rates. The dynamic cerebral autoregulation has been found to be significantly impaired in some patients with severe ICA stenoses [3, 7, 54]. After CAS, the immediate restoration of perfusion pressure does not guarantee an immediate sufficient restoration of cerebral autoregulation. In fact, hyperperfusion or impaired cerebrovascular reserve has been identified several days after CAS [3, 7, 55], indicating that cerebral autoregulation may remain insufficient for a fairly long time after CAS in some patients. In this sense, under in vivo conditions, the values of post-CAS hyperperfusion rates should be time-dependent, changing in close association with the post-CAS restoration of cerebral autoregulation. So far, the post-CAS restoring process of cerebral autoregulation remains not fully understood and seems to be strongly patient-specific [55], preventing us from developing a general model for describing post-CAS cerebral autoregulation. This limitation might be overcome if sufficient experimental data would be reported in the future.

Despite these limitations, some model-based findings regarding the factors tending to increase the risk of cerebral hyperperfusion (e.g., the existence of contralateral ICA stenosis, high-grade ICA stenosis) and the sensitivity of cerebral arterial territory to post-CAS hyperperfusion are in agreement with previous clinical findings [2, 7, 51]. Other model-based findings regarding the influences of BA stenosis and the diameters of the cerebral communicating arteries on the risk of post-CAS cerebral hyperperfusion and the change in capillary transmural pressure after CAS are reported for the first time. These findings, though awaiting further experiment-based confirmation, are of potential significance in the assessment and treatment of cerebral hyperperfusion.

Conclusions

Using a computational model, this study demonstrated the importance of considering the anatomy of the CoW in assessing the risk of post-CAS cerebral hyperperfusion. Particularly, the finding that the anatomy of the CoW and the distribution pattern of stenosis in the cerebral feeding arteries jointly determine the risk and localization of post-CAS cerebral hyperperfusion suggests that a patient-specific hemodynamic analysis aimed to help physicians identify patients at high risk of post-operative cerebral hyperperfusion should account for the combined effect of the anatomy of cerebral arteries and cerebral feeding artery stenoses on cerebral hemodynamics.

References

  1. Mas JL, Chatellier G, Beyssen B, Branchereau A, Moulin T, Becquemin JP, Larrue V, Lievre M, Leys D, Bonneville JF, Watelet J, Pruvo JP, Albucher JF, Viguier A, Piquet P, Garnier P, Viader F, Touze E, Giroud M, Hosseini H, Pillet JC, Favrole P, Neau JP, Ducrocq X, EVA-3S Investigators: Endarterectomy versus stenting in patients with symptomatic severe carotid stenosis. N Engl J Med 2006, 355: 1660–1671. 10.1056/NEJMoa061752

    Article  Google Scholar 

  2. Adhiyaman V, Alexander S: Cerebral hyperperfusion syndrome following carotid endarterectomy. QJM 2007, 100: 239–244. 10.1093/qjmed/hcm009

    Article  Google Scholar 

  3. Hosoda K, Kawaguchi T, Shibata Y, Kamei M, Kidoguchi K, Koyama J, Fujita S, Tamaki N: Cerebral vasoreactivity and internal carotid artery flow help to identify patients at risk for hyperperfusion after carotid endarterectomy. Stroke 2001, 32: 1567–1573. 10.1161/01.STR.32.7.1567

    Article  Google Scholar 

  4. Piepgras DG, Morgan MK, Sundt TM Jr, Yanagihara T, Mussman LM: Intracerebral hemorrhage after carotid endarterectomy. J Neurosurg 1988, 68: 532–536. 10.3171/jns.1988.68.4.0532

    Article  Google Scholar 

  5. Dalman JE, Beenakkers IC, Moll FL, Leusink JA, Ackerstaff RG: Transcranial Doppler monitoring during carotid endarterectomy helps to identify patients at risk of postoperative hyperperfusion. Eur J Vasc Endovasc Surg 1999, 18: 222–227. 10.1053/ejvs.1999.0846

    Article  Google Scholar 

  6. Yoshimoto T, Shirasaka T, Yoshizumi T, Fujimoto S, Kaneko S, Kashiwaba T: Evaluation of carotid distal pressure for prevention of hyperperfusion after carotid endarterectomy. Surg Neurol 2005, 63: 554–557. 10.1016/j.surneu.2004.06.016

    Article  Google Scholar 

  7. Kaku Y, Yoshimura S, Kokuzawa J: Factors predictive of cerebral hyperperfusion after carotid angioplasty and stent placement. AJNR Am J Neuroradiol 2004, 25: 1403–1408.

    Google Scholar 

  8. Kablak-Ziembicka A, Przewlocki T, Pieniazek P, Musialek P, Tekieli L, Rosławiecka A, Motyl R, Zmudka K, Tracz W, Podolec P: Predictors of cerebral reperfusion injury after carotid stenting: the role of transcranial color-coded Doppler ultrasonography. J Endovasc Ther 2010, 17: 556–563. 10.1583/09-2980.1

    Article  Google Scholar 

  9. Liebeskind DS: Collateral circulation. Stroke 2003, 34: 2279–2284. 10.1161/01.STR.0000086465.41263.06

    Article  Google Scholar 

  10. Maas MB, Lev MH, Ay H, Singhal AB, Greer DM, Smith WS, Harris GJ, Halpern E, Kemmling A, Koroshetz WJ, Furie KL: Collateral vessels on CT angiography predict outcome in acute ischemic stroke. Stroke 2009, 40: 3001–3005. 10.1161/STROKEAHA.109.552513

    Article  Google Scholar 

  11. Vernieri F, Pasqualetti P, Matteis M, Passarelli F, Troisi E, Rossini PM, Caltagirone C, Silvestrini M: Effect of collateral blood flow and cerebral vasomotor reactivity on the outcome of carotid artery occlusion. Stroke 2001, 32: 1552–1558. 10.1161/01.STR.32.7.1552

    Article  Google Scholar 

  12. Hartkamp MJ, van Der Grond J, van Everdingen KJ, Hillen B, Mali WP: Circle of Willis collateral flow investigated by magnetic resonance angiography. Stroke 1999, 30: 2671–2678. 10.1161/01.STR.30.12.2671

    Article  Google Scholar 

  13. Schneider UC, Schilling L, Schroeck H, Nebe CT, Vajkoczy P, Woitzik J: Granulocyte-macrophage colony-stimulating factor-induced vessel growth restores cerebral blood supply after bilateral carotid artery occlusion. Stroke 2007, 38: 1320–1328. 10.1161/01.STR.0000259707.43496.71

    Article  Google Scholar 

  14. Todo K, Kitagawa K, Sasaki T, Omura-Matsuoka E, Terasaki Y, Oyama N, Yagita Y, Hori M: Granulocyte-macrophage colony-stimulating factor enhances leptomeningeal collateral growth induced by common carotid artery occlusion. Stroke 2008, 39: 1875–1882. 10.1161/STROKEAHA.107.503433

    Article  Google Scholar 

  15. Hofmeijer J, Klijn CJ, Kappelle LJ, Van Huffelen AC, Van Gijn J: Collateral circulation via the ophthalmic artery or leptomeningeal vessels is associated with impaired cerebral vasoreactivity in patients with symptomatic carotid artery occlusion. Cerebrovasc Dis 2002, 14: 22–26. 10.1159/000063719

    Article  Google Scholar 

  16. Alpers BJ, Berry RG, Paddison RM: Anatomical studies of the circle of willis in normal brain. AMA Arch Neurol Psychiatry 1959, 81: 409–418.

    Article  Google Scholar 

  17. Alastruey J, Parker KH, Peiro J, Byrd SM, Sherwin SJ: Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows. J Biomech 2007, 40: 1794–1805. 10.1016/j.jbiomech.2006.07.008

    Article  Google Scholar 

  18. Liang FY, Takagi S, Himeno R, Liu H: Biomechanical characterization of ventricular-arterial coupling during aging: a multi-scale model study. J Biomech 2009, 42: 692–704. 10.1016/j.jbiomech.2009.01.010

    Article  Google Scholar 

  19. Liang FY, Takagi S, Himeno R, Liu H: Multi-scale modeling of the human cardiovascular system with applications to aortic valvular and arterial stenoses. Med Biol Eng Comp 2009, 47: 743–755. 10.1007/s11517-009-0449-9

    Article  Google Scholar 

  20. Mynard J, Nithiarasu P: A 1D arterial blood flow model incorporating ventricular pressure, aortic valve and regional coronary flow using the locally conservative Galerkin (LCG) method. Communications in Numerical Methods in Engineering 2008, 24: 367–417. 10.1002/cnm.1117

    Article  MathSciNet  Google Scholar 

  21. Formaggia L, Lamponi D, Tuveri M, Veneziani A: Numerical modeling of 1-D arterial networks coupled with a lumped parameters description of the heart. Comput Meth Biomech Biomed Eng 2006, 9: 273–288. 10.1080/10255840600857767

    Article  Google Scholar 

  22. Sherwin SJ, Franke V, Peiró J, Parker K: One-dimensional modelling of a vascular network in space-time variables. J Eng Math 2003, 47: 217–250.

    Article  Google Scholar 

  23. Mabotuwana T DS, Cheng LK, Pullan AJ: A model of blood flow in the mesenteric arterial system. Biomed Eng Online 2007, 6: 17. 10.1186/1475-925X-6-17

    Article  Google Scholar 

  24. Reymond P, Merenda F, Perren F, Rüfenacht D, Stergiopulos N: Validation of a one-dimensional model of the systemic arterial tree. Am J Physiol 2009, 297: 208–222.

    Google Scholar 

  25. Bessems D, Rutten M, van de Vosse Frans : A wave propagation model of blood flow in large vessels using an approximate velocity profile function. Journal of Fluid Mechanics 2007, 580: 145–168.

    Article  MathSciNet  Google Scholar 

  26. Olufsen MS, Peskin CS, Kim WY, Pedersen EM, Nadim A, Larsen J: Numerical simulation and experimental validation of blood flow in arteries with structured-tree outflow conditions. Ann Biomed Eng 2000, 28: 1281–1299.

    Article  Google Scholar 

  27. Devault K, Gremaud PA, Novak V, Olufsen MS, Vernieres G, Zhao P: Blood flow in the circle of Willis: Modeling and calibration. Multiscale Model Simul 2008, 7: 888–909. 10.1137/07070231X

    Article  MathSciNet  Google Scholar 

  28. Joseph C, Greenfield JR, Douglas M, Griggs JR: Relation between pressure and diameter in main pulmonary artery of man. Journal of Applied Physiology 1963, 18: 557–559.

    Google Scholar 

  29. Patel DJ, Janicki JS: Static elastic properties of the left coronary circumflex artery and the common carotid artery in dogs. Circ Res 1970, 27: 149–158.

    Article  Google Scholar 

  30. Drzewiecki G, Pilla JJ: Noninvasive measurement of the human brachial artery pressure-area relation in collapse and hypertension. Ann Biomed Eng 1998, 26: 965–974.

    Article  Google Scholar 

  31. Young DF, Tsai F: Flow characteristics in models of arterial stenoses. II. Unsteady flow. J Biomech 1973, 6: 547–559. 10.1016/0021-9290(73)90012-2

    Article  Google Scholar 

  32. Liang FY, Liu H: Simulation of hemodynamic responses to the Valsalva maneuver: an integrative computational model of the cardiovascular system and the autonomic nervous system. J Physiol Sci 2006, 56: 45–65. 10.2170/physiolsci.RP001305

    Article  Google Scholar 

  33. Shi YB, Lawford P, Hose R: Review of Zero-D and 1-D Models of Blood Flow in the Cardiovascular System. Biomed Eng Online 2011, 10: 33. 10.1186/1475-925X-10-33

    Article  Google Scholar 

  34. Lu K, Clark JW Jr, Ghorbel FH, Robertson CS, Ware DL, Zwischenberger JB, Bidani A: Cerebral autoregulation and gas exchange studied using a human cardiopulmonary model. Am J Physiol Heart Circ Physiol 2004, 286: 584–601. 10.1152/ajpheart.00594.2003

    Article  Google Scholar 

  35. Ursino M: A mathematical study of human intracranial hydrodynamics. Part I. The cerebrospinal fluid pulse pressure. Ann Biomed Eng 1988, 16: 379–401. 10.1007/BF02364625

    Article  Google Scholar 

  36. Auer L, Ishiyama N, Hodde K, Kleinert R, Pucher R: Effect of intracranial pressure on bridging veins in rats. J Neurosurg 1987, 67: 263–268. 10.3171/jns.1987.67.2.0263

    Article  Google Scholar 

  37. David T, Brown M, Ferrandez A: Auto-regulation and blood flow in the cerebral circulation. International Journal for Numerical Methods in Fluids 2003, 43: 701–713. 10.1002/fld.495

    Article  Google Scholar 

  38. Olufsen M, Tran H, Ottesen J: Modeling cerebral blood flow control during posture change from sitting to standing. Cardiovascular Engineering 2004, 4: 47–58.

    Article  Google Scholar 

  39. Panerai RB: Cerebral autoregulation: from models to clinical applications. Cardiovascular Engineering 2008, 8: 42–59. 10.1007/s10558-007-9044-6

    Article  Google Scholar 

  40. MacKenzie ET, Farrar JK, Fitch W, Graham DI, Gregory PC, Harper AM: Effects of hemorrhagic hypotension on the cerebral circulation. I. Cerebral blood flow and pial arteriolar caliber. Stroke 1979, 10: 711–718. 10.1161/01.STR.10.6.711

    Article  Google Scholar 

  41. Harper SL, Bohlen HG, Rubin MJ: Arterial and microvascular contributions to cerebral cortical autoregulation in rats. Am J Physiol 1984, 246: 17–24.

    Google Scholar 

  42. Dirnagl U, Pulsinelli W: Autoregulation of cerebral blood flow in experimental focal brain ischemia. J Cereb Blood Flow Metab 1990, 10: 327–336. 10.1038/jcbfm.1990.61

    Article  Google Scholar 

  43. Tanaka H, Fujita N, Enoki T, Matsumoto K, Watanabe Y, Murase K, Nakamura H: Relationship between variations in the circle of Willis and flow rates in internal carotid and basilar arteries determined by means of magnetic resonance imaging with semiautomated lumen segmentation: reference data from 125 healthy volunteers. AJNR Am J Neuroradiol 2006, 27: 1770–1775.

    Google Scholar 

  44. Papantchev V, Hristov S, Todorova D, Naydenov E, Paloff A, Nikolov D, Tschirkov A, Ovtscharoff W: Some variations of the circle of Willis, important for cerebral protection in aortic surgery--a study in Eastern Europeans. Eur J Cardiothorac Surg 2007, 31: 982–989. 10.1016/j.ejcts.2007.03.020

    Article  Google Scholar 

  45. Moore SM, Moorhead KT, Chase JG, David T, Fink J: One-dimensional and three-dimensional models of cerebrovascular flow. J Biomech Eng 2005, 127: 440–449. 10.1115/1.1894350

    Article  Google Scholar 

  46. Hendrikse J, Hartkamp MJ, Hillen B, Mali WP, van der Grond J: Collateral ability of the circle of Willis in patients with unilateral internal carotid artery occlusion: border zone infarcts and clinical symptoms. Stroke 2001, 32: 2768–2773. 10.1161/hs1201.099892

    Article  Google Scholar 

  47. Hoksbergen AW, Majoie CB, Hulsmans FJ, Legemate DA: Assessment of the collateral function of the circle of Willis: three-dimensional time-of-flight MR angiography compared with transcranial color-coded duplex sonography. AJNR Am J Neuroradiol 2003, 24: 456–462.

    Google Scholar 

  48. Murakami T, Ogasawara K, Yoshioka Y, Ishigaki D, Sasaki M, Kudo K, Aso K, Nishimoto H, Kobayashi M, Yoshida K, Ogawa A: Brain temperature measured by using proton MR spectroscopy predicts cerebral hyperperfusion after carotid endarterectomy. Radiology 2010, 256: 924–931. 10.1148/radiol.10090930

    Article  Google Scholar 

  49. Li JK-J, Wang TJ, Zhang HJ: Rapid noninvasive continuous monitoring of oxygenation in cerebral ischemia and hypoxia. Cardiovascular Engineering 2010, 10: 213–217. 10.1007/s10558-010-9116-x

    Article  Google Scholar 

  50. Grinberg L, Cheever E, Anor T, Madsen JR, Karniadakis GE: Modeling blood flow circulation in intracranial arterial networks: a comparative 3D/1D simulation study. Ann Biomed Eng 2011, 39: 297–309. 10.1007/s10439-010-0132-1

    Article  Google Scholar 

  51. Maltezos CK, Papanas N, Papas TT, Georgiadis GS, Dragoumanis CK, Marakis J, Maltezos E, Lazarides MK: Changes in blood flow of anterior and middle cerebral arteries following carotid endarterectomy: a transcranial Doppler study. Vasc Endovascular Surg 2007, 41: 389–396. 10.1177/1538574407302850

    Article  Google Scholar 

  52. Diana JN, Laughlin MH: Effect of ischemia on capillary pressure and equivalent pore radius in capillaries of the isolated dog hind limb. Circ Res 1974, 35: 77–101.

    Article  Google Scholar 

  53. Karapanayiotides T, Meuli R, Devuyst G, Piechowski-Jozwiak B, Dewarrat A, Ruchat P, Von Segesser L, Bogousslavsky J: Postcarotid endarterectomy hyperperfusion or reperfusion syndrome. Stroke 2005, 36: 21–26. 10.1161/01.STR.0000149946.86087.e5

    Article  Google Scholar 

  54. Reinhard M, Roth M, Muller T, Guschlbauer B, Timmer J, Czosnyka M, Hetzel A: Effect of carotid endarterectomy or stenting on impairment of dynamic cerebral autoregulation. Stroke 2004, 35: 1381–1387. 10.1161/01.STR.0000127533.46914.31

    Article  Google Scholar 

  55. Jorgensen LG, Schroeder TV: Defective cerebrovascular autoregulation after carotid endarterectomy. Eur J Vasc Surg 1993, 7: 370–379. 10.1016/S0950-821X(05)80252-X

    Article  Google Scholar 

Download references

Acknowledgements

This study was supported by Research and Development of the Next-Generation Integrated Simulation of Living Matter, a part of the Development and Use of the Next-Generation Supercomputer Project of the MEXT, Japan.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fuyou Liang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

FYL designed the study, constructed the computational model and drafted the manuscript. KF put forward the issue from a clinical perspective and was involved in drafting the manuscript. HL analyzed the computed results and corrected the manuscript. ST supervised the whole process of the study and corrected the manuscript. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Liang, F., Fukasaku, K., Liu, H. et al. A computational model study of the influence of the anatomy of the circle of willis on cerebral hyperperfusion following carotid artery surgery. BioMed Eng OnLine 10, 84 (2011). https://doi.org/10.1186/1475-925X-10-84

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-925X-10-84

Keywords