Modeling and simulation of the electric activity of the heart using graphic processing units

Mathematical modelling and simulation of the electric activity of the heart (cardiac electrophysiology) offers and ideal framework to combine clinical and experimental data in order to help understanding the underlying mechanisms behind the observed respond under physiological and pathological conditions. In this regard, solving the electric activity of the heart possess a big challenge, not only because of the structural complexities inherent to the heart tissue, but also because of the complex electric behaviour of the cardiac cells. The multi- scale nature of the electrophysiology problem makes difficult its numerical solution, requiring temporal and spatial resolutions of 0.1 ms and 0.2 mm respectively for accurate simulations, leading to models with millions degrees of freedom that need to be solved for thousand time steps. Solution of this problem requires the use of algorithms with higher level of parallelism in multi-core platforms. In this regard the newer programmable graphic processing units (GPU) has become a valid alternative due to their tremendous computational horsepower. This thesis develops around the implementation of an electrophysiology simulation software entirely developed in Compute Unified Device Architecture (CUDA) for GPU computing. The software implements fully explicit and semi-implicit solvers for the monodomain model, using operator splitting and the finite element method for space discretization. Performance is compared with classical multi-core MPI based solvers operating on dedicated high-performance computer clusters. Results obtained with the GPU based solver show enormous potential for this technology with accelerations over 50× for three-dimensional problems when using an implicit scheme for the parabolic equation, whereas accelerations reach values up to 100× for the explicit implementation. The implemented solver has been applied to study pro-arrhythmic mechanisms during acute ischemia. In particular, we investigate on how hyperkalemia affects the vulnerability window to reentry and the reentry patterns in the heterogeneous substrate caused by acute regional ischemia using an anatomically and biophysically detailed human biventricular model. A three dimensional geometrically and anatomically accurate regionally ischemic human heart model was created. The ischemic region was located in the inferolateral and posterior side of the left ventricle mimicking the occlusion of the circumflex artery, and the presence of a washed-out zone not affected by ischemia at the endocardium has been incorporated. Realistic heterogeneity and fi er anisotropy has also been considered in the model. A highly electrophysiological detailed action potential model for human has been adapted to make it suitable for modeling ischemic conditions (hyperkalemia, hipoxia, and acidic conditions) by introducing a formulation of the ATP-sensitive K+ current. The model predicts the generation of sustained re-entrant activity in the form single and double circus around a blocked area within the ischemic zone for K+ concentrations bellow 9mM, with the reentrant activity associated with ventricular tachycardia in all cases. Results suggest the washed-out zone as a potential pro-arrhythmic substrate factor helping on establishing sustained ventricular tachycardia.Colli-Franzone P, Pavarino L. A parallel solver for reaction-diffusion systems in computational electrocardiology, Math. Models Methods Appl. Sci. 14 (06):883-911, 2004.Colli-Franzone P, Deu hard P, Erdmann B, Lang J, Pavarino L F. Adaptivity in space and time for reaction-diffusion systems in electrocardiology, SIAM J. Sci. Comput. 28 (3):942-962, 2006.Ferrero J M(Jr), Saiz J, Ferrero J M, Thakor N V. Simulation of action potentials from metabolically impaired cardiac myocytes: Role of atp-sensitive K+ current. Circ Res, 79(2):208-221, 1996.Ferrero J M (Jr), Trenor B. Rodriguez B, Saiz J. Electrical acticvity and reentry during acute regional myocardial ischemia: Insights from simulations.Int J Bif Chaos, 13:3703-3715, 2003.Heidenreich E, Ferrero J M, Doblare M, Rodriguez J F. Adaptive macro finite elements for the numerical solution of monodomain equations in cardiac electrophysiology, Ann. Biomed. Eng. 38 (7):2331-2345, 2010.Janse M J, Kleber A G. Electrophysiological changes and ventricular arrhythmias in the early phase of regional myocardial ischemia. Circ. Res. 49:1069-1081, 1981.ten Tusscher K HWJ, Panlov A V. Alternans and spiral breakup in a human ventricular tissue model. Am. J.Physiol. Heart Circ. Physiol. 291(3):1088-1100, 2006.<br /

lifex-ep: a robust and efficient software for cardiac electrophysiology simulations

