A Spectral Deferred Correction Method for Solving Cardiac Models

dc.contributor.advisor

Schaeffer, David G

dc.contributor.author

Bowen, Matthew M.

dc.date.accessioned

2011-05-20T19:36:04Z

dc.date.available

2011-05-20T19:36:04Z

dc.date.issued

2011

dc.department

Mathematics

dc.description.abstract

Many numerical approaches exist to solving models of electrical activity in the heart. These models consist of a system of stiff nonlinear ordinary differential equations for the voltage and other variables governing channels, with the voltage coupled to a diffusion term. In this work, we propose a new algorithm that uses two common discretization methods, operator splitting and finite elements. Additionally, we incorporate a temporal integration process known as spectral deferred correction. Using these approaches,

we construct a numerical method that can achieve arbitrarily high order in both space and time in order to resolve important features of the models, while gaining accuracy and efficiency over lower order schemes.

Our algorithm employs an operator splitting technique, dividing the reaction-diffusion systems from the models into their constituent parts.

We integrate both the reaction and diffusion pieces via an implicit Euler method. We reduce the temporal and splitting errors by using a spectral deferred correction method, raising the temporal order and accuracy of the scheme with each correction iteration.

Our algorithm also uses continuous piecewise polynomials of high order on rectangular elements as our finite element approximation. This approximation improves the spatial discretization error over the piecewise linear polynomials typically used, especially when the spatial mesh is refined.

As part of these thesis work, we also present numerical simulations using our algorithm of one of the cardiac models mentioned, the Two-Current Model. We demonstrate the efficiency, accuracy and convergence rates of our numerical scheme by using mesh refinement studies and comparison of accuracy versus computational time. We conclude with a discussion of how our algorithm can be applied to more realistic models of cardiac electrical activity.

dc.identifier.uri

https://hdl.handle.net/10161/3883

dc.subject

Mathematics

dc.subject

cardiac models

dc.subject

spectral deferred correction

dc.title

A Spectral Deferred Correction Method for Solving Cardiac Models

dc.type

Dissertation

Files

Original bundle

Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
Bowen_duke_0066D_10943.pdf
Size:
1.08 MB
Format:
Adobe Portable Document Format

Collections