1
|
Straccia A, Barbour MC, Chassagne F, Bass D, Barros G, Leotta D, Sheehan F, Sharma D, Levitt MR, Aliseda A. Numerical Modeling of Flow in the Cerebral Vasculature: Understanding Changes in Collateral Flow Directions in the Circle of Willis for a Cohort of Vasospasm Patients Through Image-Based Computational Fluid Dynamics. Ann Biomed Eng 2024:10.1007/s10439-024-03533-w. [PMID: 38758460 DOI: 10.1007/s10439-024-03533-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2024] [Accepted: 04/30/2024] [Indexed: 05/18/2024]
Abstract
The Circle of Willis (CoW) is a ring-like network of blood vessels that perfuses the brain. Flow in the collateral pathways that connect major arterial inputs in the CoW change dynamically in response to vessel narrowing or occlusion. Vasospasm is an involuntary constriction of blood vessels following subarachnoid hemorrhage (SAH), which can lead to stroke. This study investigated interactions between localization of vasospasm in the CoW, vasospasm severity, anatomical variations, and changes in collateral flow directions. Patient-specific computational fluid dynamics (CFD) simulations were created for 25 vasospasm patients. Computed tomographic angiography scans were segmented capturing the anatomical variation and stenosis due to vasospasm. Transcranial Doppler ultrasound measurements of velocity were used to define boundary conditions. Digital subtraction angiography was analyzed to determine the directions and magnitudes of collateral flows as well as vasospasm severity in each vessel. Percent changes in resistance and viscous dissipation were analyzed to quantify vasospasm severity and localization of vasospasm in a specific region of the CoW. Angiographic severity correlated well with percent changes in resistance and viscous dissipation across all cerebral vessels. Changes in flow direction were observed in collateral pathways of some patients with localized vasospasm, while no significant changes in flow direction were observed in others. CFD simulations can be leveraged to quantify the localization and severity of vasospasm in SAH patients. These factors as well as anatomical variation may lead to changes in collateral flow directions. Future work could relate localization and vasospasm severity to clinical outcomes like the development of infarct.
Collapse
Affiliation(s)
- Angela Straccia
- Department of Mechanical Engineering, University of Washington, Seattle, WA, USA.
| | - Michael C Barbour
- Department of Mechanical Engineering, University of Washington, Seattle, WA, USA
| | | | - David Bass
- Department of Neurological Surgery, University of Washington, Seattle, WA, USA
| | - Guilherme Barros
- Department of Neurological Surgery, University of Washington, Seattle, WA, USA
| | - Daniel Leotta
- Applied Physics Laboratory, University of Washington, Seattle, WA, USA
| | - Florence Sheehan
- Department of Medicine, University of Washington, Seattle, WA, USA
| | - Deepak Sharma
- Department of Anesthesiology and Pain Medicine, University of Washington, Seattle, WA, USA
| | - Michael R Levitt
- Department of Neurological Surgery, University of Washington, Seattle, WA, USA
| | - Alberto Aliseda
- Department of Mechanical Engineering, University of Washington, Seattle, WA, USA
| |
Collapse
|
2
|
Qiao X, Zhang R, Yu J, Yan Y, Bouakaz A, Su X, Liu J, Zong Y, Wan M. Noninvasive assessment of pressure distribution and fractional flow in middle cerebral artery using microbubbles and plane wave in vitro. ULTRASONICS 2024; 138:107244. [PMID: 38237398 DOI: 10.1016/j.ultras.2024.107244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Revised: 12/20/2023] [Accepted: 01/09/2024] [Indexed: 04/02/2024]
Abstract
Fractional flow has been proposed for quantifying the degree of functional stenosis in cerebral arteries. Herein, subharmonic aided pressure estimation (SHAPE) combined with plane wave (PW) transmission was employed to noninvasively estimate the pressure distribution and fractional flow in the middle cerebral artery (MCA) in vitro. Consequently, the effects of incident sound pressure (peak negative pressures of 86-653 kPa), pulse repetition frequency (PRF), number of pulses, and blood flow rate on the subharmonic pressure relationship were investigated. The radio frequency data were stored and beamformed offline, and the subharmonic amplitude over a 0.4 MHz bandwidth was extracted using a 12-cycle PW at 4 MHz. The optimal incident sound pressure was 217 kPa without skull (sensitivity = 0.09 dB/mmHg; r2 = 0.997) and 410 kPa with skull (median sensitivity = 0.06 dB/mmHg; median r2 = 0.981). The optimal PRF was 500 Hz, as this value affords the highest sensitivity (0.09 dB/mmHg; r2 = 0.976) and temporal resolution. In addition, the blood flow rate exhibited a lesser effect on the subharmonic pressure relationship in our experimental setup. Using the optimized parameters, the blood pressure distribution and fractional flow (FFs) were measured. As such, the FFs value was in high agreement with the value measured using the pressure sensor (FFm). The mean ± standard deviations of the FF difference (FFm - FFs) were 0.03 ± 0.06 without skull and 0.01 ± 0.05 with skull.
Collapse
Affiliation(s)
- Xiaoyang Qiao
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China
| | - Ruiyan Zhang
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China
| | - Jianjun Yu
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China
| | - Yadi Yan
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China
| | - Ayache Bouakaz
- UMR 1253, iBrain, Université de Tours, Inserm, Tours, France
| | - Xiao Su
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China
| | - Jiacheng Liu
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China
| | - Yujin Zong
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China.
| | - Mingxi Wan
- Xi'an Jiaotong University, College of Life Science and Technology, Xi'an, China.
| |
Collapse
|
3
|
Moore AE, Sekhar LN, Beach KW. Ophthalmic artery reversal predicts contralateral body weakness symptoms better than carotid Doppler velocity-preliminary results. ANNALS OF TRANSLATIONAL MEDICINE 2023; 11:408. [PMID: 38213804 PMCID: PMC10777211 DOI: 10.21037/atm-23-1681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/04/2023] [Accepted: 10/11/2023] [Indexed: 01/13/2024]
Abstract
Background The measurement of blood velocity in the carotid artery has been the most popular noninvasive method of identifying and classifying carotid stenosis for half a century. Carotid stenosis is an indicator of elevated risk of stroke; anatomic revascularization reduces the chance of stroke by more than half. Controversy persists on how patients with severe carotid stenosis should be selected for anatomic revascularization. Patients with a connected circle of Willis (coW) might not benefit from anatomic revascularization; patients with two segments missing in the coW are most likely to benefit from revascularization. Methods Based on this analysis of data from carotid duplex examinations and transcranial Doppler examinations including ophthalmic artery (OA) direction in 28 patients, a refined carotid examination protocol is proposed. This refinement includes Doppler measurement of OA flow direction and documentation of internal carotid artery (ICA) bruit in addition to the adoption of an ICA peak systolic velocity (PSV) criterion exceeding 350 cm/s for identification of the patient most likely to benefit from carotid stenosis treatment. Results Sensitivity and specificity of OA direction or carotid bruit are 84.6%±5.4%, 71.4%±2.1% and for PSV >350 cm/s are 84.6%±5.4%, 59.5%±2.3% for predicting contralateral body weakness. Conclusions The proposed examination can be performed with the same duplex scanner and scan head currently used for carotid examinations with little additional time.
Collapse
Affiliation(s)
- Anne E. Moore
- Harborview Medical Center, Cerebrovascular Laboratory, University of Washington, Seattle, WA, USA
| | - Laligam N. Sekhar
- Department of Neurological Surgery, University of Washington School of Medicine, Harborview Medical Center, Seattle, WA, USA
| | - Kirk W. Beach
- Department of Surgery, University of Washington, Seattle, WA, USA
- Department of Bioengineering, University of Washington, Seattle, WA, USA
| |
Collapse
|
4
|
Kizhisseri M, Gharaie S, Boopathy SR, Lim RP, Mohammadzadeh M, Schluter J. Differential sensitivities to blood pressure variations in internal carotid and intracranial arteries: a numerical approach to stroke prediction. Sci Rep 2023; 13:22319. [PMID: 38102319 PMCID: PMC10724219 DOI: 10.1038/s41598-023-49591-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Accepted: 12/09/2023] [Indexed: 12/17/2023] Open
Abstract
Stroke remains a global health concern, necessitating early prediction for effective management. Atherosclerosis-induced internal carotid and intra cranial stenosis contributes significantly to stroke risk. This study explores the relationship between blood pressure and stroke prediction, focusing on internal carotid artery (ICA) branches: middle cerebral artery (MCA), anterior cerebral artery (ACA), and their role in hemodynamics. Computational fluid dynamics (CFD) informed by the Windkessel model were employed to simulate patient-specific ICA models with introduced stenosis. Central to our investigation is the impact of stenosis on blood pressure, flow velocity, and flow rate across these branches, incorporating Fractional Flow Reserve (FFR) analysis. Results highlight differential sensitivities to blood pressure variations, with M1 branch showing high sensitivity, ACA moderate, and M2 minimal. Comparing blood pressure fluctuations between ICA and MCA revealed heightened sensitivity to potential reverse flow compared to ICA and ACA comparisons, emphasizing MCA's role. Blood flow adjustments due to stenosis demonstrated intricate compensatory mechanisms. FFR emerged as a robust predictor of stenosis severity, particularly in the M2 branch. In conclusion, this study provides comprehensive insights into hemodynamic complexities within major intracranial arteries, elucidating the significance of blood pressure variations, flow attributes, and FFR in stenosis contexts. Subject-specific data integration enhances model reliability, aiding stroke risk assessment and advancing cerebrovascular disease understanding.
Collapse
Affiliation(s)
- Muhsin Kizhisseri
- School of Engineering, Deakin University, 75 Pigdons Rd, Waurn Ponds, VIC, 3216, Australia
| | - Saleh Gharaie
- School of Engineering, Deakin University, 75 Pigdons Rd, Waurn Ponds, VIC, 3216, Australia.
| | | | | | | | - Jorg Schluter
- School of Engineering, Deakin University, 75 Pigdons Rd, Waurn Ponds, VIC, 3216, Australia
| |
Collapse
|
5
|
Xia Y, Ravikumar N, Lassila T, Frangi AF. Virtual high-resolution MR angiography from non-angiographic multi-contrast MRIs: synthetic vascular model populations for in-silico trials. Med Image Anal 2023; 87:102814. [PMID: 37196537 DOI: 10.1016/j.media.2023.102814] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Revised: 04/04/2023] [Accepted: 04/08/2023] [Indexed: 05/19/2023]
Abstract
Despite success on multi-contrast MR image synthesis, generating specific modalities remains challenging. Those include Magnetic Resonance Angiography (MRA) that highlights details of vascular anatomy using specialised imaging sequences for emphasising inflow effect. This work proposes an end-to-end generative adversarial network that can synthesise anatomically plausible, high-resolution 3D MRA images using commonly acquired multi-contrast MR images (e.g. T1/T2/PD-weighted MR images) for the same subject whilst preserving the continuity of vascular anatomy. A reliable technique for MRA synthesis would unleash the research potential of very few population databases with imaging modalities (such as MRA) that enable quantitative characterisation of whole-brain vasculature. Our work is motivated by the need to generate digital twins and virtual patients of cerebrovascular anatomy for in-silico studies and/or in-silico trials. We propose a dedicated generator and discriminator that leverage the shared and complementary features of multi-source images. We design a composite loss function for emphasising vascular properties by minimising the statistical difference between the feature representations of the target images and the synthesised outputs in both 3D volumetric and 2D projection domains. Experimental results show that the proposed method can synthesise high-quality MRA images and outperform the state-of-the-art generative models both qualitatively and quantitatively. The importance assessment reveals that T2 and PD-weighted images are better predictors of MRA images than T1; and PD-weighted images contribute to better visibility of small vessel branches towards the peripheral regions. In addition, the proposed approach can generalise to unseen data acquired at different imaging centres with different scanners, whilst synthesising MRAs and vascular geometries that maintain vessel continuity. The results show the potential for use of the proposed approach to generating digital twin cohorts of cerebrovascular anatomy at scale from structural MR images typically acquired in population imaging initiatives.
Collapse
Affiliation(s)
- Yan Xia
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK.
| | - Nishant Ravikumar
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK
| | - Toni Lassila
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK
| | - Alejandro F Frangi
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK; Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK; Medical Imaging Research Center (MIRC), Cardiovascular Science and Electronic Engineering Departments, KU Leuven, Leuven, Belgium; Alan Turing Institute, London, UK
| |
Collapse
|
6
|
Zhou J, Li J, Qin S, Liu J, Lin Z, Xie J, Zhang Z, Chen R. High-resolution cerebral blood flow simulation with a domain decomposition method and verified by the TCD measurement. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2022; 224:107004. [PMID: 35841853 DOI: 10.1016/j.cmpb.2022.107004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 06/24/2022] [Accepted: 07/03/2022] [Indexed: 06/15/2023]
Abstract
BACKGROUND An efficient and accurate blood flow simulation can be useful for understanding many vascular diseases. Accurately resolving the blood flow velocity based on patient-specific geometries and model parameters is still a major challenge because of complex geomerty and turbulence issues. In addition, obtaining results in a short amount of computing time is important so that the simulation can be used in the clinical environment. In this work, we present a parallel scalable method for the patient-specific blood flow simulation with focuses on its parallel performance study and clinical verification. METHODS We adopt a fully implicit unstructured finite element method for a patient-specific simulation of blood flow in a full precerebral artery. The 3D artery is constructed from MRI images, and a parallel Newton-Krylov method preconditioned with a two-level domain decomposition method is adopted to solve the large nonlinear system discretized from the time-dependent 3D Navier-Stokes equations in the artery with an integral outlet boundary condition. The simulated results are verified using the clinical data measured by transcranial Doppler ultrasound, and the parallel performance of the algorithm is studied on a supercomputer. RESULTS The simulated velocity matches the clinical measured data well. Other simulated blood flow parameters, such as pressure and wall shear stress, are within reasonable ranges. The results show that the parallel algorithm scales up to 2160 processors with a 49% parallel efficiency for solving a problem with over 20 million unstructured elements on a supercomputer. For a standard cerebral blood flow simulation case with approximately 4 million finite elements, the calculation of one cardiac cycle can be finished within one hour with 1000 processors. CONCLUSION The proposed method is able to perform high-resolution 3D blood flow simulations in a patient-specific full precerebral artery within an acceptable time, and the simulated results are comparable with the clinical measured data, which demonstrates its high potential for clinical applications.
Collapse
Affiliation(s)
- Jie Zhou
- School of Mathematics and Statistics, Changsha University of Science and Technology, Changsha, China
| | - Jing Li
- School of Mathematics and Statistics, Changsha University of Science and Technology, Changsha, China
| | - Shanlin Qin
- Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.
| | - Jia Liu
- Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Zeng Lin
- Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Jian Xie
- Department of Neurology, Affiliated ZhongDa Hospital, School of Medicine, Southeast University, Nanjing, China
| | - Zhijun Zhang
- Department of Neurology, Affiliated ZhongDa Hospital, School of Medicine, Southeast University, Nanjing, China
| | - Rongliang Chen
- Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China; Shenzhen Key Laboratory for Exascale Engineering and Scientific Computing, Shenzhen, China.
| |
Collapse
|
7
|
Patient-specific brain arteries molded as a flexible phantom model using 3D printed water-soluble resin. Sci Rep 2022; 12:10172. [PMID: 35715506 PMCID: PMC9205921 DOI: 10.1038/s41598-022-14279-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Accepted: 06/03/2022] [Indexed: 11/08/2022] Open
Abstract
Visualizing medical images from patients as physical 3D models (phantom models) have many roles in the medical field, from education to preclinical preparation and clinical research. However, current phantom models are generally generic, expensive, and time-consuming to fabricate. Thus, there is a need for a cost- and time-efficient pipeline from medical imaging to patient-specific phantom models. In this work, we present a method for creating complex 3D sacrificial molds using an off-the-shelf water-soluble resin and a low-cost desktop 3D printer. This enables us to recreate parts of the cerebral arterial tree as a full-scale phantom model ([Formula: see text] cm) in transparent silicone rubber (polydimethylsiloxane, PDMS) from computed tomography angiography images (CTA). We analyzed the model with magnetic resonance imaging (MRI) and compared it with the patient data. The results show good agreement and smooth surfaces for the arteries. We also evaluate our method by looking at its capability to reproduce 1 mm channels and sharp corners. We found that round shapes are well reproduced, whereas sharp features show some divergence. Our method can fabricate a patient-specific phantom model with less than 2 h of total labor time and at a low fabrication cost.
Collapse
|
8
|
Benemerito I, Narata AP, Narracott A, Marzo A. Determining Clinically-Viable Biomarkers for Ischaemic Stroke Through a Mechanistic and Machine Learning Approach. Ann Biomed Eng 2022; 50:740-750. [PMID: 35364704 PMCID: PMC9079032 DOI: 10.1007/s10439-022-02956-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Accepted: 03/20/2022] [Indexed: 11/29/2022]
Abstract
Assessment of distal cerebral perfusion after ischaemic stroke is currently only possible through expensive and time-consuming imaging procedures which require the injection of a contrast medium. Alternative approaches that could indicate earlier the impact of blood flow occlusion on distal cerebral perfusion are currently lacking. The aim of this study was to identify novel biomarkers suitable for clinical implementation using less invasive diagnostic techniques such as Transcranial Doppler (TCD). We used 1D modelling to simulate pre- and post-stroke velocity and flow wave propagation in a typical arterial network, and Sobol’s sensitivity analysis, supported by the use of Gaussian process emulators, to identify biomarkers linked to cerebral perfusion. We showed that values of pulsatility index of the right anterior cerebral artery > 1.6 are associated with poor perfusion and may require immediate intervention. Three additional biomarkers with similar behaviour, all related to pulsatility indices, were identified. These results suggest that flow pulsatility measured at specific locations could be used to effectively estimate distal cerebral perfusion rates, and ultimately improve clinical diagnosis and management of ischaemic stroke.
Collapse
Affiliation(s)
- Ivan Benemerito
- INSIGNEO Institute for In Silico Medicine, The University of Sheffield, Sheffield, UK. .,Department of Mechanical Engineering, The University of Sheffield, Sheffield, UK.
| | - Ana Paula Narata
- Department of Neuroradiology, University Hospital of Southampton, Southampton, UK
| | - Andrew Narracott
- INSIGNEO Institute for In Silico Medicine, The University of Sheffield, Sheffield, UK.,Department of Infection, Immunity and Cardiovascular Disease, The University of Sheffield, Sheffield, UK
| | - Alberto Marzo
- INSIGNEO Institute for In Silico Medicine, The University of Sheffield, Sheffield, UK.,Department of Mechanical Engineering, The University of Sheffield, Sheffield, UK
| |
Collapse
|
9
|
Abstract
Alterations in cerebral blood flow are common in several neurological diseases among the elderly including stroke, cerebral small vessel disease, vascular dementia, and Alzheimer's disease. 4D flow magnetic resonance imaging (MRI) is a relatively new technique to investigate cerebrovascular disease, and makes it possible to obtain time-resolved blood flow measurements of the entire cerebral arterial venous vasculature and can be used to derive a repertoire of hemodynamic biomarkers indicative of cerebrovascular health. The information that can be obtained from one single 4D flow MRI scan allows both the investigation of aberrant flow patterns at a focal location in the vasculature as well as estimations of brain-wide disturbances in blood flow. Such focal and global hemodynamic biomarkers show the potential of being sensitive to impending cerebrovascular disease and disease progression and can also become useful during planning and follow-up of interventions aiming to restore a normal cerebral circulation. Here, we describe 4D flow MRI approaches for analyzing the cerebral vasculature. We then survey key hemodynamic biomarkers that can be reliably assessed using the technique. Finally, we highlight cerebrovascular diseases where one or multiple hemodynamic biomarkers are of central interest.
Collapse
Affiliation(s)
- Anders Wåhlin
- Department of Radiation Sciences, Umeå University, Umeå, Sweden.,Department of Applied Physics and Electronics, Umeå University, Umeå, Sweden.,Umeå Center for Functional Brain Imaging (UFBI), Umeå University, Umeå, Sweden
| | - Anders Eklund
- Department of Radiation Sciences, Umeå University, Umeå, Sweden.,Umeå Center for Functional Brain Imaging (UFBI), Umeå University, Umeå, Sweden
| | - Jan Malm
- Department of Clinical Science and Neurosciences, Umeå University, Umeå, Sweden
| |
Collapse
|