The main goal of our thesis is to understand qualitative properties of solutions to the Cauchy problem for the semi-linear wave model with time-dependent speed and dissipation. We greatly benefited from very precise estimates for the corresponding linear problem in order to obtain the global existence (in time) of small data solutions. This reason motivated us to introduce very carefully a complete description for classification of our models: scattering, non-effective, effective, over-damping. We have considered those separately.

In this report, we study how finite differences can be used to simulate acoustic wave propagation originating from a point source in the ocean using the Helmholtz equation. How to model sediment layers and the vast size of the ocean is studied in particular. The finite differences are implemented with summation by parts operators with boundary conditions enforced with simultaneous approximation terms and projection. The numerical solver is combined with the WaveHoltz method to improve the performance. Sediment layers are handled with interface conditions and the domain is artificially expanded using absorbing layers. The absorbing layer is implemented with an alternative approach to the super-grid method where the domain expansion is accomplished by altering the wave speed rather than with coordinate transformations. To isolate these issues, other parameters such as variations in the ocean floor are neglected. With this simplification, cylindrical coordinates are used and the angular variation is assumed to be zero. This reduces the problem to a quasi-three-dimensional system. We study how the parameters of the alternative absorbing layer approach affect its quality. The numerical solver is verified on several test cases and appears to work according to theory. Finally, a semi-realistic simulation is carried out and the solution seems correct in this setting.

This work addresses the question of the efficient numerical solution of time-domain boundary integral equations with retarded potentials arising in the problems of acoustic and electromagnetic scattering. The convolutional form of the time-domain boundary operators allows to discretize them with the help of Runge-Kutta convolution quadrature. This method combines Laplace-transform and time-stepping approaches and requires the explicit form of the fundamental solution only in the Laplace domain to be known. Recent numerical and analytical studies revealed excellent properties of Runge-Kutta convolution quadrature, e.g. high convergence order, stability, low dissipation and dispersion. As a model problem, we consider the wave scattering in three dimensions. The convolution quadrature discretization of the indirect formulation for the three-dimensional wave equation leads to the lower triangular Toeplitz system of equations. Each entry of this system is a boundary integral operator with a kernel defined by convolution quadrature. In this work we develop an efficient method of almost linear complexity for the solution of this system based on the existing recursive algorithm. The latter requires the construction of many discretizations of the Helmholtz boundary single layer operator for a wide range of complex wavenumbers. This leads to two main problems: the need to construct many dense matrices and to evaluate many singular and near-singular integrals. The first problem is overcome by the use of data-sparse techniques, namely, the high-frequency fast multipole method (HF FMM) and H-matrices. The applicability of both techniques for the discretization of the Helmholtz boundary single-layer operators with complex wavenumbers is analyzed. It is shown that the presence of decay can favorably affect the length of the fast multipole expansions and thus reduce the matrix-vector multiplication times. The performance of H-matrices and the HF FMM is compared for a range of complex wavenumbers, and the strategy to choose between two techniques is suggested. The second problem, namely, the assembly of many singular and nearly-singular integrals, is solved by the use of the Huygens principle. In this work we prove that kernels of the boundary integral operators $w_n^h(d)$ ($h$ is the time step and $t_n=nh$ is the time) exhibit exponential decay outside of the neighborhood of $d=nh$ (this is the consequence of the Huygens principle). The size of the support of these kernels for fixed $h$ increases with $n$ as $n^a,a<1$, where $a$ depends on the order of the Runge-Kutta method and is (typically) smaller for Runge-Kutta methods of higher order. Numerical experiments demonstrate that theoretically predicted values of $a$ are quite close to optimal. In the work it is shown how this property can be used in the recursive algorithm to construct only a few matrices with the near-field, while for the rest of the matrices the far-field only is assembled. The resulting method allows to solve the three-dimensional wave scattering problem with asymptotically almost linear complexity. The efficiency of the approach is confirmed by extensive numerical experiments.

