Researcher profile

Andrea Moiola

Andrea Moiola contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Interpolation of Hilbert and Sobolev Spaces: Quantitative Estimates and Counterexamples

This paper provides an overview of interpolation of Banach and Hilbert spaces, with a focus on establishing when equivalence of norms is in fact equality of norms in the key results of the theory. (In brief, our conclusion for the Hilbert space case is that, with the right normalisations, all the key results hold with equality of norms.) In the final section we apply the Hilbert space results to the Sobolev spaces $H^s(Ω)$ and $\widetilde{H}^s(Ω)$, for $s\in \mathbb{R}$ and an open $Ω\subset \mathbb{R}^n$. We exhibit examples in one and two dimensions of sets $Ω$ for which these scales of Sobolev spaces are not interpolation scales. In the cases when they are interpolation scales (in particular, if $Ω$ is Lipschitz) we exhibit examples that show that, in general, the interpolation norm does not coincide with the intrinsic Sobolev norm and, in fact, the ratio of these two norms can be arbitrarily large.

preprint2022arXiv

Spurious Quasi-Resonances in Boundary Integral Equations for the Helmholtz Transmission Problem

We consider the Helmholtz transmission problem with piecewise-constant material coefficients, and the standard associated direct boundary integral equations. For certain coefficients and geometries, the norms of the inverses of the boundary integral operators grow rapidly through an increasing sequence of frequencies, even though this is not the case for the solution operator of the transmission problem; we call this phenomenon that of spurious quasi-resonances. We give a rigorous explanation of why and when spurious quasi-resonances occur, and propose modified boundary integral equations that are not affected by them.

preprint2021arXiv

Boundary element methods for acoustic scattering by fractal screens

We study boundary element methods for time-harmonic scattering in $\mathbb{R}^n$ ($n=2,3$) by a fractal planar screen, assumed to be a non-empty bounded subset $Γ$ of the hyperplane $Γ_\infty=\mathbb{R}^{n-1}\times \{0\}$. We consider two distinct cases: (i) $Γ$ is a relatively open subset of $Γ_\infty$ with fractal boundary (e.g.\ the interior of the Koch snowflake in the case $n=3$); (ii) $Γ$ is a compact fractal subset of $Γ_\infty$ with empty interior (e.g.\ the Sierpinski triangle in the case $n=3$). In both cases our numerical simulation strategy involves approximating the fractal screen $Γ$ by a sequence of smoother "prefractal" screens, for which we compute the scattered field using boundary element methods that discretise the associated first kind boundary integral equations. We prove sufficient conditions on the mesh sizes guaranteeing convergence to the limiting fractal solution, using the framework of Mosco convergence. We also provide numerical examples illustrating our theoretical results.

preprint2021arXiv

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 $Ω\subset\mathbb R^n$, $\mathcal{D}(Ω)$ is dense in $\{u\in H^s(\mathbb R^n):{\rm supp}\, u\subset \overlineΩ\}$ whenever $\partialΩ$ has zero Lebesgue measure and $Ω$ is &#34;thick&#34; (in the sense of Triebel); and (ii) for a $d$-set $Γ\subset\mathbb R^n$ ($0<d<n$), $\{u\in H^{s_1}(\mathbb R^n):{\rm supp}\, u\subset Γ\}$ is dense in $\{u\in H^{s_2}(\mathbb R^n):{\rm supp}\, u\subset Γ\}$ 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 Γ\}=\{0\}$ for a given closed set $Γ\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.

preprint2020arXiv

A high frequency boundary element method for scattering by a class of multiple obstacles

We propose a boundary element method for problems of time-harmonic acoustic scattering by multiple obstacles in two dimensions, at least one of which is a convex polygon. By combining a Hybrid Numerical Asymptotic (HNA) approximation space on the convex polygon with standard polynomial-based approximation spaces on each of the other obstacles, we show that the number of degrees of freedom required in the HNA space to maintain a given accuracy needs to grow only logarithmically with respect to the frequency, as opposed to the (at least) linear growth required by standard polynomial-based schemes. This method is thus most effective when the convex polygon is many wavelengths in diameter and the small obstacles have a combined perimeter comparable to the problem wavelength.

preprint2020arXiv

Space-time discontinuous Galerkin approximation of acoustic waves with point singularities

