High-order spectral/hp element discretisation for reaction–diffusion problems on surfaces: Application to cardiac electrophysiology
Abstract
We present a numerical discretisation of an embedded two-dimensional manifold using high-order continuous Galerkin spectral/hp elements, which provide exponential convergence of the solution with increasing polynomial order, while retaining geometric flexibility in the representation of the domain. Our work is motivated by applications in cardiac electrophysiology where sharp gradients in the solution benefit from the high-order discretisation, while the computational cost of anatomically-realistic models can be significantly reduced through the surface representation and use of high-order methods. We describe and validate our discretisation and provide a demonstration of its application to modelling electrochemical propagation across a human left atrium.
- Authors:
- Publication Date:
- Research Org.:
- Univ. of Utah, Salt Lake City, UT (United States)
- Sponsoring Org.:
- USDOE Office of Energy Efficiency and Renewable Energy (EERE); British Heart Foundation
- OSTI Identifier:
- 1209733
- Alternate Identifier(s):
- OSTI ID: 1623691
- Grant/Contract Number:
- DOE NETL DE-EE0004449; EE0004449
- Resource Type:
- Published Article
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Name: Journal of Computational Physics Journal Volume: 257 Journal Issue: PA; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICS AND COMPUTING; Computer Science; Physics; High-order finite elements; Spectral/hp elements; Continuous Galerkin method; Surface PDE; Cardiac electrophysiology; Monodomain equation
Citation Formats
Cantwell, Chris D., Yakovlev, Sergey, Kirby, Robert M., Peters, Nicholas S., and Sherwin, Spencer J. High-order spectral/hp element discretisation for reaction–diffusion problems on surfaces: Application to cardiac electrophysiology. United States: N. p., 2014.
Web. doi:10.1016/j.jcp.2013.10.019.
Cantwell, Chris D., Yakovlev, Sergey, Kirby, Robert M., Peters, Nicholas S., & Sherwin, Spencer J. High-order spectral/hp element discretisation for reaction–diffusion problems on surfaces: Application to cardiac electrophysiology. United States. https://doi.org/10.1016/j.jcp.2013.10.019
Cantwell, Chris D., Yakovlev, Sergey, Kirby, Robert M., Peters, Nicholas S., and Sherwin, Spencer J. Wed .
"High-order spectral/hp element discretisation for reaction–diffusion problems on surfaces: Application to cardiac electrophysiology". United States. https://doi.org/10.1016/j.jcp.2013.10.019.
@article{osti_1209733,
title = {High-order spectral/hp element discretisation for reaction–diffusion problems on surfaces: Application to cardiac electrophysiology},
author = {Cantwell, Chris D. and Yakovlev, Sergey and Kirby, Robert M. and Peters, Nicholas S. and Sherwin, Spencer J.},
abstractNote = {We present a numerical discretisation of an embedded two-dimensional manifold using high-order continuous Galerkin spectral/hp elements, which provide exponential convergence of the solution with increasing polynomial order, while retaining geometric flexibility in the representation of the domain. Our work is motivated by applications in cardiac electrophysiology where sharp gradients in the solution benefit from the high-order discretisation, while the computational cost of anatomically-realistic models can be significantly reduced through the surface representation and use of high-order methods. We describe and validate our discretisation and provide a demonstration of its application to modelling electrochemical propagation across a human left atrium.},
doi = {10.1016/j.jcp.2013.10.019},
journal = {Journal of Computational Physics},
number = PA,
volume = 257,
place = {United States},
year = {Wed Jan 01 00:00:00 EST 2014},
month = {Wed Jan 01 00:00:00 EST 2014}
}
https://doi.org/10.1016/j.jcp.2013.10.019
Web of Science
Works referenced in this record:
Flows on surfaces of arbitrary topology
journal, July 2003
- Stam, Jos
- ACM Transactions on Graphics, Vol. 22, Issue 3
Simulations of (An)Isotropic Diffusion on Curved Biological Surfaces
journal, February 2006
- Sbalzarini, Ivo F.; Hayer, Arnold; Helenius, Ari
- Biophysical Journal, Vol. 90, Issue 3
Finite element approximation of elliptic partial differential equations on implicit surfaces
journal, January 2008
- Burger, Martin
- Computing and Visualization in Science, Vol. 12, Issue 3
A Finite Volume Method for Solving Parabolic Equations on Logically Cartesian Curved Surface Meshes
journal, January 2010
- Calhoun, Donna A.; Helzel, Christiane
- SIAM Journal on Scientific Computing, Vol. 31, Issue 6
Level Set Equations on Surfaces via the Closest Point Method
journal, March 2008
- Macdonald, Colin B.; Ruuth, Steven J.
- Journal of Scientific Computing, Vol. 35, Issue 2-3
The Implicit Closest Point Method for the Numerical Solution of Partial Differential Equations on Surfaces
journal, January 2010
- Macdonald, Colin B.; Ruuth, Steven J.
- SIAM Journal on Scientific Computing, Vol. 31, Issue 6
Fourth order partial differential equations on general geometries
journal, July 2006
- Greer, John B.; Bertozzi, Andrea L.; Sapiro, Guillermo
- Journal of Computational Physics, Vol. 216, Issue 1
A simple embedding method for solving partial differential equations on surfaces
journal, January 2008
- Ruuth, Steven J.; Merriman, Barry
- Journal of Computational Physics, Vol. 227, Issue 3
Eulerian finite element method for parabolic PDEs on implicit surfaces
journal, January 2008
- Dziuk, Gerhard; Elliott, Charles
- Interfaces and Free Boundaries
An h-narrow band finite-element method for elliptic equations on implicit surfaces
journal, March 2009
- Deckelnick, K.; Dziuk, G.; Elliott, C. M.
- IMA Journal of Numerical Analysis, Vol. 30, Issue 2
An Improvement of a Recent Eulerian Method for Solving PDEs on General Geometries
journal, December 2005
- Greer, John B.
- Journal of Scientific Computing, Vol. 29, Issue 3
Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model
journal, July 1998
- Courtemanche, Marc; Ramirez, Rafael J.; Nattel, Stanley
- American Journal of Physiology-Heart and Circulatory Physiology, Vol. 275, Issue 1
The Spectral Element Method for the Shallow Water Equations on the Sphere
journal, January 1997
- Taylor, Mark; Tribbia, Joseph; Iskandarani, Mohamed
- Journal of Computational Physics, Vol. 130, Issue 1
Anisotropic diffusion in vector field visualization on Euclidean domains and surfaces
journal, January 2000
- Diewald, U.; Preusser, T.; Rumpf, M.
- IEEE Transactions on Visualization and Computer Graphics, Vol. 6, Issue 2
Efficient simulation of cardiac electrical propagation using high order finite elements
journal, May 2012
- Arthurs, Christopher J.; Bishop, Martin J.; Kay, David
- Journal of Computational Physics, Vol. 231, Issue 10
TetrahedralhpFinite Elements: Algorithms and Flow Simulations
journal, March 1996
- Sherwin, S. J.; Karniadakis, G. E.
- Journal of Computational Physics, Vol. 124, Issue 1
From h to p efficiently: Implementing finite and spectral/hp element methods to achieve optimal performance for low- and high-order discretisations
journal, July 2010
- Vos, Peter E. J.; Sherwin, Spencer J.; Kirby, Robert M.
- Journal of Computational Physics, Vol. 229, Issue 13
Petascale algorithms for reactor hydrodynamics
journal, July 2008
- Fischer, P.; Lottes, J.; Pointer, D.
- Journal of Physics: Conference Series, Vol. 125
Implicit-Explicit Methods for Time-Dependent Partial Differential Equations
journal, June 1995
- Ascher, Uri M.; Ruuth, Steven J.; Wetton, Brian T. R.
- SIAM Journal on Numerical Analysis, Vol. 32, Issue 3
A generic framework for time-stepping partial differential equations (PDEs): general linear methods, object-oriented implementation and application to fluid problems
journal, March 2011
- Vos, Peter E. J.; Eskilsson, Claes; Bolis, Alessandro
- International Journal of Computational Fluid Dynamics, Vol. 25, Issue 3