1
|
Colebank MJ, Oomen PA, Witzenburg CM, Grosberg A, Beard DA, Husmeier D, Olufsen MS, Chesler NC. Guidelines for mechanistic modeling and analysis in cardiovascular research. Am J Physiol Heart Circ Physiol 2024; 327:H473-H503. [PMID: 38904851 DOI: 10.1152/ajpheart.00766.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 06/07/2024] [Accepted: 06/16/2024] [Indexed: 06/22/2024]
Abstract
Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.
Collapse
Affiliation(s)
- Mitchel J Colebank
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States
| | - Pim A Oomen
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States
| | - Colleen M Witzenburg
- Department of Biomedical Engineering, University of Wisconsin-Madison, Madison, Wisconsin, United States
| | - Anna Grosberg
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States
| | - Daniel A Beard
- Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, Michigan, United States
| | - Dirk Husmeier
- School of Mathematics and Statistics, University of Glasgow, Glasgow, United Kingdom
| | - Mette S Olufsen
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina, United States
| | - Naomi C Chesler
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States
| |
Collapse
|
2
|
Shim YD, Chen MC, Ha S, Chang HJ, Baek S, Lee EH. Multi-scaled temporal modeling of cardiovascular disease progression: An illustration of proximal arteries in pulmonary hypertension. J Biomech 2024; 168:112059. [PMID: 38631187 PMCID: PMC11096051 DOI: 10.1016/j.jbiomech.2024.112059] [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: 11/19/2023] [Revised: 03/16/2024] [Accepted: 03/21/2024] [Indexed: 04/19/2024]
Abstract
The progression of cardiovascular disease is intricately influenced by a complex interplay between physiological pathways, biochemical processes, and physical mechanisms. This study aimed to develop an in-silico physics-based approach to comprehensively model the multifaceted vascular pathophysiological adaptations. This approach focused on capturing the progression of proximal pulmonary arterial hypertension, which is significantly associated with the irreversible degradation of arterial walls and compensatory stress-induced growth and remodeling. This study incorporated critical characteristics related to the distinct time scales for the deformation, thus reflecting the impact of mean pressure on artery growth and tissue damage. The in-silico simulation of the progression of pulmonary hypertension was realized based on computational code combined with the finite element method (FEM) for the simulation of disease progression. The parametric studies further explored the consequences of these irreversible processes. This computational modeling approach may advance our understanding of pulmonary hypertension and its progression.
Collapse
Affiliation(s)
- Young-Dae Shim
- Department of Smart Fabrication Technology, Sungkyunkwan University, Suwon-si, Gyeonggi-do 16419, Republic of Korea.
| | - Mei-Cen Chen
- Department of Smart Fabrication Technology, Sungkyunkwan University, Suwon-si, Gyeonggi-do 16419, Republic of Korea.
| | - Seongmin Ha
- Biomedical Engineering, Yonsei University College of Medicine 250, Seoul, Republic of Korea.
| | - Hyuk-Jae Chang
- Division of Cardiology, Yonsei University College of Medicine 250, Seoul, Republic of Korea.
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI 48824, United States.
| | - Eun-Ho Lee
- Department of Smart Fabrication Technology, Sungkyunkwan University, Suwon-si, Gyeonggi-do 16419, Republic of Korea; School of Mechanical Engineering, Sungkyunkwan University, Suwon-si, Gyeonggi-do 16419, Republic of Korea; Department of Intelligent Robotics, Sungkyunkwan University, Suwon-si, Gyeonggi-do 16419, Republic of Korea.
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Haider MA, Pearce KJ, Chesler NC, Hill NA, Olufsen MS. Application and reduction of a nonlinear hyperelastic wall model capturing ex vivo relationships between fluid pressure, area, and wall thickness in normal and hypertensive murine left pulmonary arteries. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2024; 40:e3798. [PMID: 38214099 DOI: 10.1002/cnm.3798] [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: 11/02/2022] [Revised: 08/10/2023] [Accepted: 11/26/2023] [Indexed: 01/13/2024]
Abstract
Pulmonary hypertension is a cardiovascular disorder manifested by elevated mean arterial blood pressure (>20 mmHg) together with vessel wall stiffening and thickening due to alterations in collagen, elastin, and smooth muscle cells. Hypoxia-induced (type 3) pulmonary hypertension can be studied in animals exposed to a low oxygen environment for prolonged time periods leading to biomechanical alterations in vessel wall structure. This study introduces a novel approach to formulating a reduced order nonlinear elastic structural wall model for a large pulmonary artery. The model relating blood pressure and area is calibrated using ex vivo measurements of vessel diameter and wall thickness changes, under controlled pressure conditions, in left pulmonary arteries isolated from control and hypertensive mice. A two-layer, hyperelastic, and anisotropic model incorporating residual stresses is formulated using the Holzapfel-Gasser-Ogden model. Complex relations predicting vessel area and wall thickness with increasing blood pressure are derived and calibrated using the data. Sensitivity analysis, parameter estimation, subset selection, and physical plausibility arguments are used to systematically reduce the 16-parameter model to one in which a much smaller subset of identifiable parameters is estimated via solution of an inverse problem. Our final reduced one layer model includes a single set of three elastic moduli. Estimated ranges of these parameters demonstrate that nonlinear stiffening is dominated by elastin in the control animals and by collagen in the hypertensive animals. The pressure-area relation developed in this novel manner has potential impact on one-dimensional fluids network models of vessel wall remodeling in the presence of cardiovascular disease.
Collapse
Affiliation(s)
- Mansoor A Haider
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina, USA
| | - Katherine J Pearce
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina, USA
| | - Naomi C Chesler
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center & Department of Biomedical Engineering, University of California, Irvine (UCI), Irvine, California, USA
| | - Nicholas A Hill
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Mette S Olufsen
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina, USA
| |
Collapse
|
6
|
Szafron JM, Yang W, Feinstein JA, Rabinovitch M, Marsden AL. A computational growth and remodeling framework for adaptive and maladaptive pulmonary arterial hemodynamics. Biomech Model Mechanobiol 2023; 22:1935-1951. [PMID: 37658985 PMCID: PMC10929588 DOI: 10.1007/s10237-023-01744-z] [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: 04/19/2023] [Accepted: 07/05/2023] [Indexed: 09/05/2023]
Abstract
Hemodynamic loading is known to contribute to the development and progression of pulmonary arterial hypertension (PAH). This loading drives changes in mechanobiological stimuli that affect cellular phenotypes and lead to pulmonary vascular remodeling. Computational models have been used to simulate mechanobiological metrics of interest, such as wall shear stress, at single time points for PAH patients. However, there is a need for new approaches that simulate disease evolution to allow for prediction of long-term outcomes. In this work, we develop a framework that models the pulmonary arterial tree through adaptive and maladaptive responses to mechanical and biological perturbations. We coupled a constrained mixture theory-based growth and remodeling framework for the vessel wall with a morphometric tree representation of the pulmonary arterial vasculature. We show that non-uniform mechanical behavior is important to establish the homeostatic state of the pulmonary arterial tree, and that hemodynamic feedback is essential for simulating disease time courses. We also employed a series of maladaptive constitutive models, such as smooth muscle hyperproliferation and stiffening, to identify critical contributors to development of PAH phenotypes. Together, these simulations demonstrate an important step toward predicting changes in metrics of clinical interest for PAH patients and simulating potential treatment approaches.
Collapse
Affiliation(s)
- Jason M Szafron
- Department of Pediatrics (Cardiology), Stanford University, Stanford, CA, 94305, USA
- Cardiovascular Institute, Stanford University, Palo Alto, CA, 94305, USA
| | - Weiguang Yang
- Department of Pediatrics (Cardiology), Stanford University, Stanford, CA, 94305, USA
| | - Jeffrey A Feinstein
- Department of Pediatrics (Cardiology), Stanford University, Stanford, CA, 94305, USA
- Cardiovascular Institute, Stanford University, Palo Alto, CA, 94305, USA
| | - Marlene Rabinovitch
- Department of Pediatrics (Cardiology), Stanford University, Stanford, CA, 94305, USA
- Cardiovascular Institute, Stanford University, Palo Alto, CA, 94305, USA
| | - Alison L Marsden
- Department of Pediatrics (Cardiology), Stanford University, Stanford, CA, 94305, USA.
- Cardiovascular Institute, Stanford University, Palo Alto, CA, 94305, USA.
| |
Collapse
|
7
|
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
|
8
|
Liang L, Liu M, Elefteriades J, Sun W. PyTorch-FEA: Autograd-enabled finite element analysis methods with applications for biomechanical analysis of human aorta. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 238:107616. [PMID: 37230048 PMCID: PMC10330852 DOI: 10.1016/j.cmpb.2023.107616] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 05/15/2023] [Accepted: 05/16/2023] [Indexed: 05/27/2023]
Abstract
BACKGROUND AND OBJECTIVES Finite-element analysis (FEA) is widely used as a standard tool for stress and deformation analysis of solid structures, including human tissues and organs. For instance, FEA can be applied at a patient-specific level to assist in medical diagnosis and treatment planning, such as risk assessment of thoracic aortic aneurysm rupture/dissection. These FEA-based biomechanical assessments often involve both forward and inverse mechanics problems. Current commercial FEA software packages (e.g., Abaqus) and inverse methods exhibit performance issues in either accuracy or speed. METHODS In this study, we propose and develop a new library of FEA code and methods, named PyTorch-FEA, by taking advantage of autograd, an automatic differentiation mechanism in PyTorch. We develop a class of PyTorch-FEA functionalities to solve forward and inverse problems with improved loss functions, and we demonstrate the capability of PyTorch-FEA in a series of applications related to human aorta biomechanics. In one of the inverse methods, we combine PyTorch-FEA with deep neural networks (DNNs) to further improve performance. RESULTS We applied PyTorch-FEA in four fundamental applications for biomechanical analysis of human aorta. In the forward analysis, PyTorch-FEA achieved a significant reduction in computational time without compromising accuracy compared with Abaqus, a commercial FEA package. Compared to other inverse methods, inverse analysis with PyTorch-FEA achieves better performance in either accuracy or speed, or both if combined with DNNs. CONCLUSIONS We have presented PyTorch-FEA, a new library of FEA code and methods, representing a new approach to develop FEA methods to forward and inverse problems in solid mechanics. PyTorch-FEA eases the development of new inverse methods and enables a natural integration of FEA and DNNs, which will have numerous potential applications.
Collapse
Affiliation(s)
- Liang Liang
- Department of Computer Science, University of Miami, Coral Gables, FL, United States.
| | - Minliang Liu
- Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, United States
| | - John Elefteriades
- Aortic Institute, School of Medicine, Yale University, New Haven, CT, United States
| | - Wei Sun
- Sutra Medical Inc, Lake Forest, CA, United States
| |
Collapse
|
9
|
Gharahi H, Filonova V, Mullagura HN, Nama N, Baek S, Figueroa CA. A multiscale framework for defining homeostasis in distal vascular trees: applications to the pulmonary circulation. Biomech Model Mechanobiol 2023; 22:971-986. [PMID: 36917305 DOI: 10.1007/s10237-023-01693-7] [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: 08/10/2021] [Accepted: 01/11/2023] [Indexed: 03/16/2023]
Abstract
Pulmonary arteries constitute a low-pressure network of vessels, often characterized as a bifurcating tree with heterogeneous vessel mechanics. Understanding the vascular complexity and establishing homeostasis is important to study diseases such as pulmonary arterial hypertension (PAH). The onset and early progression of PAH can be traced to changes in the morphometry and structure of the distal vasculature. Coupling hemodynamics with vessel wall growth and remodeling (G&R) is crucial for understanding pathology at distal vasculature. Accordingly, the goal of this study is to provide a multiscale modeling framework that embeds the essential features of arterial wall constituents coupled with the hemodynamics within an arterial network characterized by an extension of Murray's law. This framework will be used to establish the homeostatic baseline characteristics of a pulmonary arterial tree, including important parameters such as vessel radius, wall thickness and shear stress. To define the vascular homeostasis and hemodynamics in the tree, we consider two timescales: a cardiac cycle and a longer period of vascular adaptations. An iterative homeostatic optimization, which integrates a metabolic cost function minimization, the stress equilibrium, and hemodynamics, is performed at the slow timescale. In the fast timescale, the pulsatile blood flow dynamics is described by a Womersley's deformable wall analytical solution. Illustrative examples for symmetric and asymmetric trees are presented that provide baseline characteristics for the normal pulmonary arterial vasculature. The results are compared with diverse literature data on morphometry, structure, and mechanics of pulmonary arteries. The developed framework demonstrates a potential for advanced parametric studies and future G&R and hemodynamics modeling of PAH.
Collapse
Affiliation(s)
- Hamidreza Gharahi
- Section of Vascular Surgery, Department of Surgery, University of Michigan, Ann Arbor, MI, USA.
| | - Vasilina Filonova
- Section of Vascular Surgery, Department of Surgery, University of Michigan, Ann Arbor, MI, USA
| | - Haritha N Mullagura
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA
| | - Nitesh Nama
- Department of Mechanical & Materials Engineering, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA
| | - C Alberto Figueroa
- Section of Vascular Surgery, Department of Surgery, University of Michigan, Ann Arbor, MI, USA
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
10
|
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
|
11
|
Tossas-Betancourt C, Li NY, Shavik SM, Afton K, Beckman B, Whiteside W, Olive MK, Lim HM, Lu JC, Phelps CM, Gajarski RJ, Lee S, Nordsletten DA, Grifka RG, Dorfman AL, Baek S, Lee LC, Figueroa CA. Data-driven computational models of ventricular-arterial hemodynamics in pediatric pulmonary arterial hypertension. Front Physiol 2022; 13:958734. [PMID: 36160862 PMCID: PMC9490558 DOI: 10.3389/fphys.2022.958734] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 08/01/2022] [Indexed: 11/13/2022] Open
Abstract
Pulmonary arterial hypertension (PAH) is a complex disease involving increased resistance in the pulmonary arteries and subsequent right ventricular (RV) remodeling. Ventricular-arterial interactions are fundamental to PAH pathophysiology but are rarely captured in computational models. It is important to identify metrics that capture and quantify these interactions to inform our understanding of this disease as well as potentially facilitate patient stratification. Towards this end, we developed and calibrated two multi-scale high-resolution closed-loop computational models using open-source software: a high-resolution arterial model implemented using CRIMSON, and a high-resolution ventricular model implemented using FEniCS. Models were constructed with clinical data including non-invasive imaging and invasive hemodynamic measurements from a cohort of pediatric PAH patients. A contribution of this work is the discussion of inconsistencies in anatomical and hemodynamic data routinely acquired in PAH patients. We proposed and implemented strategies to mitigate these inconsistencies, and subsequently use this data to inform and calibrate computational models of the ventricles and large arteries. Computational models based on adjusted clinical data were calibrated until the simulated results for the high-resolution arterial models matched within 10% of adjusted data consisting of pressure and flow, whereas the high-resolution ventricular models were calibrated until simulation results matched adjusted data of volume and pressure waveforms within 10%. A statistical analysis was performed to correlate numerous data-derived and model-derived metrics with clinically assessed disease severity. Several model-derived metrics were strongly correlated with clinically assessed disease severity, suggesting that computational models may aid in assessing PAH severity.
Collapse
Affiliation(s)
| | - Nathan Y. Li
- Department of Pharmaceutical Sciences, University of Michigan, Ann Arbor, MI, United States
| | - Sheikh M. Shavik
- Department of Mechanical Engineering, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh
| | - Katherine Afton
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Brian Beckman
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - Wendy Whiteside
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Mary K. Olive
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Heang M. Lim
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Jimmy C. Lu
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Christina M. Phelps
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - Robert J. Gajarski
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - Simon Lee
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - David A. Nordsletten
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
- Department of Surgery, University of Michigan, Ann Arbor, MI, United States
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Ronald G. Grifka
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Adam L. Dorfman
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - C. Alberto Figueroa
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
- Department of Surgery, University of Michigan, Ann Arbor, MI, United States
| |
Collapse
|
12
|
A patient-specific image-based approach to estimate pulmonary artery stiffness based on vessel constitutive model. Med Eng Phys 2022; 107:103851. [DOI: 10.1016/j.medengphy.2022.103851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 06/28/2022] [Accepted: 07/10/2022] [Indexed: 11/21/2022]
|
13
|
Shahid L, Rice J, Berhane H, Rigsby C, Robinson J, Griffin L, Markl M, Roldán-Alzate A. Enhanced 4D Flow MRI-Based CFD with Adaptive Mesh Refinement for Flow Dynamics Assessment in Coarctation of the Aorta. Ann Biomed Eng 2022; 50:1001-1016. [PMID: 35624334 PMCID: PMC11034844 DOI: 10.1007/s10439-022-02980-7] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 05/11/2022] [Indexed: 01/28/2023]
Abstract
4D Flow MRI is a diagnostic tool that can visualize and quantify patient-specific hemodynamics and help interventionalists optimize treatment strategies for repairing coarctation of the aorta (COA). Despite recent developments in 4D Flow MRI, shortcomings include phase-offset errors, limited spatiotemporal resolution, aliasing, inaccuracies due to slow aneurysmal flows, and distortion of images due to metallic artifact from vascular stents. To address these limitations, we developed a framework utilizing Computational Fluid Dynamics (CFD) with Adaptive Mesh Refinement (AMR) that enhances 4D Flow MRI visualization/quantification. We applied this framework to five pediatric patients with COA, providing in-vivo and in-silico datasets, pre- and post-intervention. These two data sets were compared and showed that CFD flow rates were within 9.6% of 4D Flow MRI, which is within a clinically acceptable range. CFD simulated slow aneurysmal flow, which MRI failed to capture due to high relative velocity encoding (Venc). CFD successfully predicted in-stent blood flow, which was not visible in the in-vivo data due to susceptibility artifact. AMR improved spatial resolution by factors of 101 to 103 and temporal resolution four-fold. This computational framework has strong potential to optimize visualization/quantification of aneurysmal and in-stent flows, improve spatiotemporal resolution, and assess hemodynamic efficiency post-COA treatment.
Collapse
Affiliation(s)
- Labib Shahid
- Department of Mechanical Engineering, University of Wisconsin-Madison, 1111 Highland Ave, Room 2476 WIMR II, Madison, WI, 53705, USA.
| | - James Rice
- Department of Mechanical Engineering, University of Wisconsin-Madison, 1111 Highland Ave, Room 2476 WIMR II, Madison, WI, 53705, USA
| | - Haben Berhane
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, USA
| | - Cynthia Rigsby
- Department of Medical Imaging, Ann & Robert H. Lurie Children's Hospital of Chicago, Chicago, IL, USA
| | - Joshua Robinson
- Department of Medical Imaging, Ann & Robert H. Lurie Children's Hospital of Chicago, Chicago, IL, USA
| | - Lindsay Griffin
- Department of Medical Imaging, Ann & Robert H. Lurie Children's Hospital of Chicago, Chicago, IL, USA
| | - Michael Markl
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, USA
| | - Alejandro Roldán-Alzate
- Department of Mechanical Engineering, University of Wisconsin-Madison, 1111 Highland Ave, Room 2476 WIMR II, Madison, WI, 53705, USA
- Department of Biomedical Engineering, University of Wisconsin-Madison, Madison, WI, USA
- Department of Radiology, University of Wisconsin-Madison, Madison, WI, USA
| |
Collapse
|
14
|
Introduction of a Novel Image-Based and Non-Invasive Method for the Estimation of Local Elastic Properties of Great Vessels. ELECTRONICS 2022. [DOI: 10.3390/electronics11132055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Background: In the context of a growing demand for the use of in silico models to meet clinical requests, image-based methods play a crucial role. In this study, we present a parametric equation able to estimate the elasticity of vessel walls, non-invasively and indirectly, from information uniquely retrievable from imaging. Methods: A custom equation was iteratively refined and tuned from the simulations of a wide range of different vessel models, leading to the definition of an indirect method able to estimate the elastic modulus E of a vessel wall. To test the effectiveness of the predictive capability to infer the E value, two models with increasing complexity were used: a U-shaped vessel and a patient-specific aorta. Results: The original formulation was demonstrated to deviate from the ground truth, with a difference of 89.6%. However, the adoption of our proposed equation was found to significantly increase the reliability of the estimated E value for a vessel wall, with a mean percentage error of 9.3% with respect to the reference values. Conclusion: This study provides a strong basis for the definition of a method able to estimate local mechanical information of vessels from data easily retrievable from imaging, thus potentially increasing the reliability of in silico cardiovascular models.
Collapse
|
15
|
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
|
16
|
Bracamonte JH, Wilson JS, Soares JS. Quantification of the heterogeneous effect of static and dynamic perivascular structures on patient-specific local aortic wall mechanics using inverse finite element modeling and DENSE MRI. J Biomech 2022; 138:111119. [PMID: 35576631 PMCID: PMC9536506 DOI: 10.1016/j.jbiomech.2022.111119] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Revised: 04/28/2022] [Accepted: 04/29/2022] [Indexed: 11/30/2022]
Abstract
Recent studies have highlighted the relevance of perivascular interactions on aortic wall mechanics. Most of the approaches assume static perivascular structures; however, the beating heart dynamically displaces the neighboring aorta. We develop a model to account for the effect of periaortic interactions due to static and dynamic structures by prescribing a moving elastic foundation boundary condition (EFBC) embedded into an inverse finite element algorithm using in vivo displacements from 2D displacement encoding with stimulated echoes (DENSE) MRI as target data. We applied this method at three different locations of interest, the distal aortic arch (DAA), descending thoracic aorta (DTA), and infrarenal abdominal aorta (IAA) for a total of 27 cases in healthy humans. The model reproduces the target diastole-to-systole deformation and bulk displacement of the aortic wall with median displacement errors below 0.5mm. The EFBC showed good agreement with the location of anatomical features and was consistent among individuals of similar characteristics. Results show that an energy source acting on the adventitia is required to reproduce the displacements measured at the vicinity of the heart, but not at the abdomen. The average adventitial load as a percentage of the luminal pulse-pressure was found to increase with age and to decrease along the descending aorta, from 61% at the DAA to 37% at the DTA, and 30% at the IAA. This approach offers a patient-specific method to estimate in vivo adventitial loads and aortic wall stiffness, which can bring a better understanding of normal and pathological in vivo aortic function.
Collapse
Affiliation(s)
- Johane H Bracamonte
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, Richmond, VA 23284, United States.
| | - John S Wilson
- Department of Biomedical Engineering & Pauley Heart Center, Virginia Commonwealth University, 601 West Main Street, Richmond, VA 23284, United States.
| | - Joao S Soares
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, Richmond, VA 23284, United States.
| |
Collapse
|
17
|
He Y, Northrup H, Le H, Cheung AK, Berceli SA, Shiu YT. Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases. Front Bioeng Biotechnol 2022; 10:855791. [PMID: 35573253 PMCID: PMC9091352 DOI: 10.3389/fbioe.2022.855791] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2022] [Accepted: 04/08/2022] [Indexed: 01/17/2023] Open
Abstract
Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.
Collapse
Affiliation(s)
- Yong He
- Division of Vascular Surgery and Endovascular Therapy, University of Florida, Gainesville, FL, United States
| | - Hannah Northrup
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, United States
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
| | - Ha Le
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
| | - Alfred K. Cheung
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
- Veterans Affairs Salt Lake City Healthcare System, Salt Lake City, UT, United States
| | - Scott A. Berceli
- Division of Vascular Surgery and Endovascular Therapy, University of Florida, Gainesville, FL, United States
- Vascular Surgery Section, Malcom Randall Veterans Affairs Medical Center, Gainesville, FL, United States
| | - Yan Tin Shiu
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
- Veterans Affairs Salt Lake City Healthcare System, Salt Lake City, UT, United States
- *Correspondence: Yan Tin Shiu,
| |
Collapse
|
18
|
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
|
19
|
Schiebler ML, Fain S, van Beek E. Vascular imaging of the lung: perspectives on current imaging methods. Br J Radiol 2022; 95:20200759. [PMID: 32795175 PMCID: PMC9153715 DOI: 10.1259/bjr.20200759] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
This commentary will discuss the use of advanced non-invasive imaging methodology for the pulmonary vascular system with special attention to a rubric for the imaging and clinical team to use for any particular clinical situation.
Collapse
Affiliation(s)
- Mark L. Schiebler
- Department of Radiology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Sean Fain
- Department of Medical Physics, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Edwin van Beek
- SINAPSE Chair of Clinical Radiology, Director Edinburgh Imaging facility QMRI, Honorary Consultant Radiologist, NHS Lothian Health Board, CO.19 Edinburgh Imaging, Queen’s Medical Research Institute, Edinburgh Bioquarter, 47 Little France Crescent, Edinburgh EH16 4TJ, United Kingdom
| |
Collapse
|
20
|
Bartolo MA, Qureshi MU, Colebank MJ, Chesler NC, Olufsen MS. Numerical predictions of shear stress and cyclic stretch in pulmonary hypertension due to left heart failure. Biomech Model Mechanobiol 2022; 21:363-381. [PMID: 35037114 PMCID: PMC10174274 DOI: 10.1007/s10237-021-01538-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Accepted: 11/24/2021] [Indexed: 12/15/2022]
Abstract
Isolated post-capillary pulmonary hypertension (Ipc-PH) occurs due to left heart failure, which contributes to 1 out of every 9 deaths in the United States. In some patients, through unknown mechanisms, Ipc-PH transitions to combined pre-/post-capillary PH (Cpc-PH) and is associated with a dramatic increase in mortality. Altered mechanical forces and subsequent biological signaling in the pulmonary vascular bed likely contribute to the transition from Ipc-PH to Cpc-PH. However, even in a healthy pulmonary circulation, the mechanical forces in the smallest vessels (the arterioles, capillary bed, and venules) have not been quantitatively defined. This study is the first to examine this question via a computational fluid dynamics model of the human pulmonary arteries, arterioles, venules, and veins. Using this model, we predict temporal and spatial dynamics of cyclic stretch and wall shear stress with healthy and diseased hemodynamics. In the normotensive case for large vessels, numerical simulations show that large arteries have higher pressure and flow than large veins, as well as more pronounced changes in area throughout the cardiac cycle. In the microvasculature, shear stress increases and cyclic stretch decreases as vessel radius decreases. When we impose an increase in left atrial pressure to simulate Ipc-PH, shear stress decreases and cyclic stretch increases as compared to the healthy case. Overall, this model predicts pressure, flow, shear stress, and cyclic stretch that providing a way to analyze and investigate hypotheses related to disease progression in the pulmonary circulation.
Collapse
Affiliation(s)
- Michelle A Bartolo
- Department of Mathematics, North Carolina State University, Raleigh, NC, 27695, USA
| | - M Umar Qureshi
- Department of Mathematics, North Carolina State University, Raleigh, NC, 27695, USA
| | - Mitchel J Colebank
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, University of California, Irvine, CA, 92697, USA
| | - Naomi C Chesler
- Edwards Lifesciences Foundation Cardiovascular Innovation and Research Center, University of California, Irvine, CA, 92697, USA
| | - Mette S Olufsen
- Department of Mathematics, North Carolina State University, Raleigh, NC, 27607, USA.
| |
Collapse
|
21
|
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
|
22
|
Delaney M, Cleveland V, Mass P, Capuano F, Mandell JG, Loke YH, Olivieri L. Right ventricular afterload in repaired D-TGA is associated with inefficient flow patterns, rather than stenosis alone. Int J Cardiovasc Imaging 2021; 38:653-662. [PMID: 34727253 DOI: 10.1007/s10554-021-02436-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/03/2021] [Accepted: 10/04/2021] [Indexed: 11/24/2022]
Abstract
Treatment of D- transposition of great arteries (DTGA) involves the Arterial Switch Operation (ASO), which can create PA branch stenosis (PABS) and alter PA blood flow energetics. This altered PA flow may contribute to elevated right ventricular (RV) afterload more significantly than stenosis alone. Our aim was to correlate RV afterload and PA flow characteristics using 4D flow cardiac magnetic resonance (CMR) imaging of a mock circulatory system (MCS) incorporating 3D printed replicas. CMR imaging and clinical characteristics were analyzed from 22 ASO patients (age 11.9 ± 8.7 years, 68% male). Segmentation was performed to create 3D printed PA replicas that were mounted in an MRI-compatible MCS. Pressure drop across the PA replica was recorded and 4D flow CMR acquisitions were analyzed for blood flow inefficiency (energy loss, vorticity). In post-ASO patients, there is no difference in RV mass (p = 0.07), nor RV systolic pressure (p = 0.26) in the presence or absence of PABS. 4D flow analysis of MCS shows energy loss is correlated to RV mass (p = 0.01, r = 0.67) and MCS pressure differential (p = 0.02, r = 0.57). Receiver operating characteristic curve shows energy loss detects elevated RV mass above 30 g/m2 (p = 0.02, AUC 0.88) while index of PA dimensions (Nakata) does not (p = 0.09, AUC 0.79). PABS alone does not account for differences in RV mass or afterload in post-ASO patients. In MCS simulations, energy loss is correlated with both RV mass and PA pressure, and can moderately detect elevated RV mass. Inefficient PA flow may be an important predictor of RV afterload in this population.
Collapse
Affiliation(s)
- Marc Delaney
- Division of Pediatrics, Children's National Medical Center, 111 Michigan Ave, NW, Washington, DC, 20010, USA.
| | - Vincent Cleveland
- Sheikh Zayed Institute for Pediatric Surgical Innovation, Children's National Medical Center, Washington, DC, USA
| | - Paige Mass
- Sheikh Zayed Institute for Pediatric Surgical Innovation, Children's National Medical Center, Washington, DC, USA
| | - Francesco Capuano
- Department of Mechanics, Mathematics and Management, Politecnico di Bari, Bari, Italy
| | - Jason G Mandell
- Division of Cardiology, Children's National Medical Center, Washington, DC, USA
| | - Yue-Hin Loke
- Division of Cardiology, Children's National Medical Center, Washington, DC, USA
| | - Laura Olivieri
- Division of Cardiology, Children's National Medical Center, Washington, DC, USA
| |
Collapse
|
23
|
Lee EH, Baek S. Plasticity and Enzymatic Degradation Coupled With Volumetric Growth in Pulmonary Hypertension Progression. J Biomech Eng 2021; 143:111012. [PMID: 34076235 PMCID: PMC8299811 DOI: 10.1115/1.4051383] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Revised: 05/27/2021] [Indexed: 12/16/2022]
Abstract
Pulmonary hypertension (PH) is one of the least understood and highly elusive cardiovascular conditions associated with elevated pulmonary arterial pressure. Although the disease mechanisms are not completely understood, evidence has accumulated from human and animal studies that irreversible processes of pulmonary arterial wall damage, compensated by stress-mediated growth, play critical roles in eliciting the mechanisms of disease progression. The aim of this study is to develop a thermodynamic modeling structure of the pulmonary artery to consider coupled plastic-degradation-growth irreversible processes to investigate the mechanical roles of the dissipative phenomena in the disease progression. The proposed model performs a model parameter study of plastic deformation and degradation processes coupled with dissipative growth subjected to elevated pulmonary arterial pressure and computationally generates in silico simulations of PH progression using the clinical features of PH, found in human morphological and mechanical data. The results show that considering plastic deformation can provide a much better fitting of the ex vivo inflation tests than a widely used pure hyperelastic model in higher pressure conditions. In addition, the parameter sensitivity study illustrates that arterial damage and growth cause the increased stiffness, and the full simulation (combining elastic-plastic-degradation-growth models) reveals a key postpathological recovery process of compensating vessel damage by vascular adaptation by reducing the rate of vessel dilation and mediating vascular wall stress. Finally, the simulation results of luminal enlargement, arterial thickening, and arterial stiffness for an anisotropic growth are found to be close to the values from the literature.
Collapse
Affiliation(s)
- Eun-Ho Lee
- School of Mechanical Engineering, Sungkyunkwan University, Suwon, Gyeonggi-do 16419, South Korea; Department of Smart Fab. Technology, Sungkyunkwan University, Suwon, Gyeonggi-do 16419, South Korea
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, 2457 Engineering Building, East Lansing, MI 488424
| |
Collapse
|
24
|
Dong ML, Lan IS, Yang W, Rabinovitch M, Feinstein JA, Marsden AL. Computational simulation-derived hemodynamic and biomechanical properties of the pulmonary arterial tree early in the course of ventricular septal defects. Biomech Model Mechanobiol 2021; 20:2471-2489. [PMID: 34585299 DOI: 10.1007/s10237-021-01519-4] [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: 02/11/2021] [Accepted: 09/12/2021] [Indexed: 01/15/2023]
Abstract
Untreated ventricular septal defects (VSDs) can lead to pulmonary arterial hypertension (PAH) characterized by elevated pulmonary artery (PA) pressure and vascular remodeling, known as PAH associated with congenital heart disease (PAH-CHD). Though previous studies have investigated hemodynamic effects on vascular mechanobiology in late-stage PAH, hemodynamics leading to PAH-CHD initiation have not been fully quantified. We hypothesize that abnormal hemodynamics from left-to-right shunting in early stage VSDs affects PA biomechanical properties leading to PAH initiation. To model PA hemodynamics in healthy, small, moderate, and large VSD conditions prior to the onset of vascular remodeling, computational fluid dynamics simulations were performed using a 3D finite element model of a healthy 1-year-old's proximal PAs and a body-surface-area-scaled 0D distal PA tree. VSD conditions were modeled with increased pulmonary blood flow to represent degrees of left-to-right shunting. In the proximal PAs, pressure, flow, strain, and wall shear stress (WSS) increased with increasing VSD size; oscillatory shear index decreased with increasing VSD size in the larger PA vessels. WSS was higher in smaller diameter vessels and increased with VSD size, with the large VSD condition exhibiting WSS >100 dyn/cm[Formula: see text], well above values typically used to study dysfunctional mechanotransduction pathways in PAH. This study is the first to estimate hemodynamic and biomechanical metrics in the entire pediatric PA tree with VSD severity at the stage leading to PAH initiation and has implications for future studies assessing effects of abnormal mechanical stimuli on endothelial cells and vascular wall mechanics that occur during PAH-CHD initiation and progression.
Collapse
Affiliation(s)
- Melody L Dong
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| | - Ingrid S Lan
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| | - Weiguang Yang
- Department of Pediatrics, Stanford University, Stanford, CA, USA
| | | | - Jeffrey A Feinstein
- Department of Pediatrics and Bioengineering, Stanford University, Stanford, CA, USA
| | - Alison L Marsden
- Department of Pediatrics and Bioengineering, Stanford University, Stanford, CA, USA.
| |
Collapse
|
25
|
A distributed lumped parameter model of blood flow with fluid-structure interaction. Biomech Model Mechanobiol 2021; 20:1659-1674. [PMID: 34076757 DOI: 10.1007/s10237-021-01468-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Accepted: 05/17/2021] [Indexed: 12/22/2022]
Abstract
A distributed lumped parameter (DLP) model of blood flow was recently developed that can be simulated in minutes while still incorporating complex sources of energy dissipation in blood vessels. The aim of this work was to extend the previous DLP modeling framework to include fluid-structure interactions (DLP-FSI). This was done by using a simple compliance term to calculate pressure that does not increase the simulation complexity of the original DLP models. Verification and validation studies found DLP-FSI simulations had good agreement compared to analytical solutions of the wave equations, experimental measurements of pulsatile flow in elastic tubes, and in vivo MRI measurements of thoracic aortic flow. This new development of DLP-FSI allows for significantly improved computational efficiency of FSI simulations compared to FSI approaches that solve the full 3D conservation of mass and momentum equations while also including the complex sources of energy dissipation occurring in cardiovascular flows that other simplified models neglect.
Collapse
|
26
|
Sevilla-Montero J, Labrousse-Arias D, Fernández-Pérez C, Fernández-Blanco L, Barreira B, Mondéjar-Parreño G, Alfaro-Arnedo E, López IP, Pérez-Rial S, Peces-Barba G, Pichel JG, Peinado VI, Cogolludo Á, Calzada MJ. Cigarette Smoke Directly Promotes Pulmonary Arterial Remodeling and Kv7.4 Channel Dysfunction. Am J Respir Crit Care Med 2021; 203:1290-1305. [PMID: 33306938 DOI: 10.1164/rccm.201911-2238oc] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 12/10/2020] [Indexed: 01/10/2023] Open
Abstract
Rationale: Cigarette smoke is considered the chief leading cause of chronic obstructive pulmonary disease (COPD). Its impact on the progressive deterioration of airways has been extensively studied, but its direct effects on the pulmonary vasculature are less known. Objectives: To prove that pulmonary arterial remodeling in patients with COPD is not just a consequence of alveolar hypoxia but also due to the direct effects of cigarette smoke on the pulmonary vascular bed. Methods: We have used different molecular and cell biology approaches, as well as traction force microscopy, wire myography, and patch-clamp techniques in human cells and freshly isolated pulmonary arteries. In addition, we relied on in vivo models and human samples to analyze the effects of cigarette smoke on pulmonary vascular tone alterations. Measurements and Main Results: Cigarette smoke extract exposure directly promoted a hypertrophic, senescent phenotype that in turn contributed, through the secretion of inflammatory molecules, to an increase in the proliferative potential of nonexposed cells. Interestingly, these effects were significantly reversed by antioxidants. Furthermore, cigarette smoke extract affected cell contractility and dysregulated the expression and activity of the voltage-gated K+ channel Kv7.4. This contributed to the impairment of vasoconstriction and vasodilation responses. Most importantly, the levels of this channel were diminished in the lungs of smoke-exposed mice, smokers, and patients with COPD. Conclusions: Cigarette smoke directly contributes to pulmonary arterial remodeling through increased cell senescence, as well as vascular tone alterations because of diminished levels and function in the Kv7.4 channel. Strategies targeting these pathways may lead to novel therapies for COPD.
Collapse
Affiliation(s)
- Javier Sevilla-Montero
- Biomedical Research Institute La Princesa Hospital, Madrid, Spain
- Department of Medicine, School of Medicine, and
- Doctoral School, Autonoma University of Madrid, Madrid, Spain
| | - David Labrousse-Arias
- Biomedical Research Institute La Princesa Hospital, Madrid, Spain
- Department of Medicine, School of Medicine, and
| | - Cintia Fernández-Pérez
- Biomedical Research Institute La Princesa Hospital, Madrid, Spain
- Department of Medicine, School of Medicine, and
| | - Laura Fernández-Blanco
- Biomedical Research Institute La Princesa Hospital, Madrid, Spain
- Department of Medicine, School of Medicine, and
| | - Bianca Barreira
- Department of Pharmacology and Toxicology, School of Medicine, Complutense University of Madrid, Madrid, Spain
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
| | - Gema Mondéjar-Parreño
- Department of Pharmacology and Toxicology, School of Medicine, Complutense University of Madrid, Madrid, Spain
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
| | - Elvira Alfaro-Arnedo
- Lung Cancer and Respiratory Diseases Unit, Biomedical Research Center of La Rioja, Rioja Salud Foundation, Logroño, Spain
| | - Icíar P López
- Lung Cancer and Respiratory Diseases Unit, Biomedical Research Center of La Rioja, Rioja Salud Foundation, Logroño, Spain
| | - Sandra Pérez-Rial
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
- Respiratory Research Unit, Biomedical Research Unit, Health Research Institute Fundación Jiménez Díaz, Madrid, Spain; and
| | - Germán Peces-Barba
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
- Respiratory Research Unit, Biomedical Research Unit, Health Research Institute Fundación Jiménez Díaz, Madrid, Spain; and
| | - José G Pichel
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
- Lung Cancer and Respiratory Diseases Unit, Biomedical Research Center of La Rioja, Rioja Salud Foundation, Logroño, Spain
| | - Víctor Ivo Peinado
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
- Department of Pulmonary Medicine, Hospital August Pi i Sunyer Biomedical Research Institute, University of Barcelona, Barcelona, Spain
| | - Ángel Cogolludo
- Department of Pharmacology and Toxicology, School of Medicine, Complutense University of Madrid, Madrid, Spain
- Biomedical Research Networking Center in Respiratory Diseases, Institute of Health Carlos III, Madrid, Spain
| | - María J Calzada
- Biomedical Research Institute La Princesa Hospital, Madrid, Spain
- Department of Medicine, School of Medicine, and
| |
Collapse
|
27
|
Pewowaruk R, Lamers L, Roldán-Alzate A. Accelerated Estimation of Pulmonary Artery Stenosis Pressure Gradients with Distributed Lumped Parameter Modeling vs. 3D CFD with Instantaneous Adaptive Mesh Refinement: Experimental Validation in Swine. Ann Biomed Eng 2021; 49:2365-2376. [PMID: 33948748 DOI: 10.1007/s10439-021-02780-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Accepted: 04/11/2021] [Indexed: 11/30/2022]
Abstract
Branch pulmonary artery stenosis (PAS) commonly occurs in congenital heart disease and the pressure gradient over a stenotic PA lesion is an important marker for re-intervention. Image based computational fluid dynamics (CFD) has shown promise for non-invasively estimating pressure gradients but one limitation of CFD is long simulation times. The goal of this study was to compare accelerated predictions of PAS pressure gradients from 3D CFD with instantaneous adaptive mesh refinement (AMR) versus a recently developed 0D distributed lumped parameter CFD model. Predictions were then experimentally validated using a swine PAS model (n = 13). 3D CFD simulations with AMR improved efficiency by 5 times compared to fixed grid CFD simulations. 0D simulations further improved efficiency by 6 times compared to the 3D simulations with AMR. Both 0D and 3D simulations underestimated the pressure gradients measured by catheterization (- 1.87 ± 4.20 and - 1.78 ± 3.70 mmHg respectively). This was partially due to simulations neglecting the effects of a catheter in the stenosis. There was good agreement between 0D and 3D simulations (ICC 0.88 [0.66-0.96]) but only moderate agreement between simulations and experimental measurements (0D ICC 0.60 [0.11-0.86] and 3D ICC 0.66 [0.21-0.88]). Uncertainty assessment indicates that this was likely due to limited medical imaging resolution causing uncertainty in the segmented stenosis diameter in addition to uncertainty in the outlet resistances. This study showed that 0D lumped parameter models and 3D CFD with instantaneous AMR both improve the efficiency of hemodynamic modeling, but uncertainty from medical imaging resolution will limit the accuracy of pressure gradient estimations.
Collapse
Affiliation(s)
- Ryan Pewowaruk
- Biomedical Engineering, University of Wisconsin, Madison, WI, USA
| | - Luke Lamers
- Pediatrics, Division of Cardiology, University of Wisconsin, Madison, WI, USA
| | - Alejandro Roldán-Alzate
- Biomedical Engineering, University of Wisconsin, Madison, WI, USA. .,Mechanical Engineering, University of Wisconsin, Madison, WI, USA. .,Radiology, University of Wisconsin, Madison, WI, USA.
| |
Collapse
|
28
|
Zambrano BA, McLean N, Zhao X, Tan JL, Zhong L, Figueroa CA, Lee LC, Baek S. Patient-Specific Computational Analysis of Hemodynamics and Wall Mechanics and Their Interactions in Pulmonary Arterial Hypertension. Front Bioeng Biotechnol 2021; 8:611149. [PMID: 33634080 PMCID: PMC7901991 DOI: 10.3389/fbioe.2020.611149] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2020] [Accepted: 12/15/2020] [Indexed: 11/13/2022] Open
Abstract
Vascular wall stiffness and hemodynamic parameters are potential biomechanical markers for detecting pulmonary arterial hypertension (PAH). Previous computational analyses, however, have not considered the interaction between blood flow and wall deformation. Here, we applied an established computational framework that utilizes patient-specific measurements of hemodynamics and wall deformation to analyze the coupled fluid-vessel wall interaction in the proximal pulmonary arteries (PA) of six PAH patients and five control subjects. Specifically, we quantified the linearized stiffness (E), relative area change (RAC), diastolic diameter (D), regurgitant flow, and time-averaged wall shear stress (TAWSS) of the proximal PA, as well as the total arterial resistance (R t ) and compliance (C t ) at the distal pulmonary vasculature. Results found that the average proximal PA was stiffer [median: 297 kPa, interquartile range (IQR): 202 kPa vs. median: 75 kPa, IQR: 5 kPa; P = 0.007] with a larger diameter (median: 32 mm, IQR: 5.25 mm vs. median: 25 mm, IQR: 2 mm; P = 0.015) and a reduced RAC (median: 0.22, IQR: 0.10 vs. median: 0.42, IQR: 0.04; P = 0.004) in PAH compared to our control group. Also, higher total resistance (R t ; median: 6.89 mmHg × min/l, IQR: 2.16 mmHg × min/l vs. median: 3.99 mmHg × min/l, IQR: 1.15 mmHg × min/l; P = 0.002) and lower total compliance (C t ; median: 0.13 ml/mmHg, IQR: 0.15 ml/mmHg vs. median: 0.85 ml/mmHg, IQR: 0.51 ml/mmHg; P = 0.041) were observed in the PAH group. Furthermore, lower TAWSS values were seen at the main PA arteries (MPAs) of PAH patients (median: 0.81 Pa, IQR: 0.47 Pa vs. median: 1.56 Pa, IQR: 0.89 Pa; P = 0.026) compared to controls. Correlation analysis within the PAH group found that E was directly correlated to the PA regurgitant flow (r = 0.84, P = 0.018) and inversely related to TAWSS (r = -0.72, P = 0.051). Results suggest that the estimated elastic modulus E may be closely related to PAH hemodynamic changes in pulmonary arteries.
Collapse
Affiliation(s)
- Byron A. Zambrano
- J. Mike Walker '66 Department of Mechanical Engineering, Texas A&M University, College Station, TX, United States
| | - Nathan McLean
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Xiaodan Zhao
- National Heart Centre Singapore, Singapore, Singapore
| | - Ju-Le Tan
- National Heart Centre Singapore, Singapore, Singapore
| | - Liang Zhong
- National Heart Centre Singapore, Singapore, Singapore
- Duke-National University of Singapore, Singapore, Singapore
| | - C. Alberto Figueroa
- Departments of Biomedical Engineering and Surgery, University of Michigan, Ann Arbor, MI, United States
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| |
Collapse
|
29
|
Pillalamarri N, Patnaik S, Piskin S, Gueldner P, Finol E. Ex Vivo Regional Mechanical Characterization of Porcine Pulmonary Arteries. EXPERIMENTAL MECHANICS 2021; 61:285-303. [PMID: 33814554 PMCID: PMC8011683 DOI: 10.1007/s11340-020-00678-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2020] [Accepted: 10/22/2020] [Indexed: 06/12/2023]
Abstract
BACKGROUND Regional mechanical characterization of pulmonary arteries can be useful in the development of computational models of pulmonary arterial mechanics. OBJECTIVE We performed a biomechanical and microstructural characterization study of porcine pulmonary arteries, inclusive of the main, left, and right pulmonary arteries (MPA, LPA, and RPA, respectively). METHODS The specimens were initially stored at -20°C and allowed to thaw for 12-24 hours prior to testing. Each artery was further subdivided into proximal, middle, and distal regions, leading to ten location-based experimental groups. Planar equibiaxial tensile testing was performed to evaluate the mechanical behavior of the specimens, from which we calculated the stress at the maximum strain (S 55), tensile modulus (TM), anisotropy index (AI), and strain energy in terms of area under the stress-strain curve (AUC). Histological quantification was performed to evaluate the area fraction of elastin and collagen content, intima-media thickness (IMT), and adventitial thickness (AT). The constitutive material behavior of each group was represented by a five-constant Holzapfel-Gasser-Ogden model. RESULTS The specimens exhibited non-linear stress-strain characteristics across all groups. The MPA exhibited the highest mean wall stress and TM in the longitudinal and circumferential directions, while the bifurcation region yielded the highest values of AI and AUC. All regions revealed a higher stiffness in the longitudinal direction compared to the circumferential direction, suggesting a degree of anisotropy that is believed to be within the margin of experimental uncertainty. Collagen content was found to be the highest in the MPA and decreased significantly at the bifurcation, LPA and RPA. Elastin content did not yield such significant differences amongst the ten groups. The MPA had the highest IMT, which decreased concomitantly to the distal LPA and RPA. No significant differences were found in the AT amongst the ten groups. CONCLUSION The mechanical properties of porcine pulmonary arteries exhibit strong regional dissimilarities, which can be used to inform future studies of high fidelity finite element models.
Collapse
Affiliation(s)
- N.R. Pillalamarri
- University of Texas at San Antonio, Department of Mechanical Engineering, San Antonio, TX
| | - S.S. Patnaik
- University of Texas at San Antonio, Department of Mechanical Engineering, San Antonio, TX
| | - S. Piskin
- University of Texas at San Antonio, Department of Mechanical Engineering, San Antonio, TX
- Istinye University, Department of Mechanical Engineering, Zeytinburnu, Istanbul, Turkey
| | - P. Gueldner
- University of Texas at San Antonio, Department of Biomedical Engineering, San Antonio, TX
| | - E.A. Finol
- University of Texas at San Antonio, Department of Mechanical Engineering, San Antonio, TX
- University of Texas at San Antonio, UTSA/UTHSA Joint Graduate Program in Biomedical Engineering, San Antonio, TX
| |
Collapse
|
30
|
Bracamonte JH, Wilson JS, Soares JS. Assessing Patient-Specific Mechanical Properties of Aortic Wall and Peri-Aortic Structures From In Vivo DENSE Magnetic Resonance Imaging Using an Inverse Finite Element Method and Elastic Foundation Boundary Conditions. J Biomech Eng 2020; 142:121011. [PMID: 32632452 DOI: 10.1115/1.4047721] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2020] [Indexed: 11/08/2022]
Abstract
The establishment of in vivo, noninvasive patient-specific, and regionally resolved techniques to quantify aortic properties is key to improving clinical risk assessment and scientific understanding of vascular growth and remodeling. A promising and novel technique to reach this goal is an inverse finite element method (FEM) approach that utilizes magnetic resonance imaging (MRI)-derived displacement fields from displacement encoding with stimulated echoes (DENSE). Previous studies using DENSE MRI suggested that the infrarenal abdominal aorta (IAA) deforms heterogeneously during the cardiac cycle. We hypothesize that this heterogeneity is driven in healthy aortas by regional adventitial tethering and interaction with perivascular tissues, which can be modeled with elastic foundation boundary conditions (EFBCs) using a collection of radially oriented springs with varying stiffness with circumferential distribution. Nine healthy IAAs were modeled using previously acquired patient-specific imaging and displacement fields from steady-state free procession (SSFP) and DENSE MRI, followed by assessment of aortic wall properties and heterogeneous EFBC parameters using inverse FEM. In contrast to traction-free boundary condition, prescription of EFBC reduced the nodal displacement error by 60% and reproduced the DENSE-derived heterogeneous strain distribution. Estimated aortic wall properties were in reasonable agreement with previously reported experimental biaxial testing data. The distribution of normalized EFBC stiffness was consistent among all patients and spatially correlated to standard peri-aortic anatomical features, suggesting that EFBC could be generalized for human adults with normal anatomy. This approach is computationally inexpensive, making it ideal for clinical research and future incorporation into cardiovascular fluid-structure analyses.
Collapse
Affiliation(s)
- Johane H Bracamonte
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, Richmond, VA 23284
| | - John S Wilson
- Department of Biomedical Engineering and Pauley Heart Center, Virginia Commonwealth University, 601 West Main Street, Richmond, VA 23284
| | - Joao S Soares
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, 401 West Main Street, Richmond, VA 23284
| |
Collapse
|
31
|
Gopalan D, Gibbs JSR. From Early Morphometrics to Machine Learning-What Future for Cardiovascular Imaging of the Pulmonary Circulation? Diagnostics (Basel) 2020; 10:diagnostics10121004. [PMID: 33255668 PMCID: PMC7760106 DOI: 10.3390/diagnostics10121004] [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: 10/24/2020] [Revised: 11/19/2020] [Accepted: 11/24/2020] [Indexed: 02/07/2023] Open
Abstract
Imaging plays a cardinal role in the diagnosis and management of diseases of the pulmonary circulation. Behind the picture itself, every digital image contains a wealth of quantitative data, which are hardly analysed in current routine clinical practice and this is now being transformed by radiomics. Mathematical analyses of these data using novel techniques, such as vascular morphometry (including vascular tortuosity and vascular volumes), blood flow imaging (including quantitative lung perfusion and computational flow dynamics), and artificial intelligence, are opening a window on the complex pathophysiology and structure-function relationships of pulmonary vascular diseases. They have the potential to make dramatic alterations to how clinicians investigate the pulmonary circulation, with the consequences of more rapid diagnosis and a reduction in the need for invasive procedures in the future. Applied to multimodality imaging, they can provide new information to improve disease characterization and increase diagnostic accuracy. These new technologies may be used as sophisticated biomarkers for risk prediction modelling of prognosis and for optimising the long-term management of pulmonary circulatory diseases. These innovative techniques will require evaluation in clinical trials and may in themselves serve as successful surrogate end points in trials in the years to come.
Collapse
Affiliation(s)
- Deepa Gopalan
- Imperial College Healthcare NHS Trust, London W12 0HS, UK
- Imperial College London, London SW7 2AZ, UK;
- Cambridge University Hospital, Cambridge CB2 0QQ, UK
- Correspondence: ; Tel.: +44-77-3000-7780
| | - J. Simon R. Gibbs
- Imperial College London, London SW7 2AZ, UK;
- National Heart & Lung Institute, Imperial College London, London SW3 6LY, UK
| |
Collapse
|
32
|
Leopold JA. Endocrine Delivery of MicroRNA-210: A Trusted Traveler That Mediates Pulmonary Hypertension. Circ Res 2020; 127:693-695. [PMID: 32790524 DOI: 10.1161/circresaha.120.317625] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Affiliation(s)
- Jane A Leopold
- From the Division of Cardiovascular Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA
| |
Collapse
|
33
|
Fanni BM, Sauvage E, Celi S, Norman W, Vignali E, Landini L, Schievano S, Positano V, Capelli C. A Proof of Concept of a Non-Invasive Image-Based Material Characterization Method for Enhanced Patient-Specific Computational Modeling. Cardiovasc Eng Technol 2020; 11:532-543. [PMID: 32748364 DOI: 10.1007/s13239-020-00479-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Accepted: 07/22/2020] [Indexed: 11/30/2022]
Abstract
PURPOSE Computational models of cardiovascular structures rely on their accurate mechanical characterization. A validated method able to infer the material properties of patient-specific large vessels is currently lacking. The aim of the present study is to present a technique starting from the flow-area (QA) method to retrieve basic material properties from magnetic resonance (MR) imaging. METHODS The proposed method was developed and tested, first, in silico and then in vitro. In silico, fluid-structure interaction (FSI) simulations of flow within a deformable pipe were run with varying elastic modules (E) between 0.5 and 32 MPa. The proposed QA-based formulation was assessed and modified based on the FSI results to retrieve E values. In vitro, a compliant phantom connected to a mock circulatory system was tested within MR scanning. Images of the phantom were acquired and post-processed according to the modified formulation to infer E of the phantom. Results of in vitro imaging assessment were verified against standard tensile test. RESULTS In silico results from FSI simulations were used to derive the correction factor to the original formulation based on the geometrical and material characteristics. In vitro, the modified QA-based equation estimated an average E = 0.51 MPa, 2% different from the E derived from tensile tests (i.e. E = 0.50 MPa). CONCLUSION This study presented promising results of an indirect and non-invasive method to establish elastic properties from solely MR images data, suggesting a potential image-based mechanical characterization of large blood vessels.
Collapse
Affiliation(s)
- B M Fanni
- BioCardioLab, Bioengineering Unit, Fondazione Toscana Gabriele Monasterio, Via Aurelia Sud, 54100, Massa, Italy.,Department of Information Engineering, University of Pisa, Via Girolamo Caruso 16, 56122, Pisa, Italy
| | - E Sauvage
- UCL Institute of Cardiovascular Science, 20c Guilford Street, London, WC1N 1DZ, UK.,Great Ormond Street Hospital for Children, NHS Foundation Trust, 30 Great Ormond Street, London, WC1N 3JH, UK
| | - S Celi
- BioCardioLab, Bioengineering Unit, Fondazione Toscana Gabriele Monasterio, Via Aurelia Sud, 54100, Massa, Italy.
| | - W Norman
- UCL Institute of Cardiovascular Science, 20c Guilford Street, London, WC1N 1DZ, UK.,Great Ormond Street Hospital for Children, NHS Foundation Trust, 30 Great Ormond Street, London, WC1N 3JH, UK
| | - E Vignali
- BioCardioLab, Bioengineering Unit, Fondazione Toscana Gabriele Monasterio, Via Aurelia Sud, 54100, Massa, Italy.,Department of Information Engineering, University of Pisa, Via Girolamo Caruso 16, 56122, Pisa, Italy
| | - L Landini
- BioCardioLab, Bioengineering Unit, Fondazione Toscana Gabriele Monasterio, Via Aurelia Sud, 54100, Massa, Italy.,Department of Information Engineering, University of Pisa, Via Girolamo Caruso 16, 56122, Pisa, Italy
| | - S Schievano
- UCL Institute of Cardiovascular Science, 20c Guilford Street, London, WC1N 1DZ, UK.,Great Ormond Street Hospital for Children, NHS Foundation Trust, 30 Great Ormond Street, London, WC1N 3JH, UK
| | - V Positano
- BioCardioLab, Bioengineering Unit, Fondazione Toscana Gabriele Monasterio, Via Aurelia Sud, 54100, Massa, Italy
| | - C Capelli
- UCL Institute of Cardiovascular Science, 20c Guilford Street, London, WC1N 1DZ, UK.,Great Ormond Street Hospital for Children, NHS Foundation Trust, 30 Great Ormond Street, London, WC1N 3JH, UK
| |
Collapse
|
34
|
Shavik SM, Tossas-Betancourt C, Figueroa CA, Baek S, Lee LC. Multiscale Modeling Framework of Ventricular-Arterial Bi-directional Interactions in the Cardiopulmonary Circulation. Front Physiol 2020; 11:2. [PMID: 32116737 PMCID: PMC7025512 DOI: 10.3389/fphys.2020.00002] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 01/03/2020] [Indexed: 02/02/2023] Open
Abstract
Ventricular-arterial coupling plays a key role in the physiologic function of the cardiovascular system. We have previously described a hybrid lumped-finite element (FE) modeling framework of the systemic circulation that couples idealized FE models of the aorta and the left ventricle (LV). Here, we describe an extension of the lumped-FE modeling framework that couples patient-specific FE models of the left and right ventricles, aorta and the large pulmonary arteries in both the systemic and pulmonary circulations. Geometries of the FE models were reconstructed from magnetic resonance (MR) images acquired in a pediatric patient diagnosed with pulmonary arterial hypertension (PAH). The modeling framework was calibrated with pressure waveforms acquired in the heart and arteries by catheterization as well as ventricular volume and arterial diameter waveforms measured from MR images. The calibrated model hemodynamic results match well with the clinically-measured waveforms (volume and pressure) in the LV and right ventricle (RV) as well as with the clinically-measured waveforms (pressure and diameter) in the aorta and main pulmonary artery. The calibrated framework was then used to simulate three cases, namely, (1) an increase in collagen in the large pulmonary arteries, (2) a decrease in RV contractility, and (3) an increase in the total pulmonary arterial resistance, all characteristics of progressive PAH. The key finding from these simulations is that hemodynamics of the pulmonary vasculature and RV wall stress are more sensitive to vasoconstriction with a 10% of reduction in the lumen diameter of the distal vessels than a 67% increase in the proximal vessel's collagen mass.
Collapse
Affiliation(s)
- Sheikh Mohammad Shavik
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States.,Department of Mechanical Engineering, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh
| | | | - C Alberto Figueroa
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States.,Department of Surgery, University of Michigan, Ann Arbor, MI, United States
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| |
Collapse
|
35
|
Abstract
Knowledge of physiologic hemodynamics is a fundamental requirement to establish pathological findings. However, little is known about the normal flow fields in the pulmonary arteries, especially for children. The purpose of this study is to characterize flow patterns in the pulmonary artery bifurcation of healthy pediatric subjects using direct numerical simulations. A realistic geometry is obtained via statistical shape modeling, by averaging five subject-specific digital models extracted from cardiovascular magnetic resonance datasets of healthy volunteers. Boundary conditions are assigned to mimic physiological conditions at rest, corresponding to a peak Reynolds number equal to 3400 and a Womersley number equal to 15. Results show that the normal bifurcation is highly hemodynamically efficient, as measured by an energy dissipation index. The curvature of the pulmonary arteries is sufficiently small to prevent flow separation along the inner walls, and no signs of a turbulent-like state are found. In line with previous imaging studies, a helical structure protruding into the right pulmonary artery is detected, and its formation mechanism is elucidated in the paper. These findings might help to identify abnormal flow features in patients with altered anatomic and physiologic states, particularly those with repaired congenital heart disease.
Collapse
|
36
|
Multiscale modeling of ventricular–vascular dysfunction in pulmonary arterial hypertension. CURRENT OPINION IN BIOMEDICAL ENGINEERING 2019. [DOI: 10.1016/j.cobme.2019.09.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
37
|
Wu D, Birukov K. Endothelial Cell Mechano-Metabolomic Coupling to Disease States in the Lung Microvasculature. Front Bioeng Biotechnol 2019; 7:172. [PMID: 31380363 PMCID: PMC6658821 DOI: 10.3389/fbioe.2019.00172] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Accepted: 07/03/2019] [Indexed: 12/15/2022] Open
Abstract
Lungs are the most vascular part of humans, accepting the totality of cardiac output in a volume much smaller than the body itself. Due to this cardiac output, the lung microvasculature is subject to mechanical forces including shear stress and cyclic stretch that vary with the cardiac and breathing cycle. Vessels are surrounded by extracellular matrix which dictates the stiffness which endothelial cells also sense and respond to. Shear stress, stiffness, and cyclic stretch are known to influence endothelial cell state. At high shear stress, endothelial cells exhibit cell quiescence marked by low inflammatory markers and high nitric oxide synthesis, whereas at low shear stress, endothelial cells are thought to "activate" into a pro-inflammatory state and have low nitric oxide. Shear stress' profound effect on vascular phenotype is most apparent in the arterial vasculature and in the pathophysiology of vascular inflammation. To conduct the flow of blood from the right heart, the lung microvasculature must be rigid yet compliant. It turns out that excessive substrate rigidity or stiffness is important in the development of pulmonary hypertension and chronic fibrosing lung diseases via excessive cell proliferation or the endothelial-mesenchymal transition. Recently, a new body of literature has evolved that couples mechanical sensing to endothelial phenotypic changes through metabolic signaling in clinically relevant contexts such as pulmonary hypertension, lung injury syndromes, as well as fibrosis, which is the focus of this review. Stretch, like flow, has profound effect on endothelial phenotype; metabolism studies due to stretch are in their infancy.
Collapse
Affiliation(s)
- David Wu
- Section of Pulmonary and Critical Care, Department of Medicine, University of Chicago, Chicago, IL, United States
| | - Konstantin Birukov
- Department of Anesthesia, University of Maryland, Baltimore, MD, United States
| |
Collapse
|
38
|
Avazmohammadi R, Mendiola E, Li D, Vanderslice P, Dixon R, Sacks M. Interactions between structural remodeling and volumetric growth in right ventricle in response to pulmonary arterial hypertension. J Biomech Eng 2019; 141:2737741. [PMID: 31260516 DOI: 10.1115/1.4044174] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2019] [Indexed: 01/22/2023]
Abstract
Pulmonary arterial hypertension (PAH) exerts substantial pressure overload on the right ventricle (RV). The associated RV free wall (RVFW) adaptation could consist of myocardial hypertrophy, augmented intrinsic contractility, collagen fibrosis, and structural remodeling in an attempt to cope with pressure overload. If RVFW adaptation cannot maintain the RV stroke volume, RV dilation will prevail as an exit mechanism which usually decompensates the RV function leading to RV failure. Our knowledge of the factors determining the transition from the upper limit of RVFW adaptation to RV decompensation and the role of fiber remodeling events in this transition remains very limited. Computational heart models that connect the growth and remodeling (G\&R) events at the fiber and tissue levels with alterations in the organ-level function are essential to predict the temporal order and the compensatory level of the underlying mechanisms. In this work, building upon our recent rodent heart models (RHM) of PAH, we integrated mathematical models that describe time-evolution volumetric growth of the RV and structural remodeling of the RVFW. Results suggest that augmentation of the intrinsic contractility of myofibers accompanied by an increase in passive stiffness of RVFW is among the first remodeling events through which the RV strives to maintain the cardiac output. Interestingly, we found that the observed reorientation of the myofibers towards the longitudinal (apex-to-base) direction was a maladaptive mechanism that impaired the contractile pattern of RVFW and advanced along with RV dilation at later stages of PAH development.
Collapse
Affiliation(s)
- Reza Avazmohammadi
- James T. Willerson Center for Cardiovascular Modeling and Simulation Oden Institute for Computational Engineering and Sciences and the Department of Biomedical Engineering
| | - Emilio Mendiola
- James T. Willerson Center for Cardiovascular Modeling and Simulation Oden Institute for Computational Engineering and Sciences and the Department of Biomedical Engineering
| | - David Li
- James T. Willerson Center for Cardiovascular Modeling and Simulation Oden Institute for Computational Engineering and Sciences and the Department of Biomedical Engineering
| | - Peter Vanderslice
- Department of Molecular Cardiology, Texas Heart Institute, Houston, TX, USA; The University of Texas at Austin, Austin, TX, USA
| | - Richard Dixon
- Department of Molecular Cardiology, Texas Heart Institute, Houston, TX, USA; The University of Texas at Austin, Austin, TX, USA
| | - Michael Sacks
- James T. Willerson Center for Cardiovascular Modeling and Simulation Oden Institute for Computational Engineering and Sciences and the Department of Biomedical Engineering
| |
Collapse
|
39
|
Pewowaruk R, Roldán-Alzate A. 4D Flow MRI Estimation of Boundary Conditions for Patient Specific Cardiovascular Simulation. Ann Biomed Eng 2019; 47:1786-1798. [PMID: 31069584 DOI: 10.1007/s10439-019-02285-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2019] [Accepted: 05/02/2019] [Indexed: 12/11/2022]
Abstract
Accurate image based cardiovascular simulations require patient specific boundary conditions (BCs) for inlets, outlets and vessel wall mechanical properties. While inlet BCs are typically determined non-invasively, invasive pressure catheterization is often used to determine patient specific outlet BCs and vessel wall mechanical properties. A method using 4D Flow MRI to non-invasively determine both patient specific outlet BCs and vessel wall mechanical properties is presented and results for both in vitro validation with a latex tube and an in vivo pulmonary artery stenosis (PAS) stent intervention are presented. For in vitro validation, acceptable agreement is found between simulation using BCs from 4D Flow MRI and benchtop measurements. For the PAS virtual intervention, simulation correctly predicts flow distribution with 9% error compared to MRI. Using 4D Flow MRI to noninvasively determine patient specific BCs increases the ability to use image based simulations as pressure catheterization is not always performed.
Collapse
Affiliation(s)
- Ryan Pewowaruk
- Biomedical Engineering, University of Wisconsin - Madison, 1111 Highland Ave, Room 2476 WIMR 2, Madison, WI, 53705, USA
| | - Alejandro Roldán-Alzate
- Biomedical Engineering, University of Wisconsin - Madison, 1111 Highland Ave, Room 2476 WIMR 2, Madison, WI, 53705, USA. .,Mechanical Engineering, University of Wisconsin - Madison, 1111 Highland Ave, Room 2476 WIMR 2, Madison, WI, 53705, USA. .,Department of Radiology, University of Wisconsin - Madison, 1111 Highland Ave, Room 2476 WIMR 2, Madison, WI, 53705, USA.
| |
Collapse
|
40
|
Capuano F, Loke YH, Cronin I, Olivieri LJ, Balaras E. Computational Study of Pulmonary Flow Patterns After Repair of Transposition of Great Arteries. J Biomech Eng 2019; 141:2727821. [DOI: 10.1115/1.4043034] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Indexed: 11/08/2022]
Abstract
Patients that undergo the arterial switch operation (ASO) to repair transposition of great arteries (TGA) can develop abnormal pulmonary trunk morphology with significant long-term complications. In this study, cardiovascular magnetic resonance was combined with computational fluid dynamics to investigate the impact of the postoperative layout on the pulmonary flow patterns. Three ASO patients were analyzed and compared to a volunteer control. Results showed the presence of anomalous shear layer instabilities, vortical and helical structures, and turbulent-like states in all patients, particularly as a consequence of the unnatural curvature of the pulmonary bifurcation. Streamlined, mostly laminar flow was instead found in the healthy subject. These findings shed light on the correlation between the post-ASO anatomy and the presence of altered flow features, and may be useful to improve surgical planning as well as the long-term care of TGA patients.
Collapse
Affiliation(s)
- Francesco Capuano
- Department of Industrial Engineering, Università di Napoli Federico II, Napoli 80125, Italy e-mail:
| | - Yue-Hin Loke
- Division of Cardiology, Children's National Health System, Washington, DC 20010 e-mail:
| | - Ileen Cronin
- Division of Cardiology, Children's National Health System, Washington, DC 20010 e-mail:
| | - Laura J. Olivieri
- Division of Cardiology, The Sheikh Zayed Institute for Pediatric Surgical Innovation, Children's National Health System, Washington, DC 20010 e-mail:
| | - Elias Balaras
- Department of Mechanical and Aerospace Engineering, George Washington University, Washington, DC 20052 e-mail:
| |
Collapse
|