Search

We consider frequency-domain acoustic scattering at a homogeneous star-shaped penetrable obstacle, whose shape is uncertain and modelled via a radial spectral parameterization with random coefficients. Using recent results on the stability of Helmholtz transmission problems with piecewise constant coefficients from [A. Moiola and E. A. Spence, Acoustic transmission problems: wavenumber-explicit bounds and resonance-free regions, Mathematical Models and Methods in Applied Sciences, 29 (2019), pp. 317-354] we obtain frequency-explicit statements on the holomorphic dependence of the scattered field and the far-field pattern on the stochastic shape parameters. This paves the way for applying general results on the efficient construction of high-dimensional surrogate models. We also take into account the effect of domain truncation by means of perfectly matched layers (PML). In addition, spatial regularity estimates which are explicit in terms of the wavenumber $k$ permit us to quantify the impact of finite-element Galerkin discretization using high-order Lagrangian finite-element spaces.

The equations governing incompressible Stokes flow are reformulated such that the velocity is sought in the space H(curl). This relaxed regularity assumption leads to conforming finite element methods using spaces common to discretizations of Maxwell's equations. A drawback of this approach, however, is that it is not immediately clear how to enforce Navier-slip boundary conditions. By recognizing the slip condition as a Robin boundary condition, we show that the continuous system is well-posed, propose finite element methods, and analyze the discrete system by deriving stability and a priori error estimates. Numerical experiments in 2D confirm the derived, optimal convergence rates., Comment: 29 pages, 5 figures

Using the mathematical theory of shape calculus as tool, we perform a rigorous investigation of the virtual work principle for the computation of electromagnetic forces in static settings. The main goal is to shed light on the widely held belief that the virtual work principle entails passively advecting the fields with the virtual displacement when computing deformation-dependent field energies gradients. The adjoint approach to differentiation of functionals under variational constraints provides a mathematical justification for this belief. However, it also shows that passively advecting the actual electromagnetic fields is sufficient only in the case of linear materials. In the general case, also fields arising as solutions of adjoint variational problems have to be taken into account.

We study frequency domain electromagnetic scattering at a bounded, penetrable, and inhomogeneous obstacle $ \Omega \subset \mathbb{R}^3 $. From the Stratton-Chu integral representation, we derive a new representation formula when constant reference coefficients are given for the interior domain. The resulting integral representation contains the usual layer potentials, but also volume potentials on $\Omega$. Then it is possible to follow a single-trace approach to obtain boundary integral equations perturbed by traces of compact volume integral operators with weakly singular kernels. The coupled boundary and volume integral equations are discretized with a Galerkin approach with usual Curl-conforming and Div-conforming finite elements on the boundary and in the volume. Compression techniques and special quadrature rules for singular integrands are required for an efficient and accurate method. Numerical experiments provide evidence that our new formulation enjoys promising properties.

We present a preconditioning method for the linear systems arising from the boundary element discretization of the Laplace hypersingular equation on a $2$-dimensional triangulated surface $\Gamma$ in $\mathbb{R}^3$. We allow $\Gamma$ to belong to a large class of geometries that we call polygonal multiscreens, which can be non-manifold. After introducing a new, simple conforming Galerkin discretization, we analyze a substructuring domain-decomposition preconditioner based on ideas originally developed for the Finite Element Method. The surface $\Gamma$ is subdivided into non-overlapping regions, and the application of the preconditioner is obtained via the solution of the hypersingular equation on each patch, plus a coarse subspace correction. We prove that the condition number of the preconditioned linear system grows poly-logarithmically with $H/h$, the ratio of the coarse mesh and fine mesh size, and our numerical results indicate that this bound is sharp. This domain-decomposition algorithm therefore guarantees significant speedups for iterative solvers, even when a large number of subdomains is used.

In this paper we will derive an non-local (``integral'') equation which transforms a three-dimensional acoustic transmission problem with \emph{variable} coefficients, non-zero absorption, and mixed boundary conditions to a non-local equation on a ``skeleton'' of the domain $\Omega\subset\mathbb{R}^{3}$, where ``skeleton'' stands for the union of the interfaces and boundaries of a Lipschitz partition of $\Omega$. To that end, we introduce and analyze abstract layer potentials as solutions of auxiliary coercive full space variational problems and derive jump conditions across domain interfaces. This allows us to formulate the non-local skeleton equation as a \emph{direct method} for the unknown Cauchy data of the solution of the original partial differential equation. We establish coercivity and continuity of the variational form of the skeleton equation based on auxiliary full space variational problems. Explicit expressions for Green's functions is not required and all our estimates are \emph{explicit} in the complex wave number.

