# Simulation of Contrast Agent Transport in Arteries with Multilayer Arterial Wall: Impact of Arterial Transmural Transport on the Bolus Delay and Dispersion

Authors: Min Xu; Xiao Liu; Ang Li; Yubo Fan; Anqiang Sun; Xiaoyan Deng; Deyu Li
Journal: The Scientific World Journal
 (2014)

---

## Abstract

One assumption of DSC-MRI is that the injected contrast agent is kept totally intravascular and the arterial wall is impermeable to contrast agent. The assumption is unreal for such small contrast agent as Gd-DTPA can leak into the arterial wall. To investigate whether the unreal assumption is valid for the estimation of the delay and dispersion of the contrast agent bolus, we simulated flow and Gd-DTPA transport in a model with multilayer arterial wall and analyzed the bolus delay and dispersion qualified by mean vascular transit time (MVTT) and the variance of the vascular transport function. Factors that may affect Gd-DTPA transport hence the delay and dispersion were further investigated, such as integrity of endothelium and disturbed flow. The results revealed that arterial transmural transport would slightly affect MVTT and moderately increase the variance. In addition, although the integrity of endothelium can significantly affect the accumulation of contrast agent in the arterial wall, it had small effects on the bolus delay and dispersion. However, the disturbed flow would significantly increase both MVTT and the variance. In conclusion, arterial transmural transport may have a small effect on the bolus delay and dispersion when compared to the flow pattern in the artery.


---

## Body

## 1. Introduction

Dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) has been shown to be a powerful technique to qualify cerebral blood flow and is playing an increasing role in diagnosis of acute ischemic stroke [1, 2]. However, the accuracy of perfusion parameters obtained from DSC-MRI is challenged by the delay and dispersion of the contrast agent bolus caused by the estimated arterial input function (AIF) [3, 4].Cerebral blood flow quantification requires knowledge of the arterial input function that is the concentration of contrast agent in the feeding vessel to the tissue of interest. In theory the AIF should be measured at the feeding vessel close to the tissue of interest [5, 6]. However, due to technical difficulties, it is commonly estimated from a major artery (e.g., the middle cerebral artery) for analysis of cerebral perfusion [5]. Delay and dispersion of the contrast agent may occur between the position of the AIF estimated and the tissue of interest, which would lead to underestimation of the blood flow [5].Another implicit assumption of the technique of DSC-MRI related to the estimation of AIF is that no or a negligible amount contrast agent penetrates the arterial wall, due to the fact that, in general, the vascular wall is considered irrelevant for the contrast mechanism in DSC-MRI. The assumption is tenable for the intact blood brain barrier (BBB) [7]. However, for a damaged BBB which usually occurs in the stroke and normal arterial wall except in the brain, the assumption is unreal. The commonly used contrast agent gadolinium diethylenetriominepentax-acetic acid (Gd-DTPA) with hydrodynamic diameter of less than 2 nm can easily leak into the arterial wall, since it has been demonstrated that the hydrophilic solute with the diameter less than 7.0 nm can be transported through the intercellular junction of endothelial cell into arterial wall [8]. Whether the unreal assumption is valid for the estimation of the delay and dispersion of the contrast agent bolus has never been verified. It seems that the leakage of contrast agent into the arterial wall may affect contrast agent transport, which might cause distortion of the contrast concentration-time curve and hence the delay and dispersion of the contrast agent bolus.To investigate the delay and dispersion effects of a bolus of contrast agent, two approaches were usually applied. One that has been commonly used was convolving the estimated AIF with a vascular transport function (VTF). The VTFs were quite differently assumed, ranging from the simple model of a single-exponential to the more sophisticated model of a feeding artery in series with small parallel vessels [3, 9]. The disadvantage of the method is that accuracy of the different VTF models is hard to assess. The other was the simulation of the blood flow and contrast agent transport using computational fluid dynamics (CFD) to describe the delay and dispersion errors of the contrast agent bolus [10–12]. The approach has been used to investigate the effect of carotid artery stenosis on the cerebral blood flow quantification and to simulate the dispersion along a simplified coronary artery with different stenosis under steady and unsteady flow condition. Although these studies indicate that CFD simulations were useful to investigate the delay and dispersion of the contrast agent bolus, all the studies were based on the unverified assumption of impermeable arterial wall.In the present study, to investigate whether the unreal basic assumption of DSC-MRI is valid for the estimation of the delay and dispersion of the contrast agent bolus, we formulated the arterial wall as a five-layer model and numerically simulated the flow and the transport of contrast agent in the model using CFD. This five-layer model included the endothelial glycocalyx layer (EGL), the endothelium, the intima, the internal elastic lamina (IEL), and the media, which were all treated as macroscopically homogeneous porous media. The effects of different factors that may affect contrast agent transport such as the integrity of endothelium and the disturbed flow after the stenosis of artery on the delay and dispersion of contrast agent were further analyzed.


