How patient specific are patient-specific computational models of cerebral aneurysms? An overview of sources of error and variability

Computational methods for biofabrication in tissue engineering and regenerative medicine - a literature review

This literature review rigorously examines the growing scientific interest in computational methods for Tissue Engineering and Regenerative Medicine biofabrication, a leading-edge area in biomedical innovation, emphasizing the need for accurate, multi-stage, and multi-component biofabrication process models. The paper presents a comprehensive bibliometric and contextual analysis, followed by a literature review, to shed light on the vast potential of computational methods in this domain. It reveals that most existing methods focus on single biofabrication process stages and components, and there is a significant gap in approaches that utilize accurate models encompassing both biological and technological aspects. This analysis underscores the indispensable role of these methods in understanding and effectively manipulating complex biological systems and the necessity for developing computational methods that span multiple stages and components. The review concludes that such comprehensive computational methods are essential for developing innovative and efficient Tissue Engineering and Regenerative Medicine biofabrication solutions, driving forward advancements in this dynamic and evolving field.

Toward a physiological model of vascular wall vibrations in the arteriovenous fistula

The mechanism behind hemodialysis arteriovenous fistula (AVF) failure remains poorly understood, despite previous efforts to correlate altered hemodynamics with vascular remodeling. We have recently demonstrated that transitional flow induces high-frequency vibrations in the AVF wall, albeit with a simplified model. This study addresses the key limitations of our original fluid-structure interaction (FSI) approach, aiming to evaluate the vibration response using a more realistic model. A 3D AVF geometry was generated from contrast-free MRI and high-fidelity FSI simulations were performed. Patient-specific inflow and pressure were incorporated, and a three-term Mooney-Rivlin model was fitted using experimental data. The viscoelastic effect of perivascular tissue was modeled with Robin boundary conditions. Prescribing pulsatile inflow and pressure resulted in a substantial increase in vein displacement ( + 400 %) and strain ( + 317 %), with a higher maximum spectral frequency becoming visible above -42 dB (from 200 to 500 Hz). Transitioning from Saint Venant-Kirchhoff to Mooney-Rivlin model led to displacement amplitudes exceeding 10 micrometers and had a substantial impact on strain ( + 116 %). Robin boundary conditions significantly damped high-frequency displacement ( - 60 %). Incorporating venous tissue properties increased vibrations by 91%, extending up to 700 Hz, with a maximum strain of 0.158. Notably, our results show localized, high levels of vibration at the inner curvature of the vein, a site known for experiencing pronounced remodeling. Our findings, consistent with experimental and clinical reports of bruits and thrills, underscore the significance of incorporating physiologically plausible modeling approaches to investigate the role of wall vibrations in AVF remodeling and failure.

A comprehensive MRI-based computational model of blood flow in compliant aorta using radial basis function interpolation

BACKGROUND

Properly understanding the origin and progression of the thoracic aortic aneurysm (TAA) can help prevent its growth and rupture. For a better understanding of this pathogenesis, the aortic blood flow has to be studied and interpreted in great detail. We can obtain detailed aortic blood flow information using magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) with a prescribed motion of the aortic wall.

METHODS

We performed two different types of simulations-static (rigid wall) and dynamic (moving wall) for healthy control and a patient with a TAA. For the latter, we have developed a novel morphing approach based on the radial basis function (RBF) interpolation of the segmented 4D-flow MRI geometries at different time instants. Additionally, we have applied reconstructed 4D-flow MRI velocity profiles at the inlet with an automatic registration protocol.

RESULTS

The simulated RBF-based movement of the aorta matched well with the original 4D-flow MRI geometries. The wall movement was most dominant in the ascending aorta, accompanied by the highest variation of the blood flow patterns. The resulting data indicated significant differences between the dynamic and static simulations, with a relative difference for the patient of 7.47±14.18% in time-averaged wall shear stress and 15.97±43.32% in the oscillatory shear index (for the whole domain).

CONCLUSIONS

In conclusion, the RBF-based morphing approach proved to be numerically accurate and computationally efficient in capturing complex kinematics of the aorta, as validated by 4D-flow MRI. We recommend this approach for future use in MRI-based CFD simulations in broad population studies. Performing these would bring a better understanding of the onset and growth of TAA.

Characterization of Maximum Wall Shear Stress Points in Unruptured Cerebral Aneurysms Using Four-dimensional Flow Magnetic Resonance Imaging

BACKGROUND

Maximum wall shear stress (maxWSS) points of unruptured cerebral aneurysms (UCAs) may cause wall remodeling leading to rupture. We characterized maxWSS points and their inherent intra-aneurysmal flow structures in a sizable cohort of saccular UCAs using four-dimensional (4D) flow magnetic resonance imaging (MRI).

METHODS

After contrast administration, 50 saccular UCAs were subjected to 4D flow MRI using a 1.5 T MRI scanner. Post-processing of obtained data was performed using commercially available software. The maxWSS points and maxWSS values were evaluated. The maxWSS values were statistically compared between aneurysm groups.

RESULTS

The maxWSS point was located on the aneurysm apex in 9 (18.0%), body in 2 (4.0%), and neck in 39 (78.0%) UCAs. The inherent intra-aneurysmal flow structure of the maxWSS point was an inflow zone in 34 (68.0%) UCAs, an inflow jet in 8 (16.0%), and an impingement zone in 8 (16.0%). The maxWSS point on the neck had significantly higher maxWSS values than those points on the other wall areas (P = 0.008). The maxWSS values of the maxWSS points on the apex and on the impingement zone were not significantly different compared with those of the other maxWSS points.

CONCLUSION

The maxWSS points existed preferentially on the aneurysmal neck adjacent to the inflow zone with higher maxWSS values. The maxWSS points existed occasionally on the aneurysmal apex adjacent to the impingement zone. 4D flow MRI may be helpful to discriminate saccular UCAs with higher-risk maxWSS points that can cause wall remodeling leading to rupture.

Turbulent blood flow in a cerebral artery with an aneurysm

