The Current Limitations and Advanced Analysis of Hemodynamic Study of Cerebral Aneurysms
Article information
Abstract
Among the various perspectives on cerebrovascular diseases, hemodynamic analysis—which has recently garnered interest—is of great help in understanding cerebrovascular diseases. Computational fluid dynamics (CFD) analysis has been the primary hemodynamic analysis method, and studies on cerebral aneurysms have been actively conducted. However, owing to the intrinsic limitations of the analysis method, the role of wall shear stress (WSS), the most representative parameter, remains controversial. High WSS affects the formation of cerebral aneurysms; however, no consensus has been reached on the role of WSS in the growth and rupture of cerebral aneurysms. Therefore, this review aimed to briefly introduce the up-to-date results and limitations made through CFD analysis and to inform the need for a new hemodynamic analysis method.
INTRODUCTION
Because the rupture of a cerebral aneurysm has severe consequences, there is a strong tendency to treat it before rupture. However, the diagnosis of unruptured cerebral aneurysms has increased with the development of diagnostic modalities, prompting the question of whether all unruptured cerebral aneurysms should be treated. As the adverse effects of treatment cannot be disregarded, it is important to select unruptured cerebral aneurysms that require treatment on a reasonable basis.
Therefore, selecting and treating unruptured cerebral aneurysms at a high risk of rupture is critical. Several studies have been conducted to predict the risk of cerebral aneurysm rupture from various perspectives. Morphological, pathophysiological, and clinical factors have been suggested as the basis for judgment; however, explanations in the field of hemodynamics have only recently been of focus [1-12]. Computational fluid dynamics (CFD) is the most commonly used method in hemodynamic studies.
With the advent of CFD analysis in clinical practice, hemodynamic studies of cerebral aneurysms have become active. CFD analysis is widely used in mechanical engineering; however, because the shape of the cerebral aneurysms was precisely implemented with digital subtraction angiography (DSA), it can be introduced in the hemodynamic field of cerebral aneurysms.
COMPUTATIONAL FLUID DYNAMICS ANALYSIS
Computational Fluid Dynamics Process
The fundamental principle of CFD analysis involves discretizing a fluid domain into control volumes and solving the governing equations that describe fluid motion. These governing equations encompass the continuity equation and the Navier-Stokes equation, which are derived from the conservation laws of mass and momentum. The continuity equation and the incompressible Navier-Stokes equation are expressed as follows:
Continuity equation:
Incompressible Navier-Stokes equation:
where
General process of computed fluid dynamics (CFD). Following the reconstruction of the 3-dimensional (3D) model of a cerebral artery, meshing and boundary condition applications are performed. Hemodynamic parameters are then computed.
CFD analysis proves to be an invaluable tool for quantitatively assessing blood flow within cerebral arteries and aneurysms. The velocity obtained through CFD analysis can be utilized in diverse applications pertaining to cerebral aneurysms. Fig. 2 demonstrates an example of velocity evaluation in the circle of Willis and a cerebral aneurysm.
Examples of computed fluid dynamics (CFD) results. Velocity evaluation in a circle of Willis (A) and a cerebral aneurysm (B). The velocity obtained through CFD analysis is color-coded and used for calculating hemodynamic parameters.
Assumption as Limitation
Several assumptions were premised for the CFD analysis introduced for the hemodynamics of the intracranial vessels. Blood flow was assumed to be Newtonian for the convenience of simulation in the CFD analysis. In addition, the blood vessel was assumed to be a rigid wall in the CFD analysis. Therefore, the effects of the wall thickness and material properties of the blood vessels on the formation, growth, and rupture of cerebral aneurysms were excluded, and only the effect of blood flow was considered. However, since a cerebral aneurysm is a morphological change in blood vessels, it cannot be assumed that blood vessels do not change.
In addition to these intrinsic limitations of CFD analysis, it is not patient-specific. To analyze natural phenomena, the material properties of the target must be identified. However, the material properties of blood vessels used in previous studies were not patient-specific and were obtained from a cohort of healthy patients. Each patient had different material properties and vessel thicknesses. Furthermore, the velocity, viscosity, and blood pressure corresponding to the boundary conditions for the CFD analysis had different values for each patient. To date, studies have assumed that the boundary conditions are uniformly the same value for each study.
CFD analysis is based on several assumptions, which can be an intrinsic limitation of this analysis method. Nevertheless, CFD analysis is the underlying method for hemodynamic analysis, and we need to study the most important parameters in the CFD analysis of cerebral aneurysms.
Hemodynamic Parameter: Wall Shear Stress
WSS is among the most important hemodynamic parameters in CFD analysis and is a frictional force from the blood flow tangential to the arterial lumen. The relative difference in velocity between two parallel objects creates shear stress [13]. In the normal range (1.5–2.5 Pa) of WSS, endothelial function is regulated [14]. However, if the WSS is outside the normal range, histological changes related to the aneurysm may occur. Although WSS is conventionally associated with the natural history (formation, growth, and rupture) of cerebral aneurysms, the effect of the WSS on each natural history remains controversial.
Formation of Aneurysm: High Wall Shear Stress
Several studies utilizing animal models have identified a high WSS as an important parameter in aneurysm formation [15-21]. In addition, Can and Du [22] reported a strong positive correlation between elevated WSS and the location of aneurysm formation in their systematic review and meta-analysis. The authors included 19 studies that investigated WSS using CFD for geometrical models of intracranial aneurysms and found that high WSS was associated with formation and low WSS was associated with rupture of intracranial aneurysms [22].
Growth and Rupture of Aneurysms: Controversy
Unlike the consensus that aneurysm formation occurs in regions with high WSS, the exact role of WSS in the growth and rupture of intracranial aneurysms is controversial. It is unclear whether a high or low WSS plays a principal role in growth and rupture. There are “high” and “low” WSS theories, both of which explain why the hemodynamic environment within the aneurysm interacts with the cellular elements of the aneurysm wall, resulting in further weakening. However, the differences revolve around the mechanisms that cause the weakening of the weakening [21].
(1) High wall shear stress theory
According to the high WSS theory, high WSS causes endothelial injury and initiates processes of wall remodeling and degeneration [13]. Castro et al. [23] investigated anterior communicating artery aneurysms and suggested that aneurysms with a high WSS were more likely to rupture than those with a low WSS. Cebral et al. [24] reported that a high WSS is related to aneurysmal rupture. Hassan et al. [25] suggested that aneurysms rupture because blood flow impingement on the aneurysmal wall produced a high WSS.
(2) Low wall shear stress theory
In contrast, the low WSS theory suggests that stagnation of blood within the aneurysm leads to red blood cell aggregation and build-up of platelets and leukocytes [21]. These changes cause intimal damage and inflammatory cell infiltration of the vessel wall, leading to wall degeneration and reduced ability to withstand physiological hemodynamic forces [21]. Several studies have found that aneurysm growth and rupture are more likely to occur in regions with abnormally low WSS [2,21,26,27].
(3) Unifying the wall shear stress theories
To unify the controversial results of previous studies on WSS, Meng et al. [28] published a review article. They indicated that the leading controversial reports on WSS depended on differences in the natural history of each analyzed aneurysm [28]. They hypothesized that the mural cell-mediated pathway associated with high WSS may be responsible for the growth and rupture of small or secondary bleb aneurysm phenotypes, whereas the inflammation-mediated pathway associated with low WSS may be responsible for the growth and rupture of large atherosclerotic and thrombotic cerebral aneurysm phenotypes [28]. Furthermore, although the exact mechanism involved is unknown, Cebral et al. [29] recently supported the idea that different hemodynamic conditions are associated with different aneurysm phenotypes.
However, some researchers have reported no significant difference in the WSS values between stable and unstable aneurysms. In their investigation of 178 patients with 198 unruptured cerebral aneurysms, Ramachandran et al. [30] found that low WSS did not statistically discriminate between stable and unstable aneurysms. Based on their study of 33 unruptured aneurysms, Sforza et al. [31] suggested that concentrated inflow streams blown into complex intrasaccular flow patterns—causing nonuniform WSS distributions—may represent the characteristics of a hemodynamic environment that predisposes the aneurysm to grow.
(4) Inconsistency of wall shear stress
As mentioned previously, high WSS experienced during aneurysm formation may lead to a low WSS as the geometry of the aneurysm changes [28]. It is unknown whether this inconsistency is a result of intrinsic study limitations or inherent complexity and heterogeneity of the biological mechanisms involved in aneurysm evolution [32,33]. Meng et al. [28] explained that both high and low WSS can independently lead to intracranial aneurysm growth and rupture, albeit through different biological mechanisms. However, despite the explanation, it is not yet possible to predict the state in which a high WSS will cause growth or rupture, and that in which a low WSS will cause growth and rupture. Clinicians must predict where growth or rupture will occur in the same patient, either in the high- or low-WSS region.
If the intrinsic limitations of CFD analysis lead to inconsistencies and controversial results, it may be helpful to apply other analytical methods.
Other Hemodynamic Parameters
There are several hemodynamic parameters in CFD analysis. Previous hemodynamic studies define several hemodynamic parameters associated with natural history of cerebral aneurysms. Among the various hemodynamic parameters, oscillatory shear index (OSI) and WSS gradient (WSSG) have demonstrated their role in aneurysm rupture risk analysis.
OSI is a dimensionless parameter that indicates how the direction of WSS changes at a specific location during a cardiac cycle.
Previous studies reported that a higher OSI was observed in ruptured than in unruptured aneurysms, or high OSI corresponded to the rupture point [27]. They suggested that low WSS and high OSI predict rupture risk of cerebral aneurysms [27, 28].
WSSG is defined as the spatial derivative of WSS along the direction of flow. It can be thought of as the change in WSS along the length of the vessel [14]. Accelerating flow creates a positive WSSG, while decelerating flow creates a negative WSSG. Meng et al. [28] explained that the formation of cerebral aneurysms occurs in regions exposed to high WSS with a positive WSSG in the mural cell-mediated pathway theory.
As shown in the defintions presented above, both OSI and WSSG are secondary parameters derived from the concept of WSS. Cho et al. [2] defined a novel parameter, the so-called combined hemodynamic parameter, associated with aneurysmal rupture. Similarly, it was derived from the perception of WSS.
FLUID–STRUCTURE INTERACTION ANALYSIS
Another Hemodynamic Study
A cerebral aneurysm is a structural change that occurs when the shape of blood vessels changes due to the influence of blood flow. However, because CFD analysis assumes rigid blood vessel walls, the effects of the blood vessels are ignored. In contrast, fluid-structure interaction (FSI) analysis assumes that the blood vessel is not rigid but deformable; therefore, the effects of the wall thickness and mechanical properties of the aneurysms can be evaluated in a more realistic manner. In FSI analysis, using the pressure induced by the blood flow calculated from CFD analysis, the deformation of a cerebral artery can be evaluated by calculating the strain (Fig. 3). Strain represents how much a cerebral artery is stretched. Higher strain indicates that the blood vessel is more severely stretched. FSI analysis involves complex computational techniques, and the results depend on the thickness or elastic modulus of the vessel, which is difficult to obtain using common diagnostic imaging devices [34]. Nevertheless, several hemodynamic studies have used FSI analysis [35-37].
Explanation of fluid-structure interaction (FSI) process. In FSI analysis, the strain used to evaluate the deformation of the cerebral artery is calculated using the pressure on the blood vessels induced by blood flow. The values of pressure in computed fluid dynamics analysis and equivalent strain in FSI analysis are color-coded.
Cho et al. [35] suggested in FSI analysis that strain was more effective than WSS in predicting the rupture risk of cerebral aneurysms. In addition, unlike WSS, which shows inconsistency as it pregresses from formation to growth and rupture, high strain continues to be involved in the natural history of cerebral aneurysms. Kim et al. [36] reported that strain has an important role in the formation of cerebral aneurysms.
Limitations of Fluid-Structure Interaction
Despite overcoming several intrinsic limitations of CFD analysis, FSI analysis still has some limitations. The material properties and wall thickness of the cerebral artery are greatly influenced by the deformation of the cerebral artery. However, because material properties and wall thicknesses vary from person to person, it is difficult to determine the material properties and wall thicknesses of individuals. Although there have been efforts to study the material properties and wall thickness of the cerebral artery, there is still a need for further research in this area.
DIFFERENCES IN SOURCE DATA
Both CFD and FSI analyses use computational simulations, and 3D angiography should be used as the source data. Therefore, the quality of source data is very important. 3D rotational angiography (3DRA) from DSA provides the highest resolution (approximately 0.2–0.3 mm) but is difficult to perform repeatedly as an invasive modality [38]. Modern multidetector computed tomographic angiography (CTA) provides resolutions of approximately 0.4–0.5 mm, but has a critical disadvantage involving artifacts from nearby bone [38]. Time-of-flight magnetic resonance angiography (MRA) offers a poorer resolution of approximately 0.6–0.8 mm and has artifacts due to slowly recirculating flows that may be present in some aneurysms [38].
In a recent study comparing 3DRA and CTA, Lauric et al. [39] reported significant differences in low WSS between the 2 modalities. In addition, unlike CTA and MRA, 3DRA is performed by injecting contrast agents directly into the arteries, which may affect their tortuosity.
These can be sources of variability among modalities, which can affect the results of hemodynamic analysis. Therefore, research on the differences in the results of hemodynamic studies according to differences in modalities needs to be conducted carefully.
LIMITATIONS OF HEMODYNAMIC ANALYSIS
The main limitation of hemodynamic analysis, which is being studied worldwide, is that standardized research methodologies have not been established. Standardized methodologies have not been presented at all stages, such as segmentation, meshing, and solvers, which are the basic stages of hemodynamic analysis or the selection of image modalities used as source data. In addition, the conditions necessary for analysis, such as the thickness and material properties of blood vessels, and inlet conditions are not yet strictly patient-specific. A consensus on these limitations must be attained through discussions with researchers worldwide.
Acknowledgements
We thank Hyeondong Yang for his assistance in this study.
Notes
Fund
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. 2021R1F1A1049435).
Ethics Statement
The approval of the institutional ethics committee and informed consents were waived since its nature lies on literature review. We do not use patient information that may identify.
Conflicts of Interest
The author has no conflicts to disclose.