## 2. Methods

### 2.1. Geometry of the Model

As a basic geometrical model, the arterial segment concerned was simplified as a straight axisymmetric cylinder. The inner diameter of model(D) and the wall thickness of the artery (t) are chosen to be D=3.7 mm and t=0.34 mm [13, 14]. The wall of the arterial segment was modeled as a five-layer structure and the thickness of each wall layer is shown in Figure 1(b) [15, 16], where the ratio of intima and media thickness was chosen to be approximately 0.75 [17].Schematic illustration of the computational geometry and boundary conditions. (a) The computational geometry of the stenosed model; (b) five-layer arterial wall with the thickness of each layer illustrated in the parentheses; (c) flow waveform at the inlet. EGL: endothelial glycocalyx layer; IEL: internal elastic lamina;D: inner diameter.
(a)
(b)
(c)For the stenosed model, the variation of the vessel radius along the stenosis was described using a cosine function and the reduction in the cross-sectional area of the lumen was set as 75% [18].


### 2.2. Governing Equations

#### 2.2.1. Lumen

The flow simulation in the lumen of the arterial segment was based on the incompressible Navier-Stokes and continuity equations:(1)ρ∂u∂t+u·∇u+∇p-μΔu=0,∇·u=0,
where u and p represent, respectively, the fluid velocity vector and the pressure. ρ and μ are the density and viscosity of blood (ρ=1050 kg·m−3, μ=3.45×10-3 kg·m−1 s−1).The mass transport of contrast agent Gd-DTPA in the flowing blood can be described by(2)∂c∂t+u·∇c-DGΔc=0,
where c is the concentration of Gd-DTPA, and DG is the diffusion coefficient of Gd-DTPA in blood. At 37°C, when assuming that the radius of Gd-DTPA (r) is 1 nm [19], DG can be calculated as 6.5847 × 10−11 m2 s−1 from the following Stokes-Einstein equation [20]:
(3)DG=kT6πμr,
where k is Boltzmann constant, and T is the absolute temperature.


#### 2.2.2. Arterial Wall Layers

The transmural flow across the arterial wall can be described by the Brinkman equation as follows [16]:
(4)ρεl∂ul∂t+μlKlul=-∇pl+μlεl∇2ul,(5)∇·ul=0,
where ul and pl represent, respectively, the superficial velocity vector and the pressure based on the volume averaged method. εl and Kl are the porosity and the hydraulic permeability of the wall layer concerned. For all layers, the viscosity of plasma μl was assumed to be 0.72 × 10−3 kg m−1 s−1 [16].The transport of Gd-DTPA across the arterial layers was modeled by the following equation [16]:
(6)∂cl∂t+1-σlul·∇cl-DlΔcl=0,
where cl is the superficial concentration of Gd-DTPA; Dl is the effective diffusivity of Gd-DTPA in the arterial layers; σl is the filtration reflection coefficient.


### 2.3. Parameters

In order to solve (4)~(6), we have to acquire the values of the 4 parameters in the equations, namely, the porosity (εl), the hydraulic permeability (Kl), the effective diffusivity of Gd-DTPA (Dl) and the filtration reflection coefficient (σl) for each layer of the arterial wall. However, because it is hard to measure the parameters, the transport parameters are scarce. In the present study, most of the parameters were obtained from theoretical models. The calculation of the parameters is presented in Appendix. The parameters used are summarized in Table 1.Table 1
Model parameters for each layer.
EGL
Endothelium
Intima
IEL
Media
Hydraulic permeability (K, (m2))
6.0383 × 10−18
1.7383 × 10−20
4.2 × 10−17
8.4 × 10−20
6.09 × 10−19
Effective diffusivity (De, (m2 s−1))
1.0128 × 10−10
3.0276 × 10−13
1.2219 × 10−10
2.4094 × 10−13
9.3392 × 10−14
Reflection coefficient (σ)
0.0555
0.1212
0.2514
0.2514
0.3617
Porosity (ε)
0.6735
0.0005
0.8025
0.002
0.258


### 2.4. Boundary Conditions

As shown in Figure1(a), flow transport equations ((1), (4)~(5)) are subject to the following boundary conditions.BC-A at the inlet of the lumen of the arterial segment, a time-dependent fully developed (parabolic) inlet flow velocity profile was used for the pulsatile flow simulation (Figure 1(c)) [21]. The mean velocity was chosen to be 0.24 m s−1 [22].BC-B the pressure at the outlet boundary of the artery lumen was set at 100 mm Hg.BC-C at the media-adventitia interface, a constant pressure boundary condition with 30 mm Hg was employed.BC-D symmetric conditions were set at the axis of symmetry.BC-E no viscous flow was set on both the axial ends of the arterial wall.The boundary conditions for the mass transport equations ((2) and (6)) are as follows.BC-1: the injection of the contrast agent at the inlet was described by a gamma-variate function [12]:
(7)ct,z=0=at-t0be-gt-t0,t>t00,t≤t0,wherea=1.013×10-3, b=2.142, and g=0.454 s−1. t0 was set to 3 s, which meant that the bolus of contrast was injected after the periodicity of the flow field was completely developed.BC-2: for other boundaries, zero concentration gradient was assumed.