Unruptured intracranial aneurysms are common in the general population, and many uncertainties remain when predicting rupture risks and treatment outcomes. One of the cutting-edge tools used to investigate this condition is computational fluid dynamics (CFD). However, CFD is not yet mature enough to guide the clinical management of this disease. In addition, recent studies have reported significant flow instabilities when refined numerical methods are used. Questions remain as to how to properly simulate and evaluate this flow, and whether these instabilities are really turbulence. The purpose of the present study is to evaluate the impact of the simulation setup on the results and investigate the occurrence of turbulence in a cerebral artery with an aneurysm. For this purpose, direct numerical simulations were performed with up to 200 cardiac cycles and with data sampling rates of up to 100,000 times per cardiac cycle. Through phase-averaging or triple decomposition, the contributions of turbulence and of laminar pulsatile waves to the velocity, pressure and wall shear stress fluctuations were distinguished. For example, the commonly used oscillatory shear index was found to be closely related to the laminar waves introduced at the inlet, rather than turbulence. The turbulence energy cascade was evaluated through energy spectrum estimates, revealing that, despite the low flow rates and Reynolds number, the flow is turbulent near the aneurysm. Phase-averaging was shown to be an approach that can help researchers better understand this flow, although the results are highly dependent on simulation setup and post-processing choices.

Identifying hemodynamic factors associated with the rupture of anterior communicating artery aneurysms based on global modeling of blood flow in the cerebral artery network

Anterior communicating artery (ACoA) aneurysms are more prone to rupture compared to aneurysms present in other cerebral arteries. We hypothesize that systemic blood flow in the cerebral artery network plays an important role in shaping intra-aneurysmal hemodynamic environment thereby affecting the rupture risk of ACoA aneurysms. The majority of existing numerical studies in this field employed local modeling methods where the physical boundaries of a model are confined to the aneurysm region, which, though having the benefit of reducing computational cost, may compromise the physiological fidelity of numerical results due to insufficient account of systemic cerebral arterial hemodynamics. In the present study, we firstly carried out numerical experiments to address the difference between the outcomes of local and global modeling methods, demonstrating that local modeling confined to the aneurysm region results in inaccurate predictions of hemodynamic parameters compared with global modeling of the ACoA aneurysm as part of the cerebral artery network. Motivated by this finding, we built global hemodynamic models for 40 ACoA aneurysms (including 20 ruptured and 20 unruptured ones) based on medical image data. Statistical analysis of the computed hemodynamic data revealed that maximum wall shear stress (WSS), minimum WSS divergence, and maximum WSS gradient differed significantly between the ruptured and unruptured ACoA aneurysms. Optimal threshold values of high/low WSS metrics were determined through a series of statistical tests. In the meantime, some morphological parameters of aneurysms, such as large nonsphericity index, aspect ratio, and bottleneck factor, were found to be associated closely with aneurysm rupture. Furthermore, multivariate logistic regression analyses were performed to derive models combining hemodynamic and morphological parameters for discriminating the rupture status of aneurysms. The capability of the models in rupture status discrimination was high, with the area under the receiver operating characteristic curve reaching up to 0.9. The findings of the study suggest that global modeling of the cerebral artery network is essential for reliable quantification of hemodynamics in ACoA aneurysms, disturbed WSS and irregular aneurysm morphology are associated closely with aneurysm rupture, and multivariate models integrating hemodynamic and morphological parameters have high potential for assessing the rupture risk of ACoA aneurysms.

Impact of left atrial wall motion assumptions in fluid simulations on proposed predictors of thrombus formation

Atrial fibrillation (AF) poses a significant risk of stroke due to thrombus formation, which primarily occurs in the left atrial appendage (LAA). Medical image-based computational fluid dynamics (CFD) simulations can provide valuable insight into patient-specific hemodynamics and could potentially enhance personalized assessment of thrombus risk. However, the importance of accurately representing the left atrial (LA) wall dynamics has not been fully resolved. In this study, we compared four modeling scenarios; rigid walls, a generic wall motion based on a reference motion, a semi-generic wall motion based on patient-specific motion, and patient-specific wall motion based on medical images. We considered a LA geometry acquired from 4D computed tomography during AF, systematically performed convergence tests to assess the numerical accuracy of our solution strategy, and quantified the differences between the four approaches. The results revealed that wall motion had no discernible impact on LA cavity hemodynamics, nor on the markers that indicate thrombus formation. However, the flow patterns within the LAA deviated significantly in the rigid model, indicating that the assumption of rigid walls may lead to errors in the estimated risk factors. In contrast, the generic, semi-generic, and patient-specific cases were qualitatively similar. The results highlight the crucial role of wall motion on hemodynamics and predictors of thrombus formation, and also demonstrate the potential of using a generic motion model as a surrogate for the more complex patient-specific motion. While the present study considered a single case, the employed CFD framework is entirely open-source and designed for adaptability, allowing for integration of additional models and generic motions.

Using DFT on ultrasound measurements to determine patient-specific blood flow boundary conditions for computational hemodynamics of intracranial aneurysms

Boundary conditions (BCs) is one pivotal factor influencing the accuracy of hemodynamic predictions on intracranial aneurysms (IAs) using computational fluid dynamics (CFD) modeling. Unfortunately, a standard procedure to secure accurate BCs for hemodynamic modeling does not exist. To bridge such a knowledge gap, two representative patient-specific IA models (Case-I and Case-II) were reconstructed and their blood flow velocity waveforms in the internal carotid artery (ICA) were measured by ultrasonic techniques and modeled by discrete Fourier transform (DFT). Then, numerical investigations were conducted to explore the appropriate number of samples (N) for DFT modeling to secure the accurate BC by comparing a series of hemodynamic parameters using in-vitro validated CFD modeling. Subsequently, a comprehensive comparison in hemodynamic characteristics under patient-specific BCs and a generalized BC based on a one-dimensional (1D) model was conducted to reinforce the understanding that a patient-specific BC is pivotal for accurate hemodynamic risk evaluations on IA pathophysiology. In addition, the influence of the variance of heart rate/cardiac pulsatile period on hemodynamic characteristics in IA models was studied preliminarily. The results showed that N ≥ 16 for DFT model is a decent choice to secure the proper BC profile to calculate time-averaged hemodynamic parameters, while more data points such as N ≥ 36 can ensure the accuracy of instantaneous hemodynamic predictions. In addition, results revealed the generalized BC could overestimate or underestimate the hemodynamic risks on IAs significantly; thus, patient-specific BCs are highly recommended for hemodynamic modeling for IA risk evaluation. Furthermore, this study discovered the variance of heart rate has rare influences on hemodynamic characteristics in both instantaneous and time-averaged parameters under the assumption of an identical blood flow rate.

