1
|
Fan L, Wang H, Kassab GS, Lee LC. Review of cardiac-coronary interaction and insights from mathematical modeling. WIREs Mech Dis 2024; 16:e1642. [PMID: 38316634 PMCID: PMC11081852 DOI: 10.1002/wsbm.1642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 12/10/2023] [Accepted: 01/08/2024] [Indexed: 02/07/2024]
Abstract
Cardiac-coronary interaction is fundamental to the function of the heart. As one of the highest metabolic organs in the body, the cardiac oxygen demand is met by blood perfusion through the coronary vasculature. The coronary vasculature is largely embedded within the myocardial tissue which is continually contracting and hence squeezing the blood vessels. The myocardium-coronary vessel interaction is two-ways and complex. Here, we review the different types of cardiac-coronary interactions with a focus on insights gained from mathematical models. Specifically, we will consider the following: (1) myocardial-vessel mechanical interaction; (2) metabolic-flow interaction and regulation; (3) perfusion-contraction matching, and (4) chronic interactions between the myocardium and coronary vasculature. We also provide a discussion of the relevant experimental and clinical studies of different types of cardiac-coronary interactions. Finally, we highlight knowledge gaps, key challenges, and limitations of existing mathematical models along with future research directions to understand the unique myocardium-coronary coupling in the heart. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.
Collapse
Affiliation(s)
- Lei Fan
- Joint Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Haifeng Wang
- Department of Mechanical Engineering, Michigan State University, East Lansing, Michigan, USA
| | - Ghassan S Kassab
- California Medical Innovations Institute, San Diego, California, USA
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, Michigan, USA
| |
Collapse
|
2
|
Kachabi A, Colebank MJ, Chesler NC. Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension. Biomech Model Mechanobiol 2024; 23:469-483. [PMID: 38017302 PMCID: PMC10963496 DOI: 10.1007/s10237-023-01786-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 10/21/2023] [Indexed: 11/30/2023]
Abstract
Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstructs flow and increases pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the small pulmonary arteries due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we applied a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) approach to investigate the impact of CTEPH on pulmonary artery stiffening, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) in extralobar (main, right, and left) pulmonary arteries and intralobar (distal to the extralobar) arteries. Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar and intralobar arteries. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI, quantified by the novel parameter φ , which is related to thrombogenicity. Our analysis reveals a strong positive correlation between increases in mean pulmonary artery pressure (mPAP) and φ from baseline to CTEPH in individual subjects, which supports the suggestion that increased φ drives disease severity. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.
Collapse
Affiliation(s)
- Amirreza Kachabi
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, USA
| | - Mitchel J Colebank
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, USA
| | - Naomi C Chesler
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, USA.
| |
Collapse
|
3
|
Ding CCA, Dokos S, Bakir AA, Zamberi NJ, Liew YM, Chan BT, Md Sari NA, Avolio A, Lim E. Simulating impaired left ventricular-arterial coupling in aging and disease: a systematic review. Biomed Eng Online 2024; 23:24. [PMID: 38388416 PMCID: PMC10885508 DOI: 10.1186/s12938-024-01206-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Accepted: 01/11/2024] [Indexed: 02/24/2024] Open
Abstract
Aortic stenosis, hypertension, and left ventricular hypertrophy often coexist in the elderly, causing a detrimental mismatch in coupling between the heart and vasculature known as ventricular-vascular (VA) coupling. Impaired left VA coupling, a critical aspect of cardiovascular dysfunction in aging and disease, poses significant challenges for optimal cardiovascular performance. This systematic review aims to assess the impact of simulating and studying this coupling through computational models. By conducting a comprehensive analysis of 34 relevant articles obtained from esteemed databases such as Web of Science, Scopus, and PubMed until July 14, 2022, we explore various modeling techniques and simulation approaches employed to unravel the complex mechanisms underlying this impairment. Our review highlights the essential role of computational models in providing detailed insights beyond clinical observations, enabling a deeper understanding of the cardiovascular system. By elucidating the existing models of the heart (3D, 2D, and 0D), cardiac valves, and blood vessels (3D, 1D, and 0D), as well as discussing mechanical boundary conditions, model parameterization and validation, coupling approaches, computer resources and diverse applications, we establish a comprehensive overview of the field. The descriptions as well as the pros and cons on the choices of different dimensionality in heart, valve, and circulation are provided. Crucially, we emphasize the significance of evaluating heart-vessel interaction in pathological conditions and propose future research directions, such as the development of fully coupled personalized multidimensional models, integration of deep learning techniques, and comprehensive assessment of confounding effects on biomarkers.
Collapse
Affiliation(s)
- Corina Cheng Ai Ding
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, 50603, Kuala Lumpur, Malaysia
- Graduate School of Biomedical Engineering, Faculty of Engineering, University of New South Wales, Sydney, NSW, 2052, Australia
| | - Socrates Dokos
- Graduate School of Biomedical Engineering, Faculty of Engineering, University of New South Wales, Sydney, NSW, 2052, Australia
| | - Azam Ahmad Bakir
- University of Southampton Malaysia Campus, 79200, Iskandar Puteri, Johor, Malaysia
| | - Nurul Jannah Zamberi
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, 50603, Kuala Lumpur, Malaysia
- Graduate School of Biomedical Engineering, Faculty of Engineering, University of New South Wales, Sydney, NSW, 2052, Australia
| | - Yih Miin Liew
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, 50603, Kuala Lumpur, Malaysia
| | - Bee Ting Chan
- Department of Mechanical, Materials and Manufacturing Engineering, Faculty of Science and Engineering, University of Nottingham Malaysia, 43500, Selangor, Malaysia
| | - Nor Ashikin Md Sari
- Department of Medicine, Faculty of Medicine, University of Malaya, 50603, Kuala Lumpur, Malaysia
| | - Alberto Avolio
- Macquarie Medical School, Faculty of Medicine, Health and Human Sciences, Macquarie University, Sydney, NSW, 2109, Australia
| | - Einly Lim
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, 50603, Kuala Lumpur, Malaysia.
| |
Collapse
|
4
|
Kachabi A, Colebank MJ, Chesler N. Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension. RESEARCH SQUARE 2023:rs.3.rs-3214385. [PMID: 37577616 PMCID: PMC10418554 DOI: 10.21203/rs.3.rs-3214385/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstruct flow and increase pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the pulmonary vasculature due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we developed a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) models to investigate the impact of CTEPH on pulmonary artery stiffening, time averaged wall shear stress (TAWSS), and oscillatory shear index (OSI). Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar (main, right and left pulmonary arteries) and intralobar vessels. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression, including the combination of low mean wall shear stress with high oscillation. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.
Collapse
|
5
|
Uncertainty Quantification in the In Vivo Image-Based Estimation of Local Elastic Properties of Vascular Walls. J Cardiovasc Dev Dis 2023; 10:jcdd10030109. [PMID: 36975873 PMCID: PMC10058982 DOI: 10.3390/jcdd10030109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2023] [Revised: 02/15/2023] [Accepted: 03/02/2023] [Indexed: 03/08/2023] Open
Abstract
Introduction: Patient-specific computational models are a powerful tool for planning cardiovascular interventions. However, the in vivo patient-specific mechanical properties of vessels represent a major source of uncertainty. In this study, we investigated the effect of uncertainty in the elastic module (E) on a Fluid–Structure Interaction (FSI) model of a patient-specific aorta. Methods: The image-based χ-method was used to compute the initial E value of the vascular wall. The uncertainty quantification was carried out using the generalized Polynomial Chaos (gPC) expansion technique. The stochastic analysis was based on four deterministic simulations considering four quadrature points. A deviation of about ±20% on the estimation of the E value was assumed. Results: The influence of the uncertain E parameter was evaluated along the cardiac cycle on area and flow variations extracted from five cross-sections of the aortic FSI model. Results of stochastic analysis showed the impact of E in the ascending aorta while an insignificant effect was observed in the descending tract. Conclusions: This study demonstrated the importance of the image-based methodology for inferring E, highlighting the feasibility of retrieving useful additional data and enhancing the reliability of in silico models in clinical practice.
Collapse
|
6
|
Marcinno’ F, Zingaro A, Fumagalli I, Dede’ L, Vergara C. A Computational Study of Blood Flow Dynamics in the Pulmonary Arteries. VIETNAM JOURNAL OF MATHEMATICS 2022; 51:127-149. [PMID: 36536831 PMCID: PMC9750052 DOI: 10.1007/s10013-022-00595-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/28/2021] [Accepted: 08/29/2022] [Indexed: 06/17/2023]
Abstract
In this work we study the blood dynamics in the pulmonary arteries by means of a 3D-0D geometric multiscale approach, where a detailed 3D model for the pulmonary arteries is coupled with a lumped parameters (0D) model of the cardiovascular system. We propose to investigate three strategies for the numerical solution of the 3D-0D coupled problem: the Splitting-Explicit and Implicit algorithms, where information are exchanged between 3D and 0D models at each time step at the interfaces, and the One-Way algorithm, where the 0D is solved first off-line. In our numerical experiments performed in a realistic patient-specific 3D domain with a physiologically calibrated 0D model, we discuss first the issue on instabilities that may arise when not suitable connections are considered between 3D and 0D models; second we compare the performance and accuracy of the three proposed numerical strategies. Finally, we report a comparison between a healthy and a hypertensive case, providing a preliminary result highlighting how our method could be used in future for clinical purposes.
Collapse
Affiliation(s)
- Fabio Marcinno’
- Dipartimento di Chimica, Materiali e Ingegneria Chimica “Giulio Natta”, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milan, 20133 Italy
| | - Alberto Zingaro
- MOX, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milan, 20133 Italy
| | - Ivan Fumagalli
- MOX, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milan, 20133 Italy
| | - Luca Dede’
- MOX, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milan, 20133 Italy
| | - Christian Vergara
- Dipartimento di Chimica, Materiali e Ingegneria Chimica “Giulio Natta”, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milan, 20133 Italy
| |
Collapse
|
7
|
Kamada H, Nakamura M, Ota H, Higuchi S, Takase K. Blood flow analysis with computational fluid dynamics and 4D-flow MRI for vascular diseases. J Cardiol 2022; 80:386-396. [PMID: 35718672 DOI: 10.1016/j.jjcc.2022.05.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Accepted: 05/16/2022] [Indexed: 10/31/2022]
Abstract
Both computational fluid dynamics (CFD) and time-resolved, three-dimensional, phase-contrast, magnetic resonance imaging (4D-flow MRI) enable visualization of time-varying blood flow structures and quantification of blood flow in vascular diseases. However, they are totally different. CFD is a method to calculate blood flow by solving the governing equations of fluid mechanics, so the obtained flow field is somewhat virtual. On the other hand, 4D-flow MRI measures blood flow in vivo, thus the flow is real. Recently, with the development and enhancement of computers, medical imaging techniques, and related software, blood flow analysis has become more accessible to clinicians and its usefulness in vascular diseases has been demonstrated. In this review, we have outlined the methods and characteristics of CFD and 4D-flow MRI, respectively. We have discussed the differences in the characteristics between both methods; reviewed the milestones achieved by blood flow analysis in various vascular diseases; and discussed the usefulness, challenges, and limitations of blood flow analysis. We have discussed the difficulties and limitations of current blood flow analysis. We have also discussed our views on future directions.
Collapse
Affiliation(s)
- Hiroki Kamada
- Department of Diagnostic Radiology, Tohoku University Hospital, Sendai, Japan.
| | - Masanori Nakamura
- Department of Electrical and Mechanical Engineering, Nagoya Institute of Technology, Nagoya, Japan
| | - Hideki Ota
- Department of Diagnostic Radiology, Tohoku University Hospital, Sendai, Japan
| | - Satoshi Higuchi
- Department of Diagnostic Radiology, Tohoku University Hospital, Sendai, Japan
| | - Kei Takase
- Department of Diagnostic Radiology, Tohoku University Hospital, Sendai, Japan
| |
Collapse
|
8
|
Bracamonte JH, Saunders SK, Wilson JS, Truong UT, Soares JS. Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications. APPLIED SCIENCES-BASEL 2022; 12:3954. [PMID: 36911244 PMCID: PMC10004130 DOI: 10.3390/app12083954] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.
Collapse
Affiliation(s)
- Johane H. Bracamonte
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA
| | - Sarah K. Saunders
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA
| | - John S. Wilson
- Department of Biomedical Engineering and Pauley Heart Center, Virginia Commonwealth University, Richmond, VA 23219, USA
| | - Uyen T. Truong
- Department of Pediatrics, School of Medicine, Children’s Hospital of Richmond at Virginia Commonwealth University, Richmond, VA 23219, USA
| | - Joao S. Soares
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA
- Correspondence:
| |
Collapse
|
9
|
Pillalamarri NR, Piskin S, Patnaik SS, Murali S, Finol EA. Patient-Specific Computational Analysis of Hemodynamics in Adult Pulmonary Hypertension. Ann Biomed Eng 2021; 49:3465-3480. [PMID: 34799807 DOI: 10.1007/s10439-021-02884-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Accepted: 10/29/2021] [Indexed: 11/25/2022]
Abstract
Pulmonary hypertension (PH) is a progressive disease characterized by elevated pressure and vascular resistance in the pulmonary arteries. Nearly 250,000 hospitalizations occur annually in the US with PH as the primary or secondary condition. A definitive diagnosis of PH requires right heart catheterization (RHC) in addition to a chest computed tomography, a walking test, and others. While RHC is the gold standard for diagnosing PH, it is invasive and posseses inherent risks and contraindications. In this work, we characterized the patient-specific pulmonary hemodynamics in silico for diverse PH WHO groups. We grouped patients on the basis of mean pulmonary arterial pressure (mPAP) into three disease severity groups: at-risk ([Formula: see text], denoted with A), mild ([Formula: see text], denoted with M), and severe ([Formula: see text], denoted with S). The pulsatile flow hemodynamics was simulated by evaluating the three-dimensional Navier-Stokes system of equations using a flow solver developed by customizing OpenFOAM libraries (v5.0, The OpenFOAM Foundation). Quasi patient-specific boundary conditions were implemented using a Womersley inlet velocity profile and transient resistance outflow conditions. Hemodynamic indices such as spatially averaged wall shear stress ([Formula: see text]), wall shear stress gradient ([Formula: see text]), time-averaged wall shear stress ([Formula: see text]), oscillatory shear index ([Formula: see text]), and relative residence time ([Formula: see text]), were evaluated along with the clinical metrics pulmonary vascular resistance ([Formula: see text]), stroke volume ([Formula: see text]) and compliance ([Formula: see text]), to assess possible spatiotemporal correlations. We observed statistically significant decreases in [Formula: see text], [Formula: see text], and [Formula: see text], and increases in [Formula: see text] and [Formula: see text] with disease severity. [Formula: see text] was moderately correlated with [Formula: see text] and [Formula: see text] at the mid-notch stage of the cardiac cycle when these indices were computed using the global pulmonary arterial geometry. These results are promising in the context of a long-term goal of identifying computational biomarkers that can serve as surrogates for invasive diagnostic protocols of PH.
Collapse
Affiliation(s)
- Narasimha R Pillalamarri
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
| | - Senol Piskin
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
- Department of Mechanical Engineering, Istinye University, Istanbul, Turkey
| | - Sourav S Patnaik
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
- Department of Bioengineering, University of Texas at Dallas, Richardson, TX, USA
| | - Srinivas Murali
- Division of Cardiovascular Medicine, Allegheny Health Network, Pittsburgh, PA, USA
| | - Ender A Finol
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA.
- UTSA/UTHSA Joint Graduate Program in Biomedical Engineering, University of Texas at San Antonio, San Antonio, TX, USA.
| |
Collapse
|