High-order finite element methods for cardiac monodomain simulations

Efficient time splitting schemes for the monodomain equation in cardiac electrophysiology

Approximating the fast dynamics of depolarization waves in the human heart described by the monodomain model is numerically challenging. Splitting methods for the PDE-ODE coupling enable the computation with very fine space and time discretizations. Here, we compare different splitting approaches regarding convergence, accuracy, and efficiency. Simulations were performed for a benchmark problem with the Beeler-Reuter cell model on a truncated ellipsoid approximating the left ventricle including a localized stimulation. For this configuration, we provide a reference solution for the transmembrane potential. We found a semi-implicit approach with state variable interpolation to be the most efficient scheme. The results are transferred to a more physiological setup using a bi-ventricular domain with a complex external stimulation pattern to evaluate the accuracy of the activation time for different resolutions in space and time.

Forward-Solution Noninvasive Computational Arrhythmia Mapping: The VMAP Study

BACKGROUND

The accuracy of noninvasive arrhythmia source localization using a forward-solution computational mapping system has not yet been evaluated in blinded, multicenter analysis. This study tested the hypothesis that a computational mapping system incorporating a comprehensive arrhythmia simulation library would provide accurate localization of the site-of-origin for atrial and ventricular arrhythmias and pacing using 12-lead ECG data when compared with the gold standard of invasive electrophysiology study and ablation.

METHODS

The VMAP study (Vectorcardiographic Mapping of Arrhythmogenic Probability) was a blinded, multicenter evaluation with final data analysis performed by an independent core laboratory. Eligible episodes included atrial and ventricular: tachycardia, fibrillation, pacing, premature atrial and ventricular complexes, and orthodromic atrioventricular reentrant tachycardia. Mapping system results were compared with the gold standard site of successful ablation or pacing during electrophysiology study and ablation. Mapping time was assessed from time-stamped logs. Prespecified performance goals were used for statistical comparisons.

RESULTS

A total of 255 episodes from 225 patients were enrolled from 4 centers. Regional accuracy for ventricular tachycardia and premature ventricular complexes in patients without significant structural heart disease (n=75, primary end point) was 98.7% (95% CI, 96.0%-100%; P<0.001 to reject predefined H0 <0.80). Regional accuracy for all episodes (secondary end point 1) was 96.9% (95% CI, 94.7%-99.0%; P<0.001 to reject predefined H0 <0.75). Accuracy for the exact or neighboring segment for all episodes (secondary end point 2) was 97.3% (95% CI, 95.2%-99.3%; P<0.001 to reject predefined H0 <0.70). Median spatial accuracy was 15 mm (n=255, interquartile range, 7-25 mm). The mapping process was completed in a median of 0.8 minutes (interquartile range, 0.4-1.4 minutes).

CONCLUSIONS

Computational ECG mapping using a forward-solution approach exceeded prespecified accuracy goals for arrhythmia and pacing localization. Spatial accuracy analysis demonstrated clinically actionable results. This rapid, noninvasive mapping technology may facilitate catheter-based and noninvasive targeted arrhythmia therapies.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT04559061.

Balancing conduction velocity error in cardiac electrophysiology using a modified quadrature approach

Conduction velocity error is often the main culprit behind the need for very fine spatial discretizations and high computational effort in cardiac electrophysiology problems. In light of this, a novel approach for simulating an accurate conduction velocity in coarse meshes with linear elements is suggested based on a modified quadrature approach. In this approach, the quadrature points are placed at arbitrary offsets of the isoparametric coordinates. A numerical study illustrates the dependence of the conduction velocity on the spatial discretization and the conductivity when using different quadrature rules and calculation approaches. Additionally, examples using the modified quadrature in coarse meshes for wave propagation demonstrate the improved accuracy of the conduction velocity with this method. This novel approach possesses great potential in reducing the computational effort required but remains limited to specific linear elements and experiences a reduction in accuracy for irregular meshes and heterogeneous conductivities. Further research can focus on developing an adaptive quadrature and extending the approach to other element formulations in order to make the approach more generally applicable.

A computationally efficient dynamic model of human epicardial tissue

We present a new phenomenological model of human ventricular epicardial cells and we test its reentry dynamics. The model is derived from the Rogers-McCulloch formulation of the FitzHugh-Nagumo equations and represents the total ionic current divided into three contributions corresponding to the excitatory, recovery and transient outward currents. Our model reproduces the main characteristics of human epicardial tissue, including action potential amplitude and morphology, upstroke velocity, and action potential duration and conduction velocity restitution curves. The reentry dynamics is stable, and the dominant period is about 270 ms, which is comparable to clinical values. The proposed model is the first phenomenological model able to accurately resemble human experimental data by using only 3 state variables and 17 parameters. Indeed, it is more computationally efficient than existing models (i.e., almost two times faster than the minimal ventricular model). Beyond the computational efficiency, the low number of parameters facilitates the process of fitting the model to the experimental data.