On the importance of fundamental computational fluid dynamics toward a robust and reliable model of left atrial flows

Computational fluid dynamics (CFD) studies of left atrial flows have reached a sophisticated level, for example, revealing plausible relationships between hemodynamics and stresses with atrial fibrillation. However, little focus has been on fundamental fluid modeling of LA flows. The purpose of this study was to investigate the spatiotemporal convergence, along with the differences between high- (HR) versus normal-resolution/accuracy (NR) solution strategies, respectively. Rigid wall CFD simulations were conducted on 12 patient-specific left atrial geometries obtained from computed tomography scans, utilizing a second-order accurate and space/time-centered solver. The convergence studies showed an average variability of around 30% and 55% for time averaged wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT), and endothelial cell activation potential (ECAP), even between intermediate spatial and temporal resolutions, in the left atrium (LA) and left atrial appendage (LAA), respectively. The comparison between HR and NR simulations showed good correlation in the LA for WSS, RRT, and ECAP (R 2 > .9 ), but not for OSI (R 2 = .63 ). However, there were poor correlations in the LAA especially for OSI, RRT, and ECAP (R 2 = .55, .63, and .61, respectively), except for WSS (R 2 = .81 ). The errors are comparable to differences previously reported with disease correlations. To robustly predict atrial hemodynamics and stresses, numerical resolutions of 10 M elements (i.e., Δ x = ∼ .5 mm) and 10 k time-steps per cycle seem necessary (i.e., one order of magnitude higher than normally used in both space and time). In conclusion, attention to fundamental numerical aspects is essential toward establishing a plausible, robust, and reliable model of LA flows.

Standardized viscosity as a source of error in computational fluid dynamic simulations of cerebral aneurysms

BACKGROUND

Computational fluid dynamics (CFD) simulations are a powerful tool for studying cerebral aneurysms, capable of evaluating hemodynamics in a way that is infeasible with imaging alone. However, the difficulty of incorporating patient-specific information and inherent obstacles of in vivo validation have limited the clinical usefulness of CFD of cerebral aneurysms. In this work we investigate the effect of using standardized blood viscosity values in CFD simulations of cerebral aneurysms when compared to simulations of the same aneurysms using patient-specific viscosity values derived from hematocrit measurements.

PURPOSE

The objective of this work is to determine the level of error, on average, that is, caused by using standardized values of viscosity in CFD simulations of cerebral aneurysms. By quantifying this error, we demonstrate the need for incorporating patient-specific viscosity in future CFD investigations of cerebral aneurysms.

METHODS

CFD simulations of forty-one cerebral aneurysms were conducted using patient-specific boundary conditions. For each aneurysm two simulations were conducted, one utilizing patient-specific blood viscosity derived from hematocrit measurements and another using a standardized value for blood viscosity. Hemodynamic parameters such as wall shear stress (WSS), wall shear stress gradient (WSSG), and the oscillatory shear index (OSI) were calculated for each of the simulations for each aneurysm. Paired t-tests for differences in the time-averaged maps of these hemodynamic parameters between standardized and patient-specific viscosity simulations were conducted for each aneurysm. Bland-Altman analysis was used to examine the cohort-wide changes in the hemodynamic parameters. Subjects were broken into two groups, those with higher than standard viscosity and those with lower than standard viscosity. An unpaired t-test was used to compare the percent change in WSS, WSSG, and OSI between patient-specific and standardized viscosity simulations for the two cohorts. The percent changes in hemodynamic parameters were correlated against the direction and magnitude of percent change in viscosity, aneurysm size, and aneurysm location. For all t-tests, a Bonferroni-corrected significance level of 0.0167 was used.

RESULTS

63.2%, 41.5%, and 48.7% of aneurysms showed statistically significant differences between patient-specific and standardized viscosity simulations for WSS, WSSG, and OSI respectively. No statistically significant difference was found in the percent changes in WSS, WSSG, and OSI between the group with higher than standard viscosity and those with lower than standard viscosity, indicating an increase in viscosity can cause either an increase or decrease in each of the hemodynamic parameters. On a study-wide level no significant bias was found in either direction for WSS, WSSG, or OSI between the simulation groups due to the bidirectional effect of changing viscosity. No correlation was found between percent change of viscosity and percent change of WSS, WSSG, or OSI, meaning an after-the-fact correction for patient-specific viscosity is not feasible.

CONCLUSION

Standardizing viscosity values in CFD of cerebral aneurysms has a large and unpredictable impact on the calculated WSS, WSSG, and OSI when compared to CFD simulations of the same aneurysms using a patient-specific viscosity. We recommend implementing hematocrit-based patient-specific blood viscosity values for all CFD simulations of cerebral aneurysms.

Fluid dynamic analysis in predicting the recanalization of intracranial aneurysms after coil embolization - A study of spatiotemporal characteristics

Purpose

Hemodynamics play a key role in the management of cerebral aneurysm recanalization after coil embolization; however, the most reliable hemodynamic parameter remains unknown. Previous studies have explored the use of both spatiotemporally averaged and maximal definitions for hemodynamic parameters, based on computational fluid dynamics (CFD) analysis, to build predictive models for aneurysmal recanalization. In this study, we aimed to assess the influence of different spatiotemporal characteristics of hemodynamic parameters on predictive performance.

Methods