### 2.5. Computation Procedures

The numerical simulations were carried out using a validated finite element algorithm COMSOL Multiphysics. Three full cardiac cycles (3 s) simulation of the pulsatile flow with a time step of 4 ms were carried out to achieve a periodic flow independent of the initialization. Based on the initial velocity field, the mass transport equations were solved coupling with flow transport equations for 37 s with time step from 1 ms to 4 ms depending on the temporal variations of contrast concentration.


### 2.6. Quantification of Dispersion

The effect of the dispersion of the contrast agent bolus on the AIF can be described by convolving the estimated AIF with a vascular transport function (VTF) [4, 23]:
(8)AIFt-t0=AIFtest⊗VTFt-t0,z,
where t0 is the bolus delay and z is distance to the inlet. Using the variance of the VTF around its mean, the degree of dispersion can be quantified:
(9)σVTF2z=∫0∞t-MVTTz2·VTFt,zdt∫0∞VTFt,zdt,
where the mean vascular transit time (MVTT) is given by the ratio of the first to the zeroth moment of the VTF and can be used to quantify the bolus delay:
(10)MVTTz=∫0∞t·VTFt,zdt∫0∞VTFt,zdt.
AIFest is prescribed directly from the contrast agent concentration at the inlet ((7)  c(t,z=0)) and AIF is computed as the area weighted average of the concentration of contrast agent on several lumen cross-sections between the inlet and the outlet perpendicular to the axial vessel direction (c(t,z)). To obtain the mean delay (MVTT) and the variance (σVTF2), it is not necessary to calculate VTF firstly, since they can be calculated directly with the zeroth, first, and second integral moment of the concentration of the contrast agent [4]:
(11)σVTF2z=c2zc0z-c20c00+c10c002-c1zc0z2,MVTTz=c1zc0z-c10c00.
The moments of the concentration of the contrast agent are
(12)cnz=∫0∞tn·ct,zdt.


## 2.1. Geometry of the Model

As a basic geometrical model, the arterial segment concerned was simplified as a straight axisymmetric cylinder. The inner diameter of model(D) and the wall thickness of the artery (t) are chosen to be D=3.7 mm and t=0.34 mm [13, 14]. The wall of the arterial segment was modeled as a five-layer structure and the thickness of each wall layer is shown in Figure 1(b) [15, 16], where the ratio of intima and media thickness was chosen to be approximately 0.75 [17].Schematic illustration of the computational geometry and boundary conditions. (a) The computational geometry of the stenosed model; (b) five-layer arterial wall with the thickness of each layer illustrated in the parentheses; (c) flow waveform at the inlet. EGL: endothelial glycocalyx layer; IEL: internal elastic lamina;D: inner diameter.
(a)
(b)
(c)For the stenosed model, the variation of the vessel radius along the stenosis was described using a cosine function and the reduction in the cross-sectional area of the lumen was set as 75% [18].


## 2.2. Governing Equations

### 2.2.1. Lumen

The flow simulation in the lumen of the arterial segment was based on the incompressible Navier-Stokes and continuity equations:(1)ρ∂u∂t+u·∇u+∇p-μΔu=0,∇·u=0,
where u and p represent, respectively, the fluid velocity vector and the pressure. ρ and μ are the density and viscosity of blood (ρ=1050 kg·m−3, μ=3.45×10-3 kg·m−1 s−1).The mass transport of contrast agent Gd-DTPA in the flowing blood can be described by(2)∂c∂t+u·∇c-DGΔc=0,
where c is the concentration of Gd-DTPA, and DG is the diffusion coefficient of Gd-DTPA in blood. At 37°C, when assuming that the radius of Gd-DTPA (r) is 1 nm [19], DG can be calculated as 6.5847 × 10−11 m2 s−1 from the following Stokes-Einstein equation [20]:
(3)DG=kT6πμr,
where k is Boltzmann constant, and T is the absolute temperature.


### 2.2.2. Arterial Wall Layers

The transmural flow across the arterial wall can be described by the Brinkman equation as follows [16]:
(4)ρεl∂ul∂t+μlKlul=-∇pl+μlεl∇2ul,(5)∇·ul=0,
where ul and pl represent, respectively, the superficial velocity vector and the pressure based on the volume averaged method. εl and Kl are the porosity and the hydraulic permeability of the wall layer concerned. For all layers, the viscosity of plasma μl was assumed to be 0.72 × 10−3 kg m−1 s−1 [16].The transport of Gd-DTPA across the arterial layers was modeled by the following equation [16]:
(6)∂cl∂t+1-σlul·∇cl-DlΔcl=0,
where cl is the superficial concentration of Gd-DTPA; Dl is the effective diffusivity of Gd-DTPA in the arterial layers; σl is the filtration reflection coefficient.