Atlas-based methods for efficient characterization of patient-specific ventricular activation patterns

AIMS

Ventricular activation patterns can aid clinical decision-making directly by providing spatial information on cardiac electrical activation or indirectly through derived clinical indices. The aim of this work was to derive an atlas of the major modes of variation of ventricular activation from model-predicted 3D bi-ventricular activation time distributions and to relate these modes to corresponding vectorcardiograms (VCGs). We investigated how the resulting dimensionality reduction can improve and accelerate the estimation of activation patterns from surface electrogram measurements.

METHODS AND RESULTS

Atlases of activation time (AT) and VCGs were derived using principal component analysis on a dataset of simulated electrophysiology simulations computed on eight patient-specific bi-ventricular geometries. The atlases provided significant dimensionality reduction, and the modes of variation in the two atlases described similar features. Utility of the atlases was assessed by resolving clinical waveforms against them and the VCG atlas was able to accurately reconstruct the patient VCGs with fewer than 10 modes. A sensitivity analysis between the two atlases was performed by calculating a compact Jacobian. Finally, VCGs generated by varying AT atlas modes were compared with clinical VCGs to estimate patient-specific activation maps, and the resulting errors between the clinical and atlas-based VCGs were less than those from more computationally expensive method.

CONCLUSION

Atlases of activation and VCGs represent a new method of identifying and relating the features of these high-dimensional signals that capture the major sources of variation between patients and may aid in identifying novel clinical indices of arrhythmia risk or therapeutic outcome.

A numerical study on the effects of spatial and temporal discretization in cardiac electrophysiology

Millions of degrees of freedom are often required to accurately represent the electrophysiology of the myocardium due to the presence of discretization effects. This study seeks to explore the influence of temporal and spatial discretization on the simulation of cardiac electrophysiology in conjunction with changes in modeling choices. Several finite element analyses are performed to examine how discretization affects solution time, conduction velocity and electrical excitation. Discretization effects are considered along with changes in the electrophysiology model and solution approach. Two action potential models are considered: the Aliev-Panfilov model and the ten Tusscher-Noble-Noble-Panfilov model. The solution approaches consist of two time integration schemes and different treatments for solving the local system of ordinary differential equations. The efficiency and stability of the calculation approaches are demonstrated to be dependent on the action potential model. The dependency of the conduction velocity on the element size and time step is shown to be different for changes in material parameters. Finally, the discrepancies between the wave propagation in coarse and fine meshes are analyzed based on the temporal evolution of the transmembrane potential at a node and its neighboring Gauss points. Insight obtained from this study can be used to suggest new methods to improve the efficiency of simulations in cardiac electrophysiology.

Accelerating simulations of cardiac electrical dynamics through a multi-GPU platform and an optimized data structure

Simulations of cardiac electrophysiological models in tissue, particularly in 3D require the solutions of billions of differential equations even for just a couple of milliseconds, thus highly demanding in computational resources. In fact, even studies in small domains with very complex models may take several hours to reproduce seconds of electrical cardiac behavior. Today's Graphics Processor Units (GPUs) are becoming a way to accelerate such simulations, and give the added possibilities to run them locally without the need for supercomputers. Nevertheless, when using GPUs, bottlenecks related to global memory access caused by the spatial discretization of the large tissue domains being simulated, become a big challenge. For simulations in a single GPU, we propose a strategy to accelerate the computation of the diffusion term through a data-structure and memory access pattern designed to maximize coalescent memory transactions and minimize branch divergence, achieving results approximately 1.4 times faster than a standard GPU method. We also combine this data structure with a designed communication strategy to take advantage in the case of simulations in multi-GPU platforms. We demonstrate that, in the multi-GPU approach performs, simulations in 3D tissue can be just 4× slower than real time.

Perspectives on Sharing Models and Related Resources in Computational Biomechanics Research

