Composite Backward Differentiation Formula for the Bidomain Equations

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.

Reduced order method for finite difference modeling of cardiac propagation

Efficient numerical simulation of cardiac electrophysiology is crucial for studying the electrical properties of the heart tissue. The cardiac bidomain model is the most widely accepted representation of the electrical behaviour of the heart muscle. The bidomain model offers fast cardiac potential variation, which can lead to high computational cost due to the required large grid sizes and small time steps. In this paper, the complexity of the finite difference approximation of the bidomain equations is reduced with the model order reduction technique. Proper orthogonal decomposition, a projection-based algorithm, is used to efficiently approximate the original high fidelity cardiac bidomain model with a low-dimensional system of equations. The low-dimensional basis functions are computed first from the ‘snapshots,’ which contain the solutions of the full-order system for different temporal and spatial parameters. Galerkin projection of the original cardiac bidomain system onto the subspace of the reduced order basis functions reduces the size of the linear system. Numerical results confirm the efficiency of the proposed reduced order modeling technique, reducing the simulation time by a factor of 9.54, while maintaining an RMS error of 0.769 mV between the original full order solution and the reduced order POD solution.

The periodic version of the Da Prato–Grisvard theorem and applications to the bidomain equations with FitzHugh–Nagumo transport

In this article, the periodic version of the classical Da Prato–Grisvard theorem on maximal $${{L}}^p$$ L p -regularity in real interpolation spaces is developed, as well as its extension to semilinear evolution equations. Applying this technique to the bidomain equations subject to ionic transport described by the models of FitzHugh–Nagumo, Aliev–Panfilov, or Rogers–McCulloch, it is proved that this set of equations admits a unique, strongT-periodic solution in a neighborhood of stable equilibrium points provided it is innervated by T-periodic forces.

An evaluation of some assumptions underpinning the bidomain equations of electrophysiology

Tissue level cardiac electrophysiology is usually modelled by the bidomain equations, or the monodomain simplification of the bidomain equations. One assumption made when deriving the bidomain equations is that both the intracellular and extracellular spaces are in electrical equilibrium. This assumption neglects the disturbance of this equilibrium in thin regions close to the cell membrane known as Debye layers. We first demonstrate that the governing equations at the cell, or microscale, level may be adapted to take account of these Debye layers with little additional complexity, provided the permittivity within the Debye layers satisfies certain conditions that are believed to be satisfied for biological cells. We then homogenize the microscale equations using a technique developed for an almost periodic microstructure. Cardiac tissue is usually modelled as sheets of cardiac fibres stacked on top of one another. A common assumption is that an orthogonal coordinate system can be defined at each point of cardiac tissue, where the first axis is in the fibre direction, the second axis is orthogonal to the first axis but lies in the sheet of cardiac fibres and the third axis is orthogonal to the cardiac sheet. It is assumed further that both the intracellular and extracellular conductivity tensors are diagonal with respect to these axes and that the diagonal entries of these tensors are constant across the whole tissue. Using the homogenization technique we find that this assumption is usually valid for cardiac tissue, but highlight situations where the assumption may not be valid.