## 2.2.1. Lumen

The flow simulation in the lumen of the arterial segment was based on the incompressible Navier-Stokes and continuity equations:(1)ρ∂u∂t+u·∇u+∇p-μΔu=0,∇·u=0,
where u and p represent, respectively, the fluid velocity vector and the pressure. ρ and μ are the density and viscosity of blood (ρ=1050 kg·m−3, μ=3.45×10-3 kg·m−1 s−1).The mass transport of contrast agent Gd-DTPA in the flowing blood can be described by(2)∂c∂t+u·∇c-DGΔc=0,
where c is the concentration of Gd-DTPA, and DG is the diffusion coefficient of Gd-DTPA in blood. At 37°C, when assuming that the radius of Gd-DTPA (r) is 1 nm [19], DG can be calculated as 6.5847 × 10−11 m2 s−1 from the following Stokes-Einstein equation [20]:
(3)DG=kT6πμr,
where k is Boltzmann constant, and T is the absolute temperature.


## 2.2.2. Arterial Wall Layers

The transmural flow across the arterial wall can be described by the Brinkman equation as follows [16]:
(4)ρεl∂ul∂t+μlKlul=-∇pl+μlεl∇2ul,(5)∇·ul=0,
where ul and pl represent, respectively, the superficial velocity vector and the pressure based on the volume averaged method. εl and Kl are the porosity and the hydraulic permeability of the wall layer concerned. For all layers, the viscosity of plasma μl was assumed to be 0.72 × 10−3 kg m−1 s−1 [16].The transport of Gd-DTPA across the arterial layers was modeled by the following equation [16]:
(6)∂cl∂t+1-σlul·∇cl-DlΔcl=0,
where cl is the superficial concentration of Gd-DTPA; Dl is the effective diffusivity of Gd-DTPA in the arterial layers; σl is the filtration reflection coefficient.


## 2.3. Parameters

In order to solve (4)~(6), we have to acquire the values of the 4 parameters in the equations, namely, the porosity (εl), the hydraulic permeability (Kl), the effective diffusivity of Gd-DTPA (Dl) and the filtration reflection coefficient (σl) for each layer of the arterial wall. However, because it is hard to measure the parameters, the transport parameters are scarce. In the present study, most of the parameters were obtained from theoretical models. The calculation of the parameters is presented in Appendix. The parameters used are summarized in Table 1.Table 1
Model parameters for each layer.
EGL
Endothelium
Intima
IEL
Media
Hydraulic permeability (K, (m2))
6.0383 × 10−18
1.7383 × 10−20
4.2 × 10−17
8.4 × 10−20
6.09 × 10−19
Effective diffusivity (De, (m2 s−1))
1.0128 × 10−10
3.0276 × 10−13
1.2219 × 10−10
2.4094 × 10−13
9.3392 × 10−14
Reflection coefficient (σ)
0.0555
0.1212
0.2514
0.2514
0.3617
Porosity (ε)
0.6735
0.0005
0.8025
0.002
0.258


## 2.4. Boundary Conditions

As shown in Figure1(a), flow transport equations ((1), (4)~(5)) are subject to the following boundary conditions.BC-A at the inlet of the lumen of the arterial segment, a time-dependent fully developed (parabolic) inlet flow velocity profile was used for the pulsatile flow simulation (Figure 1(c)) [21]. The mean velocity was chosen to be 0.24 m s−1 [22].BC-B the pressure at the outlet boundary of the artery lumen was set at 100 mm Hg.BC-C at the media-adventitia interface, a constant pressure boundary condition with 30 mm Hg was employed.BC-D symmetric conditions were set at the axis of symmetry.BC-E no viscous flow was set on both the axial ends of the arterial wall.The boundary conditions for the mass transport equations ((2) and (6)) are as follows.BC-1: the injection of the contrast agent at the inlet was described by a gamma-variate function [12]:
(7)ct,z=0=at-t0be-gt-t0,t>t00,t≤t0,wherea=1.013×10-3, b=2.142, and g=0.454 s−1. t0 was set to 3 s, which meant that the bolus of contrast was injected after the periodicity of the flow field was completely developed.BC-2: for other boundaries, zero concentration gradient was assumed.


## 2.5. Computation Procedures

The numerical simulations were carried out using a validated finite element algorithm COMSOL Multiphysics. Three full cardiac cycles (3 s) simulation of the pulsatile flow with a time step of 4 ms were carried out to achieve a periodic flow independent of the initialization. Based on the initial velocity field, the mass transport equations were solved coupling with flow transport equations for 37 s with time step from 1 ms to 4 ms depending on the temporal variations of contrast concentration.