The primary focus of this thesis is to examine the detailed seismic structure of the northern Cascadia margin, including the Cascadia basin, the deformation front and the continental shelf. The results of this study are contributing towards understanding sediment deformation and tectonics on this margin. They also have important implications for exploration of hydrocarbons (oil and gas) and natural hazards (submarine landslides, earthquakes, tsunamis, and climate change). The first part of this thesis focuses on the role of gas hydrate in slope failure observed from multibeam bathymetry data on a frontal ridge near the deformation front off Vancouver Island margin using active-source ocean bottom seismometer (OBS) data collected in 2010. Volume estimates (∼ 0.33 km^3) of the slides observed on this margin indicate that these are capable of generating large (∼ 1 − 2 m) tsunamis. Velocity models from travel time inversion of wide angle reflections and refractions recorded on OBSs and vertical incidence single channel seismic (SCS) data were used to estimate gas hydrate concentrations using effective medium modeling. Results indicate a shallow high velocity hydrate layer with a velocity of 2.0 − 2.1 km/s that corresponds to a hydrate concentration of 40% at a depth of 100 m, and a bottom simulating reflector (BSR) at a depth of 265 − 275 m beneath the seafloor (mbsf). These are comparable to drilling results on an adjacent frontal ridge. Margin perpendicular normal faults that extend down to BSR depth were also observed on SCS and bathymetric data, two of which coincide with the sidewalls of the slump indicating that the lateral extent of the slump is controlled by these faults. Analysis of bathymetric data indicates, for the first time, that the glide plane occurs at the same depth as the shallow high velocity layer (100±10 mbsf). In contrast, the glide plane coincides with the depth of the BSR on an adjacent frontal ridge. In either case, our results suggest that the contrast in sediments strengthened by hydrates and overlying or underlying sediments where there is no hydrate is what causing the slope failure on this margin. The second part of this dissertation focuses on obtaining the detailed structure of the Cascadia basin and frontal ridge region using mirror imaging of few widely spaced OBS data. Using only a small airgun source (120 cu. in.), our results indicate structures that were previously not observed on the northern Cascadia margin. Specifically, OBS migration results show dual-vergence structure, which could be related to horizontal compression associated with subduction and low basal shear stress resulting from over-pressure. Understanding the physical and mechanical properties of the basal layer has important implications for understanding earthquakes on this margin. The OBS migrated image also clearly shows the continuity of reflectors which enabled the identification of thrust faults, and also shows the top of the igneous oceanic crust at 5−6 km beneath the seafloor, which were not possible to identify in single-channel and low-fold multi-channel seismic (MCS) data. The last part of this thesis focuses on obtaining detailed seismic structure of the Vancouver Island continental shelf from MCS data using frequency domain viscoacoustic full waveform inversion, which is first of its kind on this margin. Anelastic velocity and attenuation models, derived in this study to subseafloor depths of ∼ 2 km, are useful in understanding the deformation within the Tofino basin sediments, the nature of basement structures and their relationship with underlying accreted terranes such as the Crescent and the Pacific Rim terranes. Specifically, our results indicate a low-velocity zone (LVZ) with a contrast of 200 m/s within the Tofino basin sediment section at a depth 600 − 1000 mbsf over a lateral distance of 10 km. This LVZ is associated with high attenuation values (0.015 − 0.02) and could be a result of over pressured sediments or lithology changes associated with a high porosity layer in this potential hydrocarbon environment. Shallow high velocities of 4 − 5 km/s are observed in the mid-shelf region at depths > 1.5 km, which is interpreted as the shallowest occurrence of the Eocene volcanic Crescent terrane. The sediment velocities sharply increase about 10 km west of Vancouver Island, which probably corresponds to the underlying transition to the Mesozoic marine sedimentary Pacific Rim terrane. High attenuation values of 0.03 − 0.06 are observed at depths > 1 km, which probably corresponds to increased clay content and the presence of mineralized fluids. / Graduate / 0373 / 0372 / 0605 / subbarao@uvic.ca

