1
|
Maul N, Birkhold A, Wagner F, Thies M, Rohleder M, Berg P, Kowarschik M, Maier A. Simulation-informed learning for time-resolved angiographic contrast agent concentration reconstruction. Comput Biol Med 2024; 182:109178. [PMID: 39321585 DOI: 10.1016/j.compbiomed.2024.109178] [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: 12/01/2023] [Revised: 07/29/2024] [Accepted: 09/18/2024] [Indexed: 09/27/2024]
Abstract
Three-dimensional Digital Subtraction Angiography (3D-DSA) is a well-established X-ray-based technique for visualizing vascular anatomy. Recently, four-dimensional DSA (4D-DSA) reconstruction algorithms have been developed to enable the visualization of volumetric contrast flow dynamics through time-series of volumes. This reconstruction problem is ill-posed mainly due to vessel overlap in the projection direction and geometric vessel foreshortening, which leads to information loss in the recorded projection images. However, knowledge about the underlying fluid dynamics can be leveraged to constrain the solution space. In our work, we implicitly include this information in a neural network-based model that is trained on a dataset of image-based blood flow simulations. The model predicts the spatially averaged contrast agent concentration for each centerline point of the vasculature over time, lowering the overall computational demand. The trained network enables the reconstruction of relative contrast agent concentrations with a mean absolute error of 0.02±0.02 and a mean absolute percentage error of 5.31±9.25 %. Moreover, the network is robust to varying degrees of vessel overlap and vessel foreshortening. Our approach demonstrates the potential of the integration of machine learning and blood flow simulations in time-resolved angiographic contrast agent concentration reconstruction.
Collapse
Affiliation(s)
- Noah Maul
- Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Germany; Siemens Healthineers AG, Forchheim, Germany.
| | | | - Fabian Wagner
- Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Germany
| | - Mareike Thies
- Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Germany
| | - Maximilian Rohleder
- Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Germany
| | - Philipp Berg
- Research Campus STIMULATE, University of Magdeburg, Germany; Department of Medical Engineering, University of Magdeburg, Germany
| | | | - Andreas Maier
- Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Germany
| |
Collapse
|
2
|
Wéber R, Gyürki D, Paál G. First blood: An efficient, hybrid one- and zero-dimensional, modular hemodynamic solver. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3701. [PMID: 36948891 DOI: 10.1002/cnm.3701] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 01/24/2023] [Accepted: 03/11/2023] [Indexed: 05/13/2023]
Abstract
Low-dimensional (1D or 0D) models can describe the whole human blood circulation, for example, 1D distributed parameter model for the arterial network and 0D concentrated models for the heart or other organs. This paper presents a combined 1D-0D solver, called first_blood, that solves the governing equations of fluid dynamics to model low-dimensional hemodynamic effects. An extended method of characteristics is applied here to solve the momentum, and mass conservation equations and the viscoelastic wall model equation, mimicking the material properties of arterial walls. The heart and the peripheral lumped models are solved with a general zero-dimensional (0D) nonlinear solver. The model topology can be modular, that is, first_blood can solve any 1D-0D hemodynamic model. To demonstrate the applicability of first_blood, the human arterial system, the heart and the peripherals are modelled using the solver. The simulation time of a heartbeat takes around 2 s, that is, first_blood requires only twice the real-time for the simulation using an average PC, which highlights the computational efficiency. The source code is available on GitHub, that is, it is open source. The model parameters are based on the literature suggestions and on the validation of output data to obtain physiologically relevant results.
Collapse
Affiliation(s)
- Richárd Wéber
- Department of Hydrodynamic Systems, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
| | - Dániel Gyürki
- Department of Hydrodynamic Systems, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
| | - György Paál
- Department of Hydrodynamic Systems, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
| |
Collapse
|
3
|
Fonken J, Maas E, Nievergeld A, van Sambeek M, van de Vosse F, Lopata R. The Impact of a Limited Field-of-View on Computed Hemodynamics in Abdominal Aortic Aneurysms: Evaluating the Feasibility of Completing Ultrasound Segmentations with Parametric Geometries. Ann Biomed Eng 2023; 51:1296-1309. [PMID: 36709232 PMCID: PMC10172266 DOI: 10.1007/s10439-022-03133-6] [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: 09/21/2022] [Accepted: 12/25/2022] [Indexed: 01/30/2023]
Abstract
To improve abdominal aortic aneurysm (AAA) rupture risk assessment, a large, longitudinal study on AAA hemodynamics and biomechanics is necessary, using personalized fluid-structure interaction (FSI) modeling. 3-dimensional, time-resolved ultrasound (3D+t US) is the preferred image modality to obtain the patient-specific AAA geometry for such a study, since it is safe, affordable and provides temporal information. However, the 3D+t US field-of-view (FOV) is limited and therefore often fails to capture the inlet and aorto-iliac bifurcation geometry. In this study, a framework was developed to add parametric inlet and bifurcation geometries to the abdominal aortic aneurysm geometry by employing dataset statistics and parameters of the AAA geometry. The impact of replacing the patient-specific inlet and bifurcation geometries, acquired using computed tomography (CT) scans, by parametric geometries was evaluated by examining the differences in hemodynamics (systolic and time-averaged wall shear stress and oscillatory shear index) in the aneurysm region. The results show that the inlet geometry has a larger effect on the AAA hemodynamics (median differences of 7.5 to 18.8%) than the bifurcation geometry (median differences all below 1%). Therefore, it is not feasible to replace the patient-specific inlet geometry by a generic one. Future studies should investigate the possibilities of extending the proximal FOV of 3D+t US. However, this study did show the feasibility of adding a parametric bifurcation geometry to the aneurysm geometry. After extending the proximal FOV, the obtained framework can be used to extract AAA geometries from 3D+t US for FSI simulations, despite the absence of the bifurcation geometry.
Collapse
Affiliation(s)
- Judith Fonken
- Photoacoustics & Ultrasound Laboratory Eindhoven (PULS/e), Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands. .,Department of Surgery, Catharina Hospital Eindhoven, Eindhoven, The Netherlands.
| | - Esther Maas
- Photoacoustics & Ultrasound Laboratory Eindhoven (PULS/e), Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands.,Department of Surgery, Catharina Hospital Eindhoven, Eindhoven, The Netherlands
| | - Arjet Nievergeld
- Photoacoustics & Ultrasound Laboratory Eindhoven (PULS/e), Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands.,Department of Surgery, Catharina Hospital Eindhoven, Eindhoven, The Netherlands
| | - Marc van Sambeek
- Photoacoustics & Ultrasound Laboratory Eindhoven (PULS/e), Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands.,Department of Surgery, Catharina Hospital Eindhoven, Eindhoven, The Netherlands
| | - Frans van de Vosse
- Cardiovascular Biomechanics, Department of Biomechanical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Richard Lopata
- Photoacoustics & Ultrasound Laboratory Eindhoven (PULS/e), Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| |
Collapse
|
4
|
Mathematical Models for Blood Flow Quantification in Dialysis Access Using Angiography: A Comparative Study. Diagnostics (Basel) 2021; 11:diagnostics11101771. [PMID: 34679469 PMCID: PMC8534972 DOI: 10.3390/diagnostics11101771] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2021] [Revised: 09/10/2021] [Accepted: 09/17/2021] [Indexed: 11/26/2022] Open
Abstract
Blood flow rate in dialysis (vascular) access is the key parameter to examine patency and to evaluate the outcomes of various endovascular interve7ntions. While angiography is extensively used for dialysis access–salvage procedures, to date, there is no image-based blood flow measurement application commercially available in the angiography suite. We aim to calculate the blood flow rate in the dialysis access based on cine-angiographic and fluoroscopic image sequences. In this study, we discuss image-based methods to quantify access blood flow in a flow phantom model. Digital subtraction angiography (DSA) and fluoroscopy were used to acquire images at various sampling rates (DSA—3 and 6 frames/s, fluoroscopy—4 and 10 pulses/s). Flow rates were computed based on two bolus tracking algorithms, peak-to-peak and cross-correlation, and modeled with three curve-fitting functions, gamma variate, lagged normal, and polynomial, to correct errors with transit time measurement. Dye propagation distance and the cross-sectional area were calculated by analyzing the contrast enhancement in the vessel. The calculated flow rates were correlated versus an in-line flow sensor measurement. The cross-correlation algorithm with gamma-variate curve fitting had the best accuracy and least variability in both imaging modes. The absolute percent error (mean ± SEM) of flow quantification in the DSA mode at 6 frames/s was 21.4 ± 1.9%, and in the fluoroscopic mode at 10 pulses/s was 37.4 ± 3.6%. The radiation dose varied linearly with the sampling rate in both imaging modes and was substantially low to invoke any tissue reactions or stochastic effects. The cross-correlation algorithm and gamma-variate curve fitting for DSA acquisition at 6 frames/s had the best correlation with the flow sensor measurements. These findings will be helpful to develop a software-based vascular access flow measurement tool for the angiography suite and to optimize the imaging protocol amenable for computational flow applications.
Collapse
|
5
|
Heusinkveld MHG, Holtackers RJ, Adriaans BP, Op't Roodt J, Arts T, Delhaas T, Reesink KD, Huberts W. Complementing sparse vascular imaging data by physiological adaptation rules. J Appl Physiol (1985) 2021; 130:571-588. [PMID: 33119465 DOI: 10.1152/japplphysiol.00250.2019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Mathematical modeling of pressure and flow waveforms in blood vessels using pulse wave propagation (PWP) models has tremendous potential to support clinical decision making. For a personalized model outcome, measurements of all modeled vessel radii and wall thicknesses are required. In clinical practice, however, data sets are often incomplete. To overcome this problem, we hypothesized that the adaptive capacity of vessels in response to mechanical load could be utilized to fill in the gaps of incomplete patient-specific data sets. We implemented homeostatic feedback loops in a validated PWP model to allow adaptation of vessel geometry to maintain physiological values of wall stress and wall shear stress. To evaluate our approach, we gathered vascular MRI and ultrasound data sets of wall thicknesses and radii of central and arm arterial segments of 10 healthy subjects. Reference models (i.e., termed RefModel, n = 10) were simulated using complete data, whereas adapted models (AdaptModel, n = 10) used data of one carotid artery segment only, and the remaining geometries in this model were estimated using adaptation. We evaluated agreement between RefModel and AdaptModel geometries, as well as that between pressure and flow waveforms of both models. Limits of agreement (bias ± 2 SD of difference) between AdaptModel and RefModel radii and wall thicknesses were 0.2 ± 2.6 mm and -140 ± 557 µm, respectively. Pressure and flow waveform characteristics of the AdaptModel better resembled those of the RefModels as compared with the model in which the vessels were not adapted. Our adaptation-based PWP model enables personalization of vascular geometries even when not all required data are available.NEW & NOTEWORTHY To benefit personalized pulse wave propagation (PWP) modeling, we propose a novel method that, instead of relying on extensive data sets on vascular geometries, incorporates physiological adaptation rules. The developed vascular adaptation model adequately predicted arterial radius and wall thickness compared with ultrasound and MRI estimates, obtained in humans. Our approach could be used as a tool to facilitate personalized modeling, notably in case of missing data, as routinely found in clinical settings.
Collapse
Affiliation(s)
| | - Robert J Holtackers
- Department of Radiology and Nuclear Medicine, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Bouke P Adriaans
- Department of Cardiology, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Jos Op't Roodt
- Department of Internal Medicine, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Theo Arts
- Department of Biomedical Engineering, Maastricht University, Maastricht, The Netherlands
| | - Tammo Delhaas
- Department of Biomedical Engineering, Maastricht University, Maastricht, The Netherlands
| | - Koen D Reesink
- Department of Biomedical Engineering, Maastricht University, Maastricht, The Netherlands
| | - Wouter Huberts
- Department of Biomedical Engineering, Maastricht University, Maastricht, The Netherlands.,Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| |
Collapse
|
6
|
Stevens RRF, Donders WP, Quicken S, van de Vosse FN, Mess WH, Delhaas T, Huberts W. The Role of One-Dimensional Model-Generated Inter-Subject Variations in Systemic Properties on Wall Shear Stress Indices of Intracranial Aneurysms. IEEE Trans Biomed Eng 2019; 67:1030-1039. [PMID: 31329544 DOI: 10.1109/tbme.2019.2928416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Variations in systemic properties of the arterial tree, such as aging-induced vessel stiffness, can alter the shape of pressure and flow waveforms. As a consequence, the hemodynamics around a cerebral aneurysm change, and therefore, also the corresponding in- and outlet boundary conditions (BCs) used for three-dimensional (3D) calculations of hemodynamic indices. In this study, we investigate the effects of variations in systemic properties on wall shear stress (WSS) indices of a cerebral aneurysm. We created a virtual patient database by varying systemic properties within physiological ranges. BCs for 3D-CFD simulations were derived using a pulse wave propagation model for each realization of the virtual database. WSS indices were derived from the 3D simulations and their variabilities quantified. Variations in BCs, caused by changes in systemic properties, yielded variabilities in the WSS indices that were of the same order of magnitude as differences in these WSS indices between ruptured and unruptured aneurysms. Sensitivity analysis showed that the systemic properties impacted both in- and outlet BCs simultaneously and altered the WSS indices. We conclude that the influence of variations in patient-specific systemic properties on WSS indices should be evaluated when using WSS indices in multidisciplinary rupture prediction models.
Collapse
|
7
|
Rapadamnaba R, Nicoud F, Mohammadi B. Backward sensitivity analysis and reduced-order covariance estimation in noninvasive parameter identification for cerebral arteries. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3170. [PMID: 30426715 DOI: 10.1002/cnm.3170] [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: 05/03/2018] [Revised: 11/02/2018] [Accepted: 11/03/2018] [Indexed: 06/09/2023]
Abstract
Using a previously developed inversion platform for functional cerebral medical imaging with ensemble Kalman filters, this work analyzes the sensitivity of the results with respect to different parameters entering the physical model and inversion procedure, such as the inlet flow rate from the heart, the choice of the boundary conditions, and the nonsymmetry in the network terminations. It also proposes an alternative low complexity construction for the covariance matrix of the hemodynamic parameters of a network of arteries including the circle of Willis. The platform takes as input patient-specific blood flow rates extracted from magnetic resonance angiography and magnetic resonance imaging (dicom files) and is applied to several available patients data. The paper presents full analysis of the results for one of these patients, including a sensitivity study with respect to the proximal and distal boundary conditions. The results notably show that the uncertainties on the inlet flow rate led to uncertainties of the same order of magnitude in the estimated parameters (blood pressure and elastic parameters) and that three-lumped parameters boundary conditions are necessary for a correct retrieval of the target signals.
Collapse
Affiliation(s)
| | - Franck Nicoud
- IMAG, Université de Montpellier, CC051, Montpellier, France
| | | |
Collapse
|
8
|
Kovarovic B, Woo HH, Fiorella D, Lieber BB, Sadasivan C. Pressure and Flow Rate Changes During Contrast Injections in Cerebral Angiography: Correlation to Reflux Length. Cardiovasc Eng Technol 2018; 9:226-239. [PMID: 29497965 DOI: 10.1007/s13239-018-0344-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Accepted: 02/22/2018] [Indexed: 10/17/2022]
Abstract
Cerebral angiography involves the antegrade injection of contrast media through a catheter into the vasculature to visualize the region of interest under X-ray imaging. Depending on the injection and blood flow parameters, the bolus of contrast can propagate in the upstream direction and proximal to the catheter tip, at which point contrast is said to have refluxed. In this in vitro study, we investigate the relationship of fundamental hemodynamic variables to this phenomenon. Contrast injections were carried out under steady and pulsatile flow using various vessel diameters, catheter sizes, working fluid flow rates, and injection rates. The distance from the catheter tip to the proximal edge of the contrast bolus, called reflux length, was measured on the angiograms; the relation of this reflux length to different hemodynamic parameters was evaluated. Results show that contrast reflux occurs when the pressure distal to the catheter tip increases to be greater than the pressure proximal to the catheter tip. The ratio of this pressure difference to the baseline flow rate, called reflux resistance here, was linearly correlated to the normalized reflux length (reflux length/vessel diameter). Further, the ratio of blood flow to contrast fluid momentums, called the Craya-Curtet number, was correlated to the normalized reflux length via a sigmoid function. A sigmoid function was also found to be representative of the relationship between the ratio of the Reynolds numbers of blood flow to contrast and the normalized reflux length. As described by previous reports, catheter based contrast injections cause substantial increases in local flow and pressure. Contrast reflux should generally be avoided during standard antegrade angiography. Our study shows two specific correlations between contrast reflux length and baseline and intra-injection parameters that have not been published previously. Further studies need to be conducted to fully characterize the phenomena and to extract reliable indicators of clinical utility. Parameters relevant to cerebral angiography are studied here, but the essential principles are applicable to all angiographic procedures involving antegrade catheter injections.
Collapse
Affiliation(s)
- Brandon Kovarovic
- Department of Biomedical Engineering, Stony Brook University, 102 Bioengineering Building, Stony Brook, NY, 11794-5281, USA
| | - Henry H Woo
- Department of Neurological Surgery, Stony Brook University, HSC T-12, Rm 080, Stony Brook, NY, 11794-8122, USA
| | - David Fiorella
- Department of Neurological Surgery, Stony Brook University, HSC T-12, Rm 080, Stony Brook, NY, 11794-8122, USA
| | - Baruch B Lieber
- Department of Biomedical Engineering, Stony Brook University, 102 Bioengineering Building, Stony Brook, NY, 11794-5281, USA.,Department of Neurological Surgery, Stony Brook University, HSC T-12, Rm 080, Stony Brook, NY, 11794-8122, USA
| | - Chander Sadasivan
- Department of Neurological Surgery, Stony Brook University, HSC T-12, Rm 080, Stony Brook, NY, 11794-8122, USA.
| |
Collapse
|
9
|
Flow stagnation volume and abdominal aortic aneurysm growth: Insights from patient-specific computational flow dynamics of Lagrangian-coherent structures. Comput Biol Med 2018; 92:98-109. [DOI: 10.1016/j.compbiomed.2017.10.033] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2017] [Revised: 10/09/2017] [Accepted: 10/28/2017] [Indexed: 12/23/2022]
|
10
|
Lal R, Nicoud F, Bars EL, Deverdun J, Molino F, Costalat V, Mohammadi B. Non Invasive Blood Flow Features Estimation in Cerebral Arteries from Uncertain Medical Data. Ann Biomed Eng 2017; 45:2574-2591. [DOI: 10.1007/s10439-017-1904-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2016] [Accepted: 08/12/2017] [Indexed: 11/30/2022]
|
11
|
Vali A, Abla AA, Lawton MT, Saloner D, Rayz VL. Computational Fluid Dynamics modeling of contrast transport in basilar aneurysms following flow-altering surgeries. J Biomech 2016; 50:195-201. [PMID: 27890537 DOI: 10.1016/j.jbiomech.2016.11.028] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2016] [Accepted: 11/05/2016] [Indexed: 10/20/2022]
Abstract
In vivo measurement of blood velocity fields and flow descriptors remains challenging due to image artifacts and limited resolution of current imaging methods; however, in vivo imaging data can be used to inform and validate patient-specific computational fluid dynamics (CFD) models. Image-based CFD can be particularly useful for planning surgical interventions in complicated cases such as fusiform aneurysms of the basilar artery, where it is crucial to alter pathological hemodynamics while preserving flow to the distal vasculature. In this study, patient-specific CFD modeling was conducted for two basilar aneurysm patients considered for surgical treatment. In addition to velocity fields, transport of contrast agent was simulated for the preoperative and postoperative conditions using two approaches. The transport of a virtual contrast passively following the flow streamlines was simulated to predict post-surgical flow regions prone to thrombus deposition. In addition, the transport of a mixture of blood with an iodine-based contrast agent was modeled to compare and verify the CFD results with X-ray angiograms. The CFD-predicted patterns of contrast flow were qualitatively compared to in vivo X-ray angiograms acquired before and after the intervention. The results suggest that the mixture modeling approach, accounting for the flow rates and properties of the contrast injection, is in better agreement with the X-ray angiography data. The virtual contrast modeling assessed the residence time based on flow patterns unaffected by the injection procedure, which makes the virtual contrast modeling approach better suited for prediction of thrombus deposition, which is not limited to the peri-procedural state.
Collapse
Affiliation(s)
- Alireza Vali
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA
| | - Adib A Abla
- Department of Neurosurgery, University of Arkansas for Medical Science, AR, USA
| | - Michael T Lawton
- Department of Neurological Surgery, University of California, San Francisco, CA, USA
| | - David Saloner
- Department of Radiology and Biomedical Imaging University of California, San Francisco, CA, USA
| | - Vitaliy L Rayz
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA; Department of Mechanical Engineering, University of Wisconsin, Milwaukee, WI, USA.
| |
Collapse
|
12
|
Lin CJ, Hung SC, Chang FC, Guo WY, Luo CB, Kowarschik M, Chu WF, Liou AJYK. Finding the optimal deconvolution algorithm for MR perfusion in carotid stenosis: Correlations with angiographic cerebral circulation time. J Neuroradiol 2016; 43:290-6. [PMID: 27038737 DOI: 10.1016/j.neurad.2016.02.002] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2015] [Accepted: 02/08/2016] [Indexed: 11/28/2022]
Abstract
PURPOSE The aim of our study is to explore the impacts of different deconvolution algorithms on correlations between CBF, MTT, CBV, TTP, Tmax from MR perfusion (MRP) and angiography cerebral circulation time (CCT). METHODS Retrospectively, 30 patients with unilateral carotid stenosis, and available pre-stenting MRP and angiography were included for analysis. All MRPs were conducted in a 1.5-T MR scanner. Standard singular value decomposition, block-circulant, and two delay-corrected algorithms were used as the deconvolution methods. All angiographies were obtained in the same bi-plane flat-detector angiographic machine. A contrast bolus of 12mL was administrated via angiocatheter at a rate of 8mL/s. The acquisition protocols were the same for all cases. CCT was defined as the difference between time to peak from the cavernous ICA and the parietal vein in lateral view. Pearson correlations were calculated for CCT and CBF, MTT, CBV, TTP, Tmax. RESULTS The correlation between CCT and MTT was highest with Tmax (r=0.65), followed by MTT (r=0.60), CBF (r=-0.57), and TTP (r=0.33) when standard singular value decomposition was used. No correlation with CBV was noted. CONCLUSIONS MRP using a singular value decomposition algorithm confirmed the feasibility of quantifying cerebral blood flow deficit in steno-occlusive disease within the angio-room. This approach might further improve patient safety by providing immediate cerebral hemodynamics without extraradiation and iodine contrast.
Collapse
Affiliation(s)
- Chung-Jung Lin
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan
| | - Sheng-Che Hung
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan
| | - Feng-Chi Chang
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan
| | - Wan-Yuo Guo
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan.
| | - Chao-Bao Luo
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan
| | - Markus Kowarschik
- Siemens AG, healthcare sector, angiography and interventional X-Ray systems, Forchheim, Germany
| | - Wei-Fa Chu
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan
| | - Adrian J Y Kang Liou
- Department of radiology, Taipei Veterans general hospital, 201, Shipai Road, Section 2, 11217 Taipei, Taiwan; School of medicine, National Yang Ming university, Taipei, Taiwan
| |
Collapse
|
13
|
Manini S, Antiga L, Botti L, Remuzzi A. pyNS: an open-source framework for 0D haemodynamic modelling. Ann Biomed Eng 2014; 43:1461-73. [PMID: 25549775 DOI: 10.1007/s10439-014-1234-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2014] [Accepted: 12/19/2014] [Indexed: 10/24/2022]
Abstract
A number of computational approaches have been proposed for the simulation of haemodynamics and vascular wall dynamics in complex vascular networks. Among them, 0D pulse wave propagation methods allow to efficiently model flow and pressure distributions and wall displacements throughout vascular networks at low computational costs. Although several techniques are documented in literature, the availability of open-source computational tools is still limited. We here present python Network Solver, a modular solver framework for 0D problems released under a BSD license as part of the archToolkit ( http://archtk.github.com ). As an application, we describe patient-specific models of the systemic circulation and detailed upper extremity for use in the prediction of maturation after surgical creation of vascular access for haemodialysis.
Collapse
|
14
|
Groen D, Borgdorff J, Bona-Casas C, Hetherington J, Nash RW, Zasada SJ, Saverchenko I, Mamonski M, Kurowski K, Bernabeu MO, Hoekstra AG, Coveney PV. Flexible composition and execution of high performance, high fidelity multiscale biomedical simulations. Interface Focus 2014; 3:20120087. [PMID: 24427530 DOI: 10.1098/rsfs.2012.0087] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Multiscale simulations are essential in the biomedical domain to accurately model human physiology. We present a modular approach for designing, constructing and executing multiscale simulations on a wide range of resources, from laptops to petascale supercomputers, including combinations of these. Our work features two multiscale applications, in-stent restenosis and cerebrovascular bloodflow, which combine multiple existing single-scale applications to create a multiscale simulation. These applications can be efficiently coupled, deployed and executed on computers up to the largest (peta) scale, incurring a coupling overhead of 1-10% of the total execution time.
Collapse
Affiliation(s)
- D Groen
- Centre for Computational Science, University College London, UK
| | - J Borgdorff
- Section Computational Science, University of Amsterdam, The Netherlands
| | - C Bona-Casas
- Section Computational Science, University of Amsterdam, The Netherlands
| | - J Hetherington
- Centre for Computational Science, University College London, UK
| | - R W Nash
- Centre for Computational Science, University College London, UK
| | - S J Zasada
- Centre for Computational Science, University College London, UK
| | | | - M Mamonski
- Poznan Supercomputing and Networking Center, Poznan, Poland
| | - K Kurowski
- Poznan Supercomputing and Networking Center, Poznan, Poland
| | - M O Bernabeu
- Centre for Computational Science, University College London, UK
| | - A G Hoekstra
- Section Computational Science, University of Amsterdam, The Netherlands
| | - P V Coveney
- Centre for Computational Science, University College London, UK
| |
Collapse
|
15
|
Pereira VM, Ouared R, Brina O, Bonnefous O, Satwiaski J, Aerts H, Ruijters D, van Nijnatten F, Perren F, Bijlenga P, Schaller K, Lovblad KO. Quantification of internal carotid artery flow with digital subtraction angiography: validation of an optical flow approach with Doppler ultrasound. AJNR Am J Neuroradiol 2014; 35:156-63. [PMID: 23928145 PMCID: PMC7966491 DOI: 10.3174/ajnr.a3662] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2012] [Accepted: 03/28/2013] [Indexed: 11/07/2022]
Abstract
BACKGROUND AND PURPOSE Digital subtraction angiography is the reference standard technique to evaluate intracranial vascular anatomy and used on the endovascular treatment of vascular diseases. A dedicated optical flow-based algorithm was applied to DSA to measure arterial flow. The first quantification results of internal carotid artery flow validated with Doppler sonography are reported. MATERIALS AND METHODS We included 22 consecutive patients who underwent endovascular procedures. To assess the sensitivity of the algorithm to contrast agent-blood mixing dynamics, we acquired high-frame DSA series (60 images/s) with different injection rates: 1.5 mL/s (n = 19), 2.0 mL/s (n = 18), and 3.0 mL/s (n = 13). 3D rotational angiography was used to extract the centerline of the vessel and the arterial section necessary for volume flow calculation. Optical flow was used to measure flow velocities in straight parts of the ICAs; these data were further compared with Doppler sonography data. DSA mean flow rates were linearly regressed on Doppler sonography measurements, and regression slope coefficient bias from value 1 was analyzed within the 95% confidence interval. RESULTS DSA mean flow rates measured with the optical flow approach significantly matched Doppler sonography measurements (slope regression coefficient, b = 0.83 ± 0.19, P = .05) for injection rate = 2.0 mL/s and circulating volumetric blood flow <6 mL/s. For injection rate = 1.5 mL/s, volumetric blood flow <3 mL/s correlated well with Doppler sonography (b = 0.67 ± 0.33, P = .05). Injection rate = 3.0 mL/s failed to provide DSA-optical flow measurements correlating with Doppler sonography because of the lack of measurable pulsatility. CONCLUSIONS A new model-free optical flow technique was tested reliably on the ICA. DSA-based blood flow velocity measurements were essentially validated with Doppler sonography whenever the conditions of measurable pulsatility were achieved (injection rates = 1.5 and 2.0 mL/s).
Collapse
|
16
|
Fahy P, McCarthy P, Sultan S, Hynes N, Delassus P, Morris L. An experimental investigation of the hemodynamic variations due to aplastic vessels within three-dimensional phantom models of the circle of Willis. Ann Biomed Eng 2013; 42:123-38. [PMID: 24018609 DOI: 10.1007/s10439-013-0905-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2013] [Accepted: 08/29/2013] [Indexed: 11/26/2022]
Abstract
A complete circle of Willis (CoW) is found in approximately 30-50% of the population. Anatomical variations, such as absent or surgically clamped vessels, can result in undesirable flow patterns. These can affect the brain's ability to maintain cerebral perfusion and the formation of cerebral aneurysms. An experimental test system was developed to simulate cerebral physiological conditions through three flexible 3D patient-specific models of complete and incomplete CoW geometries. Flow visualizations were performed with isobaric dyes and the mapped dye streamlines were tracked throughout the models. Three to seven flow impact locations were observed for all configurations, corresponding to known sites for aneurysmal formation. Uni and bi-directional cross-flows occurred along the communicating arteries. The greatest shunting of flow occurred for a missing pre-communicating anterior (A1) and posterior (P1) cerebral arteries. The anterior cerebral arteries had the greatest reduction (15-37%) in efferent flow rates for missing either a unilateral A1 or bilateral P1 segments. The bi-directional cross-flows, with multiple afferent flow mixing, observed along the communicating arteries may explain the propensity of aneurysm formation at these sites. Reductions in efferent flow rates due to aplastic vessel configurations may affect normal brain function.
Collapse
Affiliation(s)
- Paul Fahy
- Galway Medical Technologies Centre (GMedTech), Department of Mechanical and Industrial Engineering, Galway Mayo Institute of Technology, Galway, Ireland
| | | | | | | | | | | |
Collapse
|
17
|
Pereira V, Brina O, Marcos Gonzales A, Narata A, Bijlenga P, Schaller K, Lovblad K, Ouared R. Evaluation of the influence of inlet boundary conditions on computational fluid dynamics for intracranial aneurysms: A virtual experiment. J Biomech 2013; 46:1531-9. [DOI: 10.1016/j.jbiomech.2013.03.024] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2012] [Revised: 03/14/2013] [Accepted: 03/16/2013] [Indexed: 10/26/2022]
|
18
|
Bernabeu MO, Nash RW, Groen D, Carver HB, Hetherington J, Krüger T, Coveney PV. Impact of blood rheology on wall shear stress in a model of the middle cerebral artery. Interface Focus 2013; 3:20120094. [PMID: 24427534 PMCID: PMC3638489 DOI: 10.1098/rsfs.2012.0094] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Perturbations to the homeostatic distribution of mechanical forces exerted by blood on the endothelial layer have been correlated with vascular pathologies, including intracranial aneurysms and atherosclerosis. Recent computational work suggests that, in order to correctly characterize such forces, the shear-thinning properties of blood must be taken into account. To the best of our knowledge, these findings have never been compared against experimentally observed pathological thresholds. In this work, we apply the three-band diagram (TBD) analysis due to Gizzi et al. (Gizzi et al. 2011 Three-band decomposition analysis of wall shear stress in pulsatile flows. Phys. Rev. E 83, 031902. (doi:10.1103/PhysRevE.83.031902)) to assess the impact of the choice of blood rheology model on a computational model of the right middle cerebral artery. Our results show that, in the model under study, the differences between the wall shear stress predicted by a Newtonian model and the well-known Carreau-Yasuda generalized Newtonian model are only significant if the vascular pathology under study is associated with a pathological threshold in the range 0.94-1.56 Pa, where the results of the TBD analysis of the rheology models considered differs. Otherwise, we observe no significant differences.
Collapse
Affiliation(s)
- Miguel O. Bernabeu
- Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
- CoMPLEX, University College London, Physics Building, Gower Street, London WC1E 6BT, UK
| | - Rupert W. Nash
- Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
| | - Derek Groen
- Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
| | - Hywel B. Carver
- Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
- CoMPLEX, University College London, Physics Building, Gower Street, London WC1E 6BT, UK
| | - James Hetherington
- Research Software Development Team, Research Computing and Facilitating Services, University College London, Podium Building, 1st Floor, Gower Street, London WC1E 6BT, UK
| | - Timm Krüger
- Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
| | - Peter V. Coveney
- Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
| |
Collapse
|
19
|
Merkx MAG, Huberts W, Bosboom EMH, Bode AS, Bescós JO, Tordoir JHM, Breeuwer M, van de Vosse FN. The benefit of non contrast-enhanced magnetic resonance angiography for predicting vascular access surgery outcome: a computer model perspective. PLoS One 2013; 8:e53615. [PMID: 23390490 PMCID: PMC3563627 DOI: 10.1371/journal.pone.0053615] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2012] [Accepted: 11/30/2012] [Indexed: 11/18/2022] Open
Abstract
INTRODUCTION Vascular access (VA) surgery, a prerequisite for hemodialysis treatment of end-stage renal-disease (ESRD) patients, is hampered by complication rates, which are frequently related to flow enhancement. To assist in VA surgery planning, a patient-specific computer model for postoperative flow enhancement was developed. The purpose of this study is to assess the benefit of non contrast-enhanced magnetic resonance angiography (NCE-MRA) data as patient-specific geometrical input for the model-based prediction of surgery outcome. METHODS 25 ESRD patients were included in this study. All patients received a NCE-MRA examination of the upper extremity blood vessels in addition to routine ultrasound (US). Local arterial radii were assessed from NCE-MRA and converted to model input using a linear fit per artery. Venous radii were determined with US. The effect of radius measurement uncertainty on model predictions was accounted for by performing Monte-Carlo simulations. The resulting flow prediction interval of the computer model was compared with the postoperative flow obtained from US. Patients with no overlap between model-based prediction and postoperative measurement were further analyzed to determine whether an increase in geometrical detail improved computer model prediction. RESULTS Overlap between postoperative flows and model-based predictions was obtained for 71% of patients. Detailed inspection of non-overlapping cases revealed that the geometrical details that could be assessed from NCE-MRA explained most of the differences, and moreover, upon addition of these details in the computer model the flow predictions improved. CONCLUSIONS The results demonstrate clearly that NCE-MRA does provide valuable geometrical information for VA surgery planning. Therefore, it is recommended to use this modality, at least for patients at risk for local or global narrowing of the blood vessels as well as for patients for whom an US-based model prediction would not overlap with surgical choice, as the geometrical details are crucial for obtaining accurate flow predictions.
Collapse
Affiliation(s)
- Maarten A G Merkx
- Department of Biomedical Engineering, Maastricht University Medical Center, Maastricht, The Netherlands.
| | | | | | | | | | | | | | | |
Collapse
|
20
|
Bonnefous O, Pereira VM, Ouared R, Brina O, Aerts H, Hermans R, van Nijnatten F, Stawiaski J, Ruijters D. Quantification of arterial flow using digital subtraction angiography. Med Phys 2012; 39:6264-75. [DOI: 10.1118/1.4754299] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
21
|
Spronck B, Martens EGHJ, Gommer ED, van de Vosse FN. A lumped parameter model of cerebral blood flow control combining cerebral autoregulation and neurovascular coupling. Am J Physiol Heart Circ Physiol 2012; 303:H1143-53. [PMID: 22777421 DOI: 10.1152/ajpheart.00303.2012] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Cerebral blood flow regulation is based on a variety of different mechanisms, of which the relative regulatory role remains largely unknown. The cerebral regulatory system expresses two regulatory properties: cerebral autoregulation and neurovascular coupling. Since partly the same mechanisms play a role in cerebral autoregulation and neurovascular coupling, this study aimed to develop a physiologically based mathematical model of cerebral blood flow regulation combining these properties. A lumped parameter model of the P2 segment of the posterior cerebral artery and its distal vessels was constructed. Blood flow regulation is exerted at the arteriolar level by vascular smooth muscle and implements myogenic, shear stress based, neurogenic, and metabolic mechanisms. In eight healthy subjects, cerebral autoregulation and neurovascular coupling were challenged by squat-stand maneuvers and visual stimulation using a checkerboard pattern, respectively. Cerebral blood flow velocity was measured using transcranial Doppler, whereas blood pressure was measured by finger volume clamping. In seven subjects, the model proposed fits autoregulation and neurovascular coupling measurement data well. Myogenic regulation is found to dominate the autoregulatory response. Neurogenic regulation, although only implemented as a first-order mechanism, describes neurovascular coupling responses to a great extent. It is concluded that our single, integrated model of cerebral blood flow control may be used to identify the main mechanisms affecting cerebral blood flow regulation in individual subjects.
Collapse
Affiliation(s)
- Bart Spronck
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands.
| | | | | | | |
Collapse
|
22
|
A numerical method of reduced complexity for simulating vascular hemodynamics using coupled 0D lumped and 1D wave propagation models. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2012; 2012:156094. [PMID: 22654957 PMCID: PMC3361674 DOI: 10.1155/2012/156094] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2011] [Revised: 12/11/2011] [Accepted: 01/13/2012] [Indexed: 11/18/2022]
Abstract
A computational method of reduced complexity is developed for simulating vascular hemodynamics by combination of one-dimensional (1D) wave propagation models for the blood vessels with zero-dimensional (0D) lumped models for the microcirculation. Despite the reduced dimension, current algorithms used to solve the model equations and simulate pressure and flow are rather complex, thereby limiting acceptance in the medical field. This complexity mainly arises from the methods used to combine the 1D and the 0D model equations. In this paper a numerical method is presented that no longer requires additional coupling methods and enables random combinations of 1D and 0D models using pressure as only state variable. The method is applied to a vascular tree consisting of 60 major arteries in the body and the head. Simulated results are realistic. The numerical method is stable and shows good convergence.
Collapse
|
23
|
|