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

Simon N. Chandler-Wilde

Simon N. Chandler-Wilde 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.SPmath-phmath.APmath.MPNumerical Analysismath.CV

Catalog footprint

What is connected

14works
8topics
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

14 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the half-space matching method for real wavenumber

The Half-Space Matching (HSM) method has recently been developed as a new method for the solution of 2D scattering problems with complex backgrounds, providing an alternative to Perfectly Matched Layers (PML) or other artificial boundary conditions. Based on half-plane representations for the solution, the scattering problem is rewritten as a system coupling (1) a standard finite element discretisation localised around the scatterer and (2) integral equations whose unknowns are traces of the solution on the boundaries of a finite number of overlapping half-planes contained in the domain. While satisfactory numerical results have been obtained for real wavenumbers, well-posedness and equivalence of this HSM formulation to the original scattering problem have been established only for complex wavenumbers. In the present paper we show, in the case of a homogeneous background, that the HSM formulation is equivalent to the original scattering problem also for real wavenumbers, and so is well-posed, provided the traces satisfy radiation conditions at infinity analogous to the standard Sommerfeld radiation condition. As a key component of our argument we show that, if the trace on the boundary of a half-plane satisfies our new radiation condition, then the corresponding solution to the half-plane Dirichlet problem satisfies the Sommerfeld radiation condition in a slightly smaller half-plane. We expect that this last result will be of independent interest, in particular in studies of rough surface scattering.

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
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

On Symmetries of the Feinberg-Zee Random Hopping Matrix

In this paper we study the spectrum $Σ$ of the infinite Feinberg-Zee random hopping matrix, a tridiagonal matrix with zeros on the main diagonal and random $\pm 1$'s on the first sub- and super-diagonals; the study of this non-selfadjoint random matrix was initiated in Feinberg and Zee (Phys. Rev. E 59 (1999), 6433--6443). Recently Hagger (arXiv:1412.1937, Random Matrices: Theory Appl.}, {\bf 4} 1550016 (2015)) has shown that the so-called periodic part $Σ_π$ of $Σ$, conjectured to be the whole of $Σ$ and known to include the unit disk, satisfies $p^{-1}(Σ_π) \subset Σ_π$ for an infinite class $S$ of monic polynomials $p$. In this paper we make very explicit the membership of $S$, in particular showing that it includes $P_m(λ) = λU_{m-1}(λ/2)$, for $m\geq 2$, where $U_n(x)$ is the Chebychev polynomial of the second kind of degree $n$. We also explore implications of these inverse polynomial mappings, for example showing that $Σ_π$ is the closure of its interior, and contains the filled Julia sets of infinitely many $p\in S$, including those of $P_m$, this partially answering a conjecture of the second author.

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

Coburn's Lemma and the Finite Section Method for Random Jacobi Operators

We study the spectra and pseudospectra of finite and infinite tridiagonal random matrices, in the case where each of the diagonals varies over a separate compact set, say $U,V,W\subset\mathbb{C}$. Such matrices are sometimes termed stochastic Toeplitz matrices $A_+$ in the semi-infinite case and stochastic Laurent matrices $A$ in the bi-infinite case. Their spectra, $Σ=$ spec $A$ and $Σ_+=$ spec $A_+$, are independent of $A$ and $A_+$ as long as $A$ and $A_+$ are pseudoergodic (in the sense of E.B. Davies, Commun. Math. Phys., 2001), which holds almost surely in the random case. This was shown in Davies (2001) for $A$; that the same holds for $A_+$ is one main result of this paper. We give upper and lower bounds on $Σ$ and $Σ_+$, and we explicitly compute a set $G$ that fills the gap between the two in the sense that $Σ\cup G=Σ_+$. We show that invertibility of one operator $A_+$ implies invertibility - and uniform boundedness of the inverses - of all finite square matrices with three diagonals in $U, V$ and $W$. This implies that the so-called finite section method for the approximate solution of a system $A_+x=b$ is applicable as soon as $A_+$ is invertible, and that the same method for estimating the spectrum of $A_+$ does not suffer from spectral pollution. Both results illustrate that tridiagonal stochastic Toeplitz operators share important properties of (classical) Toeplitz operators. One of our main tools is a new version of the Coburn lemma for classical Toeplitz operators, saying that a stochastic tridiagonal Toeplitz operator, if Fredholm, is always injective or surjective. In the final part we bound and compare the norms, and the norms of inverses, of bi-infinite, semi-infinite and finite tridiagonal matrices over $U$, $V$ and $W$. This allows the study of the resolvent norms, and hence the pseudospectra, of these operators and matrices.

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

Acoustic scattering: high frequency boundary element methods and unified transform methods

We describe some recent advances in the numerical solution of acoustic scattering problems. A major focus of the paper is the efficient solution of high frequency scattering problems via hybrid numerical-asymptotic boundary element methods. We also make connections to the unified transform method due to A.S. Fokas and co-authors, analysing particular instances of this method, proposed by J.A. DeSanto and co-authors, for problems of acoustic scattering by diffraction gratings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computing Fresnel Integrals via Modified Trapezium Rules