Dans cette thèse, nous nous intéressons à l'étude de quatre équations d'évolution non-locales. Les solutions de ces quatre équations peuvent exploser en temps fini. Dans la théorie des équations d'évolution non-linéaires, une solution est qualifiée de globale si elle est définie pour tout temps positif. Au contraire, si une solution existe seulement sur un intervalle de temps [0; T) borné, elle est dite locale. Dans ce dernier cas et quand le temps maximal d'existence est relié à une alternative d'explosion, on dit aussi que la solution explose en temps fini. Dans un premier travail, nous considérons l'équation de Schrödinger non-linéaire avec une puissance fractionnaire du laplacien, et nous obtenons l'explosion de la solution en temps fini Tmax > 0 pour toute condition initiale positive et non-triviale dans le cas d'exposant sous-critique. Ensuite, nous étudions une équation des ondes amorties avec un potentiel d'espace-temps et un terme non-linéaire et non-local en temps. Nous obtenons un résultat d'existence locale d'une solution dans l'espace d'énergie sous des conditions restrictives sur les données initiales, la dimension de l'espace et la croissance du terme non-linéaire. De plus, nous obtenons l'explosion de la solution en temps fini pour toute condition initiale de moyenne strictement positive. De plus, nous étudions un problème de Cauchy pour l'équation d'évolution avec un p- Laplacien avec une non linéarité non-locale en temps. Dans ce cadre, nous nous intéressons à l'étude de l'existence locale d'une solution de cette équation ainsi qu'un résultat de non-existence de solution globale. Finalement, nous étudions l'intervalle maximal d'existence des solutions de l'équation des milieux poreux avec un terme non-linéaire non-local en temps. / In this thesis, we study four non-local evolution equations. The solutions of these four equations can blow up in finite time. In the theory of nonlinear evolution equations, a solution is qualified as global if it isdefined for any time. Otherwise, if a solution exists only on a bounded interval [0; T), it is called local solution. In this case and when the maximum time of existence is related to a blow up alternative, we say that the solution blows up in finite time. First, we consider the nonlinear Schröodinger equation with a fractional power of the Laplacien operator, and we get a blow up result in finite time Tmax > 0 for any non-trivial non-negative initial condition in the case of sub-critical exponent. Next, we study a damped wave equation with a space-time potential and a non-local in time non-linear term. We obtain a result of local existence of a solution in the energy space under some restrictions on the initial data, the dimension of the space and the growth of nonlinear term. Additionally, we get a blow up result of the solution in finite time for any initial condition positive on average. In addition, we study a Cauchy problem for the evolution p-Laplacien equation with nonlinear memory. We study the local existence of a solution of this equation as well as a result of non-existence of global solution. Finally, we study the maximum interval of existence of solutions of the porous medium equation with a nonlinear non-local in time term.

Humanitarian clearance of minefields is most often carried out by hand, conventionally using a a metal detector and a probe. Detection is a very slow process, as every piece of detected metal must treated as if it were a landmine and carefully probed and excavated, while many of them are not. The process can be safely sped up by use of Ground-Penetrating Radar (GPR) to image the subsurface, to verify metal detection results and safely ignore any objects which could not possibly be a landmine. In this thesis, we explore the possibility of using Full Wave Inversion (FWI) to improve GPR imaging for landmine detection. Posing the imaging task as FWI means solving the large-scale, non-linear and ill-posed optimisation problem of determining the physical parameters of the subsurface (such as electrical permittivity) which would best reproduce the data. This thesis begins by giving an overview of all the mathematical and implementational aspects of FWI, so as to provide an informative text for both mathematicians (perhaps already familiar with other inverse problems) wanting to contribute to the mine detection problem, as well as a wider engineering audience (perhaps already working on GPR or mine detection) interested in the mathematical study of inverse problems and FWI.We present the first numerical 3D FWI results for GPR, and consider only surface measurements from small-scale arrays as these are suitable for our application. The FWI problem requires an accurate forward model to simulate GPR data, for which we use a hybrid finite-element boundary-integral solver utilising first order curl-conforming N\'d\'{e}lec (edge) elements. We present a novel `line search' type algorithm which prioritises inversion of some target parameters in a region of interest (ROI), with the update outside of the area defined implicitly as a function of the target parameters. This is particularly applicable to the mine detection problem, in which we wish to know more about some detected metallic objects, but are not interested in the surrounding medium. We may need to resolve the surrounding area though, in order to account for the target being obscured and multiple scattering in a highly cluttered subsurface. We focus particularly on spatial sensitivity of the inverse problem, using both a singular value decomposition to analyse the Jacobian matrix, as well as an asymptotic expansion involving polarization tensors describing the perturbation of electric field due to small objects. The latter allows us to extend the current theory of sensitivity in for acoustic FWI, based on the Born approximation, to better understand how polarization plays a role in the 3D electromagnetic inverse problem. Based on this asymptotic approximation, we derive a novel approximation to the diagonals of the Hessian matrix which can be used to pre-condition the GPR FWI problem.

No description available.

No description available.