Source author record

Peter Monk

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2022arXiv

A Spectral Target Signature for Thin Surfaces with Higher Order Jump Conditions

In this paper we consider the inverse problem of determining structural properties of a thin anisotropic and dissipative inhomogeneity in ${\mathbb R}^m$, $m=2,3$ from scattering data. In the asymptotic limit as the thickness goes to zero, the thin inhomogeneity is modeled by an open $m-1$ dimensional manifold (here referred to as screen), and the field inside is replaced by jump conditions on the total field involving a second order surface differential operator. We show that all the surface coefficients (possibly matrix valued and complex) are uniquely determined from far field patterns of the scattered fields due to infinitely many incident plane waves at a fixed frequency. Then we introduce a target signature characterized by a novel eigenvalue problem such that the eigenvalues can be determined from measured scattering data, adapting the approach in \cite{Screens}. Changes in the measured eigenvalues are used to identified changes in the coefficients without making use of the governing equations that model the healthy screen. In our investigation the shape of the screen is known, since it represents the object being evaluated. We present some preliminary numerical results indicating the validity of our inversion approach.

preprint2022arXiv

Impedance eigenvalues in linear elasticity

This paper is devoted to studying impedance eigenvalues (that is, eigenvalues of a particular Dirichlet-to-Neumann map) for the time harmonic linear elastic wave problem, and their potential use as target-signatures for fluid-solid interaction problems. We first consider several possible families of eigenvalues of the elasticity problem, focusing on certain impedance eigenvalues that are an analogue of Steklov eigenvalues. We show that one of these families arises naturally in inverse scattering. We also analyse their approximation from far field measurements of the scattered pressure field in the fluid, and illustrate several alternative methods of approximation in the case of an isotropic elastic disk.

preprint2022arXiv

The Time Domain Linear Sampling Method for determining the shape of multiple scatterers using electromagnetic waves

The time domain linear sampling method (TD-LSM) solves inverse scattering problems using time domain data by creating an indicator function for the support of the unknown scatterer. It involves only solving a linear integral equation called the near-field equation using different data from sampling points that probe the domain where the scatterer is located. To date, the method has been used for the acoustic wave equation and has been tested for several different types of scatterers, i.e. sound hard, impedance, and penetrable, and for waveguides. In this paper, we extend the TD-LSM to the time dependent Maxwell's system with impedance boundary conditions - a similar analysis handles the case of a perfect electric conductor (PEC). We provide an analysis that supports the use of the TD-LSM for this problem, and preliminary numerical tests of the algorithm. Our analysis relies on the Laplace transform approach previously used for the acoustic wave equation. This is the first application of the TD-LSM in electromagnetism.

preprint2021arXiv

Target signatures for thin surfaces

We investigate an inverse scattering problem for a thin inhomogeneous scatterer in R^m, m = 2,3, which we model as a m-1 dimensional open surface. The scatterer is referred to as a screen. The goal is to design target signatures that are computable from scattering data in order to detect changes in the material properties of the screen. This target signature is characterized by a mixed Steklov eigenvalue problem for a domain whose boundary contains the screen. We show that the corresponding eigenvalues can be determined from appropriately modified scattering data by using the generalized linear sampling method. A weaker justification is provided for the classical linear sampling method. Numerical experiments are presented to support our theoretical results.

preprint2020arXiv

$L^\infty$ norm error estimates for HDG methods applied to the Poisson equation with an application to the Dirichlet boundary control problem

We prove quasi-optimal $L^\infty$ norm error estimates (up to logarithmic factors) for the solution of Poisson's problem by the standard Hybridizable Discontinuous Galerkin (HDG) method. Although such estimates are available for conforming and mixed finite element methods, this is the first proof for HDG. The method of proof is motivated by known $L^\infty$ norm estimates for mixed finite elements. We show two applications: the first is to prove optimal convergence rates for boundary flux estimates, and the second is to prove that numerically observed convergence rates for the solution of a Dirichlet boundary control problem are to be expected theoretically. Numerical examples show that the predicted rates are seen in practice.

preprint2020arXiv

