BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David P. Hewett

David P. Hewett appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.FAmath.APNumerical Analysis

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

12 published item(s)

preprint2022arXiv

Accelerated Calderón preconditioning for Maxwell transmission problems

We investigate a range of techniques for the acceleration of Calderón (operator) preconditioning in the context of boundary integral equation methods for electromagnetic transmission problems. Our objective is to mitigate as far as possible the high computational cost of the barycentrically-refined meshes necessary for the stable discretisation of operator products. Our focus is on the well-known PMCHWT formulation, but the techniques we introduce can be applied generically. By using barycentric meshes only for the preconditioner and not for the original boundary integral operator, we achieve significant reductions in computational cost by (i) using "reduced" Calderón preconditioners obtained by discarding constituent boundary integral operators that are not essential for regularisation, and (ii) adopting a "bi-parametric" approach in which we use a lower quality (cheaper) $\mathcal{H}$-matrix assembly routine for the preconditioner than for the original operator, including a novel approach of discarding far-field interactions in the preconditioner. Using the boundary element software Bempp (www.bempp.com), we compare the performance of different combinations of these techniques in the context of scattering by multiple dielectric particles. Applying our accelerated implementation to 3D electromagnetic scattering by an aggregate consisting of 8 monomer ice crystals of overall diameter 1cm at 664GHz leads to a 99% reduction in memory cost and at least a 75% reduction in total computation time compared to a non-accelerated implementation.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 "thick" (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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Calderón preconditioning of PMCHWT boundary integral equations for scattering by multiple absorbing dielectric particles

We consider the simulation of electromagnetic scattering by single and multiple isotropic homogeneous dielectric particles using boundary integral equations. Galerkin discretizations of the classical Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) boundary integral equation formulation provide accurate solutions for complex particle geometries, but are well-known to lead to ill-conditioned linear systems. In this paper we carry out an experimental investigation into the performance of Calderón preconditioning techniques for single and multiple absorbing obstacles, which involve a squaring of the PMCHWT operator to produce a well-conditioned second-kind formulation. For single-particle scattering configurations we find that Calderón preconditioning is actually often outperformed by simple "mass-matrix" preconditioning, i.e. working with the strong form of the discretized PMCHWT operator. In the case of scattering by multiple particles we find that a significant saving in computational cost can be obtained by performing block-diagonal Calderón preconditioning in which only the self-interaction blocks are preconditioned. Using the boundary element software library Bempp (www.bempp.com) the numerical performance of the different methods is compared for a range of wavenumbers, particle geometries and complex refractive indices relevant to the scattering of light by atmospheric ice crystals.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Acoustic scattering by fractal screens: mathematical formulations and wavenumber-explicit continuity and coercivity estimates

We consider time-harmonic acoustic scattering by planar sound-soft (Dirichlet) and sound-hard (Neumann) screens. In contrast to previous studies, in which the domain occupied by the screen is assumed to be Lipschitz or smoother, we consider screens occupying an arbitrary bounded open set in the plane. Thus our study includes cases where the closure of the domain occupied by the screen has larger planar Lebesgue measure than the screen, as can happen, for example, when the screen has a fractal boundary. We show how to formulate well-posed boundary value problems for such scattering problems, our arguments depending on results on the coercivity of the acoustic single-layer and hypersingular boundary integral operators, and on properties of Sobolev spaces on general open sets which appear to be new. Our analysis teases out the explicit wavenumber dependence of the continuity and coercivity constants of the boundary integral operators, viewed as mappings between fractional Sobolev spaces, this in part extending previous results of Ha-Duong. We also consider the complementary problem of propagation through a bounded aperture in an infinite planar screen.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shadow boundary effects in hybrid numerical-asymptotic methods for high frequency scattering

The hybrid numerical-asymptotic (HNA) approach aims to reduce the computational cost of conventional numerical methods for high frequency wave scattering problems by enriching the numerical approximation space with oscillatory basis functions, chosen based on partial knowledge of the high frequency solution asymptotics. In this paper we propose a new methodology for the treatment of shadow boundary effects in HNA boundary element methods, using the classical geometrical theory of diffraction phase functions combined with mesh refinement. We develop our methodology in the context of scattering by a class of sound-soft nonconvex polygons, presenting a rigorous numerical analysis (supported by numerical results) which proves the effectiveness of our HNA approximation space at high frequencies. Our analysis is based on a study of certain approximation properties of the Fresnel integral and related functions, which govern the shadow boundary behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A frequency-independent boundary element method for scattering by two-dimensional screens and apertures

We propose and analyse a hybrid numerical-asymptotic $hp$ boundary element method for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation space enriched with oscillatory basis functions, chosen to capture the high frequency asymptotics of the solution. We provide a rigorous frequency-explicit error analysis which proves that the method converges exponentially as the number of degrees of freedom $N$ increases, and that to achieve any desired accuracy it is sufficient to increase $N$ in proportion to the square of the logarithm of the frequency as the frequency increases (standard boundary element methods require $N$ to increase at least linearly with frequency to retain accuracy). Our numerical results suggest that fixed accuracy can in fact be achieved at arbitrarily high frequencies with a frequency-independent computational cost, when the oscillatory integrals required for implementation are computed using Filon quadrature. We also show how our method can be applied to the complementary "breakwater" problem of propagation through an aperture in an infinite sound-hard screen.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

In this paper we propose and analyse a hybrid numerical-asymptotic boundary element method for the solution of problems of high frequency acoustic scattering by a class of sound-soft nonconvex polygons. The approximation space is enriched with carefully chosen oscillatory basis functions; these are selected via a study of the high frequency asymptotic behaviour of the solution. We demonstrate via a rigorous error analysis, supported by numerical examples, that to achieve any desired accuracy it is sufficient for the number of degrees of freedom to grow only in proportion to the logarithm of the frequency as the frequency increases, in contrast to the at least linear growth required by conventional methods. This appears to be the first such numerical analysis result for any problem of scattering by a nonconvex obstacle. Our analysis is based on new frequency-explicit bounds on the normal derivative of the solution on the boundary and on its analytic continuation into the complex plane.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diffraction by a right-angled impedance wedge: an edge source formulation

This paper concerns the frequency domain problem of diffraction of a plane wave incident on an infinite right-angled wedge on which impedance (absorbing) boundary conditions are imposed. It is demonstrated that the exact Sommerfeld-Malyuzhinets contour integral solution for the diffracted field can be transformed to a line integral over a physical variable along the diffracting edge. This integral can be interpreted as a superposition of secondary point sources (with directivity) positioned along the edge, in the spirit of the edge source formulations for rigid (sound-hard) wedges derived in [Svensson et al., Acta Acustica/Acustica 95, 2009, pp.~568-572]. However, when surface waves are present the physical interpretation of the edge source integral must be altered: it no longer represents solely the diffracted field, but rather includes surface wave contributions.

0 citations0 reviewsNo code linkedOpen access