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 cell 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.
A fast computational model for the electrophysiology of the whole human heart
Del Corso, Giulio;Verzicco, Roberto;Viola, Francesco
2022-01-01
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 cell 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.File | Dimensione | Formato | |
---|---|---|---|
PostPrint_2022_JComputPhys_457_DelCorso.pdf
accesso aperto
Descrizione: Articolo principale
Tipologia:
Documento in Post-print
Licenza:
Creative commons
Dimensione
14.25 MB
Formato
Adobe PDF
|
14.25 MB | Adobe PDF | Visualizza/Apri |
PostPrint_2022_JComputPhys_457_DelCorso_Rid.pdf
accesso aperto
Tipologia:
Documento in Post-print
Licenza:
Creative commons
Dimensione
2.74 MB
Formato
Adobe PDF
|
2.74 MB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.