We develop a mesh-based semi-Lagrangian discretization of the time-dependent incompressible Navier-Stokes equations with free boundary conditions recast as a non-linear transport problem for a momentum 1-form. A linearly implicit fully discrete version of the scheme enjoys excellent stability properties in the vanishing viscosity limit and is applicable to inviscid incompressible Euler flows. Conservation of energy and helicity are enforced separately., Comment: Submitted to Advances in Computational Mathematics (ACM), Topical Collection on Systematic Exploitation of Structural Properties in Electromagnetism

This work introduces ``generalized meshes", a type of meshes suited for the discretization of partial differential equations in non-regular geometries. Generalized meshes extend regular simplicial meshes by allowing for overlapping elements and more flexible adjacency relations. They can have several distinct ``generalized" vertices (or edges, faces) that occupy the same geometric position. These generalized facets are the natural degrees of freedom for classical conforming spaces of discrete differential forms appearing in finite and boundary element applications. Special attention is devoted to the representation of fractured domains and their boundaries. An algorithm is proposed to construct the so-called {\em virtually inflated mesh}, which correspond to a ``two-sided" mesh of a fracture. Discrete $d$-differential forms on the virtually inflated mesh are characterized as the trace space of discrete $d$-differential forms in the surrounding volume.

We study a new notion of trace operators and trace spaces for abstract Hilbert complexes. We introduce trace spaces as quotient spaces/annihilators. We characterize the kernels and images of the related trace operators and discuss duality relationships between trace spaces. We elaborate that many properties of the classical boundary traces associated with the Euclidean de Rham complex on bounded Lipschitz domains are rooted in the general structure of Hilbert complexes. We arrive at abstract trace Hilbert complexes that can be formulated using quotient spaces/annihilators. We show that, if a Hilbert complex admits stable "regular decompositions" with compact lifting operators, then the associated trace Hilbert complex is Fredholm. Incarnations of abstract concepts and results in the concrete case of the de Rham complex in three-dimensional Euclidean space will be discussed throughout., Comment: 33 pages

We consider the Helmholtz transmission problem with piecewise-constant material coefficients, and the standard associated direct boundary integral equations. For certain coefficients and geometries, the norms of the inverses of the boundary integral operators grow rapidly through an increasing sequence of frequencies, even though this is not the case for the solution operator of the transmission problem; we call this phenomenon that of spurious quasi-resonances. We give a rigorous explanation of why and when spurious quasi-resonances occur, and propose modified boundary integral equations that are not affected by them.

We develop novel first-kind boundary integral equations for Euclidean Dirac operators in 3D Lipschitz domains comprising square-integrable potentials and involving only weakly singular kernels. Generalized Garding inequalities are derived and we establish that the obtained boundary integral operators are Fredholm of index zero. Their finite dimensional kernels are characterized and we show that their dimension is equal to the number of topological invariants of the domain's boundary, in other words to the sum of its Betti numbers. This is explained by the fundamental discovery that the associated bilinear forms agree with those induced by the 2D surface Dirac operators for H-1/2 surface de Rham Hilbert complexes whose underlying inner-products are the non-local inner products defined through the classical single-layer boundary integral operators for the Laplacian. Decay conditions for well-posedness in natural energy spaces of the Dirac system in unbounded exterior domains are also presented., Comment: 27 pages, 2 figures

