Trefftz discontinuous Galerkin methods on unstructured meshes for the wave equation

We describe and analyse a space-time Trefftz discontinuous Galerkin method for the wave equation. The method is defined for unstructured meshes whose internal faces need not be aligned to the space-time axes. We show that the scheme is well-posed and dissipative, and we prove a priori error bounds for general Trefftz discrete spaces. A concrete discretisation can be obtained using piecewise polynomials that satisfy the wave equation elementwise.Comment: 8 pages, submitted to the XXIV CEDYA / XIV CMA conference, Cadiz 8-12 June 201

Plane wave approximation in linear elasticity

We consider the approximation of solutions of the time-harmonic linear elastic wave equation by linear combinations of plane waves. We prove algebraic orders of convergence both with respect to the dimension of the approximating space and to the diameter of the domain. The error is measured in Sobolev norms and the constants in the estimates explicitly depend on the problem wavenumber. The obtained estimates can be used in the h- and p-convergence analysis of wave-based finite element schemes

Density results for Sobolev, Besov and Triebel--Lizorkin spaces on rough sets

We investigate two density questions for Sobolev, Besov and Triebel--Lizorkin spaces on rough sets. Our main results, stated in the simplest Sobolev space setting, are that: (i) for an open set $\Omega\subset\mathbb R^n$, $\mathcal{D}(\Omega)$ is dense in $\{u\in H^s(\mathbb R^n):{\rm supp}\, u\subset \overline{\Omega}\}$ whenever $\partial\Omega$ has zero Lebesgue measure and $\Omega$ is "thick" (in the sense of Triebel); and (ii) for a $d$-set $\Gamma\subset\mathbb R^n$ ($0<d<n$), $\{u\in H^{s_1}(\mathbb R^n):{\rm supp}\, u\subset \Gamma\}$ is dense in $\{u\in H^{s_2}(\mathbb R^n):{\rm supp}\, u\subset \Gamma\}$ whenever $-\frac{n-d}{2}-m-1<s_{2}\leq s_{1}<-\frac{n-d}{2}-m$ for some $m\in\mathbb N_0$. For (ii), we provide concrete examples, for any $m\in\mathbb N_0$, where density fails when $s_1$ and $s_2$ are on opposite sides of $-\frac{n-d}{2}-m$. The results (i) and (ii) are related in a number of ways, including via their connection to the question of whether $\{u\in H^s(\mathbb R^n):{\rm supp}\, u\subset \Gamma\}=\{0\}$ for a given closed set $\Gamma\subset\mathbb R^n$ and $s\in \mathbb R$. They also both arise naturally in the study of boundary integral equation formulations of acoustic wave scattering by fractal screens. We additionally provide analogous results in the more general setting of Besov and Triebel--Lizorkin spaces.Comment: 38 pages, 6 figure

Can coercive formulations lead to fast and accurate solution of the Helmholtz equation?

A new, coercive formulation of the Helmholtz equation was introduced in [Moiola, Spence, SIAM Rev. 2014]. In this paper we investigate $h$-version Galerkin discretisations of this formulation, and the iterative solution of the resulting linear systems. We find that the coercive formulation behaves similarly to the standard formulation in terms of the pollution effect (i.e. to maintain accuracy as $k\to\infty$, $h$ must decrease with $k$ at the same rate as for the standard formulation). We prove $k$-explicit bounds on the number of GMRES iterations required to solve the linear system of the new formulation when it is preconditioned with a prescribed symmetric positive-definite matrix. Even though the number of iterations grows with $k$, these are the first such rigorous bounds on the number of GMRES iterations for a preconditioned formulation of the Helmholtz equation, where the preconditioner is a symmetric positive-definite matrix.Comment: 27 pages, 7 figure

Trefftz discontinuous Galerkin methods for acoustic scattering on locally refined meshes

