A first order system least squares method for the Helmholtz equation

We present a first order system least squares (FOSLS) method for the Helmholtz equation at high wave number k, which always deduces Hermitian positive definite algebraic system. By utilizing a non-trivial solution decomposition to the dual FOSLS problem which is quite different from that of standard finite element method, we give error analysis to the hp-version of the FOSLS method where the dependence on the mesh size h, the approximation order p, and the wave number k is given explicitly. In particular, under some assumption of the boundary of the domain, the L2 norm error estimate of the scalar solution from the FOSLS method is shown to be quasi optimal under the condition that kh/p is sufficiently small and the polynomial degree p is at least O(\log k). Numerical experiments are given to verify the theoretical results

Numerical Methods and Algorithms for High Frequency Wave Scattering Problems in Homogeneous and Random Media

This dissertation consists of four integral parts with a unified objective of developing efficient numerical methods for high frequency time-harmonic wave equations defined on both homogeneous and random media. The first part investigates the generalized weak coercivity of the acoustic Helmholtz, elastic Helmholtz, and time-harmonic Maxwell wave operators. We prove that such a weak coercivity holds for these wave operators on a class of more general domains called generalized star-shape domains. As a by-product, solution estimates for the corresponding Helmholtz-type problems are obtained. The second part of the dissertation develops an absolutely stable (i.e. stable in all mesh regimes) interior penalty discontinuous Galerkin (IP-DG) method for the elastic Helmholtz equations. A special mesh-dependent sesquilinear form is proposed and is shown to be weakly coercive in all mesh regimes. We prove that the proposed IP-DG method converges with optimal rate with respect to the mesh size. Numerical experiments are carried out to demonstrate the theoretical results and compare this method to the standard finite element method. The third part of the dissertation develops a Monte Carlo interior penalty discontinuous Galerkin (MCIP-DG) method for the acoustic Helmholtz equation defined on weakly random media. We prove that the solution to the random Helmholtz problem has a multi-modes expansion (i.e., a power series in a medium- related small parameter). Using this multi-modes expansion an efficient and accurate numerical method for computing moments of the solution to the random Helmholtz problem is proposed. The proposed method is also shown to converge optimally. Numerical experiments are carried out to compare the new multi-modes MCIP-DG method to a classical Monte Carlo method. The last part of the dissertation develops a theoretical framework for Schwarz pre- conditioning methods for general nonsymmetric and indefinite variational problems which are discretized by Galerkin-type discretization methods. Such a framework has been missing in the literature and is of great theoretical and practical importance for solving convection-diffusion equations and Helmholtz-type wave equations. Condition number estimates for the additive and hybrid Schwarz preconditioners are established under some structure assumptions. Numerical experiments are carried out to test the new framework

On stability of discretizations of the Helmholtz equation (extended version)

We review the stability properties of several discretizations of the Helmholtz equation at large wavenumbers. For a model problem in a polygon, a complete $k$-explicit stability (including $k$-explicit stability of the continuous problem) and convergence theory for high order finite element methods is developed. In particular, quasi-optimality is shown for a fixed number of degrees of freedom per wavelength if the mesh size $h$ and the approximation order $p$ are selected such that $kh/p$ is sufficiently small and $p = O(\log k)$, and, additionally, appropriate mesh refinement is used near the vertices. We also review the stability properties of two classes of numerical schemes that use piecewise solutions of the homogeneous Helmholtz equation, namely, Least Squares methods and Discontinuous Galerkin (DG) methods. The latter includes the Ultra Weak Variational Formulation