1
|
Martín Tempestti J, Kim S, Lindsey BD, Veneziani A. A Pseudo-Spectral Method for Wall Shear Stress Estimation from Doppler Ultrasound Imaging in Coronary Arteries. Cardiovasc Eng Technol 2024:10.1007/s13239-024-00741-2. [PMID: 39103664 DOI: 10.1007/s13239-024-00741-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 06/24/2024] [Indexed: 08/07/2024]
Abstract
PURPOSE The Wall Shear Stress (WSS) is the component tangential to the boundary of the normal stress tensor in an incompressible fluid, and it has been recognized as a quantity of primary importance in predicting possible adverse events in cardiovascular diseases, in general, and in coronary diseases, in particular. The quantification of the WSS in patient-specific settings can be achieved by performing a Computational Fluid Dynamics (CFD) analysis based on patient geometry, or it can be retrieved by a numerical approximation based on blood flow velocity data, e.g., ultrasound (US) Doppler measurements. This paper presents a novel method for WSS quantification from 2D vector Doppler measurements. METHODS Images were obtained through unfocused plane waves and transverse oscillation to acquire both in-plane velocity components. These velocity components were processed using pseudo-spectral differentiation techniques based on Fourier approximations of the derivatives to compute the WSS. RESULTS Our Pseudo-Spectral Method (PSM) is tested in two vessel phantoms, straight and stenotic, where a steady flow of 15 mL/min is applied. The method is successfully validated against CFD simulations and compared against current techniques based on the assumption of a parabolic velocity profile. The PSM accurately detected Wall Shear Stress (WSS) variations in geometries differing from straight cylinders, and is less sensitive to measurement noise. In particular, when using synthetic data (noise free, e.g., generated by CFD) on cylindrical geometries, the Poiseuille-based methods and PSM have comparable accuracy; on the contrary, when using the data retrieved from US measures, the average error of the WSS obtained with the PSM turned out to be 3 to 9 times smaller than that obtained by state-of-the-art methods. CONCLUSION The pseudo-spectral approach allows controlling the approximation errors in the presence of noisy data. This gives a more accurate alternative to the present standard and a less computationally expensive choice compared to CFD, which also requires high-quality data to reconstruct the vessel geometry.
Collapse
Affiliation(s)
| | - Saeyoung Kim
- George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, 801 Ferst Dr., Atlanta, GA, 30332, USA
- Interdisciplinary BioEngineering Graduate Program, Georgia Institute of Technology, 315 Ferst Dr., Atlanta, GA, 30332, USA
| | - Brooks D Lindsey
- George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, 801 Ferst Dr., Atlanta, GA, 30332, USA
- Interdisciplinary BioEngineering Graduate Program, Georgia Institute of Technology, 315 Ferst Dr., Atlanta, GA, 30332, USA
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, 313 Ferst Dr NW, Atlanta, GA, 30332, USA
| | - Alessandro Veneziani
- Department of Mathematics, Emory University, 400 Dowman Dr, Atlanta, 30322, GA, USA
- Department of Computer Science, Emory University, 400 Dowman Dr, Atlanta, GA, 30322, USA
| |
Collapse
|
2
|
Fandaros M, Kwok C, Wolf Z, Labropoulos N, Yin W. Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use. Cardiovasc Eng Technol 2024:10.1007/s13239-024-00731-4. [PMID: 38710896 DOI: 10.1007/s13239-024-00731-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 04/29/2024] [Indexed: 05/08/2024]
Abstract
PURPOSE Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential. METHODS This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed. RESULTS Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data. CONCLUSIONS Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.
Collapse
Affiliation(s)
- Marina Fandaros
- Department of Biomedical Engineering, Stony Brook University, Bioengineering Building, Room 109, 11794, Stony Brook, NY, USA
| | - Chloe Kwok
- Department of Biomedical Engineering, Stony Brook University, Bioengineering Building, Room 109, 11794, Stony Brook, NY, USA
| | - Zachary Wolf
- Department of Biomedical Engineering, Stony Brook University, Bioengineering Building, Room 109, 11794, Stony Brook, NY, USA
| | - Nicos Labropoulos
- Department of Surgery, Stony Brook Medicine, 11794, Stony Brook, NY, USA
| | - Wei Yin
- Department of Biomedical Engineering, Stony Brook University, Bioengineering Building, Room 109, 11794, Stony Brook, NY, USA.
| |
Collapse
|
3
|
Zhang H, Song X, Wu R, Li N, Hou Q, Xie J, Hou Y, Qiao A. A novel method for noninvasive quantification of fractional flow reserve based on the custom function. Front Bioeng Biotechnol 2023; 11:1207300. [PMID: 37711442 PMCID: PMC10498765 DOI: 10.3389/fbioe.2023.1207300] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Accepted: 08/21/2023] [Indexed: 09/16/2023] Open
Abstract
Boundary condition settings are key risk factors for the accuracy of noninvasive quantification of fractional flow reserve (FFR) based on computed tomography angiography (i.e., FFRCT). However, transient numerical simulation-based FFRCT often ignores the three-dimensional (3D) model of coronary artery and clinical statistics of hyperemia state set by boundary conditions, resulting in insufficient computational accuracy and high computational cost. Therefore, it is necessary to develop the custom function that combines the 3D model of the coronary artery and clinical statistics of hyperemia state for boundary condition setting, to accurately and quickly quantify FFRCT under steady-state numerical simulations. The 3D model of the coronary artery was reconstructed by patient computed tomography angiography (CTA), and coronary resting flow was determined from the volume and diameter of the 3D model. Then, we developed the custom function that took into account the interaction of stenotic resistance, microcirculation resistance, inlet aortic pressure, and clinical statistics of resting to hyperemia state due to the effect of adenosine on boundary condition settings, to accurately and rapidly identify coronary blood flow for quantification of FFRCT calculation (FFRU). We tested the diagnostic accuracy of FFRU calculation by comparing it with the existing methods (CTA, coronary angiography (QCA), and diameter-flow method for calculating FFR (FFRD)) based on invasive FFR of 86 vessels in 73 patients. The average computational time for FFRU calculation was greatly reduced from 1-4 h for transient numerical simulations to 5 min per simulation, which was 2-fold less than the FFRD method. According to the results of the Bland-Altman analysis, the consistency between FFRU and invasive FFR of 86 vessels was better than that of FFRD. The area under the receiver operating characteristic curve (AUC) for CTA, QCA, FFRD and FFRU at the lesion level were 0.62 (95% CI: 0.51-0.74), 0.67 (95% CI: 0.56-0.79), 0.85 (95% CI: 0.76-0.94), and 0.93 (95% CI: 0.87-0.98), respectively. At the patient level, the AUC was 0.61 (95% CI: 0.48-0.74) for CTA, 0.65 (95% CI: 0.53-0.77) for QCA, 0.83 (95% CI: 0.74-0.92) for FFRD, and 0.92 (95% CI: 0.89-0.96) for FFRU. The proposed novel method might accurately and rapidly identify coronary blood flow, significantly improve the accuracy of FFRCT calculation, and support its wide application as a diagnostic indicator in clinical practice.
Collapse
Affiliation(s)
- Honghui Zhang
- Key Laboratory of Intelligent Manufacturing Technology, College of Engineering, Inner Mongolia Minzu University, Tongliao, China
- Faculty of Environment and Life, Beijing University of Technology, Beijing, China
| | - Xiaorui Song
- School of Radiology, Shandong First Medical University and Shandong Academy of Medical Sciences, Tai’an, China
| | - Rile Wu
- Department of Neurology, Tong Liao City Hospital, Tongliao, China
| | - Na Li
- School of Radiology, Shandong First Medical University and Shandong Academy of Medical Sciences, Tai’an, China
| | - Qianwen Hou
- Faculty of Environment and Life, Beijing University of Technology, Beijing, China
| | - Jinjie Xie
- Beijing Anzhen Hospital, Capital Medical University, Beijing, China
- Department of Echocardiography, Jiahui International Hospital, Shanghai, China
| | - Yang Hou
- Shengjing Hospital, China Medical University, Shenyang, China
| | - Aike Qiao
- Faculty of Environment and Life, Beijing University of Technology, Beijing, China
| |
Collapse
|
4
|
Kim S, Jing B, Lindsey BD. Forward-viewing estimation of 3D blood flow velocity fields by intravascular ultrasound: Influence of the catheter on velocity estimation in stenoses. ULTRASONICS 2021; 117:106558. [PMID: 34461527 PMCID: PMC8448960 DOI: 10.1016/j.ultras.2021.106558] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Revised: 07/02/2021] [Accepted: 08/17/2021] [Indexed: 06/13/2023]
Abstract
Coronary artery disease is the most common type of cardiovascular disease, affecting > 18 million adults, and is responsible for > 365 k deaths per year in the U.S. alone. Wall shear stress (WSS) is an emerging indicator of likelihood of plaque rupture in coronary artery disease, however, non-invasive estimation of 3-D blood flow velocity and WSS is challenging due to the requirement for high spatial resolution at deep penetration depths in the presence of significant cardiac motion. Thus we propose minimally-invasive imaging with a catheter-based, 3-D intravascular forward-viewing ultrasound (FV US) transducer and present experiments to quantify the effect of the catheter on flow disturbance in stenotic vessel phantoms with realistic velocities and luminal diameters for both peripheral (6.33 mm) and coronary (4.74 mm) arteries. An external linear array ultrasound transducer was used to quantify 2-D velocity fields in vessel phantoms under various conditions of catheter geometry, luminal diameter, and position of the catheter relative to the stenosis at a frame rate of 5000 frames per second via a particle imaging velocimetry (PIV) approach. While a solid catheter introduced an underestimation of velocity measurement by > 20% relative to the case without a catheter, the hollow catheter introduced < 10% velocity overestimation, indicating that a hollow catheter design allowing internal blood flow reduces hemodynamic disturbance. In addition, for both peripheral and coronary arteries, the hollow catheter introduced < 3% deviation in flow velocity at the minimum luminal area compared to the control case. Finally, an initial comparison was made between velocity measurements acquired using a low frequency, catheter-based, 3-D intravascular FV US transducer and external linear array measurements, with relative error < 12% throughout the region of interest for a flow rate of 150 mL/min. While further system development is required, results suggest intravascular ultrasound characterization of blood flow velocity fields in stenotic vessels could be feasible with appropriate catheter design.
Collapse
Affiliation(s)
- Saeyoung Kim
- Georgia Institute of Technology, George W. Woodruff School of Mechanical Engineering, 801 Ferst Dr., Atlanta, GA 30332, USA; Georgia Institute of Technology, Interdisciplinary BioEngineering Graduate Program, 315 Ferst Dr., Atlanta, GA 30332, USA
| | - Bowen Jing
- Georgia Institute of Technology and Emory University, Wallace H. Coulter Department of Biomedical Engineering, 313 Ferst Dr NW, Atlanta, GA 30332, USA
| | - Brooks D Lindsey
- Georgia Institute of Technology, Interdisciplinary BioEngineering Graduate Program, 315 Ferst Dr., Atlanta, GA 30332, USA; Georgia Institute of Technology and Emory University, Wallace H. Coulter Department of Biomedical Engineering, 313 Ferst Dr NW, Atlanta, GA 30332, USA.
| |
Collapse
|
5
|
Abstract
Atherosclerosis is one of the main causes of cardiovascular events, namely, myocardium infarction and cerebral stroke, responsible for a great number of deaths every year worldwide. This pathology is caused by the progressive accumulation of low-density lipoproteins, cholesterol, and other substances on the arterial wall, narrowing its lumen. To date, many hemodynamic studies have been conducted experimentally and/or numerically; however, this disease is not yet fully understood. For this reason, the research of this pathology is still ongoing, mainly, resorting to computational methods. These have been increasingly used in biomedical research of atherosclerosis because of their high-performance hardware and software. Taking into account the attempts that have been made in computational techniques to simulate realistic conditions of blood flow in both diseased and healthy arteries, the present review aims to give an overview of the most recent numerical studies focused on coronary arteries, by addressing the blood viscosity models, and applied physiological flow conditions. In general, regardless of the boundary conditions, numerical studies have been contributed to a better understanding of the development of this disease, its diagnosis, and its treatment.
Collapse
|
6
|
Nguyen TD, Kadri OE, Voronov RS. An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine. Front Bioeng Biotechnol 2020; 8:529365. [PMID: 33102452 PMCID: PMC7546862 DOI: 10.3389/fbioe.2020.529365] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Accepted: 08/31/2020] [Indexed: 02/05/2023] Open
Abstract
Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.
Collapse
Affiliation(s)
- Thanh Danh Nguyen
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States
| | - Olufemi E Kadri
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States.,UC-P&G Simulation Center, University of Cincinnati, Cincinnati, OH, United States
| | - Roman S Voronov
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States.,Department of Biomedical Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States
| |
Collapse
|
7
|
Liu H, Ou S, Liu P, Xu Y, Gong Y, Xia L, Leng X, Leung TWH, Shi L, Zheng D. Dataset on blood flow and instantaneous wave-free ratio in normal and stenosed coronary arteries. Data Brief 2020; 32:106011. [PMID: 32793771 PMCID: PMC7408341 DOI: 10.1016/j.dib.2020.106011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2020] [Accepted: 07/08/2020] [Indexed: 11/29/2022] Open
Abstract
Instantaneous wave-free ratio (iFR) has been proposed as a hemodynamic parameter that can reliably reflect the blood flow in stenosed coronary arteries. Currently, there are few investigations on the quantitative analysis of iFR in the patients regarding the variation of microcirculatory resistance (MR). The data aim to provide geometric (cross-section area of branches) and hemodynamic (flow rate and iFR of branches) parameters of normal and stenosed coronary arteries derived from CFD simulation. The CFD simulation was performed on the three-dimensional artery models reconstructed from computed tomography (CT) images of four subjects. The hemodynamic parameters were obtained in six situations of MR to simulate coronary microvascular dysfunction (CMD). This dataset could be used as the reference to estimate the iFR and flow rate in patients with CMD and stenosis in coronary arteries. The geometric parameters could be used in the modelling of coronary arteries.
Collapse
Affiliation(s)
- Haipeng Liu
- Research Centre of Intelligent Healthcare, Faculty of Health and Life Science, Coventry University, Coventry CV1 5FB, UK
| | - Shanxing Ou
- Department of Radiology, General Hospital of Southern Theater Command, PLA, Guangzhou, China
| | - Panli Liu
- Department of Radiology, Guangzhou First People's Hospital, Nansha Hospital, Guangzhou, China
| | - Yuhang Xu
- Research Centre of Intelligent Healthcare, Faculty of Health and Life Science, Coventry University, Coventry CV1 5FB, UK
| | - Yinglan Gong
- Department of Biomedical Engineering, Zhejiang University, Hangzhou, China
| | - Ling Xia
- Department of Biomedical Engineering, Zhejiang University, Hangzhou, China
| | - Xinyi Leng
- Department of Medicine and Therapeutics, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, Hong Kong, China
| | - Thomas Wai Hong Leung
- Department of Medicine and Therapeutics, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, Hong Kong, China
| | - Lin Shi
- Department of Imaging and Interventional Radiology, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, Hong Kong, China
- Corresponding authors.
| | - Dingchang Zheng
- Research Centre of Intelligent Healthcare, Faculty of Health and Life Science, Coventry University, Coventry CV1 5FB, UK
- Corresponding authors.
| |
Collapse
|
8
|
Accelerating massively parallel hemodynamic models of coarctation of the aorta using neural networks. Sci Rep 2020; 10:9508. [PMID: 32528104 PMCID: PMC7289812 DOI: 10.1038/s41598-020-66225-0] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2020] [Accepted: 05/18/2020] [Indexed: 11/09/2022] Open
Abstract
Comorbidities such as anemia or hypertension and physiological factors related to exertion can influence a patient’s hemodynamics and increase the severity of many cardiovascular diseases. Observing and quantifying associations between these factors and hemodynamics can be difficult due to the multitude of co-existing conditions and blood flow parameters in real patient data. Machine learning-driven, physics-based simulations provide a means to understand how potentially correlated conditions may affect a particular patient. Here, we use a combination of machine learning and massively parallel computing to predict the effects of physiological factors on hemodynamics in patients with coarctation of the aorta. We first validated blood flow simulations against in vitro measurements in 3D-printed phantoms representing the patient’s vasculature. We then investigated the effects of varying the degree of stenosis, blood flow rate, and viscosity on two diagnostic metrics – pressure gradient across the stenosis (ΔP) and wall shear stress (WSS) - by performing the largest simulation study to date of coarctation of the aorta (over 70 million compute hours). Using machine learning models trained on data from the simulations and validated on two independent datasets, we developed a framework to identify the minimal training set required to build a predictive model on a per-patient basis. We then used this model to accurately predict ΔP (mean absolute error within 1.18 mmHg) and WSS (mean absolute error within 0.99 Pa) for patients with this disease.
Collapse
|
9
|
Xu L, Chen X, Cui M, Ren C, Yu H, Gao W, Li D, Zhao W. The improvement of the shear stress and oscillatory shear index of coronary arteries during Enhanced External Counterpulsation in patients with coronary heart disease. PLoS One 2020; 15:e0230144. [PMID: 32191730 PMCID: PMC7082042 DOI: 10.1371/journal.pone.0230144] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Accepted: 02/23/2020] [Indexed: 12/25/2022] Open
Abstract
Background Enhanced External Counterpulsation (EECP) can chronically relieve ischemic chest pain and improve the prognosis of coronary heart disease (CHD). Despite its role in mitigating heart complications, EECP and the mechanisms behind its therapeutic nature, such as its effects on blood flow hemodynamics, are still not fully understood. This study aims to elucidate the effect of EECP on significant hemodynamic parameters in the coronary arterial tree. Methods A finite volume method was used in conjunction with the inlet pressure wave (surrogated by the measured aortic pressure) before and during EECP and outlet flow resistance, assuming the blood as newtonian fluid. The time-average wall shear stress (TAWSS) and oscillatory shear index (OSI) were determined from the flow field. Results Regardless of the degree of vascular stenosis, hemodynamic conditions and flow patterns could be improved during EECP. In comparison with the original tree, the tree with a severe stenosis (75% area stenosis) demonstrated more significant improvement in hemodynamic conditions and flow patterns during EECP, with surface area ratio of TAWSS risk area (SAR-TAWSS) reduced from 12.3% to 6.7% (vs. SAR-TAWSS reduced from 7.2% to 5.6% in the original tree) and surface area ratio of OSI risk area (SAR-OSI) reduced from 6.8% to 2.5% (vs. SAR-OSI of both 0% before and during EECP in the original tree because of mild stenosis). Moreover, it was also shown that small ratio of diastolic pressure (D) and systolic pressure (S) (D/S) could only improve the hemodynamic condition mildly. The SAR-TAWSS reduction ratio significantly increased as D/S became larger. Conclusions A key finding of the study was that the improvement of hemodynamic conditions along the LMCA trees during EECP became more significant with the increase of D/S and the severity degree of stenoses at the bifurcation site. These findings have important implications on EECP as adjuvant therapy before or after percutaneous coronary intervention (PCI) in patients with diffuse atherosclerosis.
Collapse
Affiliation(s)
- Ling Xu
- NHC Key Laboratory of Cardiovascular Molecular Biology and Regulatory Peptides, Department of Cardiology, Peking University Third Hospital, Beijing, China
| | - Xi Chen
- School of Biomedical Engineering, Capital Medical University, Beijing, China
| | - Ming Cui
- NHC Key Laboratory of Cardiovascular Molecular Biology and Regulatory Peptides, Department of Cardiology, Peking University Third Hospital, Beijing, China
| | - Chuan Ren
- NHC Key Laboratory of Cardiovascular Molecular Biology and Regulatory Peptides, Department of Cardiology, Peking University Third Hospital, Beijing, China
| | - Haiyi Yu
- NHC Key Laboratory of Cardiovascular Molecular Biology and Regulatory Peptides, Department of Cardiology, Peking University Third Hospital, Beijing, China
| | - Wei Gao
- NHC Key Laboratory of Cardiovascular Molecular Biology and Regulatory Peptides, Department of Cardiology, Peking University Third Hospital, Beijing, China
| | - Dongguo Li
- School of Biomedical Engineering, Capital Medical University, Beijing, China
- * E-mail: (DGL); (WZ)
| | - Wei Zhao
- NHC Key Laboratory of Cardiovascular Molecular Biology and Regulatory Peptides, Department of Cardiology, Peking University Third Hospital, Beijing, China
- * E-mail: (DGL); (WZ)
| |
Collapse
|
10
|
Wan M, Ma L, Zhao X, Leng S, Zhang JM, Tan RS, Zhong L. Automatic Segmentation of Coronary Artery Lumen via Anisotropic Graph-cuts .. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020; 2019:4871-4874. [PMID: 31946952 DOI: 10.1109/embc.2019.8856353] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Coronary fractional flow reserve (FFR) is the ratio of distal to proximal pressure in a stenotic lesion and defines its presence of ischemia (i.e. FFR≤0.8). However, FFR requires pressure to be measured invasively across the lesion, which is barrier to wider use of the methods. Computational fluid dynamics (CFD) applied to computed tomography coronary angiography (CTCA) enable the computation of FFR (i.e. CT-FFR) noninvasively. The accurate coronary artery geometrical model is the prerequisite for CT-FFR. We aimed to develop a method to automatically segment the coronary artery lumen based on anisotrophic graph-cuts. 5 patients with coronary artery disease underwent CTCA. CTCA images were enhanced by Hessian Matrix filter and then to derive vesselness images and vessel directions. Centrelines were extracted using the shortest path algorithm. We modified the edge energy term based on the dissimilarity of neighboring voxels intensities. The inner product of vessel directions of two neighboring voxels was multiplied as well. All cases were automatically segmented successfully with average computational time of 79.5 seconds. Our anisotropic graph-cuts method reduced unnatural bumps on the lumen than classic graph-cuts methods. Our anisotrophic graph-cuts based method improved coronary lumen segmentation and downstream coronary geometry modelling, will play an important role for real-time CT-FFR computation in clinical practice.
Collapse
|
11
|
Acuna A, Berman AG, Damen FW, Meyers BA, Adelsperger AR, Bayer KC, Brindise MC, Bungart B, Kiel AM, Morrison RA, Muskat JC, Wasilczuk KM, Wen Y, Zhang J, Zito P, Goergen CJ. Computational Fluid Dynamics of Vascular Disease in Animal Models. J Biomech Eng 2019; 140:2676341. [PMID: 29570754 DOI: 10.1115/1.4039678] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2017] [Indexed: 12/19/2022]
Abstract
Recent applications of computational fluid dynamics (CFD) applied to the cardiovascular system have demonstrated its power in investigating the impact of hemodynamics on disease initiation, progression, and treatment outcomes. Flow metrics such as pressure distributions, wall shear stresses (WSS), and blood velocity profiles can be quantified to provide insight into observed pathologies, assist with surgical planning, or even predict disease progression. While numerous studies have performed simulations on clinical human patient data, it often lacks prediagnosis information and can be subject to large intersubject variability, limiting the generalizability of findings. Thus, animal models are often used to identify and manipulate specific factors contributing to vascular disease because they provide a more controlled environment. In this review, we explore the use of CFD in animal models in recent studies to investigate the initiating mechanisms, progression, and intervention effects of various vascular diseases. The first section provides a brief overview of the CFD theory and tools that are commonly used to study blood flow. The following sections are separated by anatomical region, with the abdominal, thoracic, and cerebral areas specifically highlighted. We discuss the associated benefits and obstacles to performing CFD modeling in each location. Finally, we highlight animal CFD studies focusing on common surgical treatments, including arteriovenous fistulas (AVF) and pulmonary artery grafts. The studies included in this review demonstrate the value of combining CFD with animal imaging and should encourage further research to optimize and expand upon these techniques for the study of vascular disease.
Collapse
Affiliation(s)
- Andrea Acuna
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Alycia G Berman
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Frederick W Damen
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Brett A Meyers
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, IN 47907 e-mail:
| | - Amelia R Adelsperger
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Kelsey C Bayer
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Melissa C Brindise
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, IN 47907 e-mail:
| | - Brittani Bungart
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Alexander M Kiel
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Rachel A Morrison
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Joseph C Muskat
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Kelsey M Wasilczuk
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Yi Wen
- Department of Agricultural and Biological Engineering, Purdue University, 225 South University Street, West Lafayette, IN 47907 e-mail:
| | - Jiacheng Zhang
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, IN 47907 e-mail:
| | - Patrick Zito
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Craig J Goergen
- ASME Membership Bioengineering Division, Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| |
Collapse
|
12
|
Monroy-Gonzalez AG, Alexanderson-Rosas E, Prakken NHJ, Juarez-Orozco LE, Walls-Laguarda L, Berrios-Barcenas EA, Meave-Gonzalez A, Groot JC, Slart RHJA, Tio RA. Myocardial bridging of the left anterior descending coronary artery is associated with reduced myocardial perfusion reserve: a 13N-ammonia PET study. Int J Cardiovasc Imaging 2019; 35:375-382. [PMID: 30267168 PMCID: PMC6428791 DOI: 10.1007/s10554-018-1460-8] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Accepted: 07/20/2018] [Indexed: 12/30/2022]
Abstract
Myocardial Bridging (MB) refers to the band of myocardium that abnormally overlies a segment of a coronary artery. This paper quantitatively evaluates the influence of MB of the left anterior descending artery (LAD) on myocardial perfusion of the entire left ventricle. We studied 131 consecutive patients who underwent hybrid rest/stress 13N-ammonia positron emission tomography (PET) and coronary computed tomography angiography (CCTA) due to suspected myocardial ischemia. Patients with previous myocardial infarction and/or significant coronary artery disease (≥ 50% stenosis) were excluded. Myocardial perfusion measurements were compared between patients with and without LAD-MB. Additionally, we evaluated the relationship between anatomical characteristics (length and depth) of LAD-MB and myocardial perfusion measurements. 17 (13%) patients presented a single LAD-MB. Global myocardial perfusion reserve (MPR) was lower in patients with LAD-MB than in patients without LAD-MB (1.9 ± 0.5 vs. 2.3 ± 0.6, p < 0.01). Global stress myocardial blood flow (MBF) was similar in patients with and without LAD-MB (2.2 ± 0.4 vs. 2.3 ± 0.7 ml/g/min, p = 0.40). Global rest MBF was higher in patients with LAD-MB than in patients without LAD-MB (1.2 ± 0.3 vs. 1.0 ± 0.2 ml/g/min, p < 0.01). Global rest MBF, stress MBF, and MPR quantifications were similar in patients with superficial and deep LAD-MB (all p = NS). We did not find any correlation between length and global rest MBF, stress MBF nor MPR (r = - 0.14, p = 0.59; r = 0.44, p = 0.07; and r = 0.45, p = 0.07 respectively). Quantitative myocardial perfusion suggests that LAD-MB may be related to impaired perfusion reserve, an indicator of microvascular dysfunction. Anatomical characteristics of LAD-MB were not related to changes in myocardial perfusion.
Collapse
Affiliation(s)
- Andrea G. Monroy-Gonzalez
- Medical Imaging Center, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | - Erick Alexanderson-Rosas
- National Institute of Cardiology Ignacio Chavez, Mexico City, Mexico
- Department of Physiology, National Autonomous University of Mexico, Mexico City, Mexico
| | - Niek H. J. Prakken
- Medical Imaging Center, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | - Luis E. Juarez-Orozco
- Medical Imaging Center, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | | | | | | | - Jan C. Groot
- Medical Imaging Center, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | - Riemer H. J. A. Slart
- Medical Imaging Center, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
- Biomedical Photonic Imaging, Faculty of Science and Technology, University of Twente, Enschede, The Netherlands
| | - Rene A. Tio
- Department of Cardiology, Catharina Hospital, Eindhoven, The Netherlands
| |
Collapse
|
13
|
Kharche SR, So A, Salerno F, Lee TY, Ellis C, Goldman D, McIntyre CW. Computational Assessment of Blood Flow Heterogeneity in Peritoneal Dialysis Patients' Cardiac Ventricles. Front Physiol 2018; 9:511. [PMID: 29867555 PMCID: PMC5968396 DOI: 10.3389/fphys.2018.00511] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Accepted: 04/20/2018] [Indexed: 01/28/2023] Open
Abstract
Dialysis prolongs life but augments cardiovascular mortality. Imaging data suggests that dialysis increases myocardial blood flow (BF) heterogeneity, but its causes remain poorly understood. A biophysical model of human coronary vasculature was used to explain the imaging observations, and highlight causes of coronary BF heterogeneity. Post-dialysis CT images from patients under control, pharmacological stress (adenosine), therapy (cooled dialysate), and adenosine and cooled dialysate conditions were obtained. The data presented disparate phenotypes. To dissect vascular mechanisms, a 3D human vasculature model based on known experimental coronary morphometry and a space filling algorithm was implemented. Steady state simulations were performed to investigate the effects of altered aortic pressure and blood vessel diameters on myocardial BF heterogeneity. Imaging showed that stress and therapy potentially increased mean and total BF, while reducing heterogeneity. BF histograms of one patient showed multi-modality. Using the model, it was found that total coronary BF increased as coronary perfusion pressure was increased. BF heterogeneity was differentially affected by large or small vessel blocking. BF heterogeneity was found to be inversely related to small blood vessel diameters. Simulation of large artery stenosis indicates that BF became heterogeneous (increase relative dispersion) and gave multi-modal histograms. The total transmural BF as well as transmural BF heterogeneity reduced due to large artery stenosis, generating large patches of very low BF regions downstream. Blocking of arteries at various orders showed that blocking larger arteries results in multi-modal BF histograms and large patches of low BF, whereas smaller artery blocking results in augmented relative dispersion and fractal dimension. Transmural heterogeneity was also affected. Finally, the effects of augmented aortic pressure in the presence of blood vessel blocking shows differential effects on BF heterogeneity as well as transmural BF. Improved aortic blood pressure may improve total BF. Stress and therapy may be effective if they dilate small vessels. A potential cause for the observed complex BF distributions (multi-modal BF histograms) may indicate existing large vessel stenosis. The intuitive BF heterogeneity methods used can be readily used in clinical studies. Further development of the model and methods will permit personalized assessment of patient BF status.
Collapse
Affiliation(s)
- Sanjay R Kharche
- Kidney Clinical Research Unit, Lawson's Health Research Institute, Victoria Hospital, London, ON, Canada.,Department of Medical Biophysics, Schulich School of Medicine and Dentistry, University of Western Ontario, London, ON, Canada
| | - Aaron So
- Department of Medical Biophysics, Schulich School of Medicine and Dentistry, University of Western Ontario, London, ON, Canada.,Robarts Research Institute, University of Western Ontario, London, ON, Canada
| | - Fabio Salerno
- Kidney Clinical Research Unit, Lawson's Health Research Institute, Victoria Hospital, London, ON, Canada
| | - Ting-Yim Lee
- Robarts Research Institute, University of Western Ontario, London, ON, Canada
| | - Chris Ellis
- Department of Medical Biophysics, Schulich School of Medicine and Dentistry, University of Western Ontario, London, ON, Canada
| | - Daniel Goldman
- Department of Medical Biophysics, Schulich School of Medicine and Dentistry, University of Western Ontario, London, ON, Canada
| | - Christopher W McIntyre
- Kidney Clinical Research Unit, Lawson's Health Research Institute, Victoria Hospital, London, ON, Canada.,Department of Medical Biophysics, Schulich School of Medicine and Dentistry, University of Western Ontario, London, ON, Canada
| |
Collapse
|
14
|
Liu H, Gong Y, Leng X, Xia L, Wong KS, Ou S, Leung TW, Wang D, Shi L. Estimating current and long-term risks of coronary artery in silico by fractional flow reserve, wall shear stress and low-density lipoprotein filtration rate. Biomed Phys Eng Express 2018. [DOI: 10.1088/2057-1976/aa9a09] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
15
|
Pagiatakis C, Tardif JC, L'Allier PL, Mongrain R. Effect of stenosis eccentricity on the functionality of coronary bifurcation lesions-a numerical study. Med Biol Eng Comput 2017; 55:2079-2095. [PMID: 28500478 DOI: 10.1007/s11517-017-1653-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Accepted: 04/25/2017] [Indexed: 01/10/2023]
Abstract
Interventional cardiologists still rely heavily on angiography for the evaluation of coronary lesion severity, despite its poor correlation with the presence of ischemia. In order to improve the accuracy of the current diagnostic procedures, an understanding of the relative influence of geometric characteristics on the induction of ischemia is required. This idea is especially important for coronary bifurcation lesions (CBLs), whose treatment is complex and is associated with high rates of peri- and post-procedural clinical events. Overall, it is unclear which geometric and morphological parameters of CBLs influence the onset of ischemia. More specifically, the effect of stenosis eccentricity is unknown. Computational fluid dynamic simulations, under a geometric multiscale framework, were executed for seven CBL configurations within the left main coronary artery bifurcation. Both concentric and eccentric stenosis profiles of mild to severe constriction were considered. By using a geometric multiscale framework, the fractional flow reserve, which is the gold-standard clinical diagnostic index, could be calculated and was compared between the eccentric and concentric profiles for each case. The results suggested that for configurations where the supplying vessel is stenosed, eccentricity could have a notable effect on and therefore be an important factor that influences configuration functionality.
Collapse
Affiliation(s)
- Catherine Pagiatakis
- Department of Mechanical Engineering, McGill University, 817 Sherbrooke Street West, Montreal, Quebec, H3A 0C3, Canada. .,Montreal Heart Institute, 5000 Belanger Street, Montreal, Quebec, H1T 1C8, Canada.
| | - Jean-Claude Tardif
- Montreal Heart Institute, 5000 Belanger Street, Montreal, Quebec, H1T 1C8, Canada.,Faculty of Medicine, Université de Montréal - Pavillon Roger-Gaudry, 2900 Edouard-Montpetit Boulevard, Montreal, Quebec, H3T 1J4, Canada
| | - Philippe L L'Allier
- Montreal Heart Institute, 5000 Belanger Street, Montreal, Quebec, H1T 1C8, Canada.,Faculty of Medicine, Université de Montréal - Pavillon Roger-Gaudry, 2900 Edouard-Montpetit Boulevard, Montreal, Quebec, H3T 1J4, Canada
| | - Rosaire Mongrain
- Department of Mechanical Engineering, McGill University, 817 Sherbrooke Street West, Montreal, Quebec, H3A 0C3, Canada.,Montreal Heart Institute, 5000 Belanger Street, Montreal, Quebec, H1T 1C8, Canada
| |
Collapse
|
16
|
Zhang JM, Luo T, Tan SY, Lomarda AM, Wong ASL, Keng FYJ, Allen JC, Huo Y, Su B, Zhao X, Wan M, Kassab GS, Tan RS, Zhong L. Hemodynamic analysis of patient-specific coronary artery tree. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2015; 31:e02708. [PMID: 25630671 DOI: 10.1002/cnm.2708] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2014] [Revised: 12/04/2014] [Accepted: 12/05/2014] [Indexed: 05/28/2023]
Abstract
Local hemodynamic parameters, such as wall shear stress (WSS), oscillatory shear index and relative resident time (RRT), have been linked to coronary plaque initiation and progression. In this study, a left coronary artery tree model was reconstructed from computed tomography angiography images of a patient with multiple stenoses. The geometry of the coronary artery tree model was virtually restored by eliminating the lesions, essentially re-creating the virtually healthy artery anatomy. Using numerical simulations, flow characteristics and hemodynamic parameter distributions in the stenosed and virtually healthy models were investigated. In the virtually healthy artery model, disturbed flows were found at four locations, prone to initialization of plaque formation. Low WSS and high RRT were exhibited in three of the four locations, and high WSS and low RRT were exhibited in the fourth. These findings suggest that coronary plaque is more likely to form in locations with disturbed flow conditions characterized by low WSS and high RRT or high WSS and low RRT. In addition, clinical index of fractional flow reserve was found to significantly correlate with blood flow rate, rather than anatomic parameters, such as diameter stenosis, which implied the importance of hemodynamic environment in stenosis formation.
Collapse
Affiliation(s)
- Jun-Mei Zhang
- National Heart Center Singapore, 5 Hospital Drive, 169609, Singapore; Duke-NUS Graduate Medical School, Singapore, 8 College Road, 169857, Singapore
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|