A fast computational model for the electrophysiology of the whole human heart

Giulio Del Corso,Roberto Verzicco,Francesco Viola
DOI: https://doi.org/10.1016/j.jcp.2022.111084
2021-12-24
Abstract:In this study we present a novel computational model for unprecedented simulations of the whole cardiac electrophysiology. According to the heterogeneous electrophysiologic properties of the heart, the whole cardiac geometry is decomposed into a set of coupled conductive media having different topology and electrical conductivities: (i) a network of slender bundles comprising a fast conduction atrial network, the AV-node and the ventricular bundles; (ii) the Purkinje network; and (iii) the atrial and ventricular myocardium. The propagation of the action potential in these conductive media is governed by the bidomain/monodomain equations, which are discretized in space using an in-house finite volume method and coupled to three different cellular models, the Courtemanche model [1] for the atrial myocytes, the Stewart model [2] for the Purkinje Network and the ten Tusscher-Panfilov model [3] for the ventricular myocytes. The developed numerical model correctly reproduces the cardiac electrophysiology of the whole human heart in healthy and pathologic conditions and it can be tailored to study and optimize resynchronization therapies or invasive surgical procedures. Importantly, the whole solver is GPU-accelerated using CUDA Fortran providing an unprecedented speedup, thus opening the way for systematic parametric studies and uncertainty quantification analyses.
Medical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an efficient and accurate computational model for simulating the electrophysiology of the entire human heart. Specifically, the paper proposes a new computational framework that can simultaneously address the electrophysiology of the heart's rapid conduction network and the myocardium of the four heart chambers. This framework aims to: 1. **Accurately simulate cardiac electrophysiology**: By decomposing the entire heart geometry into a collection of coupled conduction media with different topologies and conductivities, including slender bundle networks (such as the rapid - conducting atrial network, AV node, and ventricular bundles), Purkinje network, and atrial and ventricular myocardium, to accurately reproduce the electrophysiology of the human heart under healthy and pathological conditions. 2. **Improve computational efficiency**: Utilize GPU acceleration technology (CUDA Fortran) to achieve unprecedented acceleration, thus paving the way for systematic parameter studies and uncertainty quantification analysis. 3. **Support clinical applications**: This model can be customized for researching and optimizing cardiac resynchronization therapy or invasive surgical procedures, especially for applications in dealing with pathological conditions such as arrhythmia and myocardial infarction. The core contribution of the paper lies in providing a comprehensive computational tool that can not only simulate the normal electrical activity of the heart but also conduct detailed electrophysiological analysis under pathological conditions, providing strong support for cardiac pathology research and clinical treatment.