In this paper we propose methods for computing Fresnel integrals based on truncated trapezium rule approximations to integrals on the real line, these trapezium rules modified to take into account poles of the integrand near the real axis. Our starting point is a method for computation of the error function of complex argument due to Matta and Reichel ({\em J. Math. Phys.} {\bf 34} (1956), 298--307) and Hunter and Regan ({\em Math. Comp.} {\bf 26} (1972), 539--541). We construct approximations which we prove are exponentially convergent as a function of $N$, the number of quadrature points, obtaining explicit error bounds which show that accuracies of $10^{-15}$ uniformly on the real line are achieved with N=12, this confirmed by computations. The approximations we obtain are attractive, additionally, in that they maintain small relative errors for small and large argument, are analytic on the real axis (echoing the analyticity of the Fresnel integrals), and are straightforward to implement.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Spectra and Pseudospectra of a Class of Non-Self-Adjoint Random Matrices and Operators

In this paper we develop and apply methods for the spectral analysis of non-self-adjoint tridiagonal infinite and finite random matrices, and for the spectral analysis of analogous deterministic matrices which are pseudo-ergodic in the sense of E.B.Davies (Commun. Math. Phys. 216 (2001), 687-704). As a major application to illustrate our methods we focus on the "hopping sign model" introduced by J.Feinberg and A.Zee (Phys. Rev. E 59 (1999), 6433-6443), in which the main objects of study are random tridiagonal matrices which have zeros on the main diagonal and random $\pm 1$'s as the other entries. We explore the relationship between spectral sets in the finite and infinite matrix cases, and between the semi-infinite and bi-infinite matrix cases, for example showing that the numerical range and $p$-norm $\eps$-pseudospectra ($\eps>0$, $p\in [1,\infty]$) of the random finite matrices converge almost surely to their infinite matrix counterparts, and that the finite matrix spectra are contained in the infinite matrix spectrum $Σ$. We also propose a sequence of inclusion sets for $Σ$ which we show is convergent to $Σ$, with the $n$th element of the sequence computable by calculating smallest singular values of (large numbers of) $n\times n$ matrices. We propose similar convergent approximations for the 2-norm $\eps$-pseudospectra of the infinite random matrices, these approximations sandwiching the infinite matrix pseudospectra from above and below.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Condition number estimates for combined potential integral operators in acoustics and their boundary element discretisation

We consider the classical coupled, combined-field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle. In recent work, we have proved lower and upper bounds on the $L^2$ condition numbers for these formulations, and also on the norms of the classical acoustic single- and double-layer potential operators. These bounds to some extent make explicit the dependence of condition numbers on the wave number $k$, the geometry of the scatterer, and the coupling parameter. For example, with the usual choice of coupling parameter they show that, while the condition number grows like $k^{1/3}$ as $k\to\infty$, when the scatterer is a circle or sphere, it can grow as fast as $k^{7/5}$ for a class of `trapping' obstacles. In this paper we prove further bounds, sharpening and extending our previous results. In particular we show that there exist trapping obstacles for which the condition numbers grow as fast as $\exp(γk)$, for some $γ>0$, as $k\to\infty$ through some sequence. This result depends on exponential localisation bounds on Laplace eigenfunctions in an ellipse that we prove in the appendix. We also clarify the correct choice of coupling parameter in 2D for low $k$. In the second part of the paper we focus on the boundary element discretisation of these operators. We discuss the extent to which the bounds on the continuous operators are also satisfied by their discrete counterparts and, via numerical experiments, we provide supporting evidence for some of the theoretical results, both quantitative and asymptotic, indicating further which of the upper and lower bounds may be sharper.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Limit Operators, Collective Compactness, and the Spectral Theory of Infinite Matrices

In the first half of this text we explore the interrelationships between the abstract theory of limit operators (see e.g. the recent monographs of Rabinovich, Roch & Silbermann and Lindner) and the concepts and results of the generalised collectively compact operator theory introduced by Chandler-Wilde and Zhang. We build up to results obtained by applying this generalised collectively compact operator theory to the set of limit operators of an operator $A$. In the second half of this text we study bounded linear operators on the generalised sequence space $\ell^p(\Z^N,U)$, where $p\in [1,\infty]$ and $U$ is some complex Banach space. We make what seems to be a more complete study than hitherto of the connections between Fredholmness, invertibility, invertibility at infinity, and invertibility or injectivity of the set of limit operators, with some emphasis on the case when the operator $A$ is a locally compact perturbation of the identity. Especially, we obtain stronger results than previously known for the subtle limiting cases of $p=1$ and $\infty$. Our tools in this study are the results from the first half of the text and an exploitation of the partial duality between $\ell^1$ and $\ell^\infty$. Results in this second half of the text include a new proof that injectivity of all limit operators (the classic Favard condition) implies invertibility for a general class of almost periodic operators, and characterisations of invertibility at infinity and Fredholmness for operators in the so-called Wiener algebra. In two final chapters our results are illustrated by and applied to concrete examples. Firstly, we study the spectra and essential spectra of discrete Schrödinger operators (both self-adjoint and non-self-adjoint), including operators with almost periodic and random potentials. In the final chapter we apply our results to integral operators on $\R^N$.

0 citations0 reviewsNo code linkedOpen access