We consider the system of equations arising from mantle dynamics introduced by McKenzie (J. Petrology, 1985). In this multi-phase model, the fluid melt velocity obeys Darcy's law while the deformable "solid" matrix is governed by a highly viscous Stokes equation. The system is then coupled through mass conservation and compaction relations. Together these equations form a coupled Darcy-Stokes system on a continuous single-domain mixture of fluid and matrix. The porosity φ, representing the relative volume of fluid melt to the bulk volume, is assumed to be much smaller than one. When coupled with solute transport and thermal evolution in a time-dependent problem, the model transitions dynamically from a non-porous single phase solid to a two-phase porous medium. Such mixture models have an advantage for numerical approximation since the free boundary between the one and two-phase regions need not be determined explicitly. The equations of mantle dynamics apply to a wide range of applications in deep earth physics such as mid-ocean ridges, subduction zones, and hot-spot volcanism, as well as to glacier dynamics and other two-phase flows in porous media. Mid-ocean ridges form when viscous corner flow of the solid mantle focuses fluid toward a central ridge. Melt is believed to migrate upward until it reaches the lithospheric "tent" where it then moves toward the ridge in a high porosity band. Simulation of this physical phenomenon required confidence in numerical methods to handle highly heterogeneous porosity as well as the single-phase to two-phase transition. In this work we present a standard mixed finite element method for the equations of mantle dynamics and investigate its limitations for vanishing porosity. While stable and optimally convergent for porosity bounded away from zero, the stability estimates we obtain suggest, and numerical results show, the method becomes unstable as porosity approaches zero. Moreover, the fluid pressure is no longer a physical variable when the fluid phase disappears and thus is not a good variable for numerical methods. Inspired by the stability estimates of the standard method, we develop a novel stable mixed method with uniqueness and existence of solutions by studying a linear degenerate elliptic sub-problem akin to the Darcy part of the full model: [mathematical equation], where a and b satisfy a(0)=b(0)=0 and are otherwise positive, and the porosity φ ≥ 0 may be zero on a set of positive measure. Using scaled variables and mild assumptions on the regularity of φ, we develop a practical mass-conservative method based on lowest order Raviart-Thomas finite elements. Finally, we adapt the numerical method for the sub-problem to the full system of equations. We show optimal convergence for sufficiently smooth solutions for a compacting column and mid-ocean ridge-like corner flow examples, and investigate accuracy and stability for less regular problems / text

In this thesis, the stability of the family of subextremal Kerr-Newman space- times is studied in the case of linear scalar perturbations. That is, nondegenerate energy bounds (NEB) and integrated local energy decay (ILED) results are proved for solutions of the wave equation on the domain of outer communications. The main obstacles to the proof of these results are superradiance, trapping and their interaction. These difficulties are surmounted by localising solutions of the wave equation in phase space and applying the vector field method. Miraculously, as in the Kerr case, superradiance and trapping occur in disjoint regions of phase space and can be dealt with individually. Trapping is a high frequency obstruction to the proof whereas superradiance occurs at both high and low frequencies. The construction of energy currents for superradiant frequencies gives rise to an unfavourable boundary term. In the high frequency regime, this boundary term is controlled by exploiting the presence of a large parameter. For low superradiant frequencies, no such parameter is available. This difficulty is overcome by proving quantitative versions of mode stability type results. The mode stability result on the real axis is then applied to prove integrated local energy decay for solutions of the wave equation restricted to a bounded frequency regime. The (ILED) statement is necessarily degenerate due to the trapping effect. This implies that a nondegenerate (ILED) statement must lose differentiability. If one uses an (ILED) result that loses differentiability to prove (NEB), this loss is passed onto the (NEB) statement as well. Here, the geometry of the subextremal Kerr-Newman background is exploited to obtain the (NEB) statement directly from the degenerate (ILED) with no loss of differentiability.