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
|
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
|
3
|
Kharvani HR, Aghanajafi C. Investigation into the two-way interaction of coronary flow and heart function in coronary artery disease predicted by a computational model of autoregulation of coronary flow. J Biomech 2024; 164:111970. [PMID: 38325193 DOI: 10.1016/j.jbiomech.2024.111970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Revised: 01/22/2024] [Accepted: 01/28/2024] [Indexed: 02/09/2024]
Abstract
This study presents a closed-loop computational model to investigate the interplay between heart function, coronary flow, and systemic circulation during exercise, with a specific focus on the impact of coronary artery stenosis. The model incorporates a lumped representation of the heart, main arteries, and coronary arteries, establishing a closed circulatory system. The simulation investigates the autoregulation of coronary flow in response to myocardial oxygen demands during physical exercise by incorporating sympathetic and parasympathetic functions. This study establishes a closed supply-demand loop and investigates the effect of coronary flow deficiency on heart function and systemic circulation in coronary artery diseases during exercise. In coronary artery diseases with low stenosis, heart function and systemic flow resemble those of a healthy person. However, as stenosis intensifies with physical exercise, an additional regulatory mechanism (reg2) is activated. This mechanism adjusts coronary flow by reducing myocardial contractility (E) and increasing heart rate (HR) while maintaining cardiac output (CO). The study results indicate that, at the highest exercise intensity for a healthy individual (HR = 150), the value of E increases from 6 to 8.65mmHg/ml. Meanwhile, for a patient with 85 % coronary artery stenosis in the same exercise intensity, the HR increases to 200, and the value of E decreases to 3.45mmHg/ml. The results also demonstrate that the initiation of the (reg2) mechanism at rest occurs at 83 % stenosis, while at the highest exercise intensity, this mechanism commences at 67 % stenosis.
Collapse
Affiliation(s)
| | - Cyrus Aghanajafi
- Faculty of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, Iran.
| |
Collapse
|
4
|
Pearce NF, Kim EJ. A new synergistic model for simulating exercise incorporating control mechanisms at cellular and organ scales. Comput Biol Med 2023; 163:107141. [PMID: 37327758 DOI: 10.1016/j.compbiomed.2023.107141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Revised: 05/24/2023] [Accepted: 06/04/2023] [Indexed: 06/18/2023]
Abstract
The physiological response of the cardio-vascular system (CVS) to physical activity is of great importance to those working in sporting research and has profound consequences for the health and well-being of people. Coronary vasodilation and the physiological mechanisms involved in exercise have frequently been the focus of numerical models for simulating exercise. This is partly achieved using the time-varying-elastance (TVE) theory, which prescribes the pressure-volume relationship of the ventricle as a periodic function of time, tuned using empirical data. The empirical foundations of the TVE method however, and its suitability for CVS modelling are frequently questioned. To overcome this challenge, we adopt a different synergistic approach in which a model for the microscale heart muscle (myofibers) activity is embedded within a macro organ-scale CVS model. We developed such a synergistic model by including the coronary flow and various control mechanisms at the circulation level through feedback and feedforward means, and at the microscale (contractile) through the regulation of ATP availability and myofiber force depending on exercise intensity or heart rate. The coronary flow produced by the model displays the well-known 2-phase character of the flow, which is preserved under exercise. The model is tested by simulating reactive hyperemia, which is a transient occlusion of the coronary flow, successfully reproducing the additional coronary flow following the block removal. On-transient exercise results reveal a rise in both cardiac output and mean ventricle pressure as expected. The stroke volume increases initially, but then declines during the latter period of HR rise, corresponding with one of the main physiological responses to exercise. The pressure-volume loop expands during exercise, as systolic pressure rises. The Myocardial oxygen demand increases during exercise and the coronary blood supply increases in response, causing an excess of oxygen supply to the heart. Off-transient exercise recovery is largely a reverse of this response, although the behaviour is slightly more varied, with sudden spikes in coronary resistance. Different levels of fitness and exercise intensity are tested and reveal that the stroke volume rises until a level of myocardial oxygen demand is reached at which point it declines. This level of demand is independent of fitness or exercise intensity. An advantage of our model is demonstrated in the correspondence between the micro and organ scale mechanics so that cellular pathologies can be traced from exercise performance with relatively little computational or experimental expense.
Collapse
Affiliation(s)
- Nicholas F Pearce
- Fluids and Complex Systems Center, Faculty of Engineering, Environment and Computing, Coventry University, Coventry, CV1 5FB, UK.
| | - Eun-Jin Kim
- Fluids and Complex Systems Center, Faculty of Engineering, Environment and Computing, Coventry University, Coventry, CV1 5FB, UK
| |
Collapse
|
5
|
Munneke AG, Lumens J, Arts T, Prinzen FW, Delhaas T. Myocardial perfusion and flow reserve in the asynchronous heart: mechanistic insight from a computational model. J Appl Physiol (1985) 2023; 135:489-499. [PMID: 37439238 PMCID: PMC10538979 DOI: 10.1152/japplphysiol.00181.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 06/08/2023] [Accepted: 06/27/2023] [Indexed: 07/14/2023] Open
Abstract
The tight coupling between myocardial oxygen demand and supply has been recognized for decades, but it remains controversial whether this coupling persists under asynchronous activation, such as during left bundle branch block (LBBB). Furthermore, it is unclear whether the amount of local cardiac wall growth, following longer-lasting asynchronous activation, can explain differences in myocardial perfusion distribution between subjects. For a better understanding of these matters, we built upon our existing modeling framework for cardiac mechanics-to-perfusion coupling by incorporating coronary autoregulation. Regional coronary flow was regulated with a vasodilator signal based on regional demand, as estimated from regional fiber stress-strain area. Volume of left ventricular wall segments was adapted with chronic asynchronous activation toward a homogeneous distribution of myocardial oxygen demand per tissue weight. Modeling results show that 1) both myocardial oxygen demand and supply are decreased in early activated regions and increased in late-activated regions; 2) but that regional hyperemic flow remains unaffected; while 3) regional myocardial flow reserve (the ratio of hyperemic to resting myocardial flow) decreases with increases in absolute regional myocardial oxygen demand as well as with decreases in wall thickness. These findings suggest that septal hypoperfusion in LBBB represents an autoregulatory response to reduced myocardial oxygen demand. Furthermore, oxygen demand-driven remodeling of wall mass can explain asymmetric hypertrophy and the related homogenization of myocardial perfusion and flow reserve. Finally, the inconsistent observations of myocardial perfusion distribution can primarily be explained by the degree of dyssynchrony, the degree of asymmetric hypertrophy, and the imaging modality used.NEW & NOTEWORTHY This versatile modeling framework couples myocardial oxygen demand to oxygen supply and myocardial growth, enabling simulation of resting and hyperemic myocardial flow during acute and chronic asynchronous ventricular activation. Model-based findings suggest that reported inconsistencies in myocardial perfusion and flow reserve responses with asynchronous ventricular activation between patients can primarily be explained by the degree of dyssynchrony and wall mass remodeling, which together determine the heterogeneity in regional oxygen demand and, hence, supply with autoregulation.
Collapse
Affiliation(s)
- Anneloes G Munneke
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Theo Arts
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Frits W Prinzen
- Department of Physiology, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Tammo Delhaas
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
6
|
Kim T, van Bakel PAJ, Nama N, Burris N, Patel HJ, Williams DM, Figueroa CA. A Computational Study of Dynamic Obstruction in Type B Aortic Dissection. J Biomech Eng 2023; 145:031008. [PMID: 36459144 PMCID: PMC10854260 DOI: 10.1115/1.4056355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Revised: 11/18/2022] [Accepted: 11/21/2022] [Indexed: 12/05/2022]
Abstract
A serious complication in aortic dissection is dynamic obstruction of the true lumen (TL). Dynamic obstruction results in malperfusion, a blockage of blood flow to a vital organ. Clinical data reveal that increases in central blood pressure promote dynamic obstruction. However, the mechanisms by which high pressures result in TL collapse are underexplored and poorly understood. Here, we developed a computational model to investigate biomechanical and hemodynamical factors involved in Dynamic obstruction. We hypothesize that relatively small pressure gradient between TL and false lumen (FL) are sufficient to displace the flap and induce obstruction. An idealized fluid-structure interaction model of type B aortic dissection was created. Simulations were performed under mean cardiac output while inducing dynamic changes in blood pressure by altering FL outflow resistance. As FL resistance increased, central aortic pressure increased from 95.7 to 115.3 mmHg. Concurrent with blood pressure increase, flap motion was observed, resulting in TL collapse, consistent with clinical findings. The maximum pressure gradient between TL and FL over the course of the dynamic obstruction was 4.5 mmHg, consistent with our hypothesis. Furthermore, the final stage of dynamic obstruction was very sudden in nature, occurring over a short time (<1 s) in our simulation, consistent with the clinical understanding of this dramatic event. Simulations also revealed sudden drops in flow and pressure in the TL in response to the flap motion, consistent with first stages of malperfusion. To our knowledge, this study represents the first computational analysis of potential mechanisms driving dynamic obstruction in aortic dissection.
Collapse
Affiliation(s)
- T Kim
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI 48105
| | - P A J van Bakel
- Department of Cardiac Surgery, University of Michigan, Ann Arbor, MI 48105
| | - N Nama
- Department of Mechanical and Materials Engineering, University of Nebraska-Lincoln, Lincoln, NE 68588
| | - N Burris
- Department of Radiology, University of Michigan, Ann Arbor, MI 48105
| | - H J Patel
- Department of Cardiac Surgery, University of Michigan, Ann Arbor, MI 48105
| | - D M Williams
- Department of Radiology, University of Michigan, Ann Arbor, MI 48105
| | - C A Figueroa
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI 48105; Department of Surgery, University of Michigan, Ann Arbor, MI 48105
| |
Collapse
|
7
|
Roy D, Mazumder O, Jaiswal D, Ghose A, Khandelwal S, Mandana K, Sinha A. In-silico cardiovascular hemodynamic model to simulate the effect of physical exercise. Biomed Signal Process Control 2023. [DOI: 10.1016/j.bspc.2022.104364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
8
|
Fan L, Sun Y, Choy JS, Kassab GS, Lee LC. Mechanism of exercise intolerance in heart diseases predicted by a computer model of myocardial demand-supply feedback system. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2022; 227:107188. [PMID: 36334525 DOI: 10.1016/j.cmpb.2022.107188] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Revised: 09/28/2022] [Accepted: 10/16/2022] [Indexed: 06/16/2023]
Abstract
BACKGROUND AND OBJECTIVE The myocardial demand-supply feedback system plays an important role in augmenting blood supply in response to exercise-induced increased myocardial demand. During this feedback process, the myocardium and coronary blood flow interact bidirectionally at many different levels. METHODS To investigate these interactions, a novel computational framework that considers the closed myocardial demand-supply feedback system was developed. In the framework coupling the systemic circulation of the left ventricle and coronary perfusion with regulation, myocardial work affects coronary perfusion via flow regulation mechanisms (e.g., metabolic regulation) and myocardial-vessel interactions, whereas coronary perfusion affects myocardial contractility in a closed feedback system. The framework was calibrated based on the measurements from healthy subjects under graded exercise conditions, and then was applied to simulate the effects of graded exercise on myocardial demand-supply under different physiological and pathological conditions. RESULTS We found that the framework can recapitulate key features found during exercise in clinical and animal studies. We showed that myocardial blood flow is increased but maximum hyperemia is reduced during exercise, which led to a reduction in coronary flow reserve. For coronary stenosis and myocardial inefficiency, the model predicts that an increase in heart rate is necessary to maintain the baseline cardiac output. Correspondingly, the resting coronary flow reserve is exhausted and the range of heart rate before exhaustion of coronary flow reserve is reduced. In the presence of metabolic regulation dysfunction, the model predicts that the metabolic vasodilator signal is higher at rest, saturates faster during exercise, and as a result, causes quicker exhaustion of coronary flow reserve. CONCLUSIONS Model predictions showed that the coronary flow reserve deteriorates faster during graded exercise, which in turn, suggests a decrease in exercise tolerance for patients with stenosis, myocardial inefficiency and metabolic flow regulation dysfunction. The findings in this study may have clinical implications in diagnosing cardiovascular diseases.
Collapse
Affiliation(s)
- Lei Fan
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA.
| | - Yuexing Sun
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA
| | - Jenny S Choy
- California Medical Innovations Institute, San Diego, CA, USA
| | | | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
9
|
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
|
10
|
Schollenberger J, Osborne NH, Hernandez-Garcia L, Figueroa CA. A Combined Computational Fluid Dynamics and Arterial Spin Labeling MRI Modeling Strategy to Quantify Patient-Specific Cerebral Hemodynamics in Cerebrovascular Occlusive Disease. Front Bioeng Biotechnol 2021; 9:722445. [PMID: 34485260 PMCID: PMC8416094 DOI: 10.3389/fbioe.2021.722445] [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: 06/08/2021] [Accepted: 07/19/2021] [Indexed: 11/13/2022] Open
Abstract
Cerebral hemodynamics in the presence of cerebrovascular occlusive disease (CVOD) are influenced by the anatomy of the intracranial arteries, the degree of stenosis, the patency of collateral pathways, and the condition of the cerebral microvasculature. Accurate characterization of cerebral hemodynamics is a challenging problem. In this work, we present a strategy to quantify cerebral hemodynamics using computational fluid dynamics (CFD) in combination with arterial spin labeling MRI (ASL). First, we calibrated patient-specific CFD outflow boundary conditions using ASL-derived flow splits in the Circle of Willis. Following, we validated the calibrated CFD model by evaluating the fractional blood supply from the main neck arteries to the vascular territories using Lagrangian particle tracking and comparing the results against vessel-selective ASL (VS-ASL). Finally, the feasibility and capability of our proposed method were demonstrated in two patients with CVOD and a healthy control subject. We showed that the calibrated CFD model accurately reproduced the fractional blood supply to the vascular territories, as obtained from VS-ASL. The two patients revealed significant differences in pressure drop over the stenosis, collateral flow, and resistance of the distal vasculature, despite similar degrees of clinical stenosis severity. Our results demonstrated the advantages of a patient-specific CFD analysis for assessing the hemodynamic impact of stenosis.
Collapse
Affiliation(s)
- Jonas Schollenberger
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
| | - Nicholas H Osborne
- Department of Surgery, University of Michigan, Ann Arbor, MI, United States
| | - Luis Hernandez-Garcia
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States.,Functional MRI Laboratory, University of Michigan, Ann Arbor, 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
|
11
|
Zhang H, Xia J, Yang Y, Yang Q, Song H, Xie J, Ma Y, Hou Y, Qiao A. Branch flow distribution approach and its application in the calculation of fractional flow reserve in stenotic coronary artery. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:5978-5994. [PMID: 34517519 DOI: 10.3934/mbe.2021299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
OBJECTIVE To calculate fractional flow reserve (FFR) based on computed tomography angiography (i.e., FFRCT) by considering the branch flow distribution in the coronary arteries. BACKGROUND FFR is the gold standard to diagnose myocardial ischemia caused by coronary stenosis. An accurate and noninvasive method for obtaining total coronary blood flow is needed for the calculation of FFRCT. METHODS A mathematical model for estimating the coronary blood flow rate and two approaches for setting the patient-specific flow boundary condition were proposed. Coronary branch flow distribution methods based on a volume-flow approach and a diameter-flow approach were employed for the numerical simulation of FFRCT. The values of simulated FFRCT for 16 patients were compared with their clinically measured FFR. RESULTS The ratio of total coronary blood flow to cardiac output and the myocardial blood flow under the condition of hyperemia were 16.97% and 4.07 mL/min/g, respectively. The errors of FFRCT compared with clinical data under the volume-flow approach and diameter-flow approach were 10.47% and 11.76%, respectively, the diagnostic accuracies of FFRCT were 65% and 85%, and the consistencies were 95% and 90%. CONCLUSIONS The mathematical model for estimating the coronary blood flow rate and the coronary branch flow distribution method can be applied to calculate the value of clinical noninvasive FFRCT.
Collapse
Affiliation(s)
- Honghui Zhang
- Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, China
- Intelligent Physiological Measurement and Clinical Translation, Beijing International Base for Scientific and Technological Cooperation, Beijing 100124, China
| | - Jun Xia
- Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, China
- Intelligent Physiological Measurement and Clinical Translation, Beijing International Base for Scientific and Technological Cooperation, Beijing 100124, China
| | - Yinlong Yang
- Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, China
- Intelligent Physiological Measurement and Clinical Translation, Beijing International Base for Scientific and Technological Cooperation, Beijing 100124, China
| | - Qingqing Yang
- Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, China
- Intelligent Physiological Measurement and Clinical Translation, Beijing International Base for Scientific and Technological Cooperation, Beijing 100124, China
| | - Hongfang Song
- School of Biomedical Engineering, Beijing Key Laboratory of Fundamental Research on Biomechanics in Clinical Application, Capital Medical University, Beijing 100069, China
| | - Jinjie Xie
- Beijing Anzhen Hospital, Capital Medical University, Beijing 100029, China
| | - Yue Ma
- Shengjing Hospital, China Medical University, Shenyang 110001, China
| | - Yang Hou
- Shengjing Hospital, China Medical University, Shenyang 110001, China
| | - Aike Qiao
- Faculty of Environment and Life, Beijing University of Technology, Beijing 100124, China
- Intelligent Physiological Measurement and Clinical Translation, Beijing International Base for Scientific and Technological Cooperation, Beijing 100124, China
| |
Collapse
|
12
|
CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation. PLoS Comput Biol 2021; 17:e1008881. [PMID: 33970900 PMCID: PMC8148362 DOI: 10.1371/journal.pcbi.1008881] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Revised: 05/25/2021] [Accepted: 03/17/2021] [Indexed: 12/16/2022] Open
Abstract
In this work, we describe the CRIMSON (CardiovasculaR Integrated Modelling and SimulatiON) software environment. CRIMSON provides a powerful, customizable and user-friendly system for performing three-dimensional and reduced-order computational haemodynamics studies via a pipeline which involves: 1) segmenting vascular structures from medical images; 2) constructing analytic arterial and venous geometric models; 3) performing finite element mesh generation; 4) designing, and 5) applying boundary conditions; 6) running incompressible Navier-Stokes simulations of blood flow with fluid-structure interaction capabilities; and 7) post-processing and visualizing the results, including velocity, pressure and wall shear stress fields. A key aim of CRIMSON is to create a software environment that makes powerful computational haemodynamics tools accessible to a wide audience, including clinicians and students, both within our research laboratories and throughout the community. The overall philosophy is to leverage best-in-class open source standards for medical image processing, parallel flow computation, geometric solid modelling, data assimilation, and mesh generation. It is actively used by researchers in Europe, North and South America, Asia, and Australia. It has been applied to numerous clinical problems; we illustrate applications of CRIMSON to real-world problems using examples ranging from pre-operative surgical planning to medical device design optimization.
Collapse
|
13
|
Sangha GS, Goergen CJ, Prior SJ, Ranadive SM, Clyne AM. Preclinical techniques to investigate exercise training in vascular pathophysiology. Am J Physiol Heart Circ Physiol 2021; 320:H1566-H1600. [PMID: 33385323 PMCID: PMC8260379 DOI: 10.1152/ajpheart.00719.2020] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Atherosclerosis is a dynamic process starting with endothelial dysfunction and inflammation and eventually leading to life-threatening arterial plaques. Exercise generally improves endothelial function in a dose-dependent manner by altering hemodynamics, specifically by increased arterial pressure, pulsatility, and shear stress. However, athletes who regularly participate in high-intensity training can develop arterial plaques, suggesting alternative mechanisms through which excessive exercise promotes vascular disease. Understanding the mechanisms that drive atherosclerosis in sedentary versus exercise states may lead to novel rehabilitative methods aimed at improving exercise compliance and physical activity. Preclinical tools, including in vitro cell assays, in vivo animal models, and in silico computational methods, broaden our capabilities to study the mechanisms through which exercise impacts atherogenesis, from molecular maladaptation to vascular remodeling. Here, we describe how preclinical research tools have and can be used to study exercise effects on atherosclerosis. We then propose how advanced bioengineering techniques can be used to address gaps in our current understanding of vascular pathophysiology, including integrating in vitro, in vivo, and in silico studies across multiple tissue systems and size scales. Improving our understanding of the antiatherogenic exercise effects will enable engaging, targeted, and individualized exercise recommendations to promote cardiovascular health rather than treating cardiovascular disease that results from a sedentary lifestyle.
Collapse
Affiliation(s)
- Gurneet S Sangha
- Fischell Department of Bioengineering, University of Maryland, College Park, Maryland
| | - Craig J Goergen
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana.,Purdue University Center for Cancer Research, Purdue University, West Lafayette, Indiana
| | - Steven J Prior
- Department of Kinesiology, University of Maryland School of Public Health, College Park, Maryland.,Baltimore Veterans Affairs Geriatric Research, Education, and Clinical Center, Baltimore, Maryland
| | - Sushant M Ranadive
- Department of Kinesiology, University of Maryland School of Public Health, College Park, Maryland
| | - Alisa M Clyne
- Fischell Department of Bioengineering, University of Maryland, College Park, Maryland
| |
Collapse
|
14
|
Arthurs CJ, Xiao N, Moireau P, Schaeffter T, Figueroa CA. A flexible framework for sequential estimation of model parameters in computational hemodynamics. ADVANCED MODELING AND SIMULATION IN ENGINEERING SCIENCES 2020; 7:48. [PMID: 33282681 PMCID: PMC7717067 DOI: 10.1186/s40323-020-00186-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2020] [Accepted: 11/06/2020] [Indexed: 06/02/2023]
Abstract
A major challenge in constructing three dimensional patient specific hemodynamic models is the calibration of model parameters to match patient data on flow, pressure, wall motion, etc. acquired in the clinic. Current workflows are manual and time-consuming. This work presents a flexible computational framework for model parameter estimation in cardiovascular flows that relies on the following fundamental contributions. (i) A Reduced-Order Unscented Kalman Filter (ROUKF) model for data assimilation for wall material and simple lumped parameter network (LPN) boundary condition model parameters. (ii) A constrained least squares augmentation (ROUKF-CLS) for more complex LPNs. (iii) A "Netlist" implementation, supporting easy filtering of parameters in such complex LPNs. The ROUKF algorithm is demonstrated using non-invasive patient-specific data on anatomy, flow and pressure from a healthy volunteer. The ROUKF-CLS algorithm is demonstrated using synthetic data on a coronary LPN. The methods described in this paper have been implemented as part of the CRIMSON hemodynamics software package.
Collapse
Affiliation(s)
| | - Nan Xiao
- Dept. of Biomedical Engineering, King’s College London, London, UK
| | - Philippe Moireau
- Inria, Inria Saclay-Ile de France, 91128 Palaiseau, France
- LMS, Ecole Polytechnique, CNRS, Institut Polytechnique de Paris, 91128 Palaiseau, France
| | - Tobias Schaeffter
- Physikalisch-Technische Bundesanstalt, Berlin, Germany
- Technical University Berlin, Berlin, Germany
| | - C. Alberto Figueroa
- Depts. of Surgery and Biomedical Engineering, University of Michigan, 2800 Plymouth Rd, Ann Arbor, MI 48109 USA
| |
Collapse
|
15
|
van Bakel TMJ, Arthurs CJ, Nauta FJH, Eagle KA, van Herwaarden JA, Moll FL, Trimarchi S, Patel HJ, Figueroa CA. Cardiac remodelling following thoracic endovascular aortic repair for descending aortic aneurysms. Eur J Cardiothorac Surg 2020; 55:1061-1070. [PMID: 30535179 DOI: 10.1093/ejcts/ezy399] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/24/2018] [Revised: 10/27/2018] [Accepted: 10/29/2018] [Indexed: 12/17/2022] Open
Abstract
OBJECTIVES Current endografts for thoracic endovascular aortic repair (TEVAR) are much stiffer than the aorta and have been shown to induce acute stiffening. In this study, we aimed to estimate the impact of TEVAR on left ventricular (LV) stroke work (SW) and mass using a non-invasive image-based workflow. METHODS The University of Michigan database was searched for patients treated with TEVAR for descending aortic pathologies (2013-2016). Patients with available pre-TEVAR and post-TEVAR computed tomography angiography and echocardiography data were selected. LV SW was estimated via patient-specific fluid-structure interaction analyses. LV remodelling was quantified through morphological measurements using echocardiography and electrocardiographic-gated computed tomography angiography data. RESULTS Eight subjects were included in this study, the mean age of the patients was 68 (73, 25) years, and 6 patients were women. All patients were prescribed antihypertensive drugs following TEVAR. The fluid-structure interaction simulations computed a 26% increase in LV SW post-TEVAR [0.94 (0.89, 0.34) J to 1.18 (1.11, 0.65) J, P = 0.012]. Morphological measurements revealed an increase in the LV mass index post-TEVAR of +26% in echocardiography [72 (73, 17) g/m2 to 91 (87, 26) g/m2, P = 0.017] and +15% in computed tomography angiography [52 (46, 29) g/m2 to 60 (57, 22) g/m2, P = 0.043]. The post- to pre-TEVAR LV mass index ratio was positively correlated with the post- to pre-TEVAR ratios of SW and the mean blood pressure (ρ = 0.690, P = 0.058 and ρ = 0.786, P = 0.021, respectively). CONCLUSIONS TEVAR was associated with increased LV SW and mass during follow-up. Medical device manufacturers should develop more compliant devices to reduce the stiffness mismatch with the aorta. Additionally, intensive antihypertensive management is needed to control blood pressure post-TEVAR.
Collapse
Affiliation(s)
- Theodorus M J van Bakel
- Department of Surgery, University of Michigan, Ann Arbor, MI, USA.,Department of Vascular Surgery, University Medical Center Utrecht, Utrecht, Netherlands.,Thoracic Aortic Research Center, Policlinico San Donato IRCCS, San Donato Milanese, Italy
| | - Christopher J Arthurs
- Division of Imaging Sciences and Biomedical Engineering, King's College London, London, UK
| | - Foeke J H Nauta
- Department of Surgery, University of Michigan, Ann Arbor, MI, USA.,Department of Vascular Surgery, University Medical Center Utrecht, Utrecht, Netherlands.,Thoracic Aortic Research Center, Policlinico San Donato IRCCS, San Donato Milanese, Italy
| | - Kim A Eagle
- Department of Cardiology, University of Michigan, Ann Arbor, MI, USA
| | | | - Frans L Moll
- Department of Vascular Surgery, University Medical Center Utrecht, Utrecht, Netherlands
| | - Santi Trimarchi
- Thoracic Aortic Research Center, Policlinico San Donato IRCCS, San Donato Milanese, Italy.,Department of Biomedical Sciences for Health, University of Milan, Milan, Italy
| | - Himanshu J Patel
- Department of Cardiac Surgery, University of Michigan, Ann Arbor, MI, USA
| | - C Alberto Figueroa
- Department of Surgery, University of Michigan, Ann Arbor, MI, USA.,Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
16
|
Tossas-Betancourt C, van Bakel TMJ, Arthurs CJ, Coleman DM, Eliason JL, Figueroa CA, Stanley JC. Computational analysis of renal artery flow characteristics by modeling aortoplasty and aortic bypass interventions for abdominal aortic coarctation. J Vasc Surg 2020; 71:505-516.e4. [PMID: 31153701 PMCID: PMC8409007 DOI: 10.1016/j.jvs.2019.02.063] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2018] [Accepted: 02/20/2019] [Indexed: 01/23/2023]
Abstract
OBJECTIVE Suprarenal abdominal aortic coarctation (SAAC) alters flow and pressure patterns to the kidneys and is often associated with severe angiotensin-mediated hypertension refractory to drug therapy. SAAC is most often treated by a thoracoabdominal bypass (TAB) or patch aortoplasty (PA). It is currently unclear what effect these interventions have on renal flow and pressure waveforms. This study, using retrospective data from a patient with SAAC subjected to a TAB, undertook computational modeling to analyze aortorenal blood flow preoperatively as well as postoperatively after a variety of TAB and PA interventions. METHODS Patient-specific anatomic models were constructed from preoperative computed tomography angiograms of a 9-year-old child with an isolated SAAC. Fluid-structure interaction (FSI) simulations of hemodynamics were performed to analyze preoperative renal flow and pressure waveforms. A parametric study was then performed to examine the hemodynamic impact of different bypass diameters and patch oversizing. RESULTS Preoperative FSI results documented diastole-dominated renal perfusion with considerable high-frequency disturbances in blood flow and pressure. The postoperative TAB right and left kidney volumes increased by 58% and 79%, respectively, reflecting the increased renal artery blood flows calculated by the FSI analysis. Postoperative increases in systolic flow accompanied decreases in high-frequency disturbances, aortic pressure, and collateral flow after all surgical interventions. In general, lesser degrees of high-frequency disturbances followed PA interventions. High-frequency disturbances were eliminated with the 0% PA in contrast to the 30% and 50% PA oversizing and TAB interventions, in which these flow disturbances remained. CONCLUSIONS Both TAB and PA dramatically improved renal artery flow and pressure waveforms, although disturbed renal waveforms remained in many of the surgical scenarios. Importantly, only the 0% PA oversizing scenario eliminated all high-frequency disturbances, resulting in nearly normal aortorenal blood flow. The study also establishes the relevance of patient-specific computational modeling in planning interventions for the midaortic syndrome.
Collapse
Affiliation(s)
| | | | - Christopher J Arthurs
- Division of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Dawn M Coleman
- Department of Surgery, University of Michigan, Ann Arbor, Mich
| | | | - C Alberto Figueroa
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, Mich; Department of Surgery, University of Michigan, Ann Arbor, Mich.
| | - James C Stanley
- Department of Surgery, University of Michigan, Ann Arbor, Mich
| |
Collapse
|
17
|
Sokolov SY, Volchkov SO, Bessonov IS, Chestukhin VV, Kurlyandskaya GV, Blyakhman FA. Estimation of Blood Flow Velocity in Coronary Arteries Based on the Movement of Radiopaque Agent. PATTERN RECOGNITION AND IMAGE ANALYSIS 2019. [DOI: 10.1134/s1054661819040163] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
18
|
Acuna A, Berman AG, Damen FW, Meyers BA, Adelsperger AR, Bayer KC, Brindise MC, Bungart B, Kiel AM, Morrison RA, Muskat JC, Wasilczuk KM, Wen Y, Zhang J, Zito P, Goergen CJ. Computational Fluid Dynamics of Vascular Disease in Animal Models. J Biomech Eng 2019; 140:2676341. [PMID: 29570754 DOI: 10.1115/1.4039678] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2017] [Indexed: 12/19/2022]
Abstract
Recent applications of computational fluid dynamics (CFD) applied to the cardiovascular system have demonstrated its power in investigating the impact of hemodynamics on disease initiation, progression, and treatment outcomes. Flow metrics such as pressure distributions, wall shear stresses (WSS), and blood velocity profiles can be quantified to provide insight into observed pathologies, assist with surgical planning, or even predict disease progression. While numerous studies have performed simulations on clinical human patient data, it often lacks prediagnosis information and can be subject to large intersubject variability, limiting the generalizability of findings. Thus, animal models are often used to identify and manipulate specific factors contributing to vascular disease because they provide a more controlled environment. In this review, we explore the use of CFD in animal models in recent studies to investigate the initiating mechanisms, progression, and intervention effects of various vascular diseases. The first section provides a brief overview of the CFD theory and tools that are commonly used to study blood flow. The following sections are separated by anatomical region, with the abdominal, thoracic, and cerebral areas specifically highlighted. We discuss the associated benefits and obstacles to performing CFD modeling in each location. Finally, we highlight animal CFD studies focusing on common surgical treatments, including arteriovenous fistulas (AVF) and pulmonary artery grafts. The studies included in this review demonstrate the value of combining CFD with animal imaging and should encourage further research to optimize and expand upon these techniques for the study of vascular disease.
Collapse
Affiliation(s)
- Andrea Acuna
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Alycia G Berman
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Frederick W Damen
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Brett A Meyers
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, IN 47907 e-mail:
| | - Amelia R Adelsperger
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Kelsey C Bayer
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Melissa C Brindise
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, IN 47907 e-mail:
| | - Brittani Bungart
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Alexander M Kiel
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Rachel A Morrison
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Joseph C Muskat
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Kelsey M Wasilczuk
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Yi Wen
- Department of Agricultural and Biological Engineering, Purdue University, 225 South University Street, West Lafayette, IN 47907 e-mail:
| | - Jiacheng Zhang
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, IN 47907 e-mail:
| | - Patrick Zito
- Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| | - Craig J Goergen
- ASME Membership Bioengineering Division, Weldon School of Biomedical Engineering, Purdue University, 206 S. Martin Jischke Drive, West Lafayette, IN 47907 e-mail:
| |
Collapse
|
19
|
Micheletti JK, Vanderlei FM, Machado AF, de Almeida AC, Nakamura FY, Netto Junior J, Pastre CM. A New Mathematical Approach to Explore the Post-exercise Recovery Process and Its Applicability in a Cold Water Immersion Protocol. J Strength Cond Res 2019; 33:1266-1275. [DOI: 10.1519/jsc.0000000000003041] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
20
|
Silva Vieira M, Arthurs CJ, Hussain T, Razavi R, Figueroa CA. Patient-specific modeling of right coronary circulation vulnerability post-liver transplant in Alagille's syndrome. PLoS One 2018; 13:e0205829. [PMID: 30408044 PMCID: PMC6224049 DOI: 10.1371/journal.pone.0205829] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2018] [Accepted: 10/02/2018] [Indexed: 01/27/2023] Open
Abstract
OBJECTIVES Cardiac output (CO) response to dobutamine can identify Alagille's syndrome (ALGS) patients at higher risk of cardiovascular complications during liver transplantation. We propose a novel patient-specific computational methodology to estimate the coronary autoregulatory responses during different hemodynamic conditions, including those experienced in a post-reperfusion syndrome (PRS), to aid cardiac risk-assessment. MATERIAL AND METHODS Data (pressure, flow, strain and ventricular volumes) from a 6-year-old ALGS patient undergoing catheter/dobutamine stress MRI (DSMRI) were used to parameterize a closed-loop coupled-multidomain (3D-0D) approach consisting of image-derived vascular models of pulmonary and systemic circulations and a series of 0D-lumped parameter networks (LPN) of the heart chambers and the distal arterial and venous circulations. A coronary microcirculation control model (CMCM) was designed to adjust the coronary resistance to match coronary blood flow (and thus oxygen delivery) with MVO2 requirements during Rest, Stress and a virtual PRS condition. RESULTS In all three simulated conditions, diastolic dominated right coronary artery (RCA) flow was observed, due to high right ventricle (RV) afterload. Despite a measured 45% increase in CO, impaired coronary flow reserve (CFR) (~1.4) at Stress was estimated by the CMCM. During modeled PRS, a marked vasodilatory response was insufficient to match RV myocardial oxygen requirements. Such exhaustion of the RCA autoregulatory response was not anticipated by the DSMRI study. CONCLUSION Impaired CFR undetected by DSMRI resulted in predicted myocardial ischemia in a computational model of PRS. This computational framework may identify ALGS patients at higher risk of complications during liver transplantation due to impaired coronary microvascular responses.
Collapse
Affiliation(s)
- Miguel Silva Vieira
- Division of Imaging Sciences & Biomedical Engineering, King's College London, London, United Kingdom
- * E-mail:
| | - Christopher J. Arthurs
- Division of Imaging Sciences & Biomedical Engineering, King's College London, London, United Kingdom
| | - Tarique Hussain
- Department of Pediatrics, University of Texas Southwestern Medical Center at Dallas, United States of America
- Pediatric Cardiology Department, Evelina Children’s Hospital London, Guy’s and St. Thomas’ NHS Foundation Trust, London, United Kingdom
| | - Reza Razavi
- Division of Imaging Sciences & Biomedical Engineering, King's College London, London, United Kingdom
- Pediatric Cardiology Department, Evelina Children’s Hospital London, Guy’s and St. Thomas’ NHS Foundation Trust, London, United Kingdom
| | - Carlos Alberto Figueroa
- Division of Imaging Sciences & Biomedical Engineering, King's College London, London, United Kingdom
- Departments of Surgery and Biomedical Engineering, University of Michigan, Michigan, United States of America
| |
Collapse
|
21
|
Ge X, Yin Z, Fan Y, Vassilevski Y, Liang F. A multi-scale model of the coronary circulation applied to investigate transmural myocardial flow. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e3123. [PMID: 29947132 DOI: 10.1002/cnm.3123] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Revised: 05/03/2018] [Accepted: 06/17/2018] [Indexed: 06/08/2023]
Abstract
Distribution of blood flow in myocardium is a key determinant of the localization and severity of myocardial ischemia under impaired coronary perfusion conditions. Previous studies have extensively demonstrated the transmural difference of ischemic vulnerability. However, it remains incompletely understood how transmural myocardial flow is regulated under in vivo conditions. In the present study, a computational model of the coronary circulation was developed to quantitatively evaluate the sensitivity of transmural flow distribution to various cardiovascular and hemodynamic factors. The model was further incorporated with the flow autoregulatory mechanism to simulate the regulation of myocardial flow in the presence of coronary artery stenosis. Numerical tests demonstrated that heart rate (HR), intramyocardial tissue pressure (Pim ), and coronary perfusion pressure (Pper ) were the major determinant factors for transmural flow distribution (evaluated by the subendocardial-to-subepicardial (endo/epi) flow ratio) and that the flow autoregulatory mechanism played an important compensatory role in preserving subendocardial perfusion against reduced Pper . Further analysis for HR variation-induced hemodynamic changes revealed that the rise in endo/epi flow ratio accompanying HR decrease was attributable not only to the prolongation of cardiac diastole relative to systole, but more predominantly to the fall in Pim . Moreover, it was found that Pim and Pper interfered with each other with respect to their influence on transmural flow distribution. These results demonstrate the interactive effects of various cardiovascular and hemodynamic factors on transmural myocardial flow, highlighting the importance of taking into account patient-specific conditions in the explanation of clinical observations.
Collapse
Affiliation(s)
- Xinyang Ge
- School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
- Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration (CISSE), Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Zhaofang Yin
- Department of Cardiology, Shanghai Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200011, China
| | - Yuqi Fan
- Department of Cardiology, Shanghai Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200011, China
| | - Yuri Vassilevski
- Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, 119333, Russia
- Moscow Institute of Physics and Technology, Dolgoprudny, 141700, Russia
- Sechenov University, Moscow, 119991, Russia
| | - Fuyou Liang
- School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
- Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration (CISSE), Shanghai Jiao Tong University, Shanghai, 200240, China
- Sechenov University, Moscow, 119991, Russia
| |
Collapse
|
22
|
Gadda G, Majka M, Zieliński P, Gambaccini M, Taibi A. A multiscale model for the simulation of cerebral and extracerebral blood flows and pressures in humans. Eur J Appl Physiol 2018; 118:2443-2454. [DOI: 10.1007/s00421-018-3971-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2018] [Accepted: 08/21/2018] [Indexed: 01/26/2023]
|
23
|
Non-invasive instantaneous wave-free ratio using coronary CT angiography: diagnostic performance for evaluation of ischaemia-causing coronary stenosis confirmed by invasive fractional flow reserve. Clin Radiol 2018; 73:983.e15-983.e22. [PMID: 30093066 DOI: 10.1016/j.crad.2018.07.098] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2017] [Accepted: 07/05/2018] [Indexed: 01/10/2023]
Abstract
AIM To determine the diagnostic performance of instantaneous wave-free ratio (iFR) derived from non-invasive coronary computed tomography angiography (CTA; iFRCT) for ischaemia-causing coronary stenosis, and to compare the diagnostic efficacy of iFRCT, CTA, and CTA plus iFRCT. METHODS AND MATERIALS Thirty-nine patients (55 vessels) with known or suspected coronary artery disease were included. All patients underwent invasive coronary angiography and fractional flow reserve (FFR) according to CTA findings and clinical indicators. The same raw data used for CTA were used to build patient-specific computed flow dynamic models and to calculate iFRCT. RESULTS On a vessel-based level, the correlation between iFRCT and FFR was moderate (r=0.65, p<0.05); the optimal iFRCT cut-off value was 0.85 based on an FFR cut-off value of 0.80, resulting in 85% sensitivity, 69% specificity, 61% positive predictive value (PPV), 89% negative predictive value (NPV), and 75% accuracy. The AUC showed significant differences between iFRCT and CTA (vessel-based: 0.84 versus 0.68; patient-based: 0.84 versus 0.62; both p<0.01). The accuracy of CTA combined with iFRCT was significantly increased compared to CTA alone for vessels with intermediate stenosis (83% versus 40%, p<0.01). CONCLUSION iFRCT showed better diagnostic performance than CTA. iFRCT may be a promising method for detection of ischaemia-causing coronary stenosis, even in vessels with intermediate stenosis.
Collapse
|
24
|
Spazzapan M, Sastry P, Dunning J, Nordsletten D, de Vecchi A. The Use of Biophysical Flow Models in the Surgical Management of Patients Affected by Chronic Thromboembolic Pulmonary Hypertension. Front Physiol 2018; 9:223. [PMID: 29593574 PMCID: PMC5859070 DOI: 10.3389/fphys.2018.00223] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Accepted: 02/28/2018] [Indexed: 11/21/2022] Open
Abstract
Introduction: Chronic Thromboembolic Pulmonary Hypertension (CTEPH) results from progressive thrombotic occlusion of the pulmonary arteries. It is treated by surgical removal of the occlusion, with success rates depending on the degree of microvascular remodeling. Surgical eligibility is influenced by the contributions of both the thrombus occlusion and microvasculature remodeling to the overall vascular resistance. Assessing this is challenging due to the high inter-individual variability in arterial morphology and physiology. We investigated the potential of patient-specific computational flow modeling to quantify pressure gradients in the pulmonary arteries of CTEPH patients to assist the decision-making process for surgical eligibility. Methods: Detailed segmentations of the pulmonary arteries were created from postoperative chest Computed Tomography scans of three CTEPH patients. A focal stenosis was included in the original geometry to compare the pre- and post-surgical hemodynamics. Three-dimensional flow simulations were performed on each morphology to quantify velocity-dependent pressure changes using a finite element solver coupled to terminal 2-element Windkessel models. In addition to transient flow simulations, a parametric modeling approach based on constant flow simulations is also proposed as faster technique to estimate relative pressure drops through the proximal pulmonary vasculature. Results: An asymmetrical flow split between left and right pulmonary arteries was observed in the stenosed models. Removing the proximal obstruction resulted in a reduction of the right-left pressure imbalance of up to 18%. Changes were also observed in the wall shear stresses and flow topology, where vortices developed in the stenosed model while the non-stenosed retained a helical flow. The predicted pressure gradients from constant flow simulations were consistent with the ones measured in the transient flow simulations. Conclusion: This study provides a proof of concept that patient-specific computational modeling can be used as a noninvasive tool for assisting surgical decisions in CTEPH based on hemodynamics metrics. Our technique enables determination of the proximal relative pressure, which could subsequently be compared to the total pressure drop to determine the degree of distal and proximal vascular resistance. In the longer term this approach has the potential to form the basis for a more quantitative classification system of CTEPH types.
Collapse
Affiliation(s)
- Martina Spazzapan
- King's College London, GKT School of Medical Education, London, United Kingdom
| | - Priya Sastry
- Cardiothoracic Surgery Unit, Papworth Hospital NHS Foundation Trust, Cambridge, United Kingdom
| | - John Dunning
- Cardiothoracic Surgery Unit, Papworth Hospital NHS Foundation Trust, Cambridge, United Kingdom
| | - David Nordsletten
- King's College London, School of Biomedical Engineering and Imaging Sciences, St. Thomas' Hospital, London, United Kingdom
| | - Adelaide de Vecchi
- King's College London, School of Biomedical Engineering and Imaging Sciences, St. Thomas' Hospital, London, United Kingdom
| |
Collapse
|
25
|
van Bakel TMJ, Lau KD, Hirsch-Romano J, Trimarchi S, Dorfman AL, Figueroa CA. Patient-Specific Modeling of Hemodynamics: Supporting Surgical Planning in a Fontan Circulation Correction. J Cardiovasc Transl Res 2018; 11:145-155. [PMID: 29313269 DOI: 10.1007/s12265-017-9781-x] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/08/2017] [Accepted: 12/26/2017] [Indexed: 01/29/2023]
Abstract
Computational fluid dynamics (CFD) is a modeling technique that enables calculation of the behavior of fluid flows in complex geometries. In cardiovascular medicine, CFD methods are being used to calculate patient-specific hemodynamics for a variety of applications, such as disease research, noninvasive diagnostics, medical device evaluation, and surgical planning. This paper provides a concise overview of the methods to perform patient-specific computational analyses using clinical data, followed by a case study where CFD-supported surgical planning is presented in a patient with Fontan circulation complicated by unilateral pulmonary arteriovenous malformations. In closing, the challenges for implementation and adoption of CFD modeling in clinical practice are discussed.
Collapse
Affiliation(s)
| | - Kevin D Lau
- Department of Surgery, University of Michigan, Ann Arbor, MI, USA
| | - Jennifer Hirsch-Romano
- University of Michigan C.S. Mott Children's Hospital Congenital Heart Center, Ann Arbor, MI, USA
| | - Santi Trimarchi
- Policlinico San Donato IRCCS, Thoracic Aortic Research Center, San Donato Milanese, Italy
| | - Adam L Dorfman
- University of Michigan C.S. Mott Children's Hospital Congenital Heart Center, Ann Arbor, MI, USA
| | - C Alberto Figueroa
- Department of Surgery, University of Michigan, Ann Arbor, MI, USA.,Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
26
|
Computational Fluid Dynamics and Additive Manufacturing to Diagnose and Treat Cardiovascular Disease. Trends Biotechnol 2017; 35:1049-1061. [PMID: 28942268 DOI: 10.1016/j.tibtech.2017.08.008] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2017] [Revised: 08/20/2017] [Accepted: 08/23/2017] [Indexed: 11/21/2022]
Abstract
Noninvasive engineering models are now being used for diagnosing and planning the treatment of cardiovascular disease. Techniques in computational modeling and additive manufacturing have matured concurrently, and results from simulations can inform and enable the design and optimization of therapeutic devices and treatment strategies. The emerging synergy between large-scale simulations and 3D printing is having a two-fold benefit: first, 3D printing can be used to validate the complex simulations, and second, the flow models can be used to improve treatment planning for cardiovascular disease. In this review, we summarize and discuss recent methods and findings for leveraging advances in both additive manufacturing and patient-specific computational modeling, with an emphasis on new directions in these fields and remaining open questions.
Collapse
|
27
|
Nauta FJH, Lau KD, Arthurs CJ, Eagle KA, Williams DM, Trimarchi S, Patel HJ, Figueroa CA. Computational Fluid Dynamics and Aortic Thrombus Formation Following Thoracic Endovascular Aortic Repair. Ann Thorac Surg 2017; 103:1914-1921. [PMID: 28063468 DOI: 10.1016/j.athoracsur.2016.09.067] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Revised: 08/28/2016] [Accepted: 09/19/2016] [Indexed: 01/10/2023]
Abstract
BACKGROUND We present the possible utility of computational fluid dynamics in the assessment of thrombus formation and virtual surgical planning illustrated in a patient with aortic thrombus in a kinked ascending aortic graft following thoracic endovascular aortic repair. METHODS A patient-specific three-dimensional model was built from computed tomography. Additionally, we modeled 3 virtual aortic interventions to assess their effect on thrombosis potential: (1) open surgical repair, (2) conformable endografting, and (3) single-branched endografting. Flow waveforms were extracted from echocardiography and used for the simulations. We used the computational index termed platelet activation potential (PLAP) representing accumulated shear rates of fluid particles within a fluid domain to assess thrombosis potential. RESULTS The baseline model revealed high PLAP in the entire arch (119.8 ± 42.5), with significantly larger PLAP at the thrombus location (125.4 ± 41.2, p < 0.001). Surgical repair showed a 37% PLAP reduction at the thrombus location (78.6 ± 25.3, p < 0.001) and a 24% reduction in the arch (91.6 ± 28.9, p < 0.001). Single-branched endografting reduced PLAP in the thrombus region by 20% (99.7 ± 24.6, p < 0.001) and by 14% in the arch (103.8 ± 26.1, p < 0.001), whereas a more conformable endograft did not have a profound effect, resulting in a modest 4% PLAP increase (130.6 ± 43.7, p < 0.001) in the thrombus region relative to the baseline case. CONCLUSIONS Regions of high PLAP were associated with aortic thrombus. Aortic repair resolved pathologic flow patterns, reducing PLAP. Branched endografting also relieved complex flow patterns reducing PLAP. Computational fluid dynamics may assist in the prediction of aortic thrombus formation in hemodynamically complex cases and help guide repair strategies.
Collapse
Affiliation(s)
- Foeke J H Nauta
- Department of Surgery, University of Michigan, Ann Arbor, Michigan, USA; Thoracic Aortic Research Center, Vascular Surgery, Policlinico San Donato IRCCS, University of Milan, Italy
| | - Kevin D Lau
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, Michigan, USA
| | - Christopher J Arthurs
- Division of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Kim A Eagle
- Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, USA
| | - David M Williams
- Department of Radiology, University of Michigan, Ann Arbor, Michigan, USA
| | - Santi Trimarchi
- Thoracic Aortic Research Center, Vascular Surgery, Policlinico San Donato IRCCS, University of Milan, Italy
| | - Himanshu J Patel
- Department of Cardiac Surgery, University of Michigan, Ann Arbor, Michigan, USA
| | - Carlos A Figueroa
- Department of Surgery, University of Michigan, Ann Arbor, Michigan, USA; Department of Biomedical Engineering, University of Michigan, Ann Arbor, Michigan, USA.
| |
Collapse
|
28
|
Arthurs CJ, Agarwal P, John AV, Dorfman AL, Grifka RG, Figueroa CA. Reproducing Patient-Specific Hemodynamics in the Blalock-Taussig Circulation Using a Flexible Multi-Domain Simulation Framework: Applications for Optimal Shunt Design. Front Pediatr 2017; 5:78. [PMID: 28491863 PMCID: PMC5405677 DOI: 10.3389/fped.2017.00078] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/11/2017] [Accepted: 03/31/2017] [Indexed: 11/13/2022] Open
Abstract
For babies born with hypoplastic left heart syndrome, several open-heart surgeries are required. During Stage I, a Norwood procedure is performed to construct an appropriate circulation to both the systemic and the pulmonary arteries. The pulmonary arteries receive flow from the systemic circulation, often using a Blalock-Taussig (BT) shunt between the innominate artery and the right pulmonary artery. This procedure causes significantly disturbed flow in the pulmonary arteries. In this study, we use computational hemodynamic simulations to demonstrate its capacity for examining the properties of the flow through and near the BT shunt. Initially, we construct a computational model which produces blood flow and pressure measurements matching the clinical magnetic resonance imaging (MRI) and catheterization data. Achieving this required us to determine the level of BT shunt occlusion; because the occlusion is below the MRI resolution, this information is difficult to recover without the aid of computational simulations. We determined that the shunt had undergone an effective diameter reduction of 22% since the time of surgery. Using the resulting geometric model, we show that we can computationally reproduce the clinical data. We, then, replace the BT shunt with a hypothetical alternative shunt design with a flare at the distal end. Investigation of the impact of the shunt design reveals that the flare can increase pulmonary pressure by as much as 7% and flow by as much as 9% in the main pulmonary branches, which may be beneficial to the pulmonary circulation.
Collapse
Affiliation(s)
- Christopher J Arthurs
- Division of Imaging Sciences and Biomedical Engineering, King's College London, King's Health Partners, St. Thomas' Hospital, London, UK
| | - Pradyumn Agarwal
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA
| | - Anna V John
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA
| | - Adam L Dorfman
- Department of Pediatric Cardiology, University of Michigan Health System, Ann Arbor, MI, USA
| | - Ronald G Grifka
- Department of Pediatric Cardiology, University of Michigan Health System, Ann Arbor, MI, USA
| | - C Alberto Figueroa
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA.,Department of Surgery, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|