Finite difference based error estimation for boundary value ODEs, 1D PDEs and 2D PDEs

Show simple item record

dc.contributor.advisor Muir, Paul
dc.creator Murtha, Thomas
dc.date.accessioned 2017-09-14T13:45:22Z
dc.date.available 2017-09-14T13:45:22Z
dc.date.issued 2017
dc.identifier.uri http://library2.smu.ca/handle/01/27076
dc.description 1 online resource (77 p.) : col. ill.
dc.description Includes abstract and appendix.
dc.description Includes bibliographical references (p. 45-46).
dc.description.abstract Differential equations are involved in many felds. Often these cannot be solved exactly, so they must be approximated. To do this, we derive a fourth order finite difference scheme for non-uniform meshes. This is paired with a second order finite difference scheme for non-uniform meshes to generate an error estimate for the second order finite difference scheme. The difference between the solutions obtained from the two schemes is used to generate an error estimate for the second order finite difference scheme. The schemes are tested in three cases for order of convergence and quality of error estimate. The first test case is boundary value ordinary differential equations, the second test case is 1D partial differential equations, and the third case is 2D partial differential equations. There were two methods to define boundary conditions in the 2D partial differential equation, either in Dirichlet or Neumann form. We found that both finite difference schemes had an experimental order of convergence consistent with the expected theoretical order of convergence in the boundary value ordinary differential equation and 1D partial differential equation (PDE) cases. This was also demonstrated in the 2D PDE case when using boundary conditions in Neumann form. Only the fourth order scheme showed experimentally and theoretically consistent convergence rates on 2D PDEs when using Dirichlet boundary conditions. The error estimate was within one order of magnitude of the true error for the second order finite difference scheme in all cases where the subintervals in the meshes were small enough to allow the Taylor's series, upon which the finite difference schemes are based, to hold. We conclude that the schemes can be used within a computational algorithm incorporating adaptive meshing and error control features. en_CA
dc.language.iso en en_CA
dc.publisher Halifax, N.S. : Saint Mary's University
dc.title Finite difference based error estimation for boundary value ODEs, 1D PDEs and 2D PDEs en_CA
dc.type Text en_CA
thesis.degree.name Bachelor of Science (Honours Mathematics)
thesis.degree.name Bachelor of Science (Honours Computing Science)
thesis.degree.level Undergraduate
thesis.degree.discipline Mathematics and Computing Science
thesis.degree.grantor Saint Mary's University (Halifax, N.S.)


Files in this item

 
 

This item appears in the following Collection(s)

Show simple item record