1
|
Kojic M, Milosevic M, Simic V, Milicevic B, Terracciano R, Filgueira CS. On the generality of the finite element modeling physical fields in biological systems by the multiscale smeared concept (Kojic transport model). Heliyon 2024; 10:e26354. [PMID: 38434281 PMCID: PMC10907537 DOI: 10.1016/j.heliyon.2024.e26354] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 02/09/2024] [Accepted: 02/12/2024] [Indexed: 03/05/2024] Open
Abstract
The biomechanical and biochemical processes in the biological systems of living organisms are extremely complex. Advances in understanding these processes are mainly achieved by laboratory and clinical investigations, but in recent decades they are supported by computational modeling. Besides enormous efforts and achievements in this modeling, there still is a need for new methods that can be used in everyday research and medical practice. In this report, we give a view of the generality of the finite element methodology introduced by the first author and supported by his collaborators. It is based on the multiscale smeared physical fields, termed as Kojic Transport Model (KTM), published in several journal papers and summarized in a recent book (Kojic et al., 2022) [1]. We review relevant literature to demonstrate the distinctions and advantages of our methodology and indicate possible further applications. We refer to our published results by a selection of a few examples which include modeling of partitioning, blood flow, molecular transport within the pancreas, multiscale-multiphysics model of coupling electrical field and ion concentration, and a model of convective-diffusive transport within the lung parenchyma. Two new examples include a model of convective-diffusive transport within a growing tumor, and drug release from nanofibers with fiber degradation.
Collapse
Affiliation(s)
- Milos Kojic
- Houston Methodist Research Institute, The Department of Nanomedicine, 6670 Bertner Ave., R7 117, Houston, TX, 77030, USA
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia
- Serbian Academy of Sciences and Arts, Knez Mihailova 35, 11000, Belgrade, Serbia
| | - Miljan Milosevic
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia
- Institute of Information Technologies, University of Kragujevac, Department of Technical- Technological Sciences, Jovana Cvijica bb, 34000, Kragujevac, Serbia
- Belgrade Metropolitan University, Tadeusa Koscuska 63, 11000, Belgrade, Serbia
| | - Vladimir Simic
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia
- Institute of Information Technologies, University of Kragujevac, Department of Technical- Technological Sciences, Jovana Cvijica bb, 34000, Kragujevac, Serbia
| | - Bogdan Milicevic
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia
- Faculty of Engineering, University of Kragujevac, Kragujevac, 34000, Serbia
| | - Rossana Terracciano
- Houston Methodist Research Institute, The Department of Nanomedicine, 6670 Bertner Ave., R7 117, Houston, TX, 77030, USA
- Department of Electronics and Telecommunications, Politecnico di Torino, Torino, Italy
| | - Carly S. Filgueira
- Houston Methodist Research Institute, The Department of Nanomedicine, 6670 Bertner Ave., R7 117, Houston, TX, 77030, USA
- Department of Cardiovascular Surgery, Houston Methodist Research Institute, Houston, TX, 77030, USA
| |
Collapse
|
2
|
Zingaro A, Vergara C, Dede' L, Regazzoni F, Quarteroni A. A comprehensive mathematical model for cardiac perfusion. Sci Rep 2023; 13:14220. [PMID: 37648701 PMCID: PMC10469210 DOI: 10.1038/s41598-023-41312-0] [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/10/2023] [Accepted: 08/24/2023] [Indexed: 09/01/2023] Open
Abstract
The aim of this paper is to introduce a new mathematical model that simulates myocardial blood perfusion that accounts for multiscale and multiphysics features. Our model incorporates cardiac electrophysiology, active and passive mechanics, hemodynamics, valve modeling, and a multicompartment Darcy model of perfusion. We consider a fully coupled electromechanical model of the left heart that provides input for a fully coupled Navier-Stokes-Darcy Model for myocardial perfusion. The fluid dynamics problem is modeled in a left heart geometry that includes large epicardial coronaries, while the multicompartment Darcy model is set in a biventricular myocardium. Using a realistic and detailed cardiac geometry, our simulations demonstrate the biophysical fidelity of our model in describing cardiac perfusion. Specifically, we successfully validate the model reliability by comparing in-silico coronary flow rates and average myocardial blood flow with clinically established values ranges reported in relevant literature. Additionally, we investigate the impact of a regurgitant aortic valve on myocardial perfusion, and our results indicate a reduction in myocardial perfusion due to blood flow taken away by the left ventricle during diastole. To the best of our knowledge, our work represents the first instance where electromechanics, hemodynamics, and perfusion are integrated into a single computational framework.
Collapse
Affiliation(s)
- Alberto Zingaro
- MOX, Laboratory of Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy.
- ELEM Biotech S.L., Pier01, Palau de Mar, Plaça Pau Vila, 1, 08003, Barcelona, Spain.
| | - Christian Vergara
- LaBS, Dipartimento di Chimica, Materiali e Ingegneria Chimica "Giulio Natta", Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
| | - Luca Dede'
- MOX, Laboratory of Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
| | - Francesco Regazzoni
- MOX, Laboratory of Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
| | - Alfio Quarteroni
- MOX, Laboratory of Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
- Institute of Mathematics, École Polytechnique Fédérale de Lausanne, Station 8, Av. Piccard, CH-1015, Lausanne, Switzerland
| |
Collapse
|
3
|
In silico trials for treatment of acute ischemic stroke: Design and implementation. Comput Biol Med 2021; 137:104802. [PMID: 34520989 DOI: 10.1016/j.compbiomed.2021.104802] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Revised: 07/30/2021] [Accepted: 08/17/2021] [Indexed: 01/21/2023]
Abstract
An in silico trial simulates a disease and its corresponding therapies on a cohort of virtual patients to support the development and evaluation of medical devices, drugs, and treatment. In silico trials have the potential to refine, reduce cost, and partially replace current in vivo studies, namely clinical trials and animal testing. We present the design and implementation of an in silico trial for treatment of acute ischemic stroke. We propose an event-based modelling approach for the simulation of a disease and injury, where changes to the state of the system (the events) are assumed to be instantaneous. Using this approach we are able to combine a diverse set of models, spanning multiple time scales, to model acute ischemic stroke, treatment, and resulting brain tissue injury. The in silico trial is designed to be modular to aid development and reproducibility. It provides a comprehensive framework for application to any potential in silico trial. A statistical population model is used to generate cohorts of virtual patients. Patient functional outcomes are also predicted with a statistical model, using treatment and injury results and the patient's clinical parameters. We demonstrate the functionality of the event-based modelling approach and trial framework by running proof of concept in silico trials. The proof of concept trials simulate the same cohort of patients twice: once with successful treatment (successful recanalisation) and once with unsuccessful treatment (unsuccessful treatment). Ways to overcome some of the challenges and difficulties in setting up such an in silico trial are discussed, such as validation and computational limitations.
Collapse
|
4
|
Józsa TI, Padmos RM, El-Bouri WK, Hoekstra AG, Payne SJ. On the Sensitivity Analysis of Porous Finite Element Models for Cerebral Perfusion Estimation. Ann Biomed Eng 2021; 49:3647-3665. [PMID: 34155569 PMCID: PMC8671295 DOI: 10.1007/s10439-021-02808-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2021] [Accepted: 06/01/2021] [Indexed: 11/08/2022]
Abstract
Computational physiological models are promising tools to enhance the design of clinical trials and to assist in decision making. Organ-scale haemodynamic models are gaining popularity to evaluate perfusion in a virtual environment both in healthy and diseased patients. Recently, the principles of verification, validation, and uncertainty quantification of such physiological models have been laid down to ensure safe applications of engineering software in the medical device industry. The present study sets out to establish guidelines for the usage of a three-dimensional steady state porous cerebral perfusion model of the human brain following principles detailed in the verification and validation (V&V 40) standard of the American Society of Mechanical Engineers. The model relies on the finite element method and has been developed specifically to estimate how brain perfusion is altered in ischaemic stroke patients before, during, and after treatments. Simulations are compared with exact analytical solutions and a thorough sensitivity analysis is presented covering every numerical and physiological model parameter. The results suggest that such porous models can approximate blood pressure and perfusion distributions reliably even on a coarse grid with first order elements. On the other hand, higher order elements are essential to mitigate errors in volumetric blood flow rate estimation through cortical surface regions. Matching the volumetric flow rate corresponding to major cerebral arteries is identified as a validation milestone. It is found that inlet velocity boundary conditions are hard to obtain and that constant pressure inlet boundary conditions are feasible alternatives. A one-dimensional model is presented which can serve as a computationally inexpensive replacement of the three-dimensional brain model to ease parameter optimisation, sensitivity analyses and uncertainty quantification. The findings of the present study can be generalised to organ-scale porous perfusion models. The results increase the applicability of computational tools regarding treatment development for stroke and other cerebrovascular conditions.
Collapse
Affiliation(s)
- T I Józsa
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK.
| | - R M Padmos
- Computational Science Laboratory, Institute for Informatics, Faculty of Science, University of Amsterdam, Science Park 904, Amsterdam, 1098 XH, The Netherlands
| | - W K El-Bouri
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK.,Liverpool Centre for Cardiovascular Science, Department of Cardiovascular and Metabolic Medicine, University of Liverpool, Thomas Drive, Liverpool, L14 3PE, UK
| | - A G Hoekstra
- Computational Science Laboratory, Institute for Informatics, Faculty of Science, University of Amsterdam, Science Park 904, Amsterdam, 1098 XH, The Netherlands
| | - S J Payne
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK
| |
Collapse
|
5
|
El-Bouri WK, MacGowan A, Józsa TI, Gounis MJ, Payne SJ. Modelling the impact of clot fragmentation on the microcirculation after thrombectomy. PLoS Comput Biol 2021; 17:e1008515. [PMID: 33711015 PMCID: PMC7990195 DOI: 10.1371/journal.pcbi.1008515] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Revised: 03/24/2021] [Accepted: 02/23/2021] [Indexed: 12/29/2022] Open
Abstract
Many ischaemic stroke patients who have a mechanical removal of their clot (thrombectomy) do not get reperfusion of tissue despite the thrombus being removed. One hypothesis for this 'no-reperfusion' phenomenon is micro-emboli fragmenting off the large clot during thrombectomy and occluding smaller blood vessels downstream of the clot location. This is impossible to observe in-vivo and so we here develop an in-silico model based on in-vitro experiments to model the effect of micro-emboli on brain tissue. Through in-vitro experiments we obtain, under a variety of clot consistencies and thrombectomy techniques, micro-emboli distributions post-thrombectomy. Blood flow through the microcirculation is modelled for statistically accurate voxels of brain microvasculature including penetrating arterioles and capillary beds. A novel micro-emboli algorithm, informed by the experimental data, is used to simulate the impact of micro-emboli successively entering the penetrating arterioles and the capillary bed. Scaled-up blood flow parameters-permeability and coupling coefficients-are calculated under various conditions. We find that capillary beds are more susceptible to occlusions than the penetrating arterioles with a 4x greater drop in permeability per volume of vessel occluded. Individual microvascular geometries determine robustness to micro-emboli. Hard clot fragmentation leads to larger micro-emboli and larger drops in blood flow for a given number of micro-emboli. Thrombectomy technique has a large impact on clot fragmentation and hence occlusions in the microvasculature. As such, in-silico modelling of mechanical thrombectomy predicts that clot specific factors, interventional technique, and microvascular geometry strongly influence reperfusion of the brain. Micro-emboli are likely contributory to the phenomenon of no-reperfusion following successful removal of a major clot.
Collapse
Affiliation(s)
- Wahbi K. El-Bouri
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, United Kingdom
- Liverpool Centre for Cardiovascular Science, Department of Cardiovascular and Metabolic Medicine, University of Liverpool, Liverpool, United Kingdom
| | - Andrew MacGowan
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, United Kingdom
| | - Tamás I. Józsa
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, United Kingdom
| | - Matthew J. Gounis
- New England Center for Stroke Research, Department of Radiology, University of Massachusetts Medical School, Worcester, Massachusetts, United States of America
| | - Stephen J. Payne
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
6
|
Józsa TI, Padmos RM, Samuels N, El-Bouri WK, Hoekstra AG, Payne SJ. A porous circulation model of the human brain for in silico clinical trials in ischaemic stroke. Interface Focus 2021; 11:20190127. [PMID: 33343874 PMCID: PMC7739914 DOI: 10.1098/rsfs.2019.0127] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/23/2020] [Indexed: 12/30/2022] Open
Abstract
The advancement of ischaemic stroke treatment relies on resource-intensive experiments and clinical trials. In order to improve ischaemic stroke treatments, such as thrombolysis and thrombectomy, we target the development of computational tools for in silico trials which can partially replace these animal and human experiments with fast simulations. This study proposes a model that will serve as part of a predictive unit within an in silico clinical trial estimating patient outcome as a function of treatment. In particular, the present work aims at the development and evaluation of an organ-scale microcirculation model of the human brain for perfusion prediction. The model relies on a three-compartment porous continuum approach. Firstly, a fast and robust method is established to compute the anisotropic permeability tensors representing arterioles and venules. Secondly, vessel encoded arterial spin labelling magnetic resonance imaging and clustering are employed to create an anatomically accurate mapping between the microcirculation and large arteries by identifying superficial perfusion territories. Thirdly, the parameter space of the problem is reduced by analysing the governing equations and experimental data. Fourthly, a parameter optimization is conducted. Finally, simulations are performed with the tuned model to obtain perfusion maps corresponding to an open and an occluded (ischaemic stroke) scenario. The perfusion map in the occluded vessel scenario shows promising qualitative agreement with computed tomography images of a patient with ischaemic stroke caused by large vessel occlusion. The results highlight that in the case of vessel occlusion (i) identifying perfusion territories is essential to capture the location and extent of underperfused regions and (ii) anisotropic permeability tensors are required to give quantitatively realistic estimation of perfusion change. In the future, the model will be thoroughly validated against experiments.
Collapse
Affiliation(s)
- T. I. Józsa
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, UK
| | - R. M. Padmos
- Computational Science Laboratory, Institute for Informatics, Faculty of Science, University of Amsterdam, Science Park 904, Amsterdam 1098 XH, The Netherlands
| | - N. Samuels
- Department of Radiology and Nuclear Medicine, Erasmus MC, University Medical Center, Rotterdam 3015 GD, The Netherlands
| | - W. K. El-Bouri
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, UK
| | - A. G. Hoekstra
- Computational Science Laboratory, Institute for Informatics, Faculty of Science, University of Amsterdam, Science Park 904, Amsterdam 1098 XH, The Netherlands
| | - S. J. Payne
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, UK
| |
Collapse
|
7
|
Myocardial Perfusion Simulation for Coronary Artery Disease: A Coupled Patient-Specific Multiscale Model. Ann Biomed Eng 2020; 49:1432-1447. [PMID: 33263155 PMCID: PMC8057976 DOI: 10.1007/s10439-020-02681-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Accepted: 10/25/2020] [Indexed: 11/26/2022]
Abstract
Patient-specific models of blood flow are being used clinically to diagnose and plan treatment for coronary artery disease. A remaining challenge is bridging scales from flow in arteries to the micro-circulation supplying the myocardium. Previously proposed models are descriptive rather than predictive and have not been applied to human data. The goal here is to develop a multiscale patient-specific model enabling blood flow simulation from large coronary arteries to myocardial tissue. Patient vasculatures are segmented from coronary computed tomography angiography data and extended from the image-based model down to the arteriole level using a space-filling forest of synthetic trees. Blood flow is modeled by coupling a 1D model of the coronary arteries to a single-compartment Darcy myocardium model. Simulated results on five patients with non-obstructive coronary artery disease compare overall well to [\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^{15}$$\end{document}15O]\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\text {H}_{{2}}$$\end{document}H2O PET exam data for both resting and hyperemic conditions. Results on a patient with severe obstructive disease link coronary artery narrowing with impaired myocardial blood flow, demonstrating the model’s ability to predict myocardial regions with perfusion deficit. This is the first report of a computational model for simulating blood flow from the epicardial coronary arteries to the left ventricle myocardium applied to and validated on human data.
Collapse
|
8
|
Köppl T, Vidotto E, Wohlmuth B. A 3D-1D coupled blood flow and oxygen transport model to generate microvascular networks. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2020; 36:e3386. [PMID: 32659047 DOI: 10.1002/cnm.3386] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2020] [Revised: 05/18/2020] [Accepted: 07/05/2020] [Indexed: 06/11/2023]
Abstract
In this work, we introduce an algorithmic approach to generate microvascular networks starting from larger vessels that can be reconstructed without noticeable segmentation errors. Contrary to larger vessels, the reconstruction of fine-scale components of microvascular networks shows significant segmentation errors, and an accurate mapping is time and cost intense. Thus there is a need for fast and reliable reconstruction algorithms yielding surrogate networks having similar stochastic properties as the original ones. The microvascular networks are constructed in a marching way by adding vessels to the outlets of the vascular tree from the previous step. To optimise the structure of the vascular trees, we use Murray's law to determine the radii of the vessels and bifurcation angles. In each step, we compute the local gradient of the partial pressure of oxygen and adapt the orientation of the new vessels to this gradient. At the same time, we use the partial pressure of oxygen to check whether the considered tissue block is supplied sufficiently with oxygen. Computing the partial pressure of oxygen, we use a 3D-1D coupled model for blood flow and oxygen transport. To decrease the complexity of a fully coupled 3D model, we reduce the blood vessel network to a 1D graph structure and use a bi-directional coupling with the tissue which is described by a 3D homogeneous porous medium. The resulting surrogate networks are analysed with respect to morphological and physiological aspects.
Collapse
Affiliation(s)
- Tobias Köppl
- Chair for Numerics, University of Technology Munich, Garching, Germany
| | - Ettore Vidotto
- Chair for Numerics, University of Technology Munich, Garching, Germany
| | - Barbara Wohlmuth
- Chair for Numerics, University of Technology Munich, Garching, Germany
- Department of Mathematics, University of Bergen, Allegaten 41, 5020 Bergen, Norway, Germany
| |
Collapse
|
9
|
Diem AK, Valen-Sendstad K. Time-Lapsing Perfusion: Proof of Concept of a Novel Method to Study Drug Delivery in Whole Organs. Biophys J 2019; 117:2316-2323. [PMID: 31648790 DOI: 10.1016/j.bpj.2019.09.029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Revised: 08/29/2019] [Accepted: 09/23/2019] [Indexed: 11/18/2022] Open
Abstract
Perfusion is one of the most important processes maintaining organ health. From a computational perspective, however, perfusion is among the least-studied physiological processes of the heart. The recent development of novel nanoparticle-based targeted cardiac therapy calls for novel simulation methods that can provide insights into the distribution patterns of therapeutic agents within the heart tissue. Additionally, resolving the distribution patterns of perfusion is crucial for gaining a full understanding of the long-term impacts of cardiovascular diseases that can lead to adverse remodeling such as myocardial ischemia and heart failure. In this study, we have developed and used a, to our knowledge, novel particle-tracking-based method to simulate the perfusion-mediated distribution of nanoparticles or other solutes. To model blood flow through perfused tissue, we follow the approach of others and treat the tissue as a porous medium in a continuum model. Classically, solutes are modeled using reaction-advection-diffusion kinetics. However, because of the discrepancy of scales between advection and diffusion in blood vessels, this method becomes practically numerically unstable. Instead, we track a bolus of solutes or nanoparticles using particle tracking based purely on advection in arteries. In capillaries, we employ diffusion kinetics, using an effective diffusion coefficient to mimic capillary blood flow. We first demonstrate the numerical validity and computational efficiency of this method on a two-dimensional benchmark problem. Finally, we demonstrate how the method is used to visualize perfusion patterns of a healthy and ischemic human left ventricle geometry. The efficiency of the method allows for nanoparticle tracking over multiple cardiac cycles using a conventional laptop, providing a framework for the simulation of experimentally relevant timeframes to advance preclinical research.
Collapse
Affiliation(s)
- Alexandra K Diem
- Department of Computational Physiology, Simula Research Laboratory, Fornebu, Norway.
| | | |
Collapse
|
10
|
Gkontra P, El‐Bouri WK, Norton K, Santos A, Popel AS, Payne SJ, Arroyo AG. Dynamic Changes in Microvascular Flow Conductivity and Perfusion After Myocardial Infarction Shown by Image-Based Modeling. J Am Heart Assoc 2019; 8:e011058. [PMID: 30897998 PMCID: PMC6509718 DOI: 10.1161/jaha.118.011058] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/25/2018] [Accepted: 02/19/2019] [Indexed: 12/26/2022]
Abstract
Background Microcirculation is a decisive factor in tissue reperfusion inadequacy following myocardial infarction ( MI ). Nonetheless, experimental assessment of blood flow in microcirculation remains a bottleneck. We sought to model blood flow properties in coronary microcirculation at different time points after MI and to compare them with healthy conditions to obtain insights into alterations in cardiac tissue perfusion. Methods and Results We developed an image-based modeling framework that permitted feeding a continuum flow model with anatomical data previously obtained from the pig coronary microvasculature to calculate physiologically meaningful permeability tensors. The tensors encompassed the microvascular conductivity and were also used to estimate the arteriole-venule drop in pressure and myocardial blood flow. Our results indicate that the tensors increased in a bimodal pattern at infarcted areas on days 1 and 7 after MI while a nonphysiological decrease in arteriole-venule drop in pressure was observed; contrary, the tensors and the arteriole-venule drop in pressure on day 3 after MI , and in remote areas, were closer to values for healthy tissue. Myocardial blood flow calculated using the condition-dependent arteriole-venule drop in pressure decreased in infarcted areas. Last, we simulated specific modes of vascular remodeling, such as vasodilation, vasoconstriction, or pruning, and quantified their distinct impact on microvascular conductivity. Conclusions Our study unravels time- and region-dependent alterations of tissue perfusion related to the structural changes occurring in the coronary microvasculature due to MI . It also paves the way for conducting simulations in new therapeutic interventions in MI and for image-based microvascular modeling by applying continuum flow models in other biomedical scenarios.
Collapse
Affiliation(s)
- Polyxeni Gkontra
- Centro Nacional de Investigaciones Cardiovasculares (CNIC)MadridSpain
- Biomedical Image Technologies (BIT), ETSI Telecomunicación,Universidad Politécnica de MadridMadridSpain
| | - Wahbi K. El‐Bouri
- Institute of Biomedical EngineeringDepartment of Engineering ScienceUniversity of OxfordUnited Kingdom
| | - Kerri‐Ann Norton
- Division of Science, Mathematics, and ComputingBard CollegeAnnandale‐on‐HudsonNY
| | - Andrés Santos
- Biomedical Image Technologies (BIT), ETSI Telecomunicación,Universidad Politécnica de MadridMadridSpain
- Centro de Investigación Biomédica en Red de BioingenieríaBiomateriales y Nanomedicina (CIBERBBN)MadridSpain
| | - Aleksander S. Popel
- Department of Biomedical EngineeringSchool of MedicineJohns Hopkins UniversityBaltimoreMD
| | - Stephen J. Payne
- Institute of Biomedical EngineeringDepartment of Engineering ScienceUniversity of OxfordUnited Kingdom
| | - Alicia G. Arroyo
- Centro Nacional de Investigaciones Cardiovasculares (CNIC)MadridSpain
- Centro de Investigaciones Biológicas (CIB‐CSIC)MadridSpain
| |
Collapse
|
11
|
Smith AF, Doyeux V, Berg M, Peyrounette M, Haft-Javaherian M, Larue AE, Slater JH, Lauwers F, Blinder P, Tsai P, Kleinfeld D, Schaffer CB, Nishimura N, Davit Y, Lorthois S. Brain Capillary Networks Across Species: A few Simple Organizational Requirements Are Sufficient to Reproduce Both Structure and Function. Front Physiol 2019; 10:233. [PMID: 30971935 PMCID: PMC6444172 DOI: 10.3389/fphys.2019.00233] [Citation(s) in RCA: 52] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Accepted: 02/22/2019] [Indexed: 02/02/2023] Open
Abstract
Despite the key role of the capillaries in neurovascular function, a thorough characterization of cerebral capillary network properties is currently lacking. Here, we define a range of metrics (geometrical, topological, flow, mass transfer, and robustness) for quantification of structural differences between brain areas, organs, species, or patient populations and, in parallel, digitally generate synthetic networks that replicate the key organizational features of anatomical networks (isotropy, connectedness, space-filling nature, convexity of tissue domains, characteristic size). To reach these objectives, we first construct a database of the defined metrics for healthy capillary networks obtained from imaging of mouse and human brains. Results show that anatomical networks are topologically equivalent between the two species and that geometrical metrics only differ in scaling. Based on these results, we then devise a method which employs constrained Voronoi diagrams to generate 3D model synthetic cerebral capillary networks that are locally randomized but homogeneous at the network-scale. With appropriate choice of scaling, these networks have equivalent properties to the anatomical data, demonstrated by comparison of the defined metrics. The ability to synthetically replicate cerebral capillary networks opens a broad range of applications, ranging from systematic computational studies of structure-function relationships in healthy capillary networks to detailed analysis of pathological structural degeneration, or even to the development of templates for fabrication of 3D biomimetic vascular networks embedded in tissue-engineered constructs.
Collapse
Affiliation(s)
- Amy F Smith
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France
| | - Vincent Doyeux
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France
| | - Maxime Berg
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France
| | - Myriam Peyrounette
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France
| | - Mohammad Haft-Javaherian
- Nancy E. and Peter C. Meinig School of Biomedical Engineering, Cornell University, Ithaca, NY, United States
| | - Anne-Edith Larue
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France
| | - John H Slater
- Department of Biomedical Engineering, University of Delaware, Newark, DE, United States
| | - Frédéric Lauwers
- Toulouse NeuroImaging Center (TONIC), Université de Toulouse, INSERM, Toulouse, France.,Department of Anatomy, LSR44, Faculty of Medicine Toulouse-Purpan, Toulouse, France
| | - Pablo Blinder
- Department of Neurobiology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel Aviv, Israel
| | - Philbert Tsai
- Department of Physics, University of California, San Diego, La Jolla, CA, United States
| | - David Kleinfeld
- Department of Physics, University of California, San Diego, La Jolla, CA, United States
| | - Chris B Schaffer
- Nancy E. and Peter C. Meinig School of Biomedical Engineering, Cornell University, Ithaca, NY, United States
| | - Nozomi Nishimura
- Nancy E. and Peter C. Meinig School of Biomedical Engineering, Cornell University, Ithaca, NY, United States
| | - Yohan Davit
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France
| | - Sylvie Lorthois
- Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France.,Nancy E. and Peter C. Meinig School of Biomedical Engineering, Cornell University, Ithaca, NY, United States
| |
Collapse
|
12
|
Rohan E, Lukeš V, Jonášová A. Modeling of the contrast-enhanced perfusion test in liver based on the multi-compartment flow in porous media. J Math Biol 2018; 77:421-454. [PMID: 29368273 DOI: 10.1007/s00285-018-1209-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2016] [Revised: 01/15/2018] [Indexed: 12/20/2022]
Abstract
The paper deals with modeling the liver perfusion intended to improve quantitative analysis of the tissue scans provided by the contrast-enhanced computed tomography (CT). For this purpose, we developed a model of dynamic transport of the contrast fluid through the hierarchies of the perfusion trees. Conceptually, computed time-space distributions of the so-called tissue density can be compared with the measured data obtained from CT; such a modeling feedback can be used for model parameter identification. The blood flow is characterized at several scales for which different models are used. Flows in upper hierarchies represented by larger branching vessels are described using simple 1D models based on the Bernoulli equation extended by correction terms to respect the local pressure losses. To describe flows in smaller vessels and in the tissue parenchyma, we propose a 3D continuum model of porous medium defined in terms of hierarchically matched compartments characterized by hydraulic permeabilities. The 1D models corresponding to the portal and hepatic veins are coupled with the 3D model through point sources, or sinks. The contrast fluid saturation is governed by transport equations adapted for the 1D and 3D flow models. The complex perfusion model has been implemented using the finite element and finite volume methods. We report numerical examples computed for anatomically relevant geometries of the liver organ and of the principal vascular trees. The simulated tissue density corresponding to the CT examination output reflects a pathology modeled as a localized permeability deficiency.
Collapse
Affiliation(s)
- Eduard Rohan
- NTIS - New Technologies for the Information Society, Faculty of Applied Sciences, University of West Bohemia, Univerzitní 8, 30614, Pilsen, Czech Republic.
| | - Vladimír Lukeš
- NTIS - New Technologies for the Information Society, Faculty of Applied Sciences, University of West Bohemia, Univerzitní 8, 30614, Pilsen, Czech Republic
| | - Alena Jonášová
- NTIS - New Technologies for the Information Society, Faculty of Applied Sciences, University of West Bohemia, Univerzitní 8, 30614, Pilsen, Czech Republic
| |
Collapse
|
13
|
Clinical Diagnostic Biomarkers from the Personalization of Computational Models of Cardiac Physiology. Ann Biomed Eng 2015; 44:46-57. [PMID: 26399986 DOI: 10.1007/s10439-015-1439-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2015] [Accepted: 08/25/2015] [Indexed: 10/23/2022]
Abstract
Computational modelling of the heart is rapidly advancing to the point of clinical utility. However, the difficulty of parameterizing and validating models from clinical data indicates that the routine application of truly predictive models remains a significant challenge. We argue there is significant value in an intermediate step towards prediction. This step is the use of biophysically based models to extract clinically useful information from existing patient data. Specifically in this paper we review methodologies for applying modelling frameworks for this goal in the areas of quantifying cardiac anatomy, estimating myocardial stiffness and optimizing measurements of coronary perfusion. Using these indicative examples of the general overarching approach, we finally discuss the value, ongoing challenges and future potential for applying biophysically based modelling in the clinical context.
Collapse
|
14
|
Sinclair MD, Lee J, Cookson AN, Rivolo S, Hyde ER, Smith NP. Measurement and modeling of coronary blood flow. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2015; 7:335-56. [PMID: 26123867 DOI: 10.1002/wsbm.1309] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2015] [Revised: 05/19/2015] [Accepted: 05/21/2015] [Indexed: 01/10/2023]
Abstract
Ischemic heart disease that comprises both coronary artery disease and microvascular disease is the single greatest cause of death globally. In this context, enhancing our understanding of the interaction of coronary structure and function is not only fundamental for advancing basic physiology but also crucial for identifying new targets for treating these diseases. A central challenge for understanding coronary blood flow is that coronary structure and function exhibit different behaviors across a range of spatial and temporal scales. While experimental studies have sought to understand this feature by isolating specific mechanisms, in tandem, computational modeling is increasingly also providing a unique framework to integrate mechanistic behaviors across different scales. In addition, clinical methods for assessing coronary disease severity are continuously being informed and updated by findings in basic physiology. Coupling these technologies, computational modeling of the coronary circulation is emerging as a bridge between the experimental and clinical domains, providing a framework to integrate imaging and measurements from multiple sources with mathematical descriptions of governing physical laws. State-of-the-art computational modeling is being used to combine mechanistic models with data to provide new insight into coronary physiology, optimization of medical technologies, and new applications to guide clinical practice.
Collapse
Affiliation(s)
- Matthew D Sinclair
- Division of Imaging Sciences and Biomedical Engineering, British Heart Foundation (BHF) Centre of Excellence, King's College London, London, UK
| | - Jack Lee
- Division of Imaging Sciences and Biomedical Engineering, British Heart Foundation (BHF) Centre of Excellence, King's College London, London, UK
| | - Andrew N Cookson
- Division of Imaging Sciences and Biomedical Engineering, British Heart Foundation (BHF) Centre of Excellence, King's College London, London, UK
| | - Simone Rivolo
- Division of Imaging Sciences and Biomedical Engineering, British Heart Foundation (BHF) Centre of Excellence, King's College London, London, UK
| | - Eoin R Hyde
- Division of Imaging Sciences and Biomedical Engineering, British Heart Foundation (BHF) Centre of Excellence, King's College London, London, UK
| | - Nicolas P Smith
- Division of Imaging Sciences and Biomedical Engineering, British Heart Foundation (BHF) Centre of Excellence, King's College London, London, UK.,Department of Engineering, University of Auckland, Auckland, New Zealand
| |
Collapse
|
15
|
Sinclair M, Lee J, Schuster A, Chiribiri A, van den Wijngaard J, van Horssen P, Siebes M, Spaan JAE, Nagel E, Smith NP. Microsphere skimming in the porcine coronary arteries: Implications for flow quantification. Microvasc Res 2015; 100:59-70. [PMID: 25963318 DOI: 10.1016/j.mvr.2015.04.005] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2015] [Revised: 03/28/2015] [Accepted: 04/17/2015] [Indexed: 11/25/2022]
Abstract
Particle skimming is a phenomenon where particles suspended in fluid flowing through vessels distribute disproportionately to bulk fluid volume at junctions. Microspheres are considered a gold standard of intra-organ perfusion measurements and are used widely in studies of flow distribution and quantification. It has previously been hypothesised that skimming at arterial junctions is responsible for a systematic over-estimation of myocardial perfusion from microspheres at the subendocardium. Our objective is to integrate coronary arterial structure and microsphere distribution, imaged at high resolution, to test the hypothesis of microsphere skimming in a porcine left coronary arterial (LCA) network. A detailed network was reconstructed from cryomicrotome imaging data and a Poiseuille flow model was used to simulate flow. A statistical approach using Clopper-Pearson confidence intervals was applied to determine the prevalence of skimming at bifurcations in the LCA. Results reveal that microsphere skimming is most prevalent at bifurcations in the larger coronary arteries, namely the epicardial and transmural arteries. Bifurcations at which skimming was identified have significantly more asymmetric branching parameters. This finding suggests that when using thin transmural segments to quantify flow from microspheres, a skimming-related deposition bias may result in underestimation of perfusion in the subepicardium, and overestimation in the subendocardium.
Collapse
Affiliation(s)
- Matthew Sinclair
- Division of Imaging Sciences and Biomedical Engineering, King's College London, British Heart Foundation (BHF) Centre of Excellence, UK; National Institute of Heath Research (NIHR) Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust, Lambeth Wing, St. Thomas' Hospital, UK; Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC) Medical Engineering Centre, Lambeth Wing, St. Thomas' Hospital, London, UK
| | - Jack Lee
- Division of Imaging Sciences and Biomedical Engineering, King's College London, British Heart Foundation (BHF) Centre of Excellence, UK; National Institute of Heath Research (NIHR) Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust, Lambeth Wing, St. Thomas' Hospital, UK; Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC) Medical Engineering Centre, Lambeth Wing, St. Thomas' Hospital, London, UK
| | - Andreas Schuster
- Division of Imaging Sciences and Biomedical Engineering, King's College London, British Heart Foundation (BHF) Centre of Excellence, UK; National Institute of Heath Research (NIHR) Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust, Lambeth Wing, St. Thomas' Hospital, UK; Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC) Medical Engineering Centre, Lambeth Wing, St. Thomas' Hospital, London, UK; Department of Cardiology and Pneumology, Georg-August-University, Göttingen, Germany; German Centre for Cardiovascular Research (DZHK, Partner Site Göttingen), Göttingen, Germany
| | - Amedeo Chiribiri
- Division of Imaging Sciences and Biomedical Engineering, King's College London, British Heart Foundation (BHF) Centre of Excellence, UK; National Institute of Heath Research (NIHR) Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust, Lambeth Wing, St. Thomas' Hospital, UK; Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC) Medical Engineering Centre, Lambeth Wing, St. Thomas' Hospital, London, UK
| | - Jeroen van den Wijngaard
- Department of Biomedical Engineering & Physics, Academic Medical Centre, Amsterdam, The Netherlands
| | - Pepijn van Horssen
- Department of Biomedical Engineering & Physics, Academic Medical Centre, Amsterdam, The Netherlands
| | - Maria Siebes
- Department of Biomedical Engineering & Physics, Academic Medical Centre, Amsterdam, The Netherlands
| | - Jos A E Spaan
- Department of Biomedical Engineering & Physics, Academic Medical Centre, Amsterdam, The Netherlands
| | - Eike Nagel
- Division of Imaging Sciences and Biomedical Engineering, King's College London, British Heart Foundation (BHF) Centre of Excellence, UK; National Institute of Heath Research (NIHR) Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust, Lambeth Wing, St. Thomas' Hospital, UK; Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC) Medical Engineering Centre, Lambeth Wing, St. Thomas' Hospital, London, UK
| | - Nicolas P Smith
- Division of Imaging Sciences and Biomedical Engineering, King's College London, British Heart Foundation (BHF) Centre of Excellence, UK; National Institute of Heath Research (NIHR) Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust, Lambeth Wing, St. Thomas' Hospital, UK; Wellcome Trust and Engineering and Physical Sciences Research Council (EPSRC) Medical Engineering Centre, Lambeth Wing, St. Thomas' Hospital, London, UK; Department of Engineering, University of Auckland, Auckland, New Zealand.
| |
Collapse
|
16
|
Pan Q, Wang R, Reglin B, Cai G, Yan J, Pries AR, Ning G. A one-dimensional mathematical model for studying the pulsatile flow in microvascular networks. J Biomech Eng 2014; 136:011009. [PMID: 24190506 DOI: 10.1115/1.4025879] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2013] [Indexed: 11/08/2022]
Abstract
Techniques that model microvascular hemodynamics have been developed for decades. While the physiological significance of pressure pulsatility is acknowledged, most of the microcirculatory models use steady flow approaches. To theoretically study the extent and transmission of pulsatility in microcirculation, dynamic models need to be developed. In this paper, we present a one-dimensional model to describe the dynamic behavior of microvascular blood flow. The model is applied to a microvascular network from a rat mesentery. Intravital microscopy was used to record the morphology and flow velocities in individual vessel segments, and boundaries are defined according to the experimental data. The system of governing equations constituting the model is solved numerically using the discontinuous Galerkin method. An implicit integration scheme is adopted to increase computing efficiency. The model allows the simulation of the dynamic properties of blood flow in microcirculatory networks, including the pressure pulsatility (quantified by a pulsatility index) and pulse wave velocity (PWV). From the main input arteriole to the main output venule, the pulsatility index decreases by 66.7%. PWV obtained along arterioles declines with decreasing diameters, with mean values of 77.16, 25.31, and 8.30 cm/s for diameters of 26.84, 17.46, and 13.33 μm, respectively. These results suggest that the 1D model developed is able to simulate the characteristics of pressure pulsatility and wave propagation in complex microvascular networks.
Collapse
|
17
|
Cookson AN, Lee J, Michler C, Chabiniok R, Hyde E, Nordsletten D, Smith NP. A spatially-distributed computational model to quantify behaviour of contrast agents in MR perfusion imaging. Med Image Anal 2014; 18:1200-16. [PMID: 25103922 PMCID: PMC4156310 DOI: 10.1016/j.media.2014.07.002] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2013] [Revised: 07/07/2014] [Accepted: 07/08/2014] [Indexed: 10/30/2022]
Abstract
Contrast agent enhanced magnetic resonance (MR) perfusion imaging provides an early, non-invasive indication of defects in the coronary circulation. However, the large variation of contrast agent properties, physiological state and imaging protocols means that optimisation of image acquisition is difficult to achieve. This situation motivates the development of a computational framework that, in turn, enables the efficient mapping of this parameter space to provide valuable information for optimisation of perfusion imaging in the clinical context. For this purpose a single-compartment porous medium model of capillary blood flow is developed which is coupled with a scalar transport model, to characterise the behaviour of both blood-pool and freely-diffusive contrast agents characterised by their ability to diffuse through the capillary wall into the extra-cellular space. A parameter space study is performed on the nondimensionalised equations using a 2D model for both healthy and diseased myocardium, examining the sensitivity of system behaviour to Peclet number, Damköhler number (Da), diffusivity ratio and fluid porosity. Assuming a linear MR signal response model, sample concentration time series data are calculated, and the sensitivity of clinically-relevant properties of these signals to the model parameters is quantified. Both upslope and peak values display significant non-monotonic behaviour with regard to the Damköhler number, with these properties showing a high degree of sensitivity in the parameter range relevant to contrast agents currently in use. However, the results suggest that signal upslope is the more robust and discerning metric for perfusion quantification, in particular for correlating with perfusion defect size. Finally, the results were examined in the context of nonlinear signal response, flow quantification via Fermi deconvolution and perfusion reserve index, which demonstrated that there is no single best set of contrast agent parameters, instead the contrast agents should be tailored to the specific imaging protocol and post-processing method to be used.
Collapse
Affiliation(s)
- A N Cookson
- Department of Biomedical Engineering, Division of Imaging Sciences & Biomedical Engineering, St. Thomas' Hospital, King's College London, London SE1 7EH, UK
| | - J Lee
- Department of Biomedical Engineering, Division of Imaging Sciences & Biomedical Engineering, St. Thomas' Hospital, King's College London, London SE1 7EH, UK
| | - C Michler
- Department of Biomedical Engineering, Division of Imaging Sciences & Biomedical Engineering, St. Thomas' Hospital, King's College London, London SE1 7EH, UK
| | - R Chabiniok
- Department of Biomedical Engineering, Division of Imaging Sciences & Biomedical Engineering, St. Thomas' Hospital, King's College London, London SE1 7EH, UK
| | - E Hyde
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - D Nordsletten
- Department of Biomedical Engineering, Division of Imaging Sciences & Biomedical Engineering, St. Thomas' Hospital, King's College London, London SE1 7EH, UK
| | - N P Smith
- Department of Biomedical Engineering, Division of Imaging Sciences & Biomedical Engineering, St. Thomas' Hospital, King's College London, London SE1 7EH, UK.
| |
Collapse
|
18
|
Smith AF, Shipley RJ, Lee J, Sands GB, LeGrice IJ, Smith NP. Transmural variation and anisotropy of microvascular flow conductivity in the rat myocardium. Ann Biomed Eng 2014; 42:1966-77. [PMID: 24866569 DOI: 10.1007/s10439-014-1028-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Accepted: 05/03/2014] [Indexed: 11/30/2022]
Abstract
Transmural variations in the relationship between structural and fluid transport properties of myocardial capillary networks are determined via continuum modeling approaches using recent three-dimensional (3D) data on the microvascular structure. Specifically, the permeability tensor, which quantifies the inverse of the blood flow resistivity of the capillary network, is computed by volume-averaging flow solutions in synthetic networks with geometrical and topological properties derived from an anatomically-detailed microvascular data set extracted from the rat myocardium. Results show that the permeability is approximately ten times higher in the principal direction of capillary alignment (the "longitudinal" direction) than perpendicular to this direction, reflecting the strong anisotropy of the microvascular network. Additionally, a 30% increase in capillary diameter from subepicardium to subendocardium is shown to translate to a 130% transmural rise in permeability in the longitudinal capillary direction. This result supports the hypothesis that perfusion is preferentially facilitated during diastole in the subendocardial microvasculature to compensate for the severely-reduced systolic perfusion in the subendocardium.
Collapse
Affiliation(s)
- Amy F Smith
- Oxford Centre for Collaborative Applied Mathematics, Mathematical Institute, University of Oxford, Woodstock Road, Oxford, OX2 6GG, UK
| | | | | | | | | | | |
Collapse
|
19
|
Hyde ER, Cookson AN, Lee J, Michler C, Goyal A, Sochi T, Chabiniok R, Sinclair M, Nordsletten DA, Spaan J, van den Wijngaard JPHM, Siebes M, Smith NP. Multi-Scale Parameterisation of a Myocardial Perfusion Model Using Whole-Organ Arterial Networks. Ann Biomed Eng 2013; 42:797-811. [DOI: 10.1007/s10439-013-0951-y] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2013] [Accepted: 11/20/2013] [Indexed: 01/13/2023]
|