HDG and CG methods for the Indefinite Time-Harmonic Maxwell's Equations under minimal regularity

We propose to use a hybridizable discontinuous Galerkin (HDG) method combined with the continuous Galerkin (CG) method to approximate Maxwell's equations. We make two contributions in this paper. First, even though there are many papers using HDG methods to approximate Maxwell's equations, to our knowledge they all assume that the coefficients are smooth (or constant). Here, we derive optimal convergence estimates for our HGD-CG approximation when the electromagnetic coefficients are {\em piecewise} smooth. This requires new techniques of analysis. Second, we use CG elements to approximate the Lagrange multiplier used to enforce the divergence condition and we obtain a discrete system in which we can decouple the discrete the Lagrange multiplier. Because we are using a continuous Lagrange multiplier space, the number of degrees of freedom devoted to this are less than for other HDG methods. We present numerical experiments to confirm our theoretical results.

preprint2016arXiv

A Plane Wave Discontinuous Galerkin Method with a Dirichlet-to-Neumann Boundary Condition for the Scattering Problem in Acoustics

We consider the numerical solution of an acoustic scattering problem by the Plane Wave Discontinuous Galerkin Method (PWDG) in the exterior of a bounded domain in $\mathbb{R}^2$. In order to apply the PWDG method, we introduce an artificial boundary to truncate the domain, and we impose a non-local Dirichlet-to-Neumann (DtN) boundary conditions on the artificial curve. To define the method, we introduce new consistent numerical fluxes that incorporate the truncated series of the DtN map. Error estimates with respect to the truncation order of the DtN map, and with respect to mesh width are derived. Numerical results suggest that the accuracy of the PWDG method for the acoustic scattering problem can be significantly improved by using DtN boundary conditions.

preprint2015arXiv

A tent pitching scheme motivated by Friedrichs theory

Certain Friedrichs systems can be posed on Hilbert spaces normed with a graph norm. Functions in such spaces arising from advective problems are found to have traces with a weak continuity property at points where the inflow and outflow boundaries meet. Motivated by this continuity property, an explicit space-time finite element scheme of the tent pitching type, with spaces that conform to the continuity property, is designed. Numerical results for a model one-dimensional wave propagation problem are presented.

preprint2015arXiv

Time Dependent Scattering from a Grating

Computing the electromagnetic field due to a periodic grating is critical for assessing the performance of thin film solar voltaic devices. In this paper we investigate the computation of these fields in the time domain (similar problems also arise in simulating antennas). Assuming a translation invariant periodic grating this reduces to solving the wave equation in a periodic domain. Materials used in practical devices have frequency dependent coefficients, and we provide a first proof of existence and uniqueness for a general class of such materials. Using Convolution Quadrature we can then prove time stepping error estimates. We end with some preliminary numerical results that demonstrate the convergence and stability of the scheme.

preprint2014arXiv

Residual based adaptivity and PWDG methods for the Helmholtz equation

We present a study of two residual a posteriori error indicators for the Plane Wave Discontinuous Galerkin (PWDG) method for the Helmholtz equation. In particular we study the h-version of PWDG in which the number of plane wave directions per element is kept fixed. First we use a slight modification of the appropriate a priori analysis to determine a residual indicator. Numerical tests show that this is reliable but pessimistic in that the ratio between the true error and the indicator increases as the mesh is refined. We therefore introduce a new analysis based on the observation that sufficiently many plane waves can approximate piecewise linear functions as the mesh is refined. Numerical results demonstrate an improvement in the efficiency of the indicators.

preprint2014arXiv

The Time Domain Lippmann-Schwinger Equation and Convolution Quadrature

We consider time domain acoustic scattering from a penetrable medium with a variable sound speed. This problem can be reduced to solving a time domain volume Lippmann-Schwinger integral equation. Using convolution quadrature in time and trigonometric collocation in space we can compute an approximate solution. We prove that the time domain Lippmann-Schwinger equation has a unique solution and prove conditional convergence and error estimates for the fully discrete solution for smooth sound speeds. Preliminary numerical results show that the method behaves well even for discontinuous sound speeds.