Lamoureux, MichaelMcDonald, Matt2012-10-032012-11-132012-10-032012http://hdl.handle.net/11023/268The numerical modelling of wave equations is a common theme in many seismic applications, and is an important tool in understanding how the physical systems of interest react in the process of a seismic experiment. In this thesis we apply state-of-the-art numerical methods based on domain-decomposition combined with local pseudospectral spatial discretization, to three physically realistic models of seismic waves, namely their propagation in acoustic, elastic, and viscoelastic media. The Galerkin formulation solves the weak form of the partial differential equation representing wave propagation and naturally includes boundary integral terms to represent free surface, rigid, and absorbing boundary effects. Stability, accuracy, and computation issues are discussed in this context along with direct comparison with finite difference methodologies. This works is an effort to bridge the gap between the development of accurate physical models to represent the real world, as seen in seismic modelling, and the implementation of modern numerical techniques for the accurate solutions of partial differential equations.engUniversity of Calgary graduate students retain copyright ownership and moral rights for their thesis. You may use this material in any way that is permitted by the Copyright Act or through licensing that has been assigned to the document. For uses that are not allowable under copyright legislation or licensing, you are required to seek permission.MathematicsGeophysicsMathematicsModellingWavesFinite-ElementsElasticityAcousticsViscoelasticityNumericalmaster thesis10.11575/PRISM/27141