Hemodynamics were simulated using CFD for 66 cerebral aneurysms from 65 patients. We evaluated 14 types of spatiotemporal definitions for two hemodynamic parameters in the pre-coiling model and five in virtual post-coiling model (VM) created by cutting the aneurysm from the pre-coiling model. A total of 91 spatiotemporal hemodynamic features were derived and utilized to develop univariate predictor (UP) and multivariate logistic regression (LR) models. The model's performance was assessed using two metrics: the area under the receiver operating characteristic curve (AUROC) and the area under the precision-recall curve (AUPRC).

Results

Different spatiotemporal hemodynamic features exhibited a wide range of AUROC values ranging from 0.224 to 0.747, with 22 feature pairs showing a significant difference in AUROC value (P-value <0.05), despite being derived from the same hemodynamic parameter. PDave,q1 was identified as the strongest UP with AUROC/AUPRC values of 0.747/0.385, yielding sensitivity and specificity value of 0.889 and 0.614 at the optimal cut-off value, respectively. The LR model further improved the prediction performance, having AUROC/AUPRC values of 0.890/0.903. At the optimal cut-off value, the LR model achieved a specificity of 0.877, sensitivity of 0.719, outperforming the UP model.

Conclusion

Our research indicated that the characteristics of hemodynamic parameters in terms of space and time had a significant impact on the development of predictive model. Our findings suggest that LR model based on spatiotemporal hemodynamic features could be clinically useful in predicting recanalization after coil embolization in patients, without the need for invasive procedures.

Sensitivity of hostile hemodynamics to aneurysm geometry via unsupervised shape interpolation

BACKGROUND AND OBJECTIVE

Vessel geometry and hemodynamics are intrinsically linked, whereby geometry determines hemodynamics, and hemodynamics influence vascular remodeling. Both have been used for testing clinical outcomes, but geometry/morphology generally has less uncertainty than hemodynamics derived from medical image-based computational fluid dynamics (CFD). To provide clinical utility, CFD-based hemodynamic parameters must be robust to modeling errors and/or uncertainties, but must also provide useful information not more-easily extracted from shape alone. The objective of this study was to methodically assess the response of hemodynamic parameters to gradual changes in shape created using an unsupervised 3D shape interpolation method.

METHODS

We trained the neural network NeuroMorph on 3 patient-derived intracranial aneurysm surfaces (labelled A, B, C), and then generated 3 distinct morph sequences (A→B, B→C, C→A) each containing 10 interpolated surfaces. From high-fidelity CFD simulation of these, we calculated a variety of common reduced hemodynamic parameters, including many previously associated with aneurysm rupture, and analyzed their responses to changes in shape, and their correlations.

RESULTS

The interpolated surfaces demonstrate complex, gradual changes in branch angles, vessel diameters, and aneurysm morphology. CFD simulation showed gradual changes in aneurysm jetting characteristics and wall-shear stress (WSS) patterns, but demonstrated a range of responses from the reduced hemodynamic parameters. Spatially and temporally averaged parameters including time-averaged WSS, time-averaged velocity, and low-shear area (LSA) showed low variation across all morph sequences, while parameters of flow complexity such as oscillatory shear, spectral broadening, and spectral bandedness indices showed high variation between slightly-altered neighboring surfaces. Correlation analysis revealed a great deal of mutual information with easier-to-measure shape-based parameters.

CONCLUSIONS

In the absence of large clinical datasets, unsupervised shape interpolation provides an ideal laboratory for exploring the delicate balance between robustness and sensitivity of nominal hemodynamic predictors of aneurysm rupture. Parameters like time-averaged WSS and LSA that are highly "robust" may, as a result, be effectively redundant to morphological predictors, whereas more sensitive parameters may be too uncertain for practical clinical use. Understanding these sensitivities may help identify parameters that are capable of providing added value to rupture risk assessment.

Near-wall hemodynamic parameters quantification in in vitro intracranial aneurysms with 7 T PC-MRI

OBJECTIVE

Wall shear stress (WSS) and its derived spatiotemporal parameters have proven to play a major role on intracranial aneurysms (IAs) growth and rupture. This study aims to demonstrate how ultra-high field (UHF) 7 T phase contrast magnetic resonance imaging (PC-MRI) coupled with advanced image acceleration techniques allows a highly resolved visualization of near-wall hemodynamic parameters patterns in in vitro IAs, paving the way for more robust risk assessment of their growth and rupture.

MATERIALS AND METHODS

We performed pulsatile flow measurements inside three in vitro models of patient-specific IAs using 7 T PC-MRI. To this end, we built an MRI-compatible test bench, which faithfully reproduced a typical physiological intracranial flow rate in the models.

RESULTS

The ultra-high field 7 T images revealed WSS patterns with high spatiotemporal resolution. Interestingly, the high oscillatory shear index values were found in the core of low WSS vortical structures and in flow stream intersecting regions. In contrast, maxima of WSS occurred around the impinging jet sites.

CONCLUSIONS

We showed that the elevated signal-to-noise ratio arising from 7 T PC-MRI enabled to resolve high and low WSS patterns with a high degree of detail.

Assessment of the flow-diverter efficacy for intracranial aneurysm treatment considering pre- and post-interventional hemodynamics

Endovascular treatment of intracranial aneurysms with flow diverters (FD) has become one of the most promising interventions. Due to its woven high-density structure they are particularly applicable for challenging lesions. Although several studies have already conducted realistic hemodynamic quantification of the FD efficacy, a comparison with morphologic post-interventional data is still missing. This study analyses the hemodynamics of ten intracranial aneurysm patients treated with a novel FD device. Based on pre- and post-interventional 3D digital subtraction angiography image data, patient-specific 3D models of both treatment states are generated applying open source threshold-based segmentation methods. Using a fast virtual stenting approach, the real stent positions available in the post-interventional data are virtually replicated and both treatment scenarios were characterized using image-based blood flow simulations. The results show FD-induced flow reductions at the ostium by a decrease in mean neck flow rate (51%), inflow concentration index (56%) and mean inflow velocity (53%). Intraluminal reductions in flow activity for time-averaged wall shear stress (47%) and kinetic energy (71%) are present as well. However, an intra-aneurysmal increase in flow pulsatility (16%) for the post-interventional cases can be observed. Patient-specific FD simulations demonstrate the desired flow redirection and activity reduction inside the aneurysm beneficial for thrombosis formation. Differences in the magnitude of hemodynamic reduction exist over the cardiac cycle which may be addressed in a clinical setting by anti-hypertensive treatment in selected cases.