We develop a convergence theory of space-time discretizations for the linear, 2nd-order wave equation in polygonal domains $Ω\subset\mathbb{R}^2$, possibly occupied by piecewise homogeneous media with different propagation speeds. Building on an unconditionally stable space-time DG formulation developed in [Moiola, Perugia 2018], we (a) prove optimal convergence rates for the space-time scheme with local isotropic corner mesh refinement on the spatial domain, and (b) demonstrate numerically optimal convergence rates of a suitable \emph{sparse} space-time version of the DG scheme. The latter scheme is based on the so-called \emph{combination formula}, in conjunction with a family of anisotropic space-time DG-discretizations. It results in optimal-order convergent schemes, also in domains with corners, with a number of degrees of freedom that scales essentially like the DG solution of one stationary elliptic problem in $Ω$ on the finest spatial grid. Numerical experiments for both smooth and singular solutions support convergence rate optimality on spatially refined meshes of the full and sparse space-time DG schemes.

preprint2018arXiv

Acoustic transmission problems: wavenumber-explicit bounds and resonance-free regions

We consider the Helmholtz transmission problem with one penetrable star-shaped Lipschitz obstacle. Under a natural assumption about the ratio of the wavenumbers, we prove bounds on the solution in terms of the data, with these bounds explicit in all parameters. In particular, the (weighted) $H^1$ norm of the solution is bounded by the $L^2$ norm of the source term, independently of the wavenumber. These bounds then imply the existence of a resonance-free strip beneath the real axis. The main novelty is that the only comparable results currently in the literature are for smooth, convex obstacles with strictly positive curvature, while here we assume only Lipschitz regularity and star-shapedness with respect to a point. Furthermore, our bounds are obtained using identities first introduced by Morawetz (essentially integration by parts), whereas the existing bounds use the much-more sophisticated technology of microlocal analysis and propagation of singularities. We also recap existing results that show that if the assumption on the wavenumbers is lifted, then no bound with polynomial dependence on the wavenumber is possible.

preprint2018arXiv

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.

preprint2017arXiv

A space-time Trefftz discontinuous Galerkin method for the acoustic wave equation in first-order formulation

We introduce a space-time Trefftz discontinuous Galerkin method for the first-order transient acoustic wave equations in arbitrary space dimensions, extending the one dimensional scheme of Kretzschmar et al. (2016, IMA J. Numer. Anal., 36, 1599-1635). Test and trial discrete functions are space-time piecewise polynomial solutions of the wave equations. We prove well-posedness and a priori error bounds in both skeleton-based and mesh-independent norms. The space-time formulation corresponds to an implicit time-stepping scheme, if posed on meshes partitioned in time slabs, or to an explicit scheme, if posed on &#34;tent-pitched&#34; meshes. We describe two Trefftz polynomial discrete spaces, introduce bases for them and prove optimal, high-order $h$-convergence bounds.

preprint2017arXiv

Sobolev spaces on non-Lipschitz subsets of $\mathbb{R}^n$ with application to boundary integral equations on fractal screens

We study properties of the classical fractional Sobolev spaces (or Bessel potential spaces) on non-Lipschitz subsets of $\mathbb{R}^n$. We investigate the extent to which the properties of these spaces, and the relations between them, that hold in the well-studied case of a Lipschitz open set, generalise to non-Lipschitz cases. Our motivation is to develop the functional analytic framework in which to formulate and analyse integral equations on non-Lipschitz sets. In particular we consider an application to boundary integral equations for wave scattering by planar screens that are non-Lipschitz, including cases where the screen is fractal or has fractal boundary.

preprint2016arXiv

A note on properties of the restriction operator on Sobolev spaces

In our companion paper (S.N. Chandler Wilde, D.P. Hewett, A. Moiola, Sobolev spaces on non-Lipschitz subsets of $\mathbb{R}^n$ with application to boundary integral equations on fractal screens, 2016) we studied a number of different Sobolev spaces on a general (non-Lipschitz) open subset $Ω$ of $\mathbb{R}^n$, defined as closed subspaces of the classical Bessel potential spaces $H^s(\mathbb{R}^n)$ for $s\in\mathbb{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.

preprint2015arXiv

A priori error analysis of space-time Trefftz discontinuous Galerkin methods for wave problems

We present and analyse a space-time discontinuous Galerkin method for wave propagation problems. The special feature of the scheme is that it is a Trefftz method, namely that trial and test functions are solution of the partial differential equation to be discretised in each element of the (space-time) mesh. The method considered is a modification of the discontinuous Galerkin schemes of Kretzschmar et al., and of Monk and Richter. For Maxwell&#39;s equations in one space dimension, we prove stability of the method, quasi-optimality, best approximation estimates for polynomial Trefftz spaces and (fully explicit) error bounds with high order in the meshwidth and in the polynomial degree. The analysis framework also applies to scalar wave problems and Maxwell&#39;s equations in higher space dimensions. Some numerical experiments demonstrate the theoretical results proved and the faster convergence compared to the non-Trefftz version of the scheme.

preprint2015arXiv

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.