Source author record

Timo Betcke

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

ResearcherUnclaimed source record

math.NA Numerical Analysis math.AP math-ph math.MP physics.med-ph

Catalog footprint

What is connected

8works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

preprint2022arXiv

Frequency-robust preconditioning of boundary integral equations for acoustic transmission

The scattering and transmission of harmonic acoustic waves at a penetrable material are commonly modelled by a set of Helmholtz equations. This system of partial differential equations can be rewritten into boundary integral equations defined at the surface of the objects and solved with the boundary element method (BEM). High frequencies or geometrical details require a fine surface mesh, which increases the number of degrees of freedom in the weak formulation. Then, matrix compression techniques need to be combined with iterative linear solvers to limit the computational footprint. Moreover, the convergence of the iterative linear solvers often depends on the frequency of the wave field and the objects' characteristic size. Here, the robust PMCHWT formulation is used to solve the acoustic transmission problem. An operator preconditioner based on on-surface radiation conditions (OSRC) is designed that yields frequency-robust convergence characteristics. Computational benchmarks compare the performance of this novel preconditioned formulation with other preconditioners and boundary integral formulations. The OSRC preconditioned PMCHWT formulation effectively simulates large-scale problems of engineering interest, such as focused ultrasound treatment of osteoid osteoma.

preprint2021arXiv

Hybrid coupling of finite element and boundary element methods using Nitsche's method and the Calderon projection

In this paper we discuss a hybridised method for FEM-BEM coupling. The coupling from both sides use a Nitsche type approach to couple to the trace variable. This leads to a formulation that is robust and flexible with respect to approximation spaces and can easily be combined as a building block with other hybridised methods. Energy error norm estimates and the convergence of Jacobi iterations are proved and the performance of the method is illustrated on some computational examples.

preprint2018arXiv

Boundary element methods with weakly imposed boundary conditions

We consider boundary element methods where the Calderón projector is used for the system matrix and boundary conditions are weakly imposed using a particular variational boundary operator designed using techniques from augmented Lagrangian methods. Regardless of the boundary conditions, both the primal trace variable and the flux are approximated. We focus on the imposition of Dirichlet, mixed Dirichlet--Neumann, and Robin conditions. A salient feature of the Robin condition is that the conditioning of the system is robust also for stiff boundary conditions. The theory is illustrated by a series of numerical examples.

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.

preprint2017arXiv

Software frameworks for integral equations in electromagnetic scattering based on Calderón identities

In recent years there have been tremendous advances in the theoretical understanding of boundary integral equations for Maxwell problems. In particular, stable dual pairing of discretisation spaces have been developed that allow robust formulations of the preconditioned electric field, magnetic field and combined field integral equations. Within the BEM++ boundary element library we have developed implementations of these frameworks that allow an intuitive formulation of the typical Maxwell boundary integral formulations within a few lines of code. The basis of these developments is an efficient and robust implementations of Calderón identities together with a product algebra that hides and automates most technicalities involved in assembling Galerkin boundary integral equations. In this paper we demonstrate this framework and use it to derive very simple and robust software formulations of the standard preconditioned electric field, magnetic field and regularised combined field integral equations for Maxwell.

preprint2014arXiv

The Factorization method for three dimensional Electrical Impedance Tomography

The use of the Factorization method for Electrical Impedance Tomography has been proved to be very promising for applications in the case where one wants to find inhomogeneous inclusions in a known background. In many situations, the inspected domain is three dimensional and is made of various materials. In this case, the main challenge in applying the Factorization method consists in computing the Neumann Green's function of the background medium. We explain how we solve this difficulty and demonstrate the capability of the Factorization method to locate inclusions in realistic inhomogeneous three dimensional background media from simulated data obtained by solving the so-called complete electrode model. We also perform a numerical study of the stability of the Factorization method with respect to various modelling errors.

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.