Mean arterial pressure-aneurysm neck ratio predicts the rupture risk of intracranial aneurysm by reflecting pressure at the dome

Background and purpose

The unruptured intracranial aneurysm (UIA) has high disability and mortality rate after rupture, it is particularly important to assess the risk of UIA and to carry out individualized treatment. The objective of this research is to introduce a novel parameter to predict the rupture risk of UIA.

Methods

A total of 649 patients with 964 intracranial aneurysms in our center were enrolled. A novel parameter named mean arterial pressure-aneurysmal neck ratio (MAPN) was defined. Ten baseline clinical features and twelve aneurysm morphological characteristics were extracted to generate the MAPN model. The discriminatory performance of the MAPN model was compared with the PHASES score and the UCAS score.

Results

In hemodynamic analysis, MAPN was positively correlated with wall shear stress and aneurysm top pressure, with Pearson correlation coefficients of 0.887 and 0.791, respectively. The MAPN was larger in the ruptured group (36.62 ± 18.96 vs. 28.38 ± 14.58, P < 0.001). The area under the curve (AUC) of the MAPN was superior than the AUC of aspect ratio (AR) and the bottleneck factor (BN), they were 0.64 (P < 0.001; 95% CI, 0.588-0.692), 0.611 (P < 0.001; 95% CI, 0.559-0.663) and 0.607 (P < 0.001; 95% CI, 0.554-0.660), respectively. The MAPN model constructed by aneurysm size, aneurysm location, presence of secondary sacs and MAPN, demonstrated good discriminatory ability. The MAPN model exhibited superior performance compared with the UCAS score and the PHASES score (the AUC values were 0.799 [P < 0.001; 95% CI, 0.756-0.840], 0.763 [P < 0.001; 95% CI,0.719-0.807] and 0.741 [P < 0.001; 95% CI, 0.695-0.787], respectively; the sensitivities were 0.849, 0.758 and 0.753, respectively).

Conclusions

Research demonstrates the potential of MAPN to augment the clinical decision-making process for assessing the rupture risk of UIAs.

Hemodynamics of anterior circulation intracranial aneurysms with daughter blebs: investigating the multidirectionality of blood flow fields

Recent advances in diagnostic neuroradiological imaging, allowed the detection of unruptured intracranial aneurysms (IAs). The shape - irregular or multilobular - of the aneurysmal dome, is considered as a possible rupture risk factor, independently of the size, the location and patient medical background. Disturbed blood flow fields in particular is thought to play a key role in IAs progression. However, there is an absence of widely-used hemodynamic indices to quantify the extent of a multi-directional disturbed flow. We simulated blood flow in twelve patient-specific anterior circulation unruptured intracranial aneurysms with daughter blebs utilizing the spectral/hp element framework Nektar++. We simulated three cardiac cycles using a volumetric flow rate waveform while we considered blood as a Newtonian fluid. To investigate the multidirectionality of the blood flow fields, besides the time-averaged wall shear stress (TAWSS), we calculated the oscillatory shear index (OSI), the relative residence time (RRT) and the time-averaged cross flow index (TACFI). Our CFD simulations suggest that in the majority of our vascular models there is a formation of complex intrasaccular flow patterns, resulting to low and highly oscillating WSS, especially in the area of the daughter blebs. The existence of disturbed multi-directional blood flow fields is also evident by the distributions of the RRT and the TACFI. These findings further support the theory that IAs with daughter blebs are linked to a potentially increased rupture risk.

Aneurysm Neck Overestimation has a Relatively Modest Impact on Simulated Hemodynamics

INTRODUCTION

Overestimation of intracranial aneurysm neck width by 3D angiography is a recognized clinical problem, and has long been a concern for image-based computational fluid dynamics (CFD). Recently, it was demonstrated that neck overestimation in 3D rotational angiography (3DRA) could be corrected via segmentation with upsampled resolution and gradient enhancement (SURGE). Our aim was to leverage this approach to determine whether and how neck overestimation actually impacts CFD-derived hemodynamics.

MATERIALS AND METHODS

A subset of 17 cases having the largest neck errors from a consecutive clinical sample of 60 was segmented from 3DRA using both standard watershed and SURGE methods. High-fidelity, pulsatile CFD was performed, and a variety of scalar hemodynamic parameters that have been associated with aneurysm growth and/or rupture status were derived.

RESULTS

With a few exceptions, flow and wall shear stress (WSS) patterns were qualitatively similar between neck-overestimated and corrected models. Sac-averaged WSS values were significantly lower after neck correction (p = 0.0005) but were highly correlated with their neck-overestimated counterparts (R² = 0.98). Jet impingement was significantly more concentrated in the neck-corrected vs. -uncorrected models (p = 0.0011), and only moderately correlated (R² = 0.61). Parameters quantifying velocity or WSS fluctuations were not significantly different after neck correction, but this reflected their poorer correlations (R² < 0.4). Nevertheless, for all hemodynamic parameters, median absolute differences were < 26%, and no parameter had more than 5/17 cases with absolute differences > 50%.

CONCLUSION

Differences in hemodynamics due to neck width overestimation were found to be at most equal to, and often less than, those reported for other sources of error/uncertainty in intracranial aneurysm CFD, such as solver settings or assumed inflow rates.

Reducing Geometric Uncertainty in Computational Hemodynamics by Deep Learning-Assisted Parallel-Chain MCMC

