Verification of cardiac tissue electrophysiology simulators using an
            <i>N</i>
            -version benchmark

Niederer, Steven A.; Kerfoot, Eric; Benson, Alan P.; Bernabéu, Miguel O.; Bernus, Olivier; Bradley, C P; Cherry, Elizabeth M.; Clayton, Richard H.; Fenton, Flavio H.; Garny, Alan; Heidenreich, Elvio; Land, Sander; Maleckar, Mary M.; Pathmanathan, Pras; Plank, Gernot; Rodrı́guez, José F.; Roy, Ian; Sachse, Frank B.; Seemann, Gunnar; Skavhaug, Ola; Smith, Nic

doi:10.1098/rsta.2011.0139

Cited by 253 publications

(363 citation statements)

References 53 publications

Supporting

Mentioning

311

Contrasting

Unclassified

Order By: Relevance

“…Similar behavior holds for the NI-SVI method [47], as in the case of GI. While we notice that a similar effect can be obtained when considering NURBS-based IGA, we remark that in this paper we are mainly focusing on the study of the effects of the continuity of the NURBS basis functions for a specific approach.…”

Section: Spatial Approximation: Iga For Surface Pdessupporting

confidence: 49%

Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium

Patelli

Dedè

Lassila

et al. 2017

Computer Methods in Applied Mechanics and Engineering

View full text Add to dashboard Cite

We consider Isogeometric Analysis in the framework of the Galerkin method for the spatial approximation of cardiac electrophysiology models defined on NURBS surfaces; specifically, we perform a numerical comparison between basis functions of degree p ≥ 1 and globally C k -continuous, with k = 0 or p − 1, to find the most accurate approximation of a propagating front with the minimal number of degrees of freedom.We show that B-spline basis functions of degree p ≥ 1, which are C p−1 -continuous capture accurately the front velocity of the transmembrane potential even with moderately refined meshes; similarly, we show that, for accurate tracking of curved fronts, high-order continuous B-spline basis functions should be used. Finally, we apply Isogeometric Analysis to an idealized human left atrial geometry described by NURBS with physiologically sound fiber directions and anisotropic conductivity tensor to demonstrate that the numerical scheme retains its favorable approximation properties also in a more realistic setting.

show abstract

Section: Spatial Approximation: Iga For Surface Pdessupporting

confidence: 49%

Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium

Patelli

Dedè

Lassila

et al. 2017

Computer Methods in Applied Mechanics and Engineering

View full text Add to dashboard Cite

show abstract

“…Similar to [7], we have also chosen a cuboid as the 3D geometry of the cardiac tissues. Zero-flux boundary conditions are assumed.…”

Section: Numerical Experiments and Discussionmentioning

confidence: 99%

“…[7]. At the same time, sophisticated cardiac cell models can involve up to one hundred state variables, thus requiring robust and computationally intensive solvers for the related systems of ordinary differential equations (ODEs).…”

Section: Introductionmentioning

confidence: 99%

Simulating Cardiac Electrophysiology in the Era of GPU-Cluster Computing

Chai

Wen

et al. 2013

IEICE Trans. Inf. & Syst.

View full text Add to dashboard Cite

SUMMARYThis paper presents a study of the applicability of clusters of GPUs to high-resolution 3D simulations of cardiac electrophysiology. By experimenting with representative cardiac cell models and ODE solvers, in association with solving the monodomain equation, we quantitatively analyze the obtainable computational capacity of GPU clusters. It is found that for a 501×501×101 3D mesh, which entails a 0.1mm spatial resolution, a 128-GPU cluster only needs a few minutes to carry out a 100,000-timestep cardiac excitation simulation that involves a four-variable cell model. Even higher spatial and temporal resolutions are achievable for such simplified mathematical models. On the other hand, our experiments also show that a dramatically larger cluster of GPUs is needed to handle a very detailed cardiac cell model.

show abstract

“…The same holds for the discrete interatrial connections at the coronary sinus, Bachmann's bundle and posteriorly. Excitation propagation was simulated using the monodomain solver acCELLerate [10,11]. Numerical field calculation in the torso was conducted on tetrahedral meshes using a reduced bidomain approach in order to obtain body surface potentials [12].…”

Section: Methodsmentioning

confidence: 99%

Left Atrial Hypertrophy Increases P:Wave Terminal Force Through Amplitude but not Duration

Loewe

Andlauer

Platonov

et al. 2016

2016 Computing in Cardiology Conference (CinC)

Self Cite

View full text Add to dashboard Cite

P-wave morphology correlates with the risk for atrial fibrillation (AF). Left atrial (LA) enlargement could explain both the higher risk for AF and higher P-wave terminal force (PTF) in ECG lead V 1 . However, PTF-V 1 has been shown to correlate poorly with LA size. We hypothesize that LA hypertrophy, i.e. a thickening of the myocardial wall, also contributes to increased PTF-V 1 and is part of the reason for the rather low specificity of increased PTF-V 1 regarding LA enlargement. To show this, atrial excitation propagation was simulated in a cohort of four anatomically individualized models including rule-based myocyte orientation and spatial electrophysiological heterogeneity using the monodomain approach. The LA wall was thickened symmetrically in steps of 0.66 mm by up to 3.96 mm. Interatrial conduction was possible via discrete connections at the coronary sinus, Bachmann's bundle and posteriorly. Body surface ECGs were computed using realistic, heterogeneous torso models. During the early P-wave stemming from sources in the RA, no changes were observed. Once the LA got activated, the voltage in V 1 tended to lower values for higher degrees of hypertrophy. Thus, the amplitude of the late positive Pwave decreased while the amplitude of the subsequent negative terminal phase increased. PTF-V 1 and LA wall thickening showed a correlation of 0.95. The P-wave duration was almost unaffected by LA wall thickening (∆ ≤2 ms). Our results show that PTF-V 1 is a sensitive marker for LA wall thickening and elucidate why it is superior to P-wave area. The interplay of LA hypertrophy and dilation might cause the poor empirical correlation of LA size and PTF-V 1 . IntroductionThe P-wave in the body surface ECG has long been used to gain insight into anatomy, function and dysfunction of the atria [1]. As a 12-lead ECG is routinely acquired non-invasively as part of a large number of examinations, ECG-derived measures represent ideal risk markers due to their availability and low associated costs [2]. These properties render ECG-based markers more attractive than alternatives like ultrasound, magnetic resonance imaging, electroanatomical mapping, or ECG imaging. Therefore, clinicians aim to stratify arrhythmia risk based on P-wave markers [3]. The assessment of morphological features of the P-wave is recommended in current guidelines for ECG interpretation [4] regarding the diagnosis of atrial abnormalities such as left or right atrial enlargement. The anatomy of the left atrium (LA) is of particular interest regarding the risk to develop atrial fibrillation (AF) as larger atria are more vulnerable to reentry in general. Therefore, left atrial enlargement (LAE) could explain both the higher risk to develop AF and higher P-wave terminal force (PTF) in ECG lead V 1 defined as the product of the duration and the amplitude of the terminal negative part of the P-wave. However, the criterion PTF-V 1 > 4 mVms correlates rather poorly with LA size in empirical studies. Truong et al. reported an accuracy of 51% using compute...

show abstract

Verification of cardiac tissue electrophysiology simulators using an N -version benchmark

Cited by 253 publications

References 53 publications

Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium

Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium

Simulating Cardiac Electrophysiology in the Era of GPU-Cluster Computing

Left Atrial Hypertrophy Increases P:Wave Terminal Force Through Amplitude but not Duration

Contact Info

Product

Resources

About