Boundary integral equation methods for the solution of scattering and transmission 2D elastodynamic problems

We introduce and analyse various regularized combined field integral equations (CFIER) formulations of time-harmonic Navier equations in media with piece-wise constant material properties. These formulations can be derived systematically starting from suitable coercive approximations of Dirichlet-to-Neumann operators (DtN), and we present a periodic pseudodifferential calculus framework within which the well posedness of CIER formulations can be established. We also use the DtN approximations to derive and analyse OS methods for the solution of elastodynamics transmission problems. The pseudodifferential calculus we develop in this paper relies on careful singularity splittings of the kernels of Navier boundary integral operators, which is also the basis of high-order Nystrom quadratures for their discretizations. Based on these high-order discretizations we investigate the rate of convergence of iterative solvers applied to CFIER and OS formulations of scattering and transmission problems. We present a variety of numerical results that illustrate that the CFIER methodology leads to important computational savings over the classical CFIE one, whenever iterative solvers are used for the solution of the ensuing discretized boundary integral equations. Finally, we show that the OS methods are competitive in the high-frequency high-contrast regime.Catalin Turc gratefully acknowledges support from National Science Foundation (NSF) through contract DMS-1614270 and DMS-1908602

Formulation and Solution of Electromagnetic Integral Equations Using Constraint-Based Helmholtz Decompositions

This dissertation develops surface integral equations using constraint-based Helmholtz decompositions for electromagnetic modeling. This new approach is applied to the electric field integral equation (EFIE), and it incorporates a Helmholtz decomposition (HD) of the current. For this reason, the new formulation is referred to as the EFIE-hd. The HD of the current is accomplished herein via appropriate surface integral constraints, and leads to a stable linear system. This strategy provides accurate solutions for the electric and magnetic fields at both high and low frequencies, it allows for the use of a locally corrected Nyström (LCN) discretization method for the resulting formulation, it is compatible with the local global solution framework, and it can be used with non-conformal meshes. To address large-scale and complex electromagnetic problems, an overlapped localizing local-global (OL-LOGOS) factorization is used to factorize the system matrix obtained from an LCN discretization of the augmented EFIE (AEFIE). The OL-LOGOS algorithm provides good asymptotic performance and error control when used with the AEFIE. This application is used to demonstrate the importance of using a well-conditioned formulation to obtain efficient performance from the factorization algorithm

Boundary integral equation methods for the elastic and thermoelastic waves in three dimensions

In this paper, we consider the boundary integral equation (BIE) method for solving the exterior Neumann boundary value problems of elastic and thermoelastic waves in three dimensions based on the Fredholm integral equations of the first kind. The innovative contribution of this work lies in the proposal of the new regularized formulations for the hyper-singular boundary integral operators (BIO) associated with the time-harmonic elastic and thermoelastic wave equations. With the help of the new regularized formulations, we only need to compute the integrals with weak singularities at most in the corresponding variational forms of the boundary integral equations. The accuracy of the regularized formulations is demonstrated through numerical examples using the Galerkin boundary element method (BEM).Comment: 24 pages, 6 figure

A LOCALLY CORRECTED NYSTRM METHOD FOR SURFACE INTEGRAL EQUATIONS: AN OBJECT ORIENTED APPROACH

Classically, researchers in Computational Physics and specifically in Computational Electromagnetics have sought to find numerical solutions to complex physical problems. Several techniques have been developed to accomplish such tasks, each of which having advantages over their counterparts. Typically, each solution method has been developed separately despite having numerous commonalities with other methods. This fact motivates a unified software tool to house each solution method to avoid duplicating previous efforts. Subsequently, these solution methods can be used alone or in conjunction with one another in a straightforward manner. The aforementioned goals can be accomplished by using an Object Oriented software approach. Thus, the goal of the presented research was to incorporate a specific solution technique, an Integral Equation Nystrm method, into a general, Object Oriented software framework

Weighted integral solvers for elastic scattering by open arcs in two dimensions

We present new methodologies for the numerical solution of problems of elastic scattering by open arcs in two dimensions. The algorithms utilize weighted versions of the classical elastic integral operators associated with Dirichlet and Neumann boundary conditions, where the integral weight accounts for (and regularizes) the singularity of the integral‐equation solutions at the open‐arc endpoints. Crucially, the method also incorporates a certain “open‐arc elastic Calderón relation” introduced in this paper, whose validity is demonstrated on the basis of numerical experiments, but whose rigorous mathematical proof is left for future work. (In fact, the aforementioned open‐arc elastic Calderón relation generalizes a corresponding elastic Calderón relation for closed surfaces, which is also introduced in this paper, and for which a rigorous proof is included.) Using the open‐surface Calderón relation in conjunction with spectrally accurate quadrature rules and the Krylov‐subspace linear algebra solver GMRES, the proposed overall open‐arc elastic solver produces results of high accuracy in small number of iterations, for both low and high frequencies. A variety of numerical examples in this paper demonstrate the accuracy and efficiency of the proposed methodology