We extend the a priori error analysis of Trefftz-discontinuous Galerkin methods for time-harmonic wave propagation problems developed in previous papers to acoustic scattering problems and locally refined meshes. To this aim, we prove refined regularity and stability results with explicit dependence of the stability constant on the wave number for non convex domains with non connected boundaries. Moreover, we devise a new choice of numerical flux parameters for which we can prove L2-error estimates in the case of locally refined meshes near the scatterer. This is the setting needed to develop a complete hp-convergence analysis

On the maximal Sobolev regularity of distributions supported by subsets of Euclidean space

This paper concerns the following question: given a subset E of Rn with empty interior and an integrability parameter 1<p<infinity, what is the maximal regularity s in R for which there exists a non-zero distribution in the Bessel potential Sobolev space Hs,p(Rn) that is supported in E? For sets of zero Lebesgue measure we apply well-known results on set capacities from potential theory to characterise the maximal regularity in terms of the Hausdorff dimension of E, sharpening previous results. Furthermore, we provide a full classification of all possible maximal regularities, as functions of p, together with the sets of values of p for which the maximal regularity is attained, and construct concrete examples for each case. Regarding sets with positive measure, for which the maximal regularity is non-negative, we present new lower bounds on the maximal Sobolev regularity supported by certain fat Cantor sets, which we obtain both by capacity-theoretic arguments, and by direct estimation of the Sobolev norms of characteristic functions. We collect several results characterising the regularity that can be achieved on certain special classes of sets, such as d-sets, boundaries of open sets, and Cartesian products, of relevance for applications in differential and integral equations

Stability results for the time-harmonic Maxwell equations with impedance boundary conditions

We consider the time-harmonic Maxwell equations with constant coefficients in a bounded, uniformly star-shaped polyhedron. We prove wavenumber-explicit norm bounds for weak solutions. This result is pivotal for convergence proofs in numerical analysis and may be a tool in the analysis of electromagnetic boundary integral operators

Vekua theory for the Helmholtz operator

Vekua operators map harmonic functions defined on domain in \mathbb R2R2 to solutions of elliptic partial differential equations on the same domain and vice versa. In this paper, following the original work of I. Vekua (Ilja Vekua (1907–1977), Soviet-Georgian mathematician), we define Vekua operators in the case of the Helmholtz equation in a completely explicit fashion, in any space dimension N ≥ 2. We prove (i) that they actually transform harmonic functions and Helmholtz solutions into each other; (ii) that they are inverse to each other; and (iii) that they are continuous in any Sobolev norm in star-shaped Lipschitz domains. Finally, we define and compute the generalized harmonic polynomials as the Vekua transforms of harmonic polynomials. These results are instrumental in proving approximation estimates for solutions of the Helmholtz equation in spaces of circular, spherical, and plane waves

A note on properties of the restriction operator on Sobolev spaces

In our companion paper [3] we studied a number of different Sobolev spaces on a general (non-Lipschitz) open subset Ω of Rn, defined as closed subspaces of the classical Bessel potential spaces Hs(Rn) for s∈R. These spaces are mapped by the restriction operator to certain spaces of distributions on Ω. In this note we make some observations about the relation between these spaces of global and local distributions. In particular, we study conditions under which the restriction operator is or is not injective, surjective and isometric between given pairs of spaces. We also provide an explicit formula for minimal norm extension (an inverse of the restriction operator in appropriate spaces) in a special case

Plane wave approximation of homogeneous Helmholtz solutions

In this paper, we study the approximation of solutions of the homogeneous Helmholtz equation Δu + ω 2 u = 0 by linear combinations of plane waves with different directions. We combine approximation estimates for homogeneous Helmholtz solutions by generalized harmonic polynomials, obtained from Vekua’s theory, with estimates for the approximation of generalized harmonic polynomials by plane waves. The latter is the focus of this paper. We establish best approximation error estimates in Sobolev norms, which are explicit in terms of the degree of the generalized polynomial to be approximated, the domain size, and the number of plane waves used in the approximations