Computational hemodynamic modeling has been widely used in cardiovascular research and healthcare. However, the reliability of model predictions is largely dependent on the uncertainties of modeling parameters and boundary conditions, which should be carefully quantified and further reduced with available measurements. In this work, we focus on propagating and reducing the uncertainty of vascular geometries within a Bayesian framework. A novel deep learning (DL)-assisted parallel Markov chain Monte Carlo (MCMC) method is presented to enable efficient Bayesian posterior sampling and geometric uncertainty reduction. A DL model is built to approximate the geometry-to-hemodynamic map, which is trained actively using online data collected from parallel MCMC chains and utilized for early rejection of unlikely proposals to facilitate convergence with less expensive full-order model evaluations. Numerical studies on two-dimensional aortic flows are conducted to demonstrate the effectiveness and merit of the proposed method.

The Role of Hemodynamics through the Circle of Willis in the Development of Intracranial Aneurysm: A Systematic Review of Numerical Models

Background: The role of regional hemodynamics in the intracranial aneurysmal formation, growth, and rupture has been widely discussed based on numerical models over the past decades. Variation of the circle of Willis (CoW), which results in hemodynamic changes, is associated with the aneurysmal formation and rupture. However, such correlation has not been further clarified yet. The aim of this systematic review is to investigate whether simulated hemodynamic indices of the CoW are relevant to the formation, growth, or rupture of intracranial aneurysm. Methods: We conducted a review of MEDLINE, Web of Science, and EMBASE for studies on the correlation between hemodynamics indices of the CoW derived from numerical models and intracranial aneurysm up to December 2020 in compliance with PRISMA guidelines. Results: Three case reports out of 1046 publications met our inclusion and exclusion criteria, reporting 13 aneurysms in six patients. Eleven aneurysms were unruptured, and the state of the other two aneurysms was unknown. Wall shear stress, oscillatory shear index, von-Mises tension, flow velocity, and flow rate were reported as hemodynamic indices. Due to limited cases and significant heterogeneity between study settings, meta-analysis could not be performed. Conclusion: Numerical models can provide comprehensive information on the cerebral blood flow as well as local flow characteristics in the intracranial aneurysm. Based on only three case reports, no firm conclusion can be drawn regarding the correlation between hemodynamic parameters in the CoW derived from numerical models and aneurysmal formation or rupture. Due to the inherent nature of numerical models, more sensitive analysis and rigorous validations are required to determine its measurement error and thus extend their application into clinical practice for personalized management. Prospero registration number: CRD42021125169.

Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases

Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.

Characterizing Intracranial Hemodynamics in Sickle Cell Anemia: Impact of Patient-Specific Viscosity

PURPOSE

Pediatric and adult patients with sickle cell anemia (SCA) are at increased risk of stroke and cerebrovascular accident. In the general adult population, there is a relationship between arterial hemodynamics and pathology; however, this relationship in SCA patients remains to be elucidated. The aim of this work was to characterize circle of Willis hemodynamics in patients with SCA and quantify the impact of viscosity choice on pathophysiologically-relevant hemodynamics measures.

METHODS

Based on measured vascular geometries, time-varying flow rates, and blood parameters, detailed patient-specific simulations of the circle of Willis were conducted for SCA patients (n = 6). Simulations quantified the impact of patient-specific and standard blood viscosities on wall shear stress (WSS).

RESULTS

These results demonstrated that use of a standard blood viscosity introduces large errors into the estimation of pathophysiologically-relevant hemodynamic parameters. Standard viscosity models overpredicted peak WSS by 55% and 49% for steady and pulsatile flow, respectively. Moreover, these results demonstrated non-uniform, spatial patterns of positive and negative WSS errors related to viscosity, and standard viscosity simulations overpredicted the time-averaged WSS by 32% (standard deviation = 7.1%). Finally, differences in shear rate demonstrated that the viscosity choice alters the simulated near-wall flow field, impacting hemodynamics measures.

CONCLUSIONS

This work presents simulations of circle of Willis arterial flow in SCA patients and demonstrates the importance and feasibility of using a patient-specific viscosity in these simulations. Accurately characterizing cerebrovascular hemodynamics in SCA populations has potential for elucidating the pathophysiology of large-vessel occlusion, aneurysms, and tissue damage in these patients.

On the prevalence of flow instabilities from high-fidelity computational fluid dynamics of intracranial bifurcation aneurysms

High-fidelity computational fluid dynamics (HF-CFD) has revealed the potential for high-frequency flow instabilities (aka "turbulent-like" flow) in intracranial aneurysms, consistent with classic in vivo and in vitro reports of bruits and/or wall vibrations. However, HF-CFD has typically been performed on limited numbers of cases, often with unphysiological inflow conditions or focused on sidewall-type aneurysms where flow instabilities may be inherently less prevalent. Here we report HF-CFD of 50 bifurcation aneurysm cases from the open-source Aneurisk model repository. These were meshed using quadratic finite elements having an average effective spatial resolution of 0.065 mm, and solved under physiologically-pulsatile flow conditions using a well-validated, minimally-dissipative solver with 20,000 time-steps per cardiac cycle Flow instability was quantified using the recently introduced spectral power index (SPI), which quantifies, from 0 to 1, the power associated with velocity fluctuations above those of the driving inflow waveform. Of the 50 cases, nearly half showed regions within the sac having SPI up to 0.5, often with non-negligible power into the 100's of Hz, and roughly 1/3 had sac-averaged SPI > 0.1. High SPI did not significantly predict rupture status in this cohort. Proper orthogonal decomposition of cases with highest SPI_avg revealed time-varying energetics consistent with transient turbulence. Our reported prevalence of high-frequency flow instabilities in HF-CFD modelling of aneurysms suggests that care must be taken to avoid routinely overlooking them if we are to understand the highly dynamic mechanical forces to which some aneurysm walls may be exposed, and their prevalence in vivo.

Baseline local hemodynamics as predictor of lumen remodeling at 1-year follow-up in stented superficial femoral arteries

