Boundary value problems (BVPs) are systems of ordinary differential equations (ODEs) with boundary conditions imposed at two or more distinct points. Such problems arise within mathematical models in a wide variety of applications. Numerically solving BVPs for ODEs generally requires the use of a series of complex numerical algorithms. Fortunately, when users are required to solve a BVP, they have a variety of BVP software packages from which to choose. However, all BVP software packages currently available implement a specific set of numerical algorithms and therefore function quite differently from each other. Users must often try multiple software packages on a BVP to find the one that solves their problem most effectively. This creates two problems for users. First, they must learn how to specify the BVP for each software package. Second, because each package solves a BVP with specific numerical algorithms, it becomes difficult to determine why one BVP package outperforms another. With that in mind, this thesis offers two contributions.
<p>
First, this thesis describes the development of the BVP component to the fully featured problem-solving environment (PSE) for the numerical solution of ODEs called pythODE. This software allows users to select between multiple numerical algorithms to solve BVPs. As a consequence, they are able to determine the numerical algorithms that are effective at each step of the solution process. Users are also able to easily add new numerical algorithms to the PSE. The effect of adding a new algorithm can be measured by making use of an automated test suite.
<p>
Second, the BVP component of pythODE is used to perform two research studies. In the first study, four known global-error estimation algorithms are compared in pythODE. These algorithms are based on the use of Richardson extrapolation, higher-order formulas, deferred corrections, and a conditioning constant. Through numerical experimentation, the algorithms based on
higher-order formulas and deferred corrections are shown to be computationally faster than Richardson extrapolation while having similar accuracy. In the second study, pythODE is used to
solve a newly developed one-dimensional model of the agglomerate in the catalyst layer of a proton exchange membrane fuel cell.
| Identifer | oai:union.ndltd.org:LACETR/oai:collectionscanada.gc.ca:SSU.etd-01182011-104957 |
| Date | 19 January 2011 |
| Creators | Boisvert, Jason J. |
| Contributors | Spiteri, Raymond J., Muir, Paul H., Ascher, Uri, Stanley, Kevin, Makaroff, Dwight |
| Publisher | University of Saskatchewan |
| Source Sets | Library and Archives Canada ETDs Repository / Centre d'archives des thèses électroniques de Bibliothèque et Archives Canada |
| Language | English |
| Detected Language | English |
| Type | text |
| Format | application/pdf |
| Source | http://library.usask.ca/theses/available/etd-01182011-104957/ |
| 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.0022 seconds