Composite Backward Differentiation Formula for the Bidomain Equations.

Journal Article (Journal Article)

The bidomain equations have been widely used to model the electrical activity of cardiac tissue. While it is well-known that implicit methods have much better stability than explicit methods, implicit methods usually require the solution of a very large nonlinear system of equations at each timestep which is computationally prohibitive. In this work, we present two fully implicit time integration methods for the bidomain equations: the backward Euler method and a second-order one-step two-stage composite backward differentiation formula (CBDF2) which is an L-stable time integration method. Using the backward Euler method as fundamental building blocks, the CBDF2 scheme is easily implementable. After solving the nonlinear system resulting from application of the above two fully implicit schemes by a nonlinear elimination method, the obtained nonlinear global system has a much smaller size, whose Jacobian is symmetric and possibly positive definite. Thus, the residual equation of the approximate Newton approach for the global system can be efficiently solved by standard optimal solvers. As an alternative, we point out that the above two implicit methods combined with operator splittings can also efficiently solve the bidomain equations. Numerical results show that the CBDF2 scheme is an efficient time integration method while achieving high stability and accuracy.

Full Text

Duke Authors

Cited Authors

  • Gao, X; Henriquez, CS; Ying, W

Published Date

  • January 2020

Published In

Volume / Issue

  • 11 /

Start / End Page

  • 591159 -

PubMed ID

  • 33381051

Pubmed Central ID

  • PMC7767930

Electronic International Standard Serial Number (EISSN)

  • 1664-042X

International Standard Serial Number (ISSN)

  • 1664-042X

Digital Object Identifier (DOI)

  • 10.3389/fphys.2020.591159


  • eng