In-stent restenosis (ISR) is the major drawback of superficial femoral artery (SFA) stenting. Abnormal hemodynamics after stent implantation seems to promote the development of ISR. Accordingly, this study aims to investigate the impact of local hemodynamics on lumen remodeling in human stented SFA lesions. Ten SFA models were reconstructed at 1-week and 1-year follow-up from computed tomography images. Patient-specific computational fluid dynamics simulations were performed to relate the local hemodynamics at 1-week, expressed in terms of time-averaged wall shear stress (TAWSS), oscillatory shear index and relative residence time, with the lumen remodeling at 1-year, quantified as the change of lumen area between 1-week and 1-year. The TAWSS was negatively associated with the lumen area change (ρ = - 0.75, p = 0.013). The surface area exposed to low TAWSS was positively correlated with the lumen area change (ρ = 0.69, p = 0.026). No significant correlations were present between the other hemodynamic descriptors and lumen area change. The low TAWSS was the best predictive marker of lumen remodeling (positive predictive value of 44.8%). Moreover, stent length and overlapping were predictor of ISR at follow-up. Despite the limited number of analyzed lesions, the overall findings suggest an association between abnormal patterns of WSS after stenting and lumen remodeling.

Hemodynamics is associated with vessel wall remodeling in patients with middle cerebral artery stenosis

OBJECTIVES

To evaluate the relationship between hemodynamics and vessel wall remodeling patterns in middle cerebral artery (MCA) stenosis based on high-resolution magnetic resonance imaging and computational fluid dynamics (CFD).

METHODS

Forty consecutive patients with recent ischemic stroke or transient ischemic attack attributed to unilateral atherosclerotic MCA stenosis (50-99%) were prospectively recruited. All patients underwent a cross-sectional scan of the stenotic MCA vessel wall. The parameters of the vessel wall, the number of patients with acute infarction, translesional wall shear stress ratio (WSSR), wall shear stress in stenosis (WSSs), and translesional pressure ratio were obtained. The patients were divided into positive remodeling (PR) and negative remodeling (NR) groups. The differences in vessel wall parameters and hemodynamics were compared. Correlations between the parameters of the vessel wall and hemodynamics were calculated.

RESULTS

Of the 40 patients, 16 had PR, 19 had NR, and the other 5 displayed non-remodeling. The PR group had a smaller lumen area (p = 0.004), larger plaque area (p < 0.001), normal wall index (p = 0.004), and higher WSSR (p = 0.004) and WSSs (p = 0.023) at the most narrowed site. The PR group had more enhanced plaques (12 vs 6, p = 0.03). The number of patients with acute stroke in the PR group was more than that in the NR group (11 vs 4, p = 0.01). The remodeling index (r = 0.376, p = 0.026) and plaque area (r = 0.407, p = 0.015) showed a positive correlation with WSSR, respectively.

CONCLUSIONS

Hemodynamics plays a role in atherosclerotic plaques and vessel wall remodeling. Individuals with greater hemodynamic values might be more prone to stroke.

KEY POINTS

• Stenotic plaques in middle cerebral artery with positive remodeling have smaller lumen area and larger resp. higher plaque area, normal wall index, translesional wall shear stress ratio, and wall shear stress than negative remodeling. • The remodeling index and plaque area are positively correlated with translesional wall shear stress ratio. • Hemodynamic may help to understand the role of positive remodeling in the development of acute stroke.

Model Verification and Error Sensitivity of Turbulence-Related Tensor Characteristics in Pulsatile Blood Flow Simulations

Model verification, validation, and uncertainty quantification are essential procedures to estimate errors within cardiovascular flow modeling, where acceptable confidence levels are needed for clinical reliability. While more turbulent-like studies are frequently observed within the biofluid community, practical modeling guidelines are scarce. Verification procedures determine the agreement between the conceptual model and its numerical solution by comparing for example, discretization and phase-averaging-related errors of specific output parameters. This computational fluid dynamics (CFD) study presents a comprehensive and practical verification approach for pulsatile turbulent-like blood flow predictions by considering the amplitude and shape of the turbulence-related tensor field using anisotropic invariant mapping. These procedures were demonstrated by investigating the Reynolds stress tensor characteristics in a patient-specific aortic coarctation model, focusing on modeling-related errors associated with the spatiotemporal resolution and phase-averaging sampling size. Findings in this work suggest that attention should also be put on reducing phase-averaging related errors, as these could easily outweigh the errors associated with the spatiotemporal resolution when including too few cardiac cycles. Also, substantially more cycles are likely needed than typically reported for these flow regimes to sufficiently converge the phase-instant tensor characteristics. Here, higher degrees of active fluctuating directions, especially of lower amplitudes, appeared to be the most sensitive turbulence characteristics.

Modeling intracranial aneurysm stability and growth: an integrative mechanobiological framework for clinical cases

We present a novel patient-specific fluid-solid-growth framework to model the mechanobiological state of clinically detected intracranial aneurysms (IAs) and their evolution. The artery and IA sac are modeled as thick-walled, non-linear elastic fiber-reinforced composites. We represent the undulation distribution of collagen fibers: the adventitia of the healthy artery is modeled as a protective sheath whereas the aneurysm sac is modeled to bear load within physiological range of pressures. Initially, we assume the detected IA is stable and then consider two flow-related mechanisms to drive enlargement: (1) low wall shear stress; (2) dysfunctional endothelium which is associated with regions of high oscillatory flow. Localized collagen degradation and remodelling gives rise to formation of secondary blebs on the aneurysm dome. Restabilization of blebs is achieved by remodelling of the homeostatic collagen fiber stretch distribution. This integrative mechanobiological modelling workflow provides a step towards a personalized risk-assessment and treatment of clinically detected IAs.

Characterization of anisotropic turbulence behavior in pulsatile blood flow

