1
|
El Moçayd N, Belhamadia Y, Seaid M. Unsupervised stochastic learning and reduced order modeling for global sensitivity analysis in cardiac electrophysiology models. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2024; 255:108311. [PMID: 39032242 DOI: 10.1016/j.cmpb.2024.108311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2024] [Revised: 06/25/2024] [Accepted: 06/26/2024] [Indexed: 07/23/2024]
Abstract
BACKGROUND AND OBJECTIVE Numerical simulations in electrocardiology are often affected by various uncertainties inherited from the lack of precise knowledge regarding input values including those related to the cardiac cell model, domain geometry, and boundary or initial conditions used in the mathematical modeling. Conventional techniques for uncertainty quantification in modeling electrical activities of the heart encounter significant challenges, primarily due to the high computational costs associated with fine temporal and spatial scales. Additionally, the need for numerous model evaluations to quantify ubiquitous uncertainties increases the computational challenges even further. METHODS In the present study, we propose a non-intrusive surrogate model to perform uncertainty quantification and global sensitivity analysis in cardiac electrophysiology models. The proposed method combines an unsupervised machine learning technique with the polynomial chaos expansion to reconstruct a surrogate model for the propagation and quantification of uncertainties in the electrical activity of the heart. The proposed methodology not only accurately quantifies uncertainties at a very low computational cost but more importantly, it captures the targeted quantity of interest as either the whole spatial field or the whole temporal period. In order to perform sensitivity analysis, aggregated Sobol indices are estimated directly from the spectral mode of the polynomial chaos expansion. RESULTS We conduct Uncertainty Quantification (UQ) and global Sensitivity Analysis (SA) considering both spatial and temporal variations, rather than limiting the analysis to specific Quantities of Interest (QoIs). To assess the comprehensive performance of our methodology in simulating cardiac electrical activity, we utilize the monodomain model. Additionally, sensitivity analysis is performed on the parameters of the Mitchell-Schaeffer cell model. CONCLUSIONS Unlike conventional techniques for uncertainty quantification in modeling electrical activities, the proposed methodology performs at a low computational cost the sensitivity analysis on the cardiac electrical activity parameters. The results are fully reproducible and easily accessible, while the proposed reduced-order model represents a significant contribution to enhancing global sensitivity analysis in cardiac electrophysiology.
Collapse
Affiliation(s)
- Nabil El Moçayd
- College of Agriculture and Environmental Sciences, University Mohammed VI Polytechnique, Ben Guerir, Morocco.
| | - Youssef Belhamadia
- Department of Mathematics and Statistics, American University of Sharjah, United Arab Emirates.
| | - Mohammed Seaid
- Department of Engineering, University of Durham, South Road, Durham DH1 3LE, United Kingdom.
| |
Collapse
|
2
|
Djoumessi R, Dongmo Vougmo I, Tadjonang Tegne J, Pelap F. Proposed cardio-pulmonary model to investigate the effects of COVID-19 on the cardiovascular system. Heliyon 2023; 9:e12908. [PMID: 36644674 PMCID: PMC9830904 DOI: 10.1016/j.heliyon.2023.e12908] [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: 05/04/2022] [Revised: 12/28/2022] [Accepted: 01/06/2023] [Indexed: 01/12/2023] Open
Abstract
In this paper, we propose a new mathematical model of cardiovascular system coupled with a respiratory system to study the effects of COVID-19 on global blood circulation parameters using the lumped parameters model. We use the fourth-order Runge-Kutta method for solving the sets of equations of motion. We validate our model by showing that the simulated flows in pulmonary and aortic valves corroborate, respectively, the results established by Smith et al. [IFAC Proceedings Volumes, 39 (2006) 453-458]. Then we examine the effects of the new coronavirus (covid-19) on the cardiopulmonary system through the impact of the high respiratory frequency and the variation of the alveoli volume. To achieve this aim, we propose a new exponential law for the time varying of the pulmonary resistance. It appears that when the respiratory frequency grows, the delay between the systemic artery flow and the flow in the pulmonary artery diminishes. Therefore, the efficiency of the cardiac pump is reduced. Moreover, our results also show that variations of the alveoli volume cause the increment of the pleural pressure in the vascular cavities that induces an exponential growth of the pulmonary resistance. Furthermore, this growth of the pulmonary resistance provokes the augmentation of pressure in some organs and its reduction in others. We found that patient with covid-19 having a prior history of cardiovascular diseases is exposed to a severe case of inflammation/damage of certain organs than those with no history of cardiovascular disease.
Collapse
|
3
|
Yogev D, Tejman-Yarden S, Feinberg O, Parmet Y, Goldberg T, Illouz S, Nagar N, Freidin D, Vazgovsky O, Chatterji S, Salem Y, Katz U, Goitein O. Proof of concept: Comparative accuracy of semiautomated VR modeling for volumetric analysis of the heart ventricles. Heliyon 2022; 8:e11250. [PMID: 36387466 PMCID: PMC9641195 DOI: 10.1016/j.heliyon.2022.e11250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2022] [Revised: 10/12/2022] [Accepted: 10/20/2022] [Indexed: 11/06/2022] Open
Abstract
Introduction Simpson's rule is generally used to estimate cardiac volumes. By contrast, modern methods such as Virtual Reality (VR) utilize mesh modeling to present the object's surface spatial structure, thus enabling intricate volumetric calculations. In this study, two types of semiautomated VR models for cardiac volumetric analysis were compared to the standard Philips dedicated cardiac imaging platform (PDP) which is based on Simpson's rule calculations. Methods This retrospective report examined the cardiac computed tomography angiography (CCTA) of twenty patients with atrial fibrillation obtained prior to a left atrial appendage occlusion procedure. We employed two VR models to evaluate each CCTA and compared them to the PDP: a VR model with Philips-similar segmentations (VR-PS) that included the trabeculae and the papillary muscles within the luminal volume, and a VR model that only included the inner blood pool (VR-IBP). Results Comparison of the VR-PS and the PDP left ventricle (LV) volumes demonstrated excellent correlation with a ρc of 0.983 (95% CI 0.96, 0.99), and a small mean difference and range. The calculated volumes of the right ventricle (RV) had a somewhat lower correlation of 0.89 (95% CI 0.781, 0.95), a small mean difference, and a broader range. The VR-IBP chamber size estimations were significantly smaller than the estimates based on the PDP. Discussion Simpson's rule and polygon summation algorithms produce similar results in normal morphological LVs. However, this correlation failed to emerge when applied to RVs and irregular chambers. Conclusions The findings suggest that the polygon summation method is preferable for RV and irregular LV volume and function calculations.
Collapse
Affiliation(s)
- David Yogev
- The Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
| | - Shai Tejman-Yarden
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
- The Edmond J. Safra International Congenital Heart Center, Sheba Medical Center, Ramat Gan, Israel
- Corresponding author.
| | - Omer Feinberg
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
| | - Yisrael Parmet
- Department of Industrial Engineering and Management, Ben Gurion University, Beer Sheva, Israel
| | - Tomer Goldberg
- The Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Shay Illouz
- The Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
| | - Netanel Nagar
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
- Industrial Design Department, Bezalel Academy of Art and Design, Jerusalem, Israel
| | - Dor Freidin
- The Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
| | - Oliana Vazgovsky
- The Engineering Medical Research Lab, Sheba Medical Center, Ramat Gan, Israel
- The Edmond J. Safra International Congenital Heart Center, Sheba Medical Center, Ramat Gan, Israel
| | - Sumit Chatterji
- The Pulmonology Unit, Sheba Medical Center, Ramat Gan, Israel
- Interventional Pulmonology Unit, Sheba Medical Center, Ramat Gan, Israel
| | - Yishay Salem
- The Edmond J. Safra International Congenital Heart Center, Sheba Medical Center, Ramat Gan, Israel
- The Leviev Heart Institute, Sheba Medical Center, Ramat Gan, Israel
| | - Uriel Katz
- The Edmond J. Safra International Congenital Heart Center, Sheba Medical Center, Ramat Gan, Israel
- The Leviev Heart Institute, Sheba Medical Center, Ramat Gan, Israel
| | - Orly Goitein
- Department of Diagnostic Imaging, Sheba Medical Center, Ramat Gan, Israel
| |
Collapse
|
4
|
Lourenço WDJ, Reis RF, Ruiz-Baier R, Rocha BM, dos Santos RW, Lobosco M. A Poroelastic Approach for Modelling Myocardial Oedema in Acute Myocarditis. Front Physiol 2022; 13:888515. [PMID: 35860652 PMCID: PMC9289286 DOI: 10.3389/fphys.2022.888515] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 05/27/2022] [Indexed: 12/28/2022] Open
Abstract
Myocarditis is a general set of mechanisms that manifest themselves into the inflammation of the heart muscle. In 2017, more than 3 million people were affected by this disease worldwide, causing about 47,000 deaths. Many aspects of the origin of this disease are well known, but several important questions regarding the disease remain open. One of them is why some patients develop a significantly localised inflammation while others develop a much more diffuse inflammation, reaching across large portions of the heart. Furthermore, the specific role of the pathogenic agent that causes inflammation as well as the interaction with the immune system in the progression of the disease are still under discussion. Providing answers to these crucial questions can have an important impact on patient treatment. In this scenario, computational methods can aid specialists to understand better the relationships between pathogens and the immune system and elucidate why some patients develop diffuse myocarditis. This paper alters a recently developed model to study the myocardial oedema formation in acute infectious myocarditis. The model describes the finite deformation regime using partial differential equations to represent tissue displacement, fluid pressure, fluid phase, and the concentrations of pathogens and leukocytes. A sensitivity analysis was performed to understand better the influence of the most relevant model parameters on the disease dynamics. The results showed that the poroelastic model could reproduce local and diffuse myocarditis dynamics in simplified and complex geometrical domains.
Collapse
Affiliation(s)
- Wesley de Jesus Lourenço
- Graduate Program on Computational Modelling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
| | - Ruy Freitas Reis
- Department of Computer Science, Institute of Exact Sciences, Federal University of Juiz de Fora, Juiz de Fora, Brazil
| | - Ricardo Ruiz-Baier
- School of Mathematics and Victorian Heart Institute, Monash University, Melbourne, VIC, Australia
- Research Core on Natural and Exact Sciences, Universidad Adventista de Chile, Chillán, Chile
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, Moscow, Russia
| | - Bernardo Martins Rocha
- Graduate Program on Computational Modelling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
- Department of Computer Science, Institute of Exact Sciences, Federal University of Juiz de Fora, Juiz de Fora, Brazil
| | - Rodrigo Weber dos Santos
- Graduate Program on Computational Modelling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
- Department of Computer Science, Institute of Exact Sciences, Federal University of Juiz de Fora, Juiz de Fora, Brazil
| | - Marcelo Lobosco
- Graduate Program on Computational Modelling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
- Department of Computer Science, Institute of Exact Sciences, Federal University of Juiz de Fora, Juiz de Fora, Brazil
- *Correspondence: Marcelo Lobosco,
| |
Collapse
|
5
|
Lazarus A, Dalton D, Husmeier D, Gao H. Sensitivity analysis and inverse uncertainty quantification for the left ventricular passive mechanics. Biomech Model Mechanobiol 2022; 21:953-982. [PMID: 35377030 PMCID: PMC9132878 DOI: 10.1007/s10237-022-01571-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Accepted: 02/28/2022] [Indexed: 01/08/2023]
Abstract
Personalized computational cardiac models are considered to be a unique and powerful tool in modern cardiology, integrating the knowledge of physiology, pathology and fundamental laws of mechanics in one framework. They have the potential to improve risk prediction in cardiac patients and assist in the development of new treatments. However, in order to use these models for clinical decision support, it is important that both the impact of model parameter perturbations on the predicted quantities of interest as well as the uncertainty of parameter estimation are properly quantified, where the first task is a priori in nature (meaning independent of any specific clinical data), while the second task is carried out a posteriori (meaning after specific clinical data have been obtained). The present study addresses these challenges for a widely used constitutive law of passive myocardium (the Holzapfel-Ogden model), using global sensitivity analysis (SA) to address the first challenge, and inverse-uncertainty quantification (I-UQ) for the second challenge. The SA is carried out on a range of different input parameters to a left ventricle (LV) model, making use of computationally efficient Gaussian process (GP) surrogate models in place of the numerical forward simulator. The results of the SA are then used to inform a low-order reparametrization of the constitutive law for passive myocardium under consideration. The quality of this parameterization in the context of an inverse problem having observed noisy experimental data is then quantified with an I-UQ study, which again makes use of GP surrogate models. The I-UQ is carried out in a Bayesian manner using Markov Chain Monte Carlo, which allows for full uncertainty quantification of the material parameter estimates. Our study reveals insights into the relation between SA and I-UQ, elucidates the dependence of parameter sensitivity and estimation uncertainty on external factors, like LV cavity pressure, and sheds new light on cardio-mechanic model formulation, with particular focus on the Holzapfel-Ogden myocardial model.
Collapse
Affiliation(s)
- Alan Lazarus
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - David Dalton
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Dirk Husmeier
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| |
Collapse
|
6
|
Borowska A, Gao H, Lazarus A, Husmeier D. Bayesian optimisation for efficient parameter inference in a cardiac mechanics model of the left ventricle. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2022; 38:e3593. [PMID: 35302293 PMCID: PMC9285944 DOI: 10.1002/cnm.3593] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Accepted: 03/12/2022] [Indexed: 06/14/2023]
Abstract
We consider parameter inference in cardio-mechanic models of the left ventricle, in particular the one based on the Holtzapfel-Ogden (HO) constitutive law, using clinical in vivo data. The equations underlying these models do not admit closed form solutions and hence need to be solved numerically. These numerical procedures are computationally expensive making computational run times associated with numerical optimisation or sampling excessive for the uptake of the models in the clinical practice. To address this issue, we adopt the framework of Bayesian optimisation (BO), which is an efficient statistical technique of global optimisation. BO seeks the optimum of an unknown black-box function by sequentially training a statistical surrogate-model and using it to select the next query point by leveraging the associated exploration-exploitation trade-off. To guarantee that the estimates based on the in vivo data are realistic also for high-pressures, unobservable in vivo, we include a penalty term based on a previously published empirical law developed using ex vivo data. Two case studies based on real data demonstrate that the proposed BO procedure outperforms the state-of-the-art inference algorithm for the HO constitutive law.
Collapse
Affiliation(s)
| | - Hao Gao
- School of Mathematics and StatisticsUniversity of GlasgowGlasgowUK
| | - Alan Lazarus
- School of Mathematics and StatisticsUniversity of GlasgowGlasgowUK
| | - Dirk Husmeier
- School of Mathematics and StatisticsUniversity of GlasgowGlasgowUK
| |
Collapse
|
7
|
Kovacheva E, Gerach T, Schuler S, Ochs M, Dössel O, Loewe A. Causes of altered ventricular mechanics in hypertrophic cardiomyopathy: an in-silico study. Biomed Eng Online 2021; 20:69. [PMID: 34294108 PMCID: PMC8296558 DOI: 10.1186/s12938-021-00900-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Accepted: 06/08/2021] [Indexed: 11/22/2022] Open
Abstract
BACKGROUND Hypertrophic cardiomyopathy (HCM) is typically caused by mutations in sarcomeric genes leading to cardiomyocyte disarray, replacement fibrosis, impaired contractility, and elevated filling pressures. These varying tissue properties are associated with certain strain patterns that may allow to establish a diagnosis by means of non-invasive imaging without the necessity of harmful myocardial biopsies or contrast agent application. With a numerical study, we aim to answer: how the variability in each of these mechanisms contributes to altered mechanics of the left ventricle (LV) and if the deformation obtained in in-silico experiments is comparable to values reported from clinical measurements. METHODS We conducted an in-silico sensitivity study on physiological and pathological mechanisms potentially underlying the clinical HCM phenotype. The deformation of the four-chamber heart models was simulated using a finite-element mechanical solver with a sliding boundary condition to mimic the tissue surrounding the heart. Furthermore, a closed-loop circulatory model delivered the pressure values acting on the endocardium. Deformation measures and mechanical behavior of the heart models were evaluated globally and regionally. RESULTS Hypertrophy of the LV affected the course of strain, strain rate, and wall thickening-the root-mean-squared difference of the wall thickening between control (mean thickness 10 mm) and hypertrophic geometries (17 mm) was >10%. A reduction of active force development by 40% led to less overall deformation: maximal radial strain reduced from 26 to 21%. A fivefold increase in tissue stiffness caused a more homogeneous distribution of the strain values among 17 heart segments. Fiber disarray led to minor changes in the circumferential and radial strain. A combination of pathological mechanisms led to reduced and slower deformation of the LV and halved the longitudinal shortening of the LA. CONCLUSIONS This study uses a computer model to determine the changes in LV deformation caused by pathological mechanisms that are presumed to underlay HCM. This knowledge can complement imaging-derived information to obtain a more accurate diagnosis of HCM.
Collapse
Affiliation(s)
- Ekaterina Kovacheva
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Kaiserstr. 12, 76131, Karlsruhe, Germany
| | - Tobias Gerach
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Kaiserstr. 12, 76131, Karlsruhe, Germany
| | - Steffen Schuler
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Kaiserstr. 12, 76131, Karlsruhe, Germany
| | - Marco Ochs
- Department of Cardiology, Theresienkrankenhaus, Academic Teaching Hospital of Heidelberg University, Bassermannstr.1, 68165, Mannheim, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Kaiserstr. 12, 76131, Karlsruhe, Germany
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Kaiserstr. 12, 76131, Karlsruhe, Germany.
| |
Collapse
|
8
|
Pagani S, Manzoni A. Enabling forward uncertainty quantification and sensitivity analysis in cardiac electrophysiology by reduced order modeling and machine learning. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3450. [PMID: 33599106 PMCID: PMC8244126 DOI: 10.1002/cnm.3450] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Revised: 02/05/2021] [Accepted: 02/07/2021] [Indexed: 06/12/2023]
Abstract
We present a new, computationally efficient framework to perform forward uncertainty quantification (UQ) in cardiac electrophysiology. We consider the monodomain model to describe the electrical activity in the cardiac tissue, coupled with the Aliev-Panfilov model to characterize the ionic activity through the cell membrane. We address a complete forward UQ pipeline, including both: (i) a variance-based global sensitivity analysis for the selection of the most relevant input parameters, and (ii) a way to perform uncertainty propagation to investigate the impact of intra-subject variability on outputs of interest depending on the cardiac potential. Both tasks exploit stochastic sampling techniques, thus implying overwhelming computational costs because of the huge amount of queries to the high-fidelity, full-order computational model obtained by approximating the coupled monodomain/Aliev-Panfilov system through the finite element method. To mitigate this computational burden, we replace the full-order model with computationally inexpensive projection-based reduced-order models (ROMs) aimed at reducing the state-space dimensionality. Resulting approximation errors on the outputs of interest are finally taken into account through artificial neural network (ANN)-based models, enhancing the accuracy of the whole UQ pipeline. Numerical results show that the proposed physics-based ROMs outperform regression-based emulators relying on ANNs built with the same amount of training data, in terms of both numerical accuracy and overall computational efficiency.
Collapse
Affiliation(s)
- Stefano Pagani
- MOX, Dipartimento di MatematicaPolitecnico di MilanoMilanItaly
| | - Andrea Manzoni
- MOX, Dipartimento di MatematicaPolitecnico di MilanoMilanItaly
| |
Collapse
|
9
|
Fedele M, Quarteroni A. Polygonal surface processing and mesh generation tools for the numerical simulation of the cardiac function. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3435. [PMID: 33415829 PMCID: PMC8244076 DOI: 10.1002/cnm.3435] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Revised: 01/01/2021] [Accepted: 01/02/2021] [Indexed: 06/05/2023]
Abstract
In order to simulate the cardiac function for a patient-specific geometry, the generation of the computational mesh is crucially important. In practice, the input is typically a set of unprocessed polygonal surfaces coming either from a template geometry or from medical images. These surfaces need ad-hoc processing to be suitable for a volumetric mesh generation. In this work we propose a set of new algorithms and tools aiming to facilitate the mesh generation process. In particular, we focus on different aspects of a cardiac mesh generation pipeline: (1) specific polygonal surface processing for cardiac geometries, like connection of different heart chambers or segmentation outputs; (2) generation of accurate boundary tags; (3) definition of mesh-size functions dependent on relevant geometric quantities; (4) processing and connecting together several volumetric meshes. The new algorithms-implemented in the open-source software vmtk-can be combined with each other allowing the creation of personalized pipelines, that can be optimized for each cardiac geometry or for each aspect of the cardiac function to be modeled. Thanks to these features, the proposed tools can significantly speed-up the mesh generation process for a large range of cardiac applications, from single-chamber single-physics simulations to multi-chambers multi-physics simulations. We detail all the proposed algorithms motivating them in the cardiac context and we highlight their flexibility by showing different examples of cardiac mesh generation pipelines.
Collapse
Affiliation(s)
- Marco Fedele
- MOX, Department of MathematicsPolitecnico di MilanoMilanItaly
| | - Alfio Quarteroni
- MOX, Department of MathematicsPolitecnico di MilanoMilanItaly
- Institute of MathematicsÉcole Polytechnique Fédérale de LausanneLausanneSwitzerland
| |
Collapse
|
10
|
Reis RF, de Melo Quintela B, de Oliveira Campos J, Gomes JM, Rocha BM, Lobosco M, Weber dos Santos R. Characterization of the COVID-19 pandemic and the impact of uncertainties, mitigation strategies, and underreporting of cases in South Korea, Italy, and Brazil. CHAOS, SOLITONS, AND FRACTALS 2020; 136:109888. [PMID: 32412556 PMCID: PMC7221372 DOI: 10.1016/j.chaos.2020.109888] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Revised: 05/06/2020] [Accepted: 05/10/2020] [Indexed: 05/09/2023]
Abstract
By April 7th, 2020, the Coronavirus disease 2019 (COVID-19) has infected one and a half million people worldwide, accounting for over 80 thousand of deaths in 209 countries and territories around the world. The new and fast dynamics of the pandemic are challenging the health systems of different countries. In the absence of vaccines or effective treatments, mitigation policies, such as social isolation and lock-down of cities, have been adopted, but the results vary among different countries. Some countries were able to control the disease at the moment, as is the case of South Korea. Others, like Italy, are now experiencing the peak of the pandemic. Finally, countries with emerging economies and social issues, like Brazil, are in the initial phase of the pandemic. In this work, we use mathematical models with time-dependent coefficients, techniques of inverse and forward uncertainty quantification, and sensitivity analysis to characterize essential aspects of the COVID-19 in the three countries mentioned above. The model parameters estimated for South Korea revealed effective social distancing and isolation policies, border control, and a high number in the percentage of reported cases. In contrast, underreporting of cases was estimated to be very high in Brazil and Italy. In addition, the model estimated a poor isolation policy at the moment in Brazil, with a reduction of contact around 40%, whereas Italy and South Korea estimated numbers for contact reduction are at 75% and 90%, respectively. This characterization of the COVID-19, in these different countries under different scenarios and phases of the pandemic, supports the importance of mitigation policies, such as social distancing. In addition, it raises serious concerns for socially and economically fragile countries, where underreporting poses additional challenges to the management of the COVID-19 pandemic by significantly increasing the uncertainties regarding its dynamics.
Collapse
Affiliation(s)
- Ruy Freitas Reis
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Brazil
| | - Bárbara de Melo Quintela
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Brazil
- Department of Industrial Engineering, Alma Mater Studiorum - University of Bologna, Italy
| | | | - Johnny Moreira Gomes
- Pós-Graduação em Modelagem Computacional, Universidade Federal de Juiz de Fora, Brazil
| | - Bernardo Martins Rocha
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Brazil
- Pós-Graduação em Modelagem Computacional, Universidade Federal de Juiz de Fora, Brazil
| | - Marcelo Lobosco
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Brazil
- Pós-Graduação em Modelagem Computacional, Universidade Federal de Juiz de Fora, Brazil
| | - Rodrigo Weber dos Santos
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Brazil
- Pós-Graduação em Modelagem Computacional, Universidade Federal de Juiz de Fora, Brazil
| |
Collapse
|
11
|
Campos JO, Sundnes J, dos Santos RW, Rocha BM. Uncertainty quantification and sensitivity analysis of left ventricular function during the full cardiac cycle. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190381. [PMID: 32448074 PMCID: PMC7287338 DOI: 10.1098/rsta.2019.0381] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/10/2020] [Indexed: 05/21/2023]
Abstract
Patient-specific computer simulations can be a powerful tool in clinical applications, helping in diagnostics and the development of new treatments. However, its practical use depends on the reliability of the models. The construction of cardiac simulations involves several steps with inherent uncertainties, including model parameters, the generation of personalized geometry and fibre orientation assignment, which are semi-manual processes subject to errors. Thus, it is important to quantify how these uncertainties impact model predictions. The present work performs uncertainty quantification and sensitivity analyses to assess the variability in important quantities of interest (QoI). Clinical quantities are analysed in terms of overall variability and to identify which parameters are the major contributors. The analyses are performed for simulations of the left ventricle function during the entire cardiac cycle. Uncertainties are incorporated in several model parameters, including regional wall thickness, fibre orientation, passive material parameters, active stress and the circulatory model. The results show that the QoI are very sensitive to active stress, wall thickness and fibre direction, where ejection fraction and ventricular torsion are the most impacted outputs. Thus, to improve the precision of models of cardiac mechanics, new methods should be considered to decrease uncertainties associated with geometrical reconstruction, estimation of active stress and of fibre orientation. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- J. O. Campos
- Centro Federal de Educação Tecnológica de Minas Gerais, Leopoldina, Brazil
- Graduate Program in Computational Modeling, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
| | - J. Sundnes
- Simula Research Laboratory, PO Box 134 1325 Lysaker, Norway
| | - R. W. dos Santos
- Graduate Program in Computational Modeling, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
| | - B. M. Rocha
- Graduate Program in Computational Modeling, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
- e-mail:
| |
Collapse
|
12
|
Niederer SA, Aboelkassem Y, Cantwell CD, Corrado C, Coveney S, Cherry EM, Delhaas T, Fenton FH, Panfilov AV, Pathmanathan P, Plank G, Riabiz M, Roney CH, dos Santos RW, Wang L. Creation and application of virtual patient cohorts of heart models. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190558. [PMID: 32448064 PMCID: PMC7287335 DOI: 10.1098/rsta.2019.0558] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/06/2020] [Indexed: 05/21/2023]
Abstract
Patient-specific cardiac models are now being used to guide therapies. The increased use of patient-specific cardiac simulations in clinical care will give rise to the development of virtual cohorts of cardiac models. These cohorts will allow cardiac simulations to capture and quantify inter-patient variability. However, the development of virtual cohorts of cardiac models will require the transformation of cardiac modelling from small numbers of bespoke models to robust and rapid workflows that can create large numbers of models. In this review, we describe the state of the art in virtual cohorts of cardiac models, the process of creating virtual cohorts of cardiac models, and how to generate the individual cohort member models, followed by a discussion of the potential and future applications of virtual cohorts of cardiac models. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
| | | | | | | | | | - E. M. Cherry
- Georgia Institute of Technology, Atlanta, GA, USA
| | - T. Delhaas
- Maastricht University, Maastricht, the Netherlands
| | - F. H. Fenton
- Georgia Institute of Technology, Atlanta, GA, USA
| | - A. V. Panfilov
- Ghent University, Gent, Belgium
- Laboratory of Computational Biology and Medicine, Ural Federal University, Ekaterinburg, Russia
| | - P. Pathmanathan
- Center for Devices and Radiological Health, U.S. Food and Administration, Rockville, MD, USA
| | - G. Plank
- Medical University of Graz, Graz, Austria
| | | | | | | | - L. Wang
- Rochester Institute of Technology, La JollaRochester, NY, USA
| |
Collapse
|
13
|
Peirlinck M, Sahli Costabal F, Sack KL, Choy JS, Kassab GS, Guccione JM, De Beule M, Segers P, Kuhl E. Using machine learning to characterize heart failure across the scales. Biomech Model Mechanobiol 2019; 18:1987-2001. [PMID: 31240511 DOI: 10.1007/s10237-019-01190-w] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Accepted: 06/16/2019] [Indexed: 12/31/2022]
Abstract
Heart failure is a progressive chronic condition in which the heart undergoes detrimental changes in structure and function across multiple scales in time and space. Multiscale models of cardiac growth can provide a patient-specific window into the progression of heart failure and guide personalized treatment planning. Yet, the predictive potential of cardiac growth models remains poorly understood. Here, we quantify predictive power of a stretch-driven growth model using a chronic porcine heart failure model, subject-specific multiscale simulation, and machine learning techniques. We combine hierarchical modeling, Bayesian inference, and Gaussian process regression to quantify the uncertainty of our experimental measurements during an 8-week long study of volume overload in six pigs. We then propagate the experimental uncertainties from the organ scale through our computational growth model and quantify the agreement between experimentally measured and computationally predicted alterations on the cellular scale. Our study suggests that stretch is the major stimulus for myocyte lengthening and demonstrates that a stretch-driven growth model alone can explain [Formula: see text] of the observed changes in myocyte morphology. We anticipate that our approach will allow us to design, calibrate, and validate a new generation of multiscale cardiac growth models to explore the interplay of various subcellular-, cellular-, and organ-level contributors to heart failure. Using machine learning in heart failure research has the potential to combine information from different sources, subjects, and scales to provide a more holistic picture of the failing heart and point toward new treatment strategies.
Collapse
Affiliation(s)
- M Peirlinck
- Biofluid, Tissue and Solid Mechanics for Medical Applications (IBiTech, bioMMeda), Ghent University, Ghent, Belgium
| | - F Sahli Costabal
- Department of Mechanical Engineering, Stanford University, Stanford, CA, USA
| | - K L Sack
- Department of Human Biology, University of Cape Town, Cape Town, South Africa
- Department of Surgery, University of California at San Francisco, San Francisco, CA, USA
| | - J S Choy
- California Medical Innovations Institute, Inc., San Diego, CA, USA
| | - G S Kassab
- California Medical Innovations Institute, Inc., San Diego, CA, USA
| | - J M Guccione
- Department of Surgery, University of California at San Francisco, San Francisco, CA, USA
| | - M De Beule
- Biofluid, Tissue and Solid Mechanics for Medical Applications (IBiTech, bioMMeda), Ghent University, Ghent, Belgium
| | - P Segers
- Biofluid, Tissue and Solid Mechanics for Medical Applications (IBiTech, bioMMeda), Ghent University, Ghent, Belgium
| | - E Kuhl
- Departments of Mechanical Engineering and Bioengineering, Stanford University, Stanford, CA, USA.
| |
Collapse
|