Background: Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both normal and pathological conditions, are still in the process of achieving full maturity within the scientific community. Results: This work introduces lifex-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both normal and pathological conditions. lifex-ep employs the monodomain equation to model the heart's electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, lifex-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within lifex-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying lifex-ep, along with comprehensive implementation details and instructions for users. lifex-ep features exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We showcase the key features of lifex-ep through various idealized and realistic simulations conducted in both normal and pathological scenarios. Furthermore, the software offers a user-friendly and flexible interface, simplifying the setup of simulations using self-documenting parameter files. Conclusions: lifex-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. lifex-ep represents a valuable tool for conducting in silico patient-specific simulations

Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units

In this paper we studied the implementation and performance of adaptive step methods for large systems of ordinary differential equations systems in graphics processing units, focusing on the simulation of three-dimensional electric cardiac activity. The Rush-Larsen method was applied in all the implemented solvers to improve efficiency. We compared the adaptive methods with the fixed step methods, and we found that the fixed step methods can be faster while the adaptive step methods are better in terms of accuracy and robustness. (c) 2013 Elsevier Ltd. All rights reserved.This work has been partially funded by Universitat Politecnica de Valencia through Programa de Apoyo a la InvestigaciOn y Desarrollo (PAID-06-11) and (PAID-05-12), by Generalitat Valenciana through projects PROMETEO/2009/013 and Ayudas para la realizacion de proyectos de I+D para grupos de investigacion emergentes GV/2012/039, and by Ministerio Espafiol de Economia y Competitividad through project TEC2012-38142-004.García Mollá, VM.; Liberos Mascarell, A.; Vidal Maciá, AM.; Guillem Sánchez, MS.; Millet Roig, J.; González Salvador, A.; Martínez Zaldívar, FJ.... (2014). Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units. Computers in Biology and Medicine. 44:15-26. https://doi.org/10.1016/j.compbiomed.2013.10.023S15264

Finite element simulations: computations and applications to aerodynamics and biomedicine.

171 p.Las ecuaciones en derivadas parciales describen muchos fenómenos de interés práctico y sus solucionessuelen necesitar correr simulaciones muy costosas en clústers de cálculo.En el ámbito de los flujos turbulentos, en particular, el coste de las simulaciones es demasiado grande sise utilizan métodos básicos, por eso es necesario modelizar el sistema.Esta tesis doctoral trata principalmente de dos temas en Cálculo Científico.Por un lado, estudiamos nuevos desarrollos en la modelización y simulación de flujos turbulentos;utilizamos un Método de Elementos Finitos adaptativo y un modelo de ¿número de Reynolds infinito¿para reducir el coste computacional de simulaciones que, sin estas modificaciones, serían demasiadocostosas.De esta manera conseguimos lograr simulaciones evolutivas de flujos turbulentos con número deReynolds muy grande, lo cual se considera uno de los mayores retos en aerodinámica.El otro pilar de esta tesis es una aplicación biomédica.Desarrollamos un modelo computacional de Ablación (Cardiaca) por Radiofrecuencia, una terapiacomún para tratar varias enfermedades, por ejemplo algunas arritmias.Nuestro modelo mejora los modelos existentes en varias maneras, y en particular en tratar de obteneruna aproximación fiel de la geometría del sistema, lo cual se descubre ser crítico para simularcorrectamente la física del fenómeno

Exploration of the Human Purkinje Network in Virtual Populations

This thesis investigates the Purkinje network (PN) and its dependency on the heart shape (HS) through cardiac simulation on virtual populations (VPs). The heart is a complex organ and essential to the wellbeing of humans; its dysfunction is responsible for more than 27% of all deaths in the UK. The PN delivers the activation impulse to the ventricles of the heart and ensures their synchronous activation. Thus, the morphology of the PN is important, but it varies between species and in vivo imaging is not feasible. However, computer simulation could provide an alternative experimental tool. In simulation of the cardiac electrophysiology, the PN is often replaced by stimulus points on the HS that are fitted to physiological measurements (heart activation times, ECG). Thus, not allowing the study of the PN morphology, nor studies of arrhythmia involving re-entry into the PN. In this thesis, three studies involving explicit models of PNs have been conducted. First, an efficient algorithm for solving electrophysiology models for the PN is introduced. These allow performing simulations of physiological activations. To minimise the time for simulations, parallelisation with CPU and GPU architectures are investigated, which is of interest for VP studies. In the second study, false tendons (FTs) are studied, which provide an additional connection from the left bundle branch (LBB) and are potentially beneficial in case of LBB block. Therefore, the reduction in activation times by FT is studied as a function of the HS. In the third study, an automatically generated VP is used to explore uncertainty in the PN morphology. The conjecture is that the PN structure adapts to the HS. The coverage of the septum and the minimum distance of the PN to the base are varied. The features of the resulting ECG are used to find the PN that gives maximally synchronised contraction