## 2.6. Quantification of Dispersion

The effect of the dispersion of the contrast agent bolus on the AIF can be described by convolving the estimated AIF with a vascular transport function (VTF) [4, 23]:
(8)AIFt-t0=AIFtest⊗VTFt-t0,z,
where t0 is the bolus delay and z is distance to the inlet. Using the variance of the VTF around its mean, the degree of dispersion can be quantified:
(9)σVTF2z=∫0∞t-MVTTz2·VTFt,zdt∫0∞VTFt,zdt,
where the mean vascular transit time (MVTT) is given by the ratio of the first to the zeroth moment of the VTF and can be used to quantify the bolus delay:
(10)MVTTz=∫0∞t·VTFt,zdt∫0∞VTFt,zdt.
AIFest is prescribed directly from the contrast agent concentration at the inlet ((7)  c(t,z=0)) and AIF is computed as the area weighted average of the concentration of contrast agent on several lumen cross-sections between the inlet and the outlet perpendicular to the axial vessel direction (c(t,z)). To obtain the mean delay (MVTT) and the variance (σVTF2), it is not necessary to calculate VTF firstly, since they can be calculated directly with the zeroth, first, and second integral moment of the concentration of the contrast agent [4]:
(11)σVTF2z=c2zc0z-c20c00+c10c002-c1zc0z2,MVTTz=c1zc0z-c10c00.
The moments of the concentration of the contrast agent are
(12)cnz=∫0∞tn·ct,zdt.


## 3. Results

### 3.1. Effects of the Arterial Transmural Transport on Contrast Agent Transport and the Bolus Dispersion

To demonstrate the deformation of the concentration-time curve along the flow direction in the unstenosed model, the area weighted average of the concentration of contrast agent on cross-sections at axial direction of10D and 20D(c(t,z)) was calculated. As shown in the Figures 2 and 3, the concentration-time curve illustrates the characteristics of bolus delay with phase shift and bolus dispersion with relatively low maximal concentration and spread width profile. For instance, in the lumen of the artery, the maximal concentration at 20D shifts about 0.5 second and decreases 2%, when compared with that at the inlet (Figure 2). The delay and dispersion are enlarged in the arterial wall, resulting in that the maximal concentration at 20D in endothelial layer is about 1.3 second later and 18% lower than that at the 10D (Figure 3). In addition to the delay and dispersion in the axial direction, as demonstrated in Figure 2, the concentration of the contrast in the radial direction is much more decreased and spread. For instance, the maximal concentration of intima layer is approximately 20 times smaller than that of the endothelial layer and the time interval between the half maximal and maximal concentration at 20D at the intima layer is approximately 4 times longer than that at the endothelial layer (Figure 3).Figure 2
The concentration-time curves at inlet,10D, and 20D in the model with arterial transmural transport.The concentration-time curves across the arterial wall at the axial direction of10D and 20D in the model with arterial transmural transport.
(a)
(b)
(c)
(d)The development of the bolus delay and dispersion in the unstenosed model are quantified by the mean vascular transit time (MVTT) and the variance of the VTF from (11). As illustrated in Figure 4(a), MVTT and the variance are increased almost linearly in relation to the distance to the inlet, which reflects consistently the concentration-time curve in Figure 2. When compared to the simulation with contrast agent transport only in the lumen, the arterial transmural transport would lead to slight increase in the MVTT and small increase in the variance (Figures 4(a)-4(b)) resulting in only increase of 0.186 s2 at 10D and 0.288 s2 at 20D. The quantitative difference in variance in Figure 4(c) further demonstrates that although the increase in the variance would increase gradually in relation to the distance to the inlet, the increase ratio would decrease sharply resulting in approximately 25% at 20D.Comparison of the quantitative parameters of the bolus delay (MVTT) and dispersion (variance) in the unstenosed artery between the model with arterial transmural transport and without. (a) The mean vascular transit time (MVTT); (b) the variance of vascular transport function (variance); (c) the difference and percentage difference in variance.
(a)
(b)
(c)


### 3.2. Effects of Endothelium on Contrast Agent Transport and the Bolus Dispersion

When the endothelium is damaged in diseased condition, the transport of the contrast agent may be affected and hence the bolus dispersion. To investigate the role of endothelium in the transport of the contrast agent, the endothelium was assumed to be totally damaged and the four transport parameters of the damaged endothelium and EGL were simplified to be similar to that of the intima. In addition, the transport parameters of other layers were consistent with that of the model with intact endothelium.Figures5(a)–5(d) demonstrate the concentration-time curve of the damaged endothelial model in the four arterial layers. As evident from these figures, due to the damage of the endothelium, the contrast agent is sharply accumulated in the arterial wall.Effects of endothelium on contrast agent transport and the bolus dispersion in the unstenosed model. (a–d) Comparison of the concentration-time curve across the arterial wall between the damaged endothelium model and the intact one at the axial direction of10D. (e-f) Comparison of the quantitative parameters of the bolus delay (MVTT) and dispersion (variance) between the damaged endothelium model and the intact one.
(a)
(b)
(c)
(d)
(e)
(f)As illustrated in Figures5(e)-5(f), the difference in the MVTT between the damaged endothelial model and intact one is slight and the difference in the variance between the two models decreases along the axial direction resulting in almost no difference at the outlet. These results indicate that the damage to the integrity of endothelium would have a small effect on the bolus delay and dispersion.


