BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

David P. Hewett

David P. Hewett contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.FAmath.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
3topics
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

7 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 &#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.

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 &#34;mass-matrix&#34; 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