Abstract

The numerical simulation of multiple scenarios easily becomes computationally prohibitive for cardiac electrophysiology (EP) problems if relying on usual high-fidelity, full order models (FOMs). Likewise, the use of traditional reduced order models (ROMs) for parametrized PDEs to speed up the solution of the aforementioned problems can be problematic. This is primarily due to the strong variability characterizing the solution set and to the nonlinear nature of the input-output maps that we intend to reconstruct numerically. To enhance ROM efficiency, we proposed a new generation of non-intrusive, nonlinear ROMs, based on deep learning (DL) algorithms, such as convolutional, feedforward, and autoencoder neural networks. In the proposed DL-ROM, both the nonlinear solution manifold and the nonlinear reduced dynamics used to model the system evolution on that manifold can be learnt in a non-intrusive way thanks to DL algorithms trained on a set of FOM snapshots. DL-ROMs were shown to be able to accurately capture complex front propagation processes, both in physiological and pathological cardiac EP, very rapidly once neural networks were trained, however, at the expense of huge training costs. In this study, we show that performing a prior dimensionality reduction on FOM snapshots through randomized proper orthogonal decomposition (POD) enables to speed up training times and to decrease networks complexity. Accuracy and efficiency of this strategy, which we refer to as POD-DL-ROM, are assessed in the context of cardiac EP on an idealized left atrium (LA) geometry and considering snapshots arising from a NURBS (non-uniform rational B-splines)-based isogeometric analysis (IGA) discretization. Once the ROMs have been trained, POD-DL-ROMs can efficiently solve both physiological and pathological cardiac EP problems, for any new scenario, in real-time, even in extremely challenging contexts such as those featuring circuit re-entries, that are among the factors triggering cardiac arrhythmias.

Highlights

  • Computational cardiac electrophysiology (EP) is built upon mathematical and numerical models that aim at simulating both physiological and pathological heart rhythm, such as, e.g., ventricular tachycardia and atrial fibrillation

  • We show that proper orthogonal decomposition (POD)-deep learning (DL)-reduced order models (ROMs) can handle parametrized problems in cardiac EP effectively and provide fast and accurate solutions to EP problems set on realistic geometries

  • The cardiac EP problems addressed in this paper fit into both (i) a multi-query context, since repetitive evaluations of the input-output map are required in order to perform multiscenario analysis, in order to deal with inter- and intra-subject variability and to consider specific pathological scenarios, and a (ii) real-time context, due to the need, in a clinical setting, to compute outputs of interest in a very limited amount of time

Read more

Summary

Introduction

Computational cardiac electrophysiology (EP) is built upon mathematical and numerical models that aim at simulating both physiological and pathological heart rhythm, such as, e.g., ventricular tachycardia and atrial fibrillation (see, e.g., Vigmond et al, 2002, 2008; Niederer et al, 2009, 2011; Trayanova, 2011; Prakosa et al, 2018; Strocchi et al, 2020). Simulating the electrical behavior of the heart, from the cellular scale to the tissue level, requires the numerical approximation of coupled nonlinear dynamical systems, such as, e.g. the Bidomain equations (see, e.g., Colli Franzone et al, 2005, 2006), coupled with suitable ionic models, such as the FitzHugh-Nagumo (FitzHugh, 1961; Nagumo et al, 1962), the Aliev-Panfilov (Aliev and Panfilov, 1996; Nash and Panfilov, 2004), the Roger-McCulloch (Rogers and McCulloch, 1994), the ten Tusscher-Panfilov (ten Tusscher and Panfilov, 2006), or the Mitchell and Schaeffer models (Mitchell and Schaeffer, 2003) Multiple solutions of these systems, corresponding to different model inputs parameters and data, such as, e.g., electrical conductivities, ionic model parameters, and applied currents, need to be computed to evaluate outputs of clinical interest, such as activation maps (ACs) and action potential (AP) duration. They yield more significant computational savings than low-fidelity models (such as, e.g., FOMs built on coarser meshes) by replacing the FOM by a reduced order model (ROM), featuring a much lower dimension, yet capable to express the physical features of the problem at hand

Objectives
Methods
Results
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call