The imaging of the interior of the Earth using ground motion data, or seismic tomography, has been a subject of great interest for over a century. The full elastic wave equations are not typically used in standard tomography codes. Instead, the elastic waves are idealized as rays and only phase velocity and travel times are considered as input data. This results in the inability to resolve features which are on the order of one wavelength in scale. To overcome this problem, models which use the full elastic wave equation and consider total seismograms as input data have recently been developed. Unfortunately, those methods are much more computationally expensive and are only in their infancy. While the finite element method is very popular in many applications in solid mechanics, it is still not the method of choice in many seismic applications due to high pollution error. The pollution effect creates an increasing ratio of discretization to best approximation error for problems with increasing wave numbers. It has been shown that standard finite element methods cannot overcome this issue. To compensate, the meshes for solving high wave number problems in seismology must be increasingly refined, and are computationally infeasible due to the large scale requirements. A new generalized least squares method was recently introduced. The main idea is to select test spaces such that the discrete problem inherits the stability of the continuous problem. In this dissertation, a discontinuous Petrov-Galerkin method with optimal test functions for 2D time-harmonic seismic tomography problems is developed. First, the abstract DPG framework and key results are reviewed. 2D DPG methods for both static and time-harmonic elasticity problems are then introduced and results indicating the low-pollution property are shown. Finally, a matrix-free inexact-Newton method for the seismic inverse problem is developed. To conclude, results obtained from both DPG and standard continuous Galerkin discretization schemes are compared and the potential effectiveness of DPG as a practical seismic inversion tool is discussed. / text

The non-hydrostatic dynamical core ENDGame (Even Newer Dynamics for the General Atmospheric Modelling of the Environment) is extended into the thermosphere to test its feasability as a whole-atmosphere dynamical core that can simulate the large scale fluid dynamics of the whole atmosphere from the surface to the top of the thermosphere at 600km. This research may have applications in the development of a Sun-to-Earth modelling system involving the Met Office Unified Model, which will be useful for space weather forecasting and chemical climate modelling. Initial attempts to raise the top boundary of ENDGame above ∼100km give rise to instabilities. To explore the potential causes of these instabilities, a one dimensional column version of ENDGame: ENDGame1D, is developed to study the effects of vertically propagating acoustic waves in the dynamical core. A 2D ray-tracing scheme is also developed, which accounts for the numerical effects on wave propagation. It is found that ENDGame’s numerics have a tendency towards the excessive focussing of wave energy towards vertical propagation, and have poor handling of large amplitude waves, also being unable to handle shocks. A key finding is that the physical processes of vertical molecular viscosity and diffusion prevent the excessive growth of wave amplitudes in the thermosphere in ENDGame, which may be crucial to improving ENDGame’s stability as it is extended upwards. Therefore, a fully implicit-in-time implementation of vertical molecular viscosity and diffusion is developed in both ENDGame1D and the full three-dimensional version of ENDGame: ENDGame3D. A new scheme is developed to deal with the viscous and diffusive terms with the dynamics terms in a fully coupled way to avoid time-splitting errors that may arise. The combination of a small amount of off-centring of ENDGame’s semi-implicit formulation and the inclusion of vertical molecular viscosity and diffusion act to make ENDGame significantly more stable, as long as the simulation is able to remain stable up to the molecularly diffused region above an altitude of ∼130km.