We present a quantum algorithm to solve systems of linear equations of the form $A\mathbf{x}=\mathbf{b}$, where $A$ is a tridiagonal Toeplitz matrix and $\mathbf{b}$ results from discretizing an analytic function, with a circuit complexity of $poly(\log(\kappa), 1/\sqrt{\epsilon}, \log(N))$, where $N$ denotes the number of equations, $\epsilon$ is the accuracy, and $\kappa$ the condition number. The \emph{repeat-until-success} algorithm has to be run $\mathcal{O}\left(\kappa/(1-\epsilon)\right)$ times to succeed, leveraging amplitude amplification, and sampled $\mathcal{O}(1/\epsilon^2)$ times. Thus, the algorithm achieves an exponential improvement with respect to $N$ over classical methods. In particular, we present efficient oracles for state preparation, Hamiltonian simulation and a set of observables together with the corresponding error and complexity analyses. As the main result of this work, we show how to use Richardson extrapolation to enhance Hamiltonian simulation, resulting in an implementation of Quantum Phase Estimation (QPE) within the algorithm with $1/\sqrt{\epsilon}$ circuit complexity instead of $1/\epsilon$ and which can be parallelized. Furthermore, we analyze necessary conditions for the overall algorithm to achieve an exponential speedup compared to classical methods. Our approach is not limited to the considered setting and can be applied to more general problems where Hamiltonian simulation is approximated via product formulae, although, our theoretical results would need to be extended accordingly. All the procedures presented are implemented with Qiskit and tested for small systems using classical simulation as well as using real quantum devices available through the IBM Quantum Experience., Comment: 25 pages, 15 figures

We develop a framework shedding light on common features of coupled variational formulations arising in electromagnetic scattering and acoustics. We show that spurious resonances haunting coupled domain-boundary formulations based on direct boundary integral equations of the first kind originate from the formal structure of their Calderon identities. Using this observation, the kernel of the coupled problem is characterized explicitly and we show that it completely vanishes under the exterior representation formula.

We couple the mixed variational problem for the generalized Hodge-Helmholtz or Hodge-Laplace equation posed on a bounded three-dimensional Lipschitz domain with the first-kind boundary integral equation arising from the latter when constant coefficients are assumed in the unbounded complement. Recently developed Calderon projectors for the relevant boundary integral operators are used to perform a symmetric coupling. We prove stability of the coupled problem away from resonant frequencies by establishing a generalized Garding inequality (T-coercivity). The resulting system of equations describes the scattering of monochromatic electromagnetic waves at a bounded inhomogeneous isotropic body possibly having a "rough" surface. The low-frequency robustness of the potential formulation of Maxwell's equations makes this model a promising starting point for Galerkin discretization., Comment: 31 pages

In this paper, we consider an acoustic wave transmission problem with mixed boundary conditions of Dirichlet, Neumann, and impedance type. The transmission interfaces may join the domain boundary in a general way independent of the location of the boundary conditions. We will derive a formulation as a \textit{direct}, \textit{space-time retarded boundary integral equation}, where both Cauchy data are kept as unknowns on the impedance part of the boundary. This requires the definition of single-trace spaces which incorporate homogeneous Dirichlet and Neumann conditions on the corresponding parts on the boundary. We prove the continuity and coercivity of the formulation by employing the technique of operational calculus in the Laplace domain., Comment: 15 pages, 1 figure

By introducing a shape manifold as a solution set to solve inverse obstacle scattering problems we allow the reconstruction of general, not necessarily star-shaped curves. The bending energy is used as a stabilizing term in Tikhonov regularization to gain independence of the parametrization. Moreover, we discuss how self-intersections can be avoided by penalization with the M\"obius energy and prove the regularizing property of our approach as well as convergence rates under variational source conditions. In the second part of the paper the discrete setting is introduced, and we describe a numerical method for finding the minimizer of the Tikhonov functional on a shape-manifold. Numerical examples demonstrate the feasibility of reconstructing non-star-shaped obstacles.

We introduce new boundary integral operators which are the exact inverses of the weakly singular and hypersingular operators for the Laplacian on flat disks. Moreover, we provide explicit closed forms for them and prove the continuity and ellipticity of their corresponding bilinear forms in the natural Sobolev trace spaces. This permit us to derive new Calder\'on-type identities that can provide the foundation for optimal operator preconditioning in Galerkin boundary element methods.