Turbulent-like hemodynamics with prominent cycle-to-cycle flow variations have received increased attention as a potential stimulus for cardiovascular diseases. These turbulent conditions are typically evaluated in a statistical sense from single scalars extracted from ensemble-averaged tensors (such as the Reynolds stress tensor), limiting the amount of information that can be used for physical interpretations and quality assessments of numerical models. In this study, barycentric anisotropy invariant mapping was used to demonstrate an efficient and comprehensive approach to characterize turbulence-related tensor fields in patient-specific cardiovascular flows, obtained from scale-resolving large eddy simulations. These techniques were also used to analyze some common modeling compromises as well as MRI turbulence measurements through an idealized constriction. The proposed method found explicit sites of elevated turbulence anisotropy, including a broad but time-varying spectrum of characteristics over the flow deceleration phase, which was different for both the steady inflow and Reynolds-averaged Navier–Stokes modeling assumptions. Qualitatively, the MRI results showed overall expected post-stenotic turbulence characteristics, however, also with apparent regions of unrealizable or conceivably physically unrealistic conditions, including the highest turbulence intensity ranges. These findings suggest that more detailed studies of MRI-measured turbulence fields are needed, which hopefully can be assisted by more comprehensive evaluation tools such as the once described herein.

Computational Hemodynamic Modeling of Arterial Aneurysms: A Mini-Review

Arterial aneurysms are pathological dilations of blood vessels, which can be of clinical concern due to thrombosis, dissection, or rupture. Aneurysms can form throughout the arterial system, including intracranial, thoracic, abdominal, visceral, peripheral, or coronary arteries. Currently, aneurysm diameter and expansion rates are the most commonly used metrics to assess rupture risk. Surgical or endovascular interventions are clinical treatment options, but are invasive and associated with risk for the patient. For aneurysms in locations where thrombosis is the primary concern, diameter is also used to determine the level of therapeutic anticoagulation, a treatment that increases the possibility of internal bleeding. Since simple diameter is often insufficient to reliably determine rupture and thrombosis risk, computational hemodynamic simulations are being developed to help assess when an intervention is warranted. Created from subject-specific data, computational models have the potential to be used to predict growth, dissection, rupture, and thrombus-formation risk based on hemodynamic parameters, including wall shear stress, oscillatory shear index, residence time, and anomalous blood flow patterns. Generally, endothelial damage and flow stagnation within aneurysms can lead to coagulation, inflammation, and the release of proteases, which alter extracellular matrix composition, increasing risk of rupture. In this review, we highlight recent work that investigates aneurysm geometry, model parameter assumptions, and other specific considerations that influence computational aneurysm simulations. By highlighting modeling validation and verification approaches, we hope to inspire future computational efforts aimed at improving our understanding of aneurysm pathology and treatment risk stratification.

Cardiovascular models for personalised medicine: Where now and where next?

The aim of this position paper is to provide a brief overview of the current status of cardiovascular modelling and of the processes required and some of the challenges to be addressed to see wider exploitation in both personal health management and clinical practice. In most branches of engineering the concept of the digital twin, informed by extensive and continuous monitoring and coupled with robust data assimilation and simulation techniques, is gaining traction: the Gartner Group listed it as one of the top ten digital trends in 2018. The cardiovascular modelling community is starting to develop a much more systematic approach to the combination of physics, mathematics, control theory, artificial intelligence, machine learning, computer science and advanced engineering methodology, as well as working more closely with the clinical community to better understand and exploit physiological measurements, and indeed to develop jointly better measurement protocols informed by model-based understanding. Developments in physiological modelling, model personalisation, model outcome uncertainty, and the role of models in clinical decision support are addressed and 'where-next' steps and challenges discussed.

Hemodynamics of Cerebral Aneurysms: Connecting Medical Imaging and Biomechanical Analysis

In the last two decades, numerous studies have conducted patient-specific computations of blood flow dynamics in cerebral aneurysms and reported correlations between various hemodynamic metrics and aneurysmal disease progression or treatment outcomes. Nevertheless, intra-aneurysmal flow analysis has not been adopted in current clinical practice, and hemodynamic factors usually are not considered in clinical decision making. This review presents the state of the art in cerebral aneurysm imaging and image-based modeling, discussing the advantages and limitations of each approach and focusing on the translational value of hemodynamic analysis. Combining imaging and modeling data obtained from different flow modalities can improve the accuracy and fidelity of resulting velocity fields and flow-derived factors that are thought to affect aneurysmal disease progression. It is expected that predictive models utilizing hemodynamic factors in combination with patient medical history and morphological data will outperform current risk scores and treatment guidelines. Possible future directions include novel approaches enabling data assimilation and multimodality analysis of cerebral aneurysm hemodynamics.

4D-CT angiography versus 3D-rotational angiography as the imaging modality for computational fluid dynamics of cerebral aneurysms

Background and purpose

Computational fluid dynamics (CFD) can provide valuable information regarding intracranial hemodynamics. Patient-specific models can be segmented from various imaging modalities, which may influence the geometric output and thus hemodynamic results. This study aims to compare CFD results from aneurysm models segmented from three-dimensional rotational angiography (3D-RA) versus novel four-dimensional CT angiography (4D-CTA).

Methods

Fourteen patients with 16 cerebral aneurysms underwent novel 4D-CTA followed by 3D-RA. Endoluminal geometries were segmented from each modality using an identical workflow, blinded to the other modality, to produce 28 'original' models. Each was then minimally edited a second time to match length of branches, producing 28 additional 'matched' models. CFD simulations were performed using estimated flow rates for 'original' models (representing real-world experience) and patient-specific flow rates from 4D-CTA for 'matched' models (to control for influence of modality alone).

Results

Overall, geometric and hemodynamic results were consistent between models segmented from 3D-RA and 4D-CTA, with correlations improving after matching to control for operator-introduced variability. Despite smaller 4D-CTA parent artery diameters (3.49±0.97 mm vs 3.78±0.92 mm for 3D-RA; p=0.005) and sac volumes (157 (37–750 mm3) vs 173 (53–770 mm3) for 3D-RA; p=0.0002), sac averages of time-averaged wall shear stress (TAWSS), oscillatory shear (OSI), and high frequency fluctuations (measured by spectral power index, SPI) were well correlated between 3D-RA and 4D-CTA 'matched' control models (TAWSS, R2=0.91; OSI, R2=0.79; SPI, R2=0.90).

Conclusions

Our study shows that CFD performed using 4D-CTA models produces reliable geometric and hemodynamic information in the intracranial circulation. 4D-CTA may be considered as a follow-up imaging tool for hemodynamic assessment of cerebral aneurysms.