Effects of size and elasticity on the relation between flow velocity and wall shear stress in side-wall aneurysms: A lattice Boltzmann-based computer simulation study
Authors:
Haifeng Wang aff001; Timm Krüger aff002; Fathollah Varnik aff001
Authors place of work:
Interdisciplinary Center for Advanced Materials Simulation (ICAMS), Ruhr-Universität Bochum, Bochum, Germany
aff001; School of Engineering, Institute for Multiscale Thermofluids, University of Edinburgh, Edinburgh, Scotland, United Kingdom
aff002
Published in the journal:
PLoS ONE 15(1)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0227770
Summary
Blood flow in an artery is a fluid-structure interaction problem. It is widely accepted that aneurysm formation, enlargement and failure are associated with wall shear stress (WSS) which is exerted by flowing blood on the aneurysmal wall. To date, the combined effect of aneurysm size and wall elasticity on intra-aneurysm (IA) flow characteristics, particularly in the case of side-wall aneurysms, is poorly understood. Here we propose a model of three-dimensional viscous flow in a compliant artery containing an aneurysm by employing the immersed boundary-lattice Boltzmann-finite element method. This model allows to adequately account for the elastic deformation of both the blood vessel and aneurysm walls. Using this model, we perform a detailed investigation of the flow through aneurysm under different conditions with a focus on the parameters which may influence the wall shear stress. Most importantly, it is shown in this work that the use of flow velocity as a proxy for wall shear stress is well justified only in those sections of the vessel which are close to the ideal cylindrical geometry. Within the aneurysm domain, however, the correlation between wall shear stress and flow velocity is largely lost due to the complexity of the geometry and the resulting flow pattern. Moreover, the correlations weaken further with the phase shift between flow velocity and transmural pressure. These findings have important implications for medical applications since wall shear stress is believed to play a crucial role in aneurysm rupture.
Keywords:
blood pressure – Arteries – Aneurysms – Deformation – Shear stresses – Fluid flow – Fluids – Mechanical energy
Introduction
Brain aneurysms lead to almost 500,000 deaths per year [1]. Physiological flows, such as arterial blood flow, are generally characterized by the transport of fluids in compliant tubes [2], which involves complex fluid-structure interaction problems. The interplay of hemodynamics and elastic arterial vessel walls (e.g. through wall shear stress (WSS) and blood pressure) is believed to play a central role in the aneurysm initiation, growth and rupture [3–7].
It is known that abnormal WSS drives degradation of the vascular wall [8–10]. Flow velocity is one of the principal factors determining the magnitude of WSS. Compared with WSS, in practice, flow velocity is easier to be measured and often acts as a proxy observable for WSS [11–13]. In this work, we study how size, elastic deformation of the aneurysm and the waveform of transmural pressure determine the connection between flow velocity and wall shear stress as an important hemodynamic factor.
Aneurysms can be classified into two major categories: fusiform and saccular aneurysms; the latter is the most common type of aneurysm and the type most prone to rupture [14]. The saccular aneurysm has two subtypes: end-wall and side-wall aneurysms [15].
Intra-aneurysm (IA) hemodynamics is sensitive to morphological factors such as shape and size of the aneurysm [16–19]. It has been reported that an increase of aneurysm size leads to a decrease of the average IA flow velocity, U¯, and wall shear stress, σw¯, in the case of side-wall-type aneurysms [20, 21]. In these studies, the aneurysm wall is assumed to be rigid [20] or modeled as a stiff hyperelastic wall [21]. To date, however, the effect of wall elasticity and the resulting deformation is poorly studied.
Another objective of this work, therefore, is to investigate the combined effect of aneurysm size and wall softness on the variations of the average IA flow velocity and WSS, centering on the side-wall-type aneurysm.
We model the interaction of hemodynamics and wall elasticity in an idealized aneurysm using the lattice Boltzmann method (LBM) for the fluid flow, the finite element method (FEM) for the thin-walled vessel dynamics and the immersed boundary method (IBM) for the fluid-structure interaction.
The paper is organized as follows. The physical model is introduced in Sect Physical model. The numerical methods are presented in Sect Numerical methods, and benchmark tests are shown in Sect Validation. Simulation results using a curved artery with a side-wall aneurysm are presented and discussed in Sect Results and discussion, with a particular focus on the effect of aneurysm softness on IA hemodynamics. Our work is concluded in Sect Conclusion.
Physical model
Physical ingredients
Aneurysms occur in different sizes and shapes, and their detailed properties depend on the subject and progression state. According to [20, 21], the trends of the variations of flow velocity and wall shear stress due to the effect of aneurysm size are not specific to any particular side-wall aneurysm geometry. Therefore, after benchmark studies of flow through a straight cylindrical channel (Fig 1a), we use a representative model of a curved artery segment including a simple side-wall aneurysmal dome (Fig 1b and 1c). Without losing generality, this approach is similar to the simplification used in studies of the end-wall aneurysms [22–24].
Vascular tissues are generally heterogeneous. Healthy arterial walls are normally strain-hardening due to the existence of collagen fiber constituents; aneurysms, however, often exhibit strain-softening behavior [27–30]. The arterial wall in our model is, therefore, divided into different regions (Fig 1b and 1c) with different constitutive laws. Similar to Charalambous et al. [31], we employ the strain-hardening Skalak (SK) model [26] for the healthy arterial wall:
The vessel wall is represented as a two-dimensional surface with zero thickness. However, one important effect closely related to a finite vessel wall thickness is resistance to bending. In order to account for this property in our membrane model, we explicitly introduce a local bending energy with a bending constant κb, which can then give rise to forces normal to the wall (for further details, please refer to [32]).
The two free parameters κα and κb are related to the Young’s modulus E, wall thickness d and Poisson ratio ν according to [33]
Thus, given the two material coefficients, E and ν and the thickness of the tissue, d, it is possible to determine the parameters of our model such that it reproduces the properties of a real vessel. This is explicitly shown below via a benchmark of the present study versus a 3D model. Obviously, the inverse mapping from a two-parameter model to a three-parameter one is not unique. This is not necessarily a drawback. Rather, it enhances the generic character of our model and allows the results obtained within these simulations to be transferred to vessels with different materials properties and wall thickness.
The transmural pressure Ptr is the pressure difference between the mean arterial pressure and the intracranial pressure. It is generally believed that, given heterogeneity of elastic properties along the arterial wall, Ptr is mainly responsible for the wall deformation and plays a crucial role in the progress of aneurysm growth as well as failure. Therefore, the transmural pressure needs to be taken into account in our model.
Although blood is non-Newtonian, it is often treated as an incompressible Newtonian fluid (e.g., [20, 21, 34–37]). Since the focus of this work is on the basic understanding of the effect of aneurysm size and softness rather than making accurate predictions, we model blood as a homogeneous Newtonian fluid with constant dynamic viscosity η.
Finally, we can assume the no-slip velocity boundary condition at the wall.
Hemodynamic observables
Wall shear stress (WSS) is a crucial factor for the development of aneurysms. The fluid stress tensor includes pressure (isotropic stress) and viscous contributions (deviatoric stress):
For the aneurysm dome Ω (Fig 1b and 1c) with corresponding surface area A and volume V, we define the spatial WSS average
Dimensionless groups
Based on the known material and flow properties (Young’s modulus E, artery radius R, transmural pressure Ptr, arterial flow velocity U and fluid viscosity η), we can define a series of dimensionless groups characterizing the deformation of the arterial wall. The effect of transmural pressure can be described by
Numerical methods
We employ a computational scheme combining the lattice Boltzmann method (LBM), finite element method (FEM) and immersed boundary method (IBM) [32] for a study of the coupled hemodynamic-artery problem outlined in Sect. Physical model. This hybrid method has previously been applied to suspensions of deformable red blood cells [40–42].
Elastic wall model and force computation
The artery and aneurysm are modeled as thin elastic walls represented by a set of vertices and flat triangular facets. There are four force contributions acting on each vertex: the in-plane elastic shear force Fs, the normal bending force Fb, the transmural-pressure force Ftp and a tether-like spring force Fsp.
The surface elastic shear energy Ws = ∫dAωs corresponding to the energy densities in Eq (1) or Eq (2) is computed via
The present model includes a bending resistance to avoid buckling of the arterial wall. The total bending energy of the artery is numerically approximated by
Given the surface elastic shear energy Ws and the bending energy Wb, the in-plane force Fsi and bending force Fbi acting on vertex i at position xi can be computed via the principle of virtual work:
The transmural pressure Ptr is known to be one of the main factors dominating the enlargement and rupture of aneurysms [44–46]. There are basically two ways of modeling the transmural pressure. The most obvious one is to treat the transmural pressure as different pressures on either side of the arterial wall. This approach, however, poses a major challenge when combined with the immersed boundary method. We follow an alternative approach, introducing a transmural-pressure-based force Ftp acting on each facet i along its normal direction:
The fourth force contribution is a spring force Fsp which is used to tether the ends of the artery:
Flow solver and fluid-structure interaction
In order to solve the incompressible Navier-Stokes equations, we employ a standard LBM on a D3Q19 lattice with the multiple-relaxation-time (MRT) collision operator for improved stability and accuracy [47–49]. Forces are included according to the Guo’s forcing scheme [50]. The viscous stress tensor is computed from the non-equilibrium populations according to [51].
Peskin’s IBM is used to couple the arterial wall and fluid dynamics [52]. Elastic forces calculated for each vertex i are spread to surrounding lattice nodes where they are used as input for the LBM. After the updated fluid velocity field has been obtained, fluid velocities are interpolated at the vertex positions via IBM. A tri-linear stencil with 23 lattice nodes for each vertex is used during spreading and interpolation. The forward-Euler scheme is used to update all vertex positions. Further details can be found in [32].
Interaction-related elastic energy
In addition to parameters such as maximum pressure difference along the blood vessel, which drives the flow, and size as well as flexibility of the aneurysm, the dynamics of the aneurysm wall may also strongly influence the IA flow field and thus lead to significant variations of the wall shear stress. In order to address this aspect, we will monitor the rate of variation of the elastic energy and investigate its correlations with WSS. The rate of elastic energy variation due to fluid-wall interactions can be calculated as
Validation
To validate the present model, we perform two different types of simulations: (i) steady and unsteady flow in a rigid tube and (ii) steady flow in an elastic tube.
Steady and unsteady flows in rigid tube
We consider three-dimensional flow in a stiff tube with circular cross-section (Fig 1a). With the analytical solutions known for steady Poiseuille and pulsatile Womersley flows, this test allows a quantitative analysis of our model. All relevant simulation parameters are given in Table 1.
A body force F along the axis of the tube is applied on all fluid nodes inside the tube to drive the flow. For Poiseuille flow, F is equivalent to the constant pressure gradient −p′ along the axis. For Womersley flow, the body force is given by F(t) = −p′ cos(ωt) where ω = 2π/T is the angular frequency and T is the sampling period.
The analytical Poiseuille solution is
The simulated flow velocities and shear stresses agree well with the corresponding analytical predictions in the two flow tests (Figs 2 and 3). However, an extrapolation-based correction scheme is necessary to reach good agreement near the wall. The velocity and stress in direct wall vicinity are affected by the diffuse nature of the immersed boundary method. The basic idea behind the correction scheme is to enforce the no-slip condition in the diffuse interface IBM [54]. Specifically, fluid velocity and stress at the fluid nodes near the wall (i.e., at |r/R*| = 11/11.035 in Figs 2 and 3) are corrected using a second-order Lagrange extrapolation from corresponding information of their extrapolation points (i.e., at |r/R*| = 10/11.035, 9/11.035 and 8/11.035). For general geometries, the extrapolation points may be off lattice lines. In such cases, fluid quantities (denoted by f) on each extrapolation node are computed by a linear interpolation of information at neighboring lattice nodes as f = (∑i fi/ℓi)/(∑i 1/ℓi) where ℓi indicates the distance between the extrapolation node and its neighboring lattice node i. For the sake of stability, when the extrapolation point is very close to a lattice node (e.g., ℓ < 0.0001), the value at that node is taken (f = fi). With this correction approach, the maximum relative L2 error of the shear stress σ is reduced from 31.2% to 2.2% in the Poiseuille flow and from 32.7% to 3.5% in the Womersley flow.
Note that the wall location as recovered from the simulation is not exactly where it is expected to be. Due to the diffuse nature of the IBM, the location of zero velocity (as obtained from a quadratic extrapolation from the interior solution) is shifted with respect to the location of the Lagrangian markers. We treated the tube radius as fit parameter to obtain the best agreement with the simulated velocity profile. For our simulation parameters the apparent tube radius is R* = 11.035Δx compared to the input value of R = 11.459Δx. Similar deviations have been reported in the literature [55–57] and are acceptably small.
Steady flow in elastic tube
In order to show that the numerical model can faithfully reproduce flow in flexible three-dimensional geometries, we simulate steady flow in an elastic tube with a circular cross section (Fig 1a). In contrast to the steady Poiseuille flow, we consider the transmural pressure and the radial deformation of the tube walls. We use the well-developed software SimVascular [58] to generate reference results. SimVascular employs a simple linear elastic model for the walls.
To compare results from two different software tools, we need to ensure that the material models and properties used are identical. SimVascular requires the Young’s modulus E, the wall thickness d, and the Poisson ratio ν as input parameters. In the present model, the input parameters are the surface elastic shear modulus κs, the area dilation modulus κα, and the bending modulus κb. Both sets of parameters are related according to Eqs (3) and (4). For the sake of stability of simulations with a certain resolution, small values of Young’s modulus E and transmural pressure Ptr are utilized. Nevertheless, it is Cp, rather than Ptr or E alone, that dominates the wall deformation. Typical values used in SimVascular are E = 0.05 MPa, d = 0.15 mm and ν = 0.45. To match these parameters we set κs = 0.025, κα = 0.025 and κb = 0.002 in simulation units. The other simulation parameters are listed in Table 1. The length-diameter aspect ratio of the tube is 7.5. In all simulations, the ends of the tube are tethered and thus the deformation is not homogeneous along the tube axis; the radial displacements are zero at the ends and maximum in the middle region where the diametric strains (ΔR/R) are measured.
The combined effect of transmural pressure and elasticity on wall deformation is expressed by the dimensionless parameter Cp, see Eq (10). We vary Cp by changing the transmural pressure Ptr in the present test case. Specifically, Ptr varies from 5.1 Pa to 153.2 Pa (i.e., 0.0001 ≤ Cp ≤ 0.003) in SimVascular simulations. Correspondingly, Ptr in our model ranges from 3.3 ⋅ 10−6 to 1 ⋅ 10−4 in simulation units. The radial strain ΔR/R of the wall is linearly proportional to Cp in SimVascular simulations due to the underlying linear elastic model. In order to investigate the strain-hardening behavior of the tube in the present model, five additional simulations with Ptr up to 6.5 ⋅ 10−4 in simulation units are carried out.
We find that simulation results using our model compare well with those obtained from SimVascular within the small strain regime (Fig 4a). In the large strain range, however, the strain-softening and strain-hardening features of the Neo-Hookean and Skalak models, respectively, become dominant and deviate from the behavior of the linear elastic model used in SimVascular.
Our simulation results also compare well with the so-called tube law (Fig 4b). In cases of wall expansion (A/A0 > 1), the tube law is [59, 60]
Results and discussion
In the following, the combined effect of aneurysm size and wall softness on the flow characteristics within the side-wall aneurysm dome is investigated, to our best knowledge, for the first time. Brief introductions to both characteristic parameters and simulation setup are given before results are reported.
Characteristic parameters
Blood and artery properties show variations among patients, but reference values based on a reasonable number of measurements have been reported. In general, the radius R of the cerebral artery ranges from 1.5 to 2.5 mm [21, 61, 62]. Young’s modulus E and Poisson’s ratio ν of the healthy cerebral artery are approximately 1–2 MPa and 0.45, respectively [21, 63]. The wall thickness is between 0.03 and 0.2 mm in most aneurysms, and is here assumed to be uniform and 0.15 mm, as adopted by [63]. There is no uniform ideal value for the transmural pressure Ptr. However, its range lies roughly between 50 and 70 mmHg (i.e., between 6.7 and 9.3 kPa). Blood density ρ and dynamic viscosity η are within the range of 1–1.06 g/cm3 and 3–4 mPa s, respectively [21, 63]. The flow velocity u in cerebral arteries ranges from 0.1 to 1.0 m/s, and those in the common carotid and middle cerebral arteries are approximately 0.4 and 0.6 m/s, respectively [64–67].
In order to better highlight the effects arising from the interplay between elasticity, complex geometry and the time-dependent forces acting on the wall and on the fluid, it is useful to express physical quantities in such a way that trivial effects such as quadratic dependence of the maximum flow velocity and the linear dependence of maximum shear stress on tube radius as well as their linear dependence on the applied pressure gradient are “divided out”. A way to achieve this goal is to divide the quantities of interest such as flow velocity, u, wall shear stress, σw, and viscous dissipation rate, ϕ˙, by their corresponding characteristic values in a straight cylindrical channel subject to a stationary Poiseuille flow, u*, σw* and ϕ˙*. Throughout the paper, we will use the following reference values to make the corresponding quantities dimensionless: u* = −p′R2/(4η) = 0.355 m/s, σw*=−p′R/2=1.42Pa, and ϕ˙*=(σw*)2/(2η)=252W/m3. In this estimate, we used the values of the mean pressure gradient, p′ = 14.2 Pa/cm (Fig 5), vessel radius, R = 2 mm and dynamic blood viscosity, η = 4 ⋅ 10−3 Pa s.
Simulation setup
We consider pulsatile flow in seven different cases:
Reference channel: straight cylindrical artery (Fig 1a) with elastic Young’s modulus E0 = 2 MPa. In the following, we will refer to this choice of E0 as the “Normal Elastic” (NE) case.
LA/stiff: large aneurysm (Fig 1b) with Ptr = 0
LA/NE: large aneurysm (Fig 1b) with normal elasticity E0
LA/RE: large aneurysm (Fig 1b) with reduced elasticity E0/3
SA/stiff: small aneurysm (Fig 1c) with Ptr = 0
SA/NE: small aneurysm (Fig 1c) with normal elasticity E0
SA/RE: small aneurysm (Fig 1c) with reduced elasticity E0/3
The initial (equilibrium) volume of the large aneurysm is approximately twice as large as that of the small aneurysm. The cases with zero transmural pressure mimic nearly rigid walls and allow to focus on the effects of size and deformability separately. All relevant parameter values used in the following simulations are given in Table 2.
Pulsatile flow is generated by applying a time-dependent body force (equivalent to the pressure gradient p′) at the inlet. The transmural pressure is applied on the entire wall. In reality, it is challenging to obtain both pressure gradient (or flow velocity) and transmural pressure waveforms from the same underlying mechanism (heart beat). Therefore, they are fed separately into the present model.
In practice, fluid velocity is often used as a proxy observable for wall shear stress [11–13]. While obvious in the case of a straight cylindrical channel, the complex shape of the aneurysm and its deformability introduce new aspects whose effects on the connection between flow velocity and wall shear stress needs a thorough investigation. Here, we address this issue via two different classes of spatially and temporally resolved simulations using the present hybrid model.
First we employ realistic waveforms (Fig 5) for both the pressure gradient and the transmural pressure. Given the chosen input parameters, the Reynolds, Womersley and Cp values inside the parent artery region with normal elasticity are approximately 350, 3 and 0.005, respectively. Within the aneurysm dome, the Reynolds number is almost 10 times smaller and Cp changes according to the variation of aneurysmal wall elasticity.
In the second class of simulations, to simplify the analysis, we assume sinusoidal waveforms for the pressure gradient, p′(t) = p′(0) sin(2πt/T), and the transmural pressure, Ptr(t) = Ptr(0) sin(θ + 2πt/T). We introduce a phase shift θ between both waveforms in order to quantify the individual contributions of both waveforms to the flow field.
Simulations are repeated for six cycles, which is sufficient to reach a time-periodic flow field. In all the simulations whose results are reported below, the open ends of the blood vessel are tethered. A body force (equivalent to pressure gradient p′ in Fig 5) is imposed at the inlet, and periodic boundary conditions are applied along the vessel axis.
Realistic waveforms: Effect of aneurysm softness
Fig 6 shows the time evolution of the space-averaged WSS and flow velocity in the straight artery (reference case). A comparison with the imposed waveforms in Fig 5 reveals that the time evolution of the hydrodynamic observables in the reference case basically follow the pressure gradient and transmural pressure. This is to be expected since the straight artery does not have complex geometrical features that could affect the flow. The average flow velocity and WSS are similar to those obtained for the corresponding steady Poiseuille flow.
The time evolution of the intra-aneurysm (IA) WSS and flow velocity for the large (LA) and small (SA) aneurysm cases are shown in Fig 7. Similar to the straight tube case, the key features of the temporal variations of WSS and flow velocity stem from the pulsatility of flow rather than the existence of the aneurysm (compare Figs 5 and 7). However, both the aneurysm size and wall softness impact the finer details of the IA hemodynamic characteristics as will be discussed in the following.
For the small-deformation cases (normal elasticity, LA/NE and SA/NE, and non-deformable, LA/stiff and SA/stiff), we find that the spatial average wall shear stress, σw¯, and flow velocity, U¯, decrease with increasing the aneurysm size. In these cases the IA flow can be approximated by flow in a rigid cavity. The larger aneurysm has more potential to result in flow regions of recirculation and stagnation which are generally accompanied by lower WSS. This is in line with reported results in [20] and [21].
In real aneurysms, arterial wall stiffness may decrease due to tissue degradation. For the softer aneurysm (cases SA/RE and LA/RE), we find that σw¯ and U¯ increase with aneurysm size (solid and dashed red lines in Fig 7). Note that the size of the initially larger aneurysm (LA) always remains larger than that of the smaller aneurysm. The variation of aneurysm volume (i.e., volumetric strain) directly follows the time evolution of the transmural pressure in Fig 5. The maximum volumetric strains are 61% and 86% in the cases SA/RE and LA/RE, respectively. Apparently, the aneurysm size no longer dominates the IA WSS and flow velocity behavior if the aneurysm is sufficiently soft. We can, therefore, conclude that there are at least two independent mechanisms controlling the IA flow properties: (i) size and (ii) deformability of the aneurysm dome.
Focusing on either the large or the small aneurysm, Fig 7 demonstrates that, the softer the aneurysm, the larger the IA WSS and the flow velocity, all averaged within the aneurysm domain. In order to gain more insight into the problem, a detailed view of the velocity field is shown in Fig 8 for three choices of a non-deformable (rigid) wall (indicated as stiff in the plot), a membrane with an intermediate (‘normal’) elastic constant, E0, and a soft aneurysm tissue with reduced elasticity. The plot shows that the flow velocity within the soft aneurysm reaches higher values than the two other cases. This observation may be rationalized as follows. Since the velocity of fluid at the wall is equal to that of the wall (stick boundary condition), a rigid and immobile vessel acts like a sink with respect to fluid momentum and thus decelerates the flow. A deformable tissue, on the other hand, can move under the action of transmural pressure and thus represents a less severe obstacle to the flow. In fact, depending on the phase of the flow with respect to the motion of the membrane, a partial enhancement of the fluid velocity is also possible. One possibility here would be that an elastic wall can store energy during one period of a cycle and release it during a different period again. This effect becomes stronger as the membrane softens further and its deformability increases.
Such correlation between wall softness and flow parameters may, however, be reduced or even violated to some extent when effects of input waveforms become significant, as will be discussed in the following section.
From the viewpoint of energy, the effect of wall softness can be indicated by the variation of elastic energy. The elastic energy follows the transmural pressure (Fig 9a) since the shear and inertial contributions are relatively negligible (at least two orders of magnitude smaller) compared with the transmural-pressure contribution. As expected, the elastic contribution disappears in the case of a stiff wall (zero transmural pressure) (Fig 9a).
In the case of side-wall aneurysms with realistic waveforms, we find a linear correlation between the rate of elastic energy variation and the wall shear stress (Fig 9b). No clear relation between the IA flow velocity and the rate of change of elastic energy has been found, though. Note, however, that it might be impossible to find a universal correlation between elastic energy rate and wall shear stress for all types of aneurysm geometries and input signals because the hemodynamics in aneurysms is sensitive to a number of conditions, including shape parameters such as inflow angle.
Naively, the IA WSS should follow the imposed pressure gradient. The deformation of the aneurysm, largely caused by the transmural pressure waveform, however, modifies the flow field in the aneurysm, which in turn affects the WSS. In order to investigate this effect, we consider sinusoidal input waveforms for the large aneurysm (Fig 1b) in the following.
Sinusoidal waveforms: Correlations between WSS and fluid velocity
Using sinusoidal waveforms for both pressure gradient p′ and transmural pressure Ptr, together with a phase shift θ for the latter, we consider the following cases:
θ = −π/2: cosine-type Ptr, i.e., p′ and Ptr are uncorrelated
θ = 0: sine-type Ptr, i.e., p′ and Ptr are positively correlated
θ = π: sine-type Ptr with inverse amplitude, i.e., p′ and Ptr are negatively correlated
We begin with an analysis of the flow in the parent vessel before moving on to the aneurysmal dome. Well away from the aneurysm, flow velocity and WSS are governed by the input pressure gradient, and the phase shift between this driving force and the transmural pressure does not play any significant role (Fig 10). The cross-sectional average values of the flow velocity and wall shear stress are linearly correlated in the parent vessel (Fig 11a). There also exists a linear correlation between maximum velocity and average WSS and between maximum velocity and maximum WSS (the corresponding data is very similar to Fig 11a and is thus not shown). This strong correlation is also reflected in a high Pearson correlation coefficient, which turns out to be close to unity outside the aneurysm domain for all the investigated cases: 0.98, 0.99, 0.99 for θ = 0, −π/2 and π, respectively. Other combinations of size and transmural pressure give similarly large correlation coefficients: 0.99 (LA/NE/θ = π), 0.96 (LA/stiff) and 0.98 (SA/stiff). For the fluid outside the aneurysm, therefore, one can safely use the flow velocity as a proxy for wall shear stress.
Within the aneurysmal domain, however, the flow has a more complex pattern so that a one-to-one correspondence between characteristic flow velocity and wall shear stress is no longer ensured. This is particularly the case for different values of the phase shift θ between pressure gradient and transmural pressure waveforms (Fig 12). From Fig 11b it is evident that the correlation between the velocity and WSS within the aneurysm is deteriorating with increasing phase shift angle θ. In terms of Pearson correlation coefficient, one obtains 0.91 (θ = 0), 0.76 (θ = −π/2), 0.67 (θ = π), 0.84 (LA/NE/θ = π), 0.72 (LA/stiff) and 0.93 (SA/stiff). One contributing factor to this deterioration is that the transmural pressure directly controls the aneurysm wall motion, which in return affects the IA flow and WSS.
Results shown in Fig 11b also reveal effects of the aneurysm size on the correlations between flow velocity and WSS. It can be seen from this data that the proxy role of flow velocity as an estimate of WSS is less justified in the case of larger and more flexible domes. This finding is of great practical importance as serious pathological cases are often associated with large aneurysm sizes and an advanced state of degradation, which leads to a weak elasticity and hence large deformability.
Conclusion
In this work, an immersed-boundary-lattice-Boltzmann-finite-element method is employed to simulate three-dimensional viscous flow in an elastic vessel containing a side-wall aneurysm, with a special focus on the relevance of the aneurysm deformability. Accuracy of the method is shown through steady Poiseuille and pulsatile Womersley flow tests in a circular rigid tube. Furthermore, the deformation behavior of an elastic pipe under the effect of transmural pressure is studied using the present model and validated by comparisons to both SimVascular and the tube law.
Hemodynamic quantities, such as flow velocity and wall shear stress (WSS), play a crucial role in the development and failure of aneurysms. Using the above described simulation tool, we investigate the combined effect of aneurysm size and wall softness on intra-aneurysm (IA) flow parameters. Under the rigid wall assumption, the IA flow velocity and WSS decrease with increasing aneurysm size. However, the situation becomes more complex when the aneurysm deformation becomes important.
The central result of our investigation concerns the connection between WSS and flow velocity, the latter often used as a proxy observable to estimate the former. While an indirect evaluation of the WSS via the flow velocity is a reliable method for regular geometries (planar or cylindrical), the present study reveals that the correlation between the two quantities deteriorates within the aneurysmal domain. The situation becomes more problematic in the case of large and easily deformable aneurysms.
In aneurysms, arterial wall stiffness can significantly decrease due to tissue degradation. This is often accompanied by an increase of the aneurysm’s size. Our results suggest that, in such serious pathological cases, fluid velocity is no longer a reliable proxy for WSS.
Supporting information
S1 Appendix [wmv]
Temporal variation of fluid velocity field within a planar cut of the aneurysm.
Zdroje
1. Bederson JB, Awad IA, Wiebers DO, Piepgras D, Haley EC Jr, Brott T, et al. Recommendations for the Management of Patients With Unruptured Intracranial Aneurysms: A Statement for Healthcare Professionals From the Stroke Council of the American Heart Association. Stroke. 2000;31:2742–2750. doi: 10.1161/01.str.31.11.2742 11062304
2. Grotberg JB, Jensen OE. Biofluid mechanics in flexible tubes. Annu Rev Fluid Mech. 2004;36:121–147. doi: 10.1146/annurev.fluid.36.050802.121918
3. Castro MA, Putman CM, Cebral JR. Patient-specific computational fluid dynamics modeling of anterior communicating artery aneurysms: a study of the sensitivity of intra-aneurysmal flow patterns to flow conditions in the carotid arteries. AJNR Am J Neuroradiol. 2006;27(10):2061–8. 17110667
4. Castro MA. Understanding the role of hemodynamics in the initiation, progression, rupture, and treatment outcome of cerebral aneurysm from medical image-based computational studies. ISRN Radiology. 2013;013:602707.
5. Furukawa K, Ishida F, Tsuji M, Miura Y, Kishimoto T, Shiba M, et al. Hemodynamic characteristics of hyperplastic remodeling lesions in cerebral aneurysms. PLoS ONE. 2018;13(1):e0191287. doi: 10.1371/journal.pone.0191287 29338059
6. Xu L, Liang F, Zhao B, Wan J, Liu H. Influence of aging-induced flow waveform variation on hemodynamics in aneurysms present at the internal carotid artery: A computational model-based study. Comput Biol Med. 2018;101:51–60. doi: 10.1016/j.compbiomed.2018.08.004 30099239
7. Diagbouga MR, Morel S, Bijlenga P, Kwak BR. Role of hemodynamics in initiation/growth of intracranial aneurysms. Eur J Clin Invest. 2018;48(9):e12992. doi: 10.1111/eci.12992 29962043
8. Humphrey JD, Taylor CA. Intracranial and abdominal aortic aneurysms: similarities, differences, and need for a new class of computational models. Annu Rev Biomed Eng. 2008;10:221–46. doi: 10.1146/annurev.bioeng.10.061807.160439 18647115
9. Kolega J, Gao L, Mandelbaum M, Mocco J, Siddiqui AH, Natarajan SK, et al. Cellular and molecular responses of the basilar terminus to hemodynamics during intracranial aneurysm initiation in a rabbit model. J Vasc Res. 2011;48(5):429–42. doi: 10.1159/000324840 21625176
10. Meng H, Tutino VM, Xiang J, Siddiqui A. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol. 2014;35(7):1254–62. doi: 10.3174/ajnr.A3558 23598838
11. Oshinski J, Ku D, Mukundan S Jr, Loth F, Pettigrew R. Determination of wall shear stress in the aorta with the use of MR phase velocity mapping. J Magn Reson Imaging. 1995;5:640–647. doi: 10.1002/jmri.1880050605 8748480
12. Carvalho J, Nielsen J, Nayak K. Feasibility of in vivo measurement of carotid wall shear rate using spiral fourier velocity encoded MRI. Magn Reson Med. 2010;63:1537–1547. doi: 10.1002/mrm.22325 20512857
13. Szajer J, Ho-Shon K. A comparison of 4D flow MRI-derived wall shear stress with computational fluid dynamics methods for intracranial aneurysms and carotid bifurcations—A review. Magn Reson Imaging. 2018;48:62–69. doi: 10.1016/j.mri.2017.12.005 29223732
14. Novitzke J. The basics of brain aneurysms: A guide for patients. J Vasc Interv Neurol. 2008;1(3):89–90. 22518230
15. Hassan T, Timofeev EV, Saito T, Shimizu H, Ezura M, Matsumoto Y, et al. A proposed parent vessel geometry-based categorization of saccular intracranial aneurysms: Computational flow dynamics analysis of the risk factors for lesion rupture. J Neurosurg. 2005;103:662–680. doi: 10.3171/jns.2005.103.4.0662 16266049
16. Ujiie H, Tachibana H, Hiramatsu O, Hazel AL, Matsumoto T, Ogasawara Y, et al. Effects of size and shape (aspect ratio) on the hemodynamics of saccular aneurysms: a possible index for surgical treatment of intracranial aneurysms. Neurosurgery. 1999;45(1):119–29. doi: 10.1097/00006123-199907000-00028 10414574
17. Baharoglu MI, Schirmer CM, Hoit DA, Gao BL, Malek AM. Aneurysm inflow-angle as a discriminant for rupture in sidewall cerebral aneurysms: morphometric and computational fluid dynamic analysis. Stroke. 2010;41(7):1423–30. doi: 10.1161/STROKEAHA.109.570770 20508183
18. Xu J, Wu Z, Yu Y, Lv N, Wang S, Karmonik C, et al. Combined Effects of Flow Diverting Strategies and Parent Artery Curvature on Aneurysmal Hemodynamics: A CFD Study. PLoS ONE. 2015;10(9):e0138648. doi: 10.1371/journal.pone.0138648 26398847
19. Duan Z, Li Y, Guan S, Ma C, Han Y, Ren X, et al. Morphological parameters and anatomical locations associated with rupture status of small intracranial aneurysms. Scientific Reports. 2018;8(1):6440. doi: 10.1038/s41598-018-24732-1 29691446
20. Tateshima S, Chien A, Sayre J, Cebral J, Viñuela F. The effect of aneurysm geometry on the intra-aneurysmal flow condition. Neuroradiology. 2010;52(12):1135–41. doi: 10.1007/s00234-010-0687-4 20373097
21. Long Y, Yu H, Zhuo Z, Zhang Y, Wang Y, Yang X, et al. A geometric scaling model for assessing the impact of aneurysm size ratio on hemodynamic characteristics. Biomed Eng Online. 2014;17:13–17.
22. Valencia Alvaro, Solis Francisco. Blood flow dynamics and arterial wall interaction in a saccular aneurysm model of the basilar artery. Computers & Structures. 2006;84(21):1326–1337. doi: 10.1016/j.compstruc.2006.03.008
23. Yamaguchi R, Tanaka G, Liu H. Effect of Elasticity on Flow Characteristics Inside Intracranial Aneurysms. Int J Neurol Neurother. 2016;3:049. doi: 10.23937/2378-3001/3/3/1049
24. Xu L, Sugawara M, Tanaka G, Ohta M, Liu H, Yamaguchi R. Effect of elasticity on wall shear stress inside cerebral aneurysm at anterior cerebral artery. Technol Health Care. 2016;24(3):349–57. doi: 10.3233/THC-161135 26835728
25. Chadwick P. The existence and uniqueness of solutions of two problems in the Mooney-Rivlin theory for rubber. Journal of Elasticity. 1972;2(2):123–128. doi: 10.1007/BF00046061
26. Skalak R, Tozeren A, Zarda RP, Chien S. Strain Energy Function of Red Blood Cell Membranes. Biophys J. 1973;13(3):245–264. doi: 10.1016/S0006-3495(73)85983-1 4697236
27. Fung YC. Biomechanics: Mechanical Properties of Living Tissues. Springer-Verlag. 1981.
28. Volokh KY, Vorp DA. A model of growth and rupture of abdominal aortic aneurysm. J Biomech. 2008;41(5):1015–1021. doi: 10.1016/j.jbiomech.2007.12.014 18255074
29. Rodríguez J, Merodio J. A new derivation of the bifurcation conditions of inflated cylindrical membranes of elastic material under axial loading: Application to aneurysm formation. Mechanics Research Communications. 2011;38(3):203–210. doi: 10.1016/j.mechrescom.2011.02.004
30. Pierce DM, Maier F, Weisbecker H, Viertler C, Verbrugghe P, Famaey N, et al. Human thoracic and abdominal aortic aneurysmal tissues: Damage experiments, statistical analysis and constitutive modeling. J Mech Behav Biomed Mater. 2015;41:92–107. doi: 10.1016/j.jmbbm.2014.10.003 25460406
31. Charalambous HP, Roussis PC, Giannakopoulos AE. The Effect of Strain Hardening on the Dynamic Response of Human Artery Segments. Open Biomed Eng J. 2017;26(11):85–110. doi: 10.2174/1874120701711010085
32. Krüger T, Varnik F, Raabe D. Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method. Comput Math with Appl. 2011;61:3485–3505. doi: 10.1016/j.camwa.2010.03.057
33. Kamm R, Grodzinsky A. 20.310J Molecular, Cellular, and Tissue Biomechanics. Spring 2015. Massachusetts Institute of Technology: MIT OpenCourseWare, https://ocw.mit.edu.
34. Krafczyk M, Cerrolaza M, Schulz M, Rank E. Analysis of 3D transient blood flow passing through an artificial aortic valve by Lattice-Boltzmann methods. J Biomech. 1998;31(5):453–462. doi: 10.1016/s0021-9290(98)00036-0 9727343
35. Cebral JR, Castro MA, Appanaboyina S, Putman CM, Millan D, Frangi AF. Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: Technique and sensitivity. IEEE Trans Med Imaging. 2005;(24):457–467. doi: 10.1109/tmi.2005.844159 15822804
36. Shishir SS, Miah MAK, Islam AKMS, Hasan ABMT. Flow Dynamics in Cerebral Aneurysm—A CFD Simulation. Procedia Engineering. 2015;105:919–927.
37. Ariane M, Allouche MH, Bussone M, Giacosa F, Bernard F, Barigou M, Alexiadis A. Discrete multi-physics: A mesh-free model of blood flow in flexible biological valve including solid aggregate formation. PLoS ONE. 2017;12(4):e0174795. doi: 10.1371/journal.pone.0174795 28384341
38. Han HC, Fung YC. Longitudinal strain of canine and porcine aortas. J Biomech. 1995;28(5):637–41. doi: 10.1016/0021-9290(94)00091-h 7775500
39. 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(11):1281–99. doi: 10.1114/1.1326031 11212947
40. Krüger T, Gross M, Raabe D, Varnik F. Crossover from tumbling to tank-treading-like motion in dense simulated suspensions of red blood cells. Soft Matter. 2013;9(37):9008–9015. doi: 10.1039/c3sm51645h 25353617
41. Gross M, Krüger T, Varnik F. Fluctuations and diffusion in sheared athermal suspensions of deformable particles. EPL (Europhysics Letters). 2014;108(6):68006. doi: 10.1209/0295-5075/108/68006
42. Gross M, Krüger T, Varnik F. Rheology of dense suspensions of elastic capsules: normal stresses, yield stress, jamming and confinement effects. Soft Matter. 2014;10(24):4360–4372. doi: 10.1039/c4sm00081a 24796957
43. Gompper G, Kroll D. Random Surface Discretizations and the Renormalization of the Bending Rigidity. J Phys I France. 1996;6(10):1305–1320. doi: 10.1051/jp1:1996246
44. Taylor CL, Yuan Z, Selman WR, Ratcheson RA, Rimm AA. Cerebral arterial neurysm formation and rupture in 20,767 elderly patients: hypertension and other risk factors. J Neurosurg. 1995;83(5):812–819. doi: 10.3171/jns.1995.83.5.0812 7472548
45. Juvela S. Prehemorrhage risk factors for fatal intracranial aneurysm rupture. Stroke. 2003;34(8):1852–7. doi: 10.1161/01.STR.0000080380.56799.DD 12829865
46. Kim MO, Adji A, O’Rourke MF, Avolio AP, Smielewski P, Pickard JD, et al. Principles of cerebral hemodynamics when intracranial pressure is raised: lessons from the peripheral circulation. J Hypertens. 2015;33(6):1233–41. doi: 10.1097/HJH.0000000000000539 25764046
47. Lallemand P, Luo LS. Theory of the Lattice Boltzmann Method: Dispersion, dissipation, isotropy, Galilean invariance, and stability. Phys Rev E. 2000;61(6):6546–62. doi: 10.1103/PhysRevE.61.6546
48. d’Humières D, Ginzburg I. Viscosity independent numerical errors for Lattice Boltzmann models: From recurrence equations to “magic” collision numbers. Comput Math Appl. 2009;58:823–840. doi: 10.1016/j.camwa.2009.02.008
49. Asinari P, Karlin I. Quasiequilibrium lattice Boltzmann models with tunable bulk viscosity for enhancing stability. Phys Rev E. 2010;81:016702. doi: 10.1103/PhysRevE.81.016702
50. Guo Z, Zheng C, Shi B. Discrete lattice effects on the forcing term in the lattice Boltzmann method. Phys Rev E. 2002;65:046308. doi: 10.1103/PhysRevE.65.046308
51. Krüger T, Varnik F, Raabe D. Shear stress in lattice Boltzmann simulations. Phys Rev E. 2009;79:046704. doi: 10.1103/PhysRevE.79.046704
52. Peskin CS. Flow patterns around heart valves: A digital computer method for solving the equations of motion. Ph.D. thesis, Albert Einstein College of Medicine. 1972.
53. Womersley JR. Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known. J Physiol. 1955;127(3):553–563. doi: 10.1113/jphysiol.1955.sp005276 14368548
54. Gilmanov A, Sotiropoulos F, Balaras E. A general reconstruction algorithm for simulating flows with complex 3D immersed boundaries on Cartesian grids. J Comput Phys. 2003;191:660–669. doi: 10.1016/S0021-9991(03)00321-8
55. Arthurs KM, Moore LC, Peskin CS, Pitman EB, Layton HE. Modeling arteriolar flow and mass transport using the immersed boundary method. J Comput Phys. 1998;147(2):402–440. doi: 10.1006/jcph.1998.6097
56. Peskin CS. The immersed boundary method. Acta Numerica. 2002;11:479–517. doi: 10.1017/S0962492902000077
57. Smith KM, Moore LC, Layton HE. Advective transport of nitric oxide in a mathematical model of the afferent arteriole. Am J Physiol Renal Physiol. 2003;284(5):F1080–96. doi: 10.1152/ajprenal.00141.2002 12712988
58. Updegrove A, Wilson NM, Merkow J, Lan H, Marsden AL, Shadden SC. SimVascular—An open source pipeline for cardiovascular simulation. Ann Biomed Eng. 2017;45(3):525–541. doi: 10.1007/s10439-016-1762-8 27933407
59. Shapiro AH. Steady flow in collapsed tubes. Trans. ASME K: J Biomech Engng. 1977;99:126–147.
60. Jensen OE, Pedley TJ. The existence of steady flow in a collapsed tube. J Fluid Mech. 1989;206:339–374. doi: 10.1017/S0022112089002326
61. Gibo H, Carver CC, Rhoton AL Jr, Lenkey C, Mitchell RJ. Microsurgical anatomy of the middle cerebral artery. J Neurosurg. 1981;54(2):151–69. doi: 10.3171/jns.1981.54.2.0151 7452329
62. Mouches P, Forkert ND. A statistical atlas of cerebral arteries generated using multi-center MRA datasets from healthy subjects. Scientific Data. 2019;6:29. doi: 10.1038/s41597-019-0034-5 30975990
63. Isaksen JG, Bazilevs Y, Kvamsdal T, Zhang Y, Kaspersen JH, Waterloo K, et al. Determination of wall tension in cerebral artery aneurysms by numerical simulation. Stroke. 2008;39(12):3172–8. doi: 10.1161/STROKEAHA.107.503698 18818402
64. Aaslid R, Markwalder TM, Nornes H. Noninvasive transcranial Doppler ultrasound recording of flow velocity in basal cerebral arteries. J Neurosurg. 1982;57(6):769–74. doi: 10.3171/jns.1982.57.6.0769 7143059
65. Krejza J, Szydlik P, Liebeskind DS, Kochanowicz J, Bronov O, Mariak Z, et al. Age and sex variability and normal reference values for the V(MCA)/V(ICA) index. AJNR Am J Neuroradiol. 2005;26(4):730–5. 15814913
66. Ogoh S, Fadel PJ, Zhang R, Selmer C, Jans Ø, Secher NH, et al. Middle cerebral artery flow velocity and pulse pressure during dynamic exercise in humans. Am J Physiol Heart Circ Physiol. 2005;288(4):H1526–31. doi: 10.1152/ajpheart.00979.2004 15591094
67. Bouvy WH, Geurts LJ, Kuijf HJ, Luijten PR, Kappelle LJ, Biessels GJ, et al. Assessment of blood flow velocity and pulsatility in cerebral perforating arteries with 7-T quantitative flow MRI. NMR Biomed. 2016;29(9):1295–304. doi: 10.1002/nbm.3306 25916399
68. Warriner RK, Johnston KW, Cobbold RS. A viscoelastic model of arterial wall motion in pulsatile flow: implications for Doppler ultrasound clutter assessment. Physiol Meas. 2008;29(2):157–79. doi: 10.1088/0967-3334/29/2/001 18256449
Článek vyšel v časopise
PLOS One
2020 Číslo 1
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Severity of misophonia symptoms is associated with worse cognitive control when exposed to misophonia trigger sounds
- Chemical analysis of snus products from the United States and northern Europe
- Calcium dobesilate reduces VEGF signaling by interfering with heparan sulfate binding site and protects from vascular complications in diabetic mice
- Effect of Lactobacillus acidophilus D2/CSL (CECT 4529) supplementation in drinking water on chicken crop and caeca microbiome