Trefftz methods are finite element-type schemes whose test and trial functions are (locally) solutions of the targeted differential equation. They are particularly popular for time-harmonic wave problems, as their trial spaces contain oscillating basis functions and may achieve better approximation properties than classical piecewise-polynomial spaces. We review the construction and properties of several Trefftz variational formulations developed for the Helmholtz equation, including least squares, discontinuous Galerkin, ultra weak variational formulation, variational theory of complex rays and wave based methods. The most common discrete Trefftz spaces used for this equation employ generalised harmonic polynomials (circular and spherical waves), plane and evanescent waves, fundamental solutions and multipoles as basis functions; we describe theoretical and computational aspects of these spaces, focusing in particular on their approximation properties. One of the most promising, but not yet well developed, features of Trefftz methods is the use of adaptivity in the choice of the propagation directions for the basis functions. The main difficulties encountered in the implementation are the assembly and the ill-conditioning of linear systems, we briefly survey some strategies that have been proposed to cope with these problems., Comment: 41 pages, 2 figures, to appear as a chapter in Springer Lecture Notes in Computational Science and Engineering. Differences from v1: added a few sentences in Sections 2.1, 2.2.2 and 2.3.1; inserted small corrections

We consider generalized linear transient convection-diffusion problems for differential forms on bounded domains in $\mathbb{R}^{n}$. These involve Lie derivatives with respect to a prescribed smooth vector field. We construct both new Eulerian and semi-Lagrangian approaches to the discretization of the Lie derivatives in the context of a Galerkin approximation based on discrete differential forms. Details of implementation are discussed as well as an application to the discretization of eddy current equations in moving media.

In this paper we prove the discrete compactness property for a wide class of p-version finite element approximations of non-elliptic variational eigenvalue problems in two and three space dimensions. In a very general framework, we find sufficient conditions for the p-version of a generalized discrete compactness property, which is formulated in the setting of discrete differential forms of any order on a d-dimensional polyhedral domain. One of the main tools for the analysis is a recently introduced smoothed Poincar\'e lifting operator [M. Costabel and A. McIntosh, On Bogovskii and regularized Poincar\'e integral operators for de Rham complexes on Lipschitz domains, Math. Z., (2010)]. For forms of order 1 our analysis shows that several widely used families of edge finite elements satisfy the discrete compactness property in p-version and hence provide convergent solutions to the Maxwell eigenvalue problem. In particular, N\'ed\'elec elements on triangles and tetrahedra (first and second kind) and on parallelograms and parallelepipeds (first kind) are covered by our theory.

We consider the variational formulation of the electric field integral equation (EFIE) on bounded polyhedral open or closed surfaces. We employ a conforming Galerkin discretization based on div-conforming Raviart-Thomas boundary elements (BEM) of locally variable polynomial degree on shape-regular surface meshes. We establish asymptotic quasi-optimality of Galerkin solutions on sufficiently fine meshes or for sufficiently high polynomial degree.

We consider H(curl)-elliptic variational problems on bounded Lipschitz polyhedra and their finite element Galerkin discretization by means of lowest order edge elements. We assume that the underlying tetrahedral mesh has been created by successive local mesh refinement, either by local uniform refinement with hanging nodes or bisection refinement. In this setting we develop a convergence theory for the the so-called local multigrid correction scheme with hybrid smoothing. We establish that its convergence rate is uniform with respect to the number of refinement steps. The proof relies on corresponding results for local multigrid in a H1-context along with local discrete Helmholtz-type decompositions of the edge element space.

We consider the first family of $\Hcurl$-conforming Ned\'el\'ec finite elements on tetrahedral meshes. Spectral approximation ($p$-version) is achieved by keeping the mesh fixed and raising the polynomial degree $p$ uniformly in all mesh cells. We prove that the associated subspaces of discretely weakly divergence free piecewise polynomial vector fields enjoy a long conjectured discrete compactness property as $p\to\infty$. This permits us to conclude asymptotic spectral correctness of spectral Galerkin finite element approximations of Maxwell eigenvalue problems., Comment: 23 pages

In computational relativity, critical behaviour near the black hole threshold has been studied numerically for several models in the last decade. In this paper we present a spatial Galerkin method, suitable for finding numerical solutions of the Einstein-Dirac equations in spherically symmetric spacetime (in polar/areal coordinates). The method features exact conservation of the total electric charge and allows for a spatial mesh adaption based on physical arclength. Numerical experiments confirm excellent robustness and convergence properties of our approach. Hence, this new algorithm is well suited for studying critical behaviour., Comment: 24 pages, 8 figures