### 3.3. Effects of Disturbed Flow on Contrast Agent Transport and the Bolus Dispersion

Until now, the numerical simulations have been carried out only for a simplified straight axisymmetric blood vessel. However, the geometry of the physiological blood vessel is much more complex with such characteristics as branching, twisting, taper, and curvature, which would lead to far more than parabolic flow profile in the simplified model but very complicated flow patterns [24, 25]. One common flow pattern developed in the arterial system is the disturbed flow, which is inclined to occur in locations such as aneurysm, stenosed artery, the inner wall of curved segments, and the outer walls of arterial bifurcations [26, 27]. To further investigate the effects of arterial transmural transport on the bolus delay and dispersion under disturbed flow, contrast agent transport in the stenosed model was calculated.Figure6(a) shows the concentration distribution in the stenosed region at different moments of the contrast injection. It can be observed that the concentration in the disturbed flow region after the stenosis is quite uneven and significantly later than that of the regions of high flow.Effects of disturbed flow on contrast agent transport. (a) Concentration distribution in the stenosed region at three time points of the contrast injection at the inlet.t1 is the moment with half maximal concentration in the increasing phase; t2 is the moment with maximal concentration; t3 is the moment with half maximal concentration in the decreasing phase. (b–e) Comparison of the concentration-time curve across the arterial wall between the unstenosed model and the stenosed one at the axial direction of 10D.
(a)
(b)
(c)
(d)
(e)Figures6(b)–6(e) demonstrate the concentration-time curve in the four arterial layers at the disturbed flow region (10D). The concentration in the disturbed flow model decreases and spreads much more significantly than that in the straight model at the same axial location.In the stenosed model the influence of the disturbed flow on the bolus delay and dispersion was observed immediately behind the stenosis (Figure7). After the initial significant increase in the MVTT and variance behind the stenosis a reduction in the MVTT and variance follows in the disturbed flow zone. The arterial transmural transport would also affect the bolus delay and dispersion by decreasing the MVTT and increasing the variance. However, the effect of the arterial transmural transport on the MVTT and the variance is much less than that of the disturbed flow caused by the stenosis.Comparison of the quantitative parameters of the bolus delay (MVTT) and dispersion (variance) in the stenosed artery between the model with arterial transmural transport and without. (a) The mean vascular transit time (MVTT); (b) the variance of vascular transport function (variance).
(a)
(b)


## 3.1. Effects of the Arterial Transmural Transport on Contrast Agent Transport and the Bolus Dispersion

To demonstrate the deformation of the concentration-time curve along the flow direction in the unstenosed model, the area weighted average of the concentration of contrast agent on cross-sections at axial direction of10D and 20D(c(t,z)) was calculated. As shown in the Figures 2 and 3, the concentration-time curve illustrates the characteristics of bolus delay with phase shift and bolus dispersion with relatively low maximal concentration and spread width profile. For instance, in the lumen of the artery, the maximal concentration at 20D shifts about 0.5 second and decreases 2%, when compared with that at the inlet (Figure 2). The delay and dispersion are enlarged in the arterial wall, resulting in that the maximal concentration at 20D in endothelial layer is about 1.3 second later and 18% lower than that at the 10D (Figure 3). In addition to the delay and dispersion in the axial direction, as demonstrated in Figure 2, the concentration of the contrast in the radial direction is much more decreased and spread. For instance, the maximal concentration of intima layer is approximately 20 times smaller than that of the endothelial layer and the time interval between the half maximal and maximal concentration at 20D at the intima layer is approximately 4 times longer than that at the endothelial layer (Figure 3).Figure 2
The concentration-time curves at inlet,10D, and 20D in the model with arterial transmural transport.The concentration-time curves across the arterial wall at the axial direction of10D and 20D in the model with arterial transmural transport.
(a)
(b)
(c)
(d)The development of the bolus delay and dispersion in the unstenosed model are quantified by the mean vascular transit time (MVTT) and the variance of the VTF from (11). As illustrated in Figure 4(a), MVTT and the variance are increased almost linearly in relation to the distance to the inlet, which reflects consistently the concentration-time curve in Figure 2. When compared to the simulation with contrast agent transport only in the lumen, the arterial transmural transport would lead to slight increase in the MVTT and small increase in the variance (Figures 4(a)-4(b)) resulting in only increase of 0.186 s2 at 10D and 0.288 s2 at 20D. The quantitative difference in variance in Figure 4(c) further demonstrates that although the increase in the variance would increase gradually in relation to the distance to the inlet, the increase ratio would decrease sharply resulting in approximately 25% at 20D.Comparison of the quantitative parameters of the bolus delay (MVTT) and dispersion (variance) in the unstenosed artery between the model with arterial transmural transport and without. (a) The mean vascular transit time (MVTT); (b) the variance of vascular transport function (variance); (c) the difference and percentage difference in variance.
(a)
(b)
(c)