The role of computational modeling for biomechanics research and related clinical care will be increasingly prominent. The biomechanics community has been developing computational models routinely for exploration of the mechanics and mechanobiology of diverse biological structures. As a result, a large array of models, data, and discipline-specific simulation software has emerged to support endeavors in computational biomechanics. Sharing computational models and related data and simulation software has first become a utilitarian interest, and now, it is a necessity. Exchange of models, in support of knowledge exchange provided by scholarly publishing, has important implications. Specifically, model sharing can facilitate assessment of reproducibility in computational biomechanics and can provide an opportunity for repurposing and reuse, and a venue for medical training. The community's desire to investigate biological and biomechanical phenomena crossing multiple systems, scales, and physical domains, also motivates sharing of modeling resources as blending of models developed by domain experts will be a required step for comprehensive simulation studies as well as the enhancement of their rigor and reproducibility. The goal of this paper is to understand current perspectives in the biomechanics community for the sharing of computational models and related resources. Opinions on opportunities, challenges, and pathways to model sharing, particularly as part of the scholarly publishing workflow, were sought. A group of journal editors and a handful of investigators active in computational biomechanics were approached to collect short opinion pieces as a part of a larger effort of the IEEE EMBS Computational Biology and the Physiome Technical Committee to address model reproducibility through publications. A synthesis of these opinion pieces indicates that the community recognizes the necessity and usefulness of model sharing. There is a strong will to facilitate model sharing, and there are corresponding initiatives by the scientific journals. Outside the publishing enterprise, infrastructure to facilitate model sharing in biomechanics exists, and simulation software developers are interested in accommodating the community's needs for sharing of modeling resources. Encouragement for the use of standardized markups, concerns related to quality assurance, acknowledgement of increased burden, and importance of stewardship of resources are noted. In the short-term, it is advisable that the community builds upon recent strategies and experiments with new pathways for continued demonstration of model sharing, its promotion, and its utility. Nonetheless, the need for a long-term strategy to unify approaches in sharing computational models and related resources is acknowledged. Development of a sustainable platform supported by a culture of open model sharing will likely evolve through continued and inclusive discussions bringing all stakeholders at the table, e.g., by possibly establishing a consortium.

A framework for biomechanics simulations using four-chamber cardiac models

Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly-fewer elements are needed to build the cardiac geometry-compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.

Tissue-Specific Optical Mapping Models of Swine Atria Informed by Optical Coherence Tomography

Computational models and experimental optical mapping of cardiac electrophysiology serve as powerful tools to investigate the underlying mechanisms of arrhythmias. Modeling can also aid the interpretation of optical mapping signals, which may have different characteristics with respect to the underlying electrophysiological signals they represent. However, despite the prevalence of atrial arrhythmias such as atrial fibrillation, models of optical electrical mapping incorporating realistic structure of the atria are lacking. Therefore, we developed image-based models of atrial tissue using structural information extracted from optical coherence tomography (OCT), which can provide volumetric tissue characteristics in high resolution. OCT volumetric data of four swine atrial tissue samples were used to develop models incorporating tissue geometry, tissue-specific myofiber orientation, and ablation lesion regions. We demonstrated the use of these models through electrophysiology and photon scattering simulations. Changes in transmural electrical conduction were observed with the inclusion of OCT-derived, depth-resolved fiber orientation. Additionally, the amplitude of optical mapping signals were not found to correspond with lesion transmurality because of lesion geometry and electrical propagation occurring beyond excitation light penetration. This work established a framework for the development of tissue-specific models of atrial tissue derived from OCT imaging data, which can be useful in future investigations of electrophysiology and optical mapping signals with respect to realistic atrial tissue structure.

Space-discretization error analysis and stabilization schemes for conduction velocity in cardiac electrophysiology

In cardiac electrophysiology, the propagation of the action potential may be described by a set of reaction-diffusion equations known as the bidomain model. The shape of the solution is determined by a balance of a strong reaction and a relatively weak diffusion, which leads to steep variations in space and time. From a numerical point of view, the sharp spatial gradients may be seen as particularly problematic, because computational grid resolution on the order of 0.1 mm or less is required, yielding considerable computational efforts on human geometries. In this paper, we discuss a number of well-known numerical schemes for the bidomain equation and show how the quality of the solution is affected by the spatial discretization. In particular, we study in detail the effect of discretization on the conduction velocity (CV), which is an important quantity from a physiological point of view. We show that commonly applied finite element techniques tend to overestimate the CV on coarse grids, while it tends to be underestimated by finite difference schemes. Furthermore, the choice of interpolation and discretization scheme for the nonlinear reaction term has a strong impact on the CV. Finally, we exploit the results of the error analysis to propose improved numerical methods, including a stabilized scheme that tends to correct the CV on coarse grids but converges to the correct solution as the grid is refined. Copyright © 2016 John Wiley & Sons, Ltd.

Biomechanics Simulations Using Cubic Hermite Meshes with Extraordinary Nodes for Isogeometric Cardiac Modeling

Cubic Hermite hexahedral finite element meshes have some well-known advantages over linear tetrahedral finite element meshes in biomechanical and anatomic modeling using isogeometric analysis. These include faster convergence rates as well as the ability to easily model rule-based anatomic features such as cardiac fiber directions. However, it is not possible to create closed complex objects with only regular nodes; these objects require the presence of extraordinary nodes (nodes with 3 or >= 5 adjacent elements in 2D) in the mesh. The presence of extraordinary nodes requires new constraints on the derivatives of adjacent elements to maintain continuity. We have developed a new method that uses an ensemble coordinate frame at the nodes and a local-to-global mapping to maintain continuity. In this paper, we make use of this mapping to create cubic Hermite models of the human ventricles and a four-chamber heart. We also extend the methods to the finite element equations to perform biomechanics simulations using these meshes. The new methods are validated using simple test models and applied to anatomically accurate ventricular meshes with valve annuli to simulate complete cardiac cycle simulations.