High-order finite element methods for cardiac monodomain simulations.
Bottom Line:
The steep gradients inherent to cardiac action potential propagation necessitate fine spatial scales and therefore a substantial computational burden.Additionally, we propose a dimensionless number, the cell Thiele modulus, as a more useful metric for determining solution convergence than element size alone.Finally, we use the cell Thiele modulus to examine convergence criteria for obtaining clinically useful activation patterns for applications such as patient-specific modeling where the total activation time is known a priori.
View Article:
PubMed Central - PubMed
Affiliation: Department of Bioengineering, University of California San Diego La Jolla, CA, USA.
ABSTRACT
Computational modeling of tissue-scale cardiac electrophysiology requires numerically converged solutions to avoid spurious artifacts. The steep gradients inherent to cardiac action potential propagation necessitate fine spatial scales and therefore a substantial computational burden. The use of high-order interpolation methods has previously been proposed for these simulations due to their theoretical convergence advantage. In this study, we compare the convergence behavior of linear Lagrange, cubic Hermite, and the newly proposed cubic Hermite-style serendipity interpolation methods for finite element simulations of the cardiac monodomain equation. The high-order methods reach converged solutions with fewer degrees of freedom and longer element edge lengths than traditional linear elements. Additionally, we propose a dimensionless number, the cell Thiele modulus, as a more useful metric for determining solution convergence than element size alone. Finally, we use the cell Thiele modulus to examine convergence criteria for obtaining clinically useful activation patterns for applications such as patient-specific modeling where the total activation time is known a priori. No MeSH data available. Related in: MedlinePlus |
Related In:
Results -
Collection
License getmorefigures.php?uid=PMC4525671&req=5
Mentions: Finally, we sought to determine convergence criteria for electrophysiology solutions when only activation patterns rather than absolute conduction times are needed. This need arises commonly in patient-specific modeling, where the total activation time is known (e.g., from a measured QRS duration or electrocardiological mapping), but the conductivity is unknown. Electrical propagation in a patient-derived human biventricular mesh was simulated at four levels of spatial refinement resulting in simulations with an average cell Thiele modulus in the primary direction of propagation of 1.0, 2.0, 4.0, and 8.0 for the four mesh refinements. For this exercise, the simulation using cubic Hermite elements with ϕc = 1.0 was considered fully converged. Activation times in the three less-converged simulations were compared to the converged activation times at node locations from the coarsest mesh (Figure 5A). Since propagation of the high cell Thiele modulus simulations was too fast, the total activation time for those simulations was scaled to give a regression line with a slope of 1.0. The Bland-Altman plot in Figure 5B comparing the scaled activation times and the fully converged activation times identifies a pattern of outliers (>2 standard deviations) that were too fast compared with the fully converged solution. These outliers were located in the basal right ventricular free wall. The root-mean-squared (RMS) error between fully converged activation pattern and the scaled unconverged activation patterns decreases as the cell Thiele modulus decreases toward 1.0 for all three interpolation methods (Figure 5C). The RMS error was less than 5 ms for all unconverged simulations and twice the RMS error was less than 5 ms for the ϕc = 2.0 and ϕc = 4.0 simulations with cubic Hermite elements, the ϕc = 2.0 and ϕc = 1.0 simulations with serendipity Hermite elements, and the ϕc = 1.0 simulation with the linear Lagrange elements. This error is comparable to the lowest uncertainly in clinically measured activation times (Gold et al., 2011; Villongco et al., 2014). Figure 5D demonstrates that the high-order methods also have smaller RMS error than the linear elements on a degree of freedom basis. |
View Article: PubMed Central - PubMed
Affiliation: Department of Bioengineering, University of California San Diego La Jolla, CA, USA.
No MeSH data available.