## 3.2. Effects of Endothelium on Contrast Agent Transport and the Bolus Dispersion

When the endothelium is damaged in diseased condition, the transport of the contrast agent may be affected and hence the bolus dispersion. To investigate the role of endothelium in the transport of the contrast agent, the endothelium was assumed to be totally damaged and the four transport parameters of the damaged endothelium and EGL were simplified to be similar to that of the intima. In addition, the transport parameters of other layers were consistent with that of the model with intact endothelium.Figures5(a)–5(d) demonstrate the concentration-time curve of the damaged endothelial model in the four arterial layers. As evident from these figures, due to the damage of the endothelium, the contrast agent is sharply accumulated in the arterial wall.Effects of endothelium on contrast agent transport and the bolus dispersion in the unstenosed model. (a–d) Comparison of the concentration-time curve across the arterial wall between the damaged endothelium model and the intact one at the axial direction of10D. (e-f) Comparison of the quantitative parameters of the bolus delay (MVTT) and dispersion (variance) between the damaged endothelium model and the intact one.
(a)
(b)
(c)
(d)
(e)
(f)As illustrated in Figures5(e)-5(f), the difference in the MVTT between the damaged endothelial model and intact one is slight and the difference in the variance between the two models decreases along the axial direction resulting in almost no difference at the outlet. These results indicate that the damage to the integrity of endothelium would have a small effect on the bolus delay and dispersion.


## 3.3. Effects of Disturbed Flow on Contrast Agent Transport and the Bolus Dispersion

Until now, the numerical simulations have been carried out only for a simplified straight axisymmetric blood vessel. However, the geometry of the physiological blood vessel is much more complex with such characteristics as branching, twisting, taper, and curvature, which would lead to far more than parabolic flow profile in the simplified model but very complicated flow patterns [24, 25]. One common flow pattern developed in the arterial system is the disturbed flow, which is inclined to occur in locations such as aneurysm, stenosed artery, the inner wall of curved segments, and the outer walls of arterial bifurcations [26, 27]. To further investigate the effects of arterial transmural transport on the bolus delay and dispersion under disturbed flow, contrast agent transport in the stenosed model was calculated.Figure6(a) shows the concentration distribution in the stenosed region at different moments of the contrast injection. It can be observed that the concentration in the disturbed flow region after the stenosis is quite uneven and significantly later than that of the regions of high flow.Effects of disturbed flow on contrast agent transport. (a) Concentration distribution in the stenosed region at three time points of the contrast injection at the inlet.t1 is the moment with half maximal concentration in the increasing phase; t2 is the moment with maximal concentration; t3 is the moment with half maximal concentration in the decreasing phase. (b–e) Comparison of the concentration-time curve across the arterial wall between the unstenosed model and the stenosed one at the axial direction of 10D.
(a)
(b)
(c)
(d)
(e)Figures6(b)–6(e) demonstrate the concentration-time curve in the four arterial layers at the disturbed flow region (10D). The concentration in the disturbed flow model decreases and spreads much more significantly than that in the straight model at the same axial location.In the stenosed model the influence of the disturbed flow on the bolus delay and dispersion was observed immediately behind the stenosis (Figure7). After the initial significant increase in the MVTT and variance behind the stenosis a reduction in the MVTT and variance follows in the disturbed flow zone. The arterial transmural transport would also affect the bolus delay and dispersion by decreasing the MVTT and increasing the variance. However, the effect of the arterial transmural transport on the MVTT and the variance is much less than that of the disturbed flow caused by the stenosis.Comparison of the quantitative parameters of the bolus delay (MVTT) and dispersion (variance) in the stenosed artery between the model with arterial transmural transport and without. (a) The mean vascular transit time (MVTT); (b) the variance of vascular transport function (variance).
(a)
(b)


## 4. Discussion

One unreal assumption of DSC-MRI is that the arteries are impermeable to contrast agent. In the present study, we numerically coupled contrast agent transport in the arterial wall with that in arterial lumen to investigate whether the unreal assumption is valid for the estimation of the delay and dispersion of the contrast agent bolus. The results obtained reveal that the arterial transmural transport would slightly affect the bolus delay qualified by MVTT and moderately increase the bolus delay qualified by the variance. The MVTT and the variance induced by the arterial transmural transport are much less than that by the disturbed flow after the stenosis. In addition, although the integrity of endothelium can significantly affect the accumulation of contrast agent in the arterial wall, it has small effects on the bolus delay and dispersion.The bolus dispersion would lead to a systematic blood flow underestimation. Graafen et al. simulated the dispersion in coronary arteries using a computational fluid dynamics approach and demonstrated that the variance between 1.0 and 2.5 s2 would lead to the myocardial blood flow underestimation between about 6% and 10% [11]. The present study demonstrated that the variance in both the stenosed and the unstenosed model is less than 2.5 s2 and the variance caused by the arterial transmural transport is less than 0.5 s2, which indicated that the myocardial blood flow error induced by the arterial transmural transport may be well below 10%.As the transport of contrast agent is mainly governed by the convective transport in the lumen and the diffusion transport in lumen and the arterial wall, the local flow pattern, the flow rate, and the arterial transmural transport would affect the degree of bolus dispersion. Among the three factors, the arterial transmural transport may be the least important one, since the present simulation indicated that the disturbed flow produced much more dispersion of the bolus than the arterial transmural transport. However, the contribution of local flow pattern and the flow rate to the dispersion is comparable. Our simulations demonstrated the disturbed flow induced by the stenosis would significantly increase the bolus dispersion. In contrast, Graafen et al. indicated that stenosis in the coronary arterial model leads to a reduction of dispersion [11]. Calamante et al. found that for the two patients with similar degree carotid stenosis, only one of them had larger dispersion, while the other was found to produce a dispersion of the bolus similar to that found in the normal subjects [10]. These studies and our results indicated that, to achieve an accurate estimation of the dispersion, the simulations with patient-specific geometrical and boundary layer model should be performed and the effects of the arterial transmural transport may be neglected in practice.This study revealed that local flow pattern caused by the stenosis would significantly affect the local contrast agent transport and hence the local bolus dispersion, which is constant with mass transport simulations based on the patient-specific model. For instance, Liu et al. demonstrated that the disturbed flow would significantly hinder the mass transport in the human aorta and computational studies indicated that the oxygen transport was significantly low in the outer wall of carotid artery where disturbed flow developed [28, 29]. Due to the fact that flow patterns are quite complex in the arterial segments such as the aneurysm, stenosed artery, the inner wall of curved segments, and the outer walls of arterial bifurcations [26, 27], the present results combining previous studies indicated that the arterial input function may be estimated beyond these regions to avoid large bolus dispersion [12].In this pilot study, the accuracy of the results may be reduced by the simplifications of the transport parameters, the blood rheological properties, and the geometries.Due to the scarcity of the transport parameters, most of the parameters were obtained from theoretical models. The estimated parameters, especially the diffusivity of the contrast agent, would affect the simulation results. The diffusivity used in the present study (6.5847 × 10−11 m2 s−1) is lower than the estimated diffusion coefficient of contrast agent in the previous studies (5.5 × 10−10 m2 s−1 and 1.5 × 10−10m2 s−1) [11, 12]. The lower diffusion coefficient in the lumen would lead to an increase of the dispersion due to a smaller exchange of contrast agent particles perpendicular to the flow direction. However, the lower diffusion coefficient in the arterial wall reduces the arterial transmural transport and hence decreases the dispersion in the lumen, which is shown in the present result that when compared to the damaged endothelial model the normal one with relatively low diffusivity leads to a small dispersion.The blood in the present study was simplified as Newtonian fluid and the non-Newtonian rheological properties such as shear thinning nature of the blood were neglected. It is reported that the shear thinning non-Newtonian nature of blood could slightly reduce oxygen flux (similar micromolecule as contrast agent) in most regions of the arteries, and this effect became much more apparent in areas with disturbed flow [29]. Therefore, it may be valid to assume the blood Newtonian fluid to estimate the bolus dispersion in most regions of the model However, the assumption of Newtonian fluid may lead to an underestimation of dispersion in the disturbed flow region.Another limitation of the present study is that all simulations were performed on simplified geometries. The geometry of the physiological blood vessel is very complex with such characteristics as branching, twisting, taper, and curvature, which would lead to complicated flow patterns. Therefore, it is necessary to use the realistic geometries to simulate the exact dispersion of the contrast agent bolus. As the main aim of the present study is not to estimate the exact dispersion errors in an individual patient, the results obtained from the typical parabolic flow and disturbed flow can still shed some light on the effects of arterial transmural transport of the contrast agent on the bolus delay and dispersion.


## 5. Conclusion

The arterial transmural transport would slightly affect the bolus delay and small increase the bolus dispersion, which may have small effects on the estimation of the blood flow estimation. In comparison with disturbed flow induced by the presence of stenosis, the arterial transmural transport plays a less important role in the bolus delay. Although the assumption of impermeable arterial wall of DSC-MRI is unreal, it may still be a good simplification for the estimation of the delay and dispersion of the contrast agent bolus.



---

*Source: 803276-2014-11-17.xml*