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

Catalin Turc

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

ResearcherUnclaimed source record
math.NAmath.APNumerical Analysismath-phmath.MPphysics.comp-ph

Catalog footprint

What is connected

10works
6topics
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

10 published item(s)

preprint2022arXiv

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

We introduce and analyze 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 analyze Optimized Schwarz (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 Nyström 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nystrom discretizations of boundary integral equations for the solution of 2D elastic scattering problems

We present three high-order Nystrom discretization strategies of various boundary integral equation formulations of the impenetrable time-harmonic Navier equations in two dimensions. One class of such formulations is based on the four classical Boundary Integral Operators (BIOs) associated with the Green's function of the Navier operator. We consider two types of Nystrom discretizations of these operators, one that relies on Kussmaul-Martensen logarithmic splittings and the other on Alpert quadratures. In addition, we consider an alternative formulation of Navier scattering problems based on Helmholtz decompositions of the elastic fields, which can be solved via a system of boundary integral equations that feature integral operators associated with the Helmholtz equation. Owing to the fact that some of the BIOs that are featured in those formulations are non-standard, we use Quadrature by Expansion (QBX) methods for their high order Nystrom discretization. Alternatively, we use Maue integration by parts techniques to recast those non-standard operators in terms of single and double layer Helmholtz BIOs whose Nystrom discretizations is amenable to the Kussmaul-Martensen methodology. We present a variety of numerical results concerning the high order accuracy that our Nystrom discretization elastic scattering solvers achieve for both smooth and Lipschitz boundaries. We also present extensive comparisons regarding the iterative behavior of solvers based on different integral equations in the high frequency regime. Finally, we illustrate how some of the Nystrom discretizations we considered can be incorporated seamlessly into the Convolution Quadrature (CQ) methodology to deliver high-order solutions of the time domain elastic scattering problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Schur complement Domain Decomposition Methods for the solution of multiple scattering problems

We present a Schur complement Domain Decomposition (DD) algorithm for the solution of frequency domain multiple scattering problems. Just as in the classical DD methods we (1) enclose the ensemble of scatterers in a domain bounded by an artificial boundary, (2) we subdivide this domain into a collection of nonoverlapping subdomains so that the boundaries of the subdomains do not intersect any of the scatterers, and (3) we connect the solutions of the subproblems via Robin boundary conditions matching on the common interfaces between subdomains. We use subdomain Robin-to-Robin maps to recast the DD problem as a sparse linear system whose unknown consists of Robin data on the interfaces between subdomains---two unknowns per interface. The Robin-to-Robin maps are computed in terms of well-conditioned boundary integral operators. Unlike classical DD, we do not reformulate the Domain Decomposition problem in the form a fixed point iteration, but rather we solve the ensuing linear system by Gaussian elimination of the unknowns corresponding to inner interfaces between subdomains via Schur complements. Once all the unknowns corresponding to inner subdomains interfaces have been eliminated, we solve a much smaller linear system involving unknowns on the inner and outer artificial boundary. We present numerical evidence that our Schur complement DD algorithm can produce accurate solutions of very large multiple scattering problems that are out of reach for other existing approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Superalgebraically Convergent Smoothly-Windowed Lattice Sums for Doubly Periodic Green Functions in Three-Dimensional Space

This paper, Part I in a two-part series, presents (i) A simple and highly efficient algorithm for evaluation of quasi-periodic Green functions, as well as (ii) An associated boundary-integral equation method for the numerical solution of problems of scattering of waves by doubly periodic arrays of scatterers in three-dimensional space. Except for certain "Wood frequencies" at which the quasi-periodic Green function ceases to exist, the proposed approach, which is based on use of smooth windowing functions, gives rise to lattice sums which converge superalgebraically fast--that is, faster than any power of the number of terms used--in sharp contrast with the extremely slow convergence exhibited by the corresponding sums in absence of smooth windowing. (The Wood-frequency problem is treated in Part II.) A proof presented in this paper establishes rigorously the superalgebraic convergence of the windowed lattice sums. A variety of numerical results demonstrate the practical efficiency of the proposed approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Well-conditioned boundary integral equation formulations and Nyström discretizations for the solution of Helmholtz problems with impedance boundary conditions in two-dimensional Lipschitz domains

We present a regularization strategy that leads to well-conditioned boundary integral equation formulations of Helmholtz equations with impedance boundary conditions in two-dimensional Lipschitz domains. We consider both the case of classical impedance boundary conditions, as well as the case of transmission impedance conditions wherein the impedances are certain coercive operators. The latter type of problems is instrumental in the speed up of the convergence of Domain Decomposition Methods for Helmholtz problems. Our regularized formulations use as unknowns the Dirichlet traces of the solution on the boundary of the domain. Taking advantage of the increased regularity of the unknowns in our formulations, we show through a variety of numerical results that a graded-mesh based Nyström discretization of these regularized formulations leads to efficient and accurate solutions of interior and exterior Helmholtz problems with impedance boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Well-posed boundary integral equation formulations and Nyström discretizations for the solution of Helmholtz transmission problems in two-dimensional Lipschitz domains

We present a comparison between the performance of solvers based on Nyström discretizations of several well-posed boundary integral equation formulations of Helmholtz transmission problems in two-dimensional Lipschitz domains. Specifically, we focus on the following four classes of boundary integral formulations of Helmholtz transmission problems (1) the classical first kind integral equations for transmission problems, (2) the classical second kind integral equations for transmission problems, (3) the {\em single} integral equation formulations, and (4) certain direct counterparts of recently introduced Generalized Combined Source Integral Equations. The former two formulations were the only formulations whose well-posedness in Lipschitz domains was rigorously established. We establish the well-posedness of the latter two formulations in appropriate functional spaces of boundary traces of solutions of transmission Helmholtz problems in Lipschitz domains. We give ample numerical evidence that Nyström solvers based on formulations (3) and (4) are computationally more advantageous than solvers based on the classical formulations (1) and (2), especially in the case of high-contrast transmission problems at high frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Windowed Green Function method for layered-media scattering

This paper introduces a new Windowed Green Function (WGF) method for the numerical integral-equation solution of problems of electromagnetic scattering by obstacles in presence of dielectric or conducting half-planes. The WGF method, which is based on use of smooth windowing functions and integral kernels that can be expressed directly in terms of the free-space Green function, does not require evaluation of expensive Sommerfeld integrals. The proposed approach is fast, accurate, flexible and easy to implement. In particular, straightforward modifications of existing (accelerated or unaccelerated) solvers suffice to incorporate the WGF capability. The mathematical basis of the method is simple: the method relies on a certain integral equation posed on the union of the boundary of the obstacle and a small flat section of the interface between the penetrable media. Numerical experiments demonstrate that both the near- and far-field errors resulting from the proposed approach decrease faster than any negative power of the window size. In the examples considered in this paper the proposed method is up to thousands of times faster, for a given accuracy, than a corresponding method based on the layer-Green-function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Well-conditioned boundary integral equation formulations for the solution of high-frequency electromagnetic scattering problems

We present several versions of Regularized Combined Field Integral Equation (CFIER) formulations for the solution of three dimensional frequency domain electromagnetic scattering problems with Perfectly Electric Conducting (PEC) boundary conditions. Just as in the Combined Field Integral Equations (CFIE), we seek the scattered fields in the form of a combined magnetic and electric dipole layer potentials that involves a composition of the latter type of boundary layers with regularizing operators. The regularizing operators are of two types: (1) modified versions of electric field integral operators with complex wavenumbers, and (2) principal symbols of those operators in the sense of pseudodifferential operators. We show that the boundary integral operators that enter these CFIER formulations are Fredholm of the second kind, and invertible with bounded inverses in the classical trace spaces of electromagnetic scattering problems. We present a spectral analysis of CFIER operators with regularizing operators that have purely imaginary wavenumbers for spherical geometries. Under certain assumptions on the coupling constants and the absolute values of the imaginary wavenumbers of the regularizing operators, we show that the ensuing CFIER operators are coercive for spherical geometries. These properties allow us to derive wavenumber explicit bounds on the condition numbers of certain CFIER operators that have been proposed in the literature. When regularizing operators with complex wavenumbers with non-zero real parts are used, we show numerical evidence that those complex wavenumbers can be selected in a manner that leads to CFIER formulations whose condition numbers can be bounded independently of frequency for spherical geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Integral equations requiring small numbers of Krylov-subspace iterations for two-dimensional penetrable scattering problems

This paper presents a class of boundary integral equations for the solution of problems of electromagnetic and acoustic scattering by two dimensional homogeneous penetrable scatterers with smooth boundaries. The new integral equations, which, as is established in this paper, are uniquely solvable Fredholm equations of the second kind, result from representations of fields as combinations of single and double layer potentials acting on appropriately chosen regularizing operators. As demonstrated in this text by means of a variety of numerical examples (that resulted from a high-order Nystrom computational implementation of the new equations), these "regularized combined equations" can give rise to important reductions in computational costs, for a given accuracy, over those resulting from previous boundary integral formulations for transmission problems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Regularized combined field integral equations for acoustic transmission problems

We present a new class of well conditioned integral equations for the solution of two and three dimensional scattering problems by homogeneous penetrable scatterers. Our novel boundary integral equations result from suitable representations of the fields inside and outside the scatterer as combinations of single and double layer potentials acting on suitably defined regularizing operators. The regularizing operators are constructed to be suitable approximations of the admittance operators that map the transmission boundary conditions to the exterior and respectively interior Cauchy data on the interface between the media. The latter operators can be expressed in terms of Dirichlet-to-Neumann operators. We refer to these regularized boundary integral equations as Generalized Combined Source Integral Equations (GCSIE). The ensuing GCSIE are shown to be integral equations of the second kind in the case when the interface of material discontinuity is a smooth curve in two dimensions and a smooth surface in three dimensions.

0 citations0 reviewsNo code linkedOpen access