Mathematical models of electric activity in cardiac tissue are becoming increasingly powerful tools in the study of cardiac arrhythmias. Considered here are mathematical models based on ordinary differential equations (ODEs) and partial differential equations (PDEs) that describe the behaviour of this electrical activity. Generating an efficient numerical solution of these models is a challenging task, and in fact the physiological accuracy of tissue-scale models is often limited by the efficiency of the numerical solution process. In this thesis, we discuss two sets of experiments that test ideas for making the numerical solution process more efficient. In the first set of experiments, we examine the numerical solution of four single cell cardiac electrophysiological models, which consist solely of ODEs. We study the efficiency of using implicit-explicit Runge-Kutta (IMEX-RK) splitting methods to solve these models. We find that variable step-size implementations of IMEX-RK methods (ARK3 and ARK5) that take advantage of Jacobian structure clearly outperform most methods commonly used in practice for two of the models, and they outperform all methods commonly used in practice for the remaining models. In the second set of experiments, we examine the solution of the bidomain model, a model consisting of both ODEs and PDEs that are typically solved separately. We focus these experiments on numerical methods for the solution of the two PDEs in the bidomain model. The most popular method for this task, the Crank-Nicolson method, produces unphysical oscillations; we propose a method based on a second-order L-stable singly diagonally implicit Runge-Kutta (SDIRK) method to eliminate these oscillations.<p>
We find that although the SDIRK method is able to eliminate these unphysical oscillations, it is only more efficient for crude error tolerances.
Identifer | oai:union.ndltd.org:USASK/oai:usask.ca:etd-02112009-124158 |
Date | 13 April 2009 |
Creators | Dean, Ryan Christopher |
Contributors | Patrick, George W., Osgood, Nathaniel, Eramian, Mark G., Spiteri, Raymond J. |
Publisher | University of Saskatchewan |
Source Sets | University of Saskatchewan Library |
Language | English |
Detected Language | English |
Type | text |
Format | application/pdf |
Source | http://library.usask.ca/theses/available/etd-02112009-124158/ |
Rights | unrestricted, I hereby certify that, if appropriate, I have obtained and attached hereto a written permission statement from the owner(s) of each third party copyrighted matter to be included in my thesis, dissertation, or project report, allowing distribution as specified below. I certify that the version I submitted is the same as that approved by my advisory committee. I hereby grant to University of Saskatchewan or its agents the non-exclusive license to archive and make accessible, under the conditions specified below, my thesis, dissertation, or project report in whole or in part in all forms of media, now or hereafter known. I retain all other ownership rights to the copyright of the thesis, dissertation or project report. I also retain the right to use in future works (such as articles or books) all or part of this thesis, dissertation, or project report. |
Page generated in 0.0697 seconds