Source author record

Victor Dominguez

Victor Dominguez 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

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

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.

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.

preprint2019arXiv

An overlapping decomposition framework for wave propagation in heterogeneous and unbounded media: Formulation, analysis, algorithm, and simulation

A natural medium for wave propagation comprises a coupled bounded heterogeneous region and an unbounded homogeneous free-space. Frequency-domain wave propagation models in the medium, such as the variable coefficient Helmholtz equation, include a faraway decay radiation condition (RC). It is desirable to develop algorithms that incorporate the full physics of the heterogeneous and unbounded medium wave propagation model, and avoid an approximation of the RC. In this work we first present and analyze an overlapping decomposition framework that is equivalent to the full-space heterogeneous-homogenous continuous model, governed by the Helmholtz equation with a spatially dependent refractive index and the RC. Our novel overlapping framework allows the user to choose two free boundaries, and gain the advantage of applying established high-order finite and boundary element methods (FEM and BEM) to simulate an equivalent coupled model. The coupled model comprises auxiliary interior bounded heterogeneous and exterior unbounded homogeneous Helmholtz problems. A smooth boundary can be chosen for simulating the exterior problem using a spectrally accurate BEM, and a simple boundary can be used to setup a high-order FEM for the interior problem. Thanks to the spectral accuracy of the exterior computational model, the resulting coupled system in the overlapping region is relatively very small. Using the decomposed equivalent framework, we develop a novel overlapping FEM-BEM algorithm for simulating the acoustic or electromagnetic wave propagation in two dimensions. Our FEM-BEM algorithm for the full-space model incorporates the RC exactly. Numerical experiments demonstrate the efficiency of the FEM-BEM approach for simulating smooth and non-smooth wave fields, with the latter induced by a complex heterogeneous medium and a discontinuous refractive index.

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.

preprint2013arXiv

Filon-Clenshaw-Curtis rules for a class of highly-oscillatory integrals with logarithmic singularities

In this work we propose and analyse a numerical method for computing a family of highly oscillatory integrals with logarithmic singularities. For these quadrature rules we derive error estimates in terms of $N$, the number of nodes, $k$ the rate of oscillations and a Sobolev-like regularity of the function. We prove that that the method is not only robust but the error even decreases, for fixed $N$, as $k$ increases. Practical issues about the implementation of the rule are also covered in this paper by: (a) writing down ready-to-implement algorithms; (b) analysing the numerical stability of the computations and (c) estimating the overall computational cost. We finish by showing some numerical experiments which illustrate the theoretical results presented in this paper.

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.

preprint2012arXiv

A fully discrete Calderon Calculus for two dimensional time harmonic waves

In this paper, we present a fully discretized Calderón Calculus for the two dimensional Helmholtz equation. This full discretization can be understood as highly non-conforming Petrov-Galerkin methods, based on two staggered grids of mesh size $h$, Dirac delta distributions substituting acoustic charge densities and piecewise constant functions for approximating acoustic dipole densities. The resulting numerical schemes from this calculus are all of order $h^2$ provided that the continuous equations are well posed. We finish by presenting some numerical experiments illustrating the performance of this discrete calculus.

preprint2012arXiv

A Nystrom method for the two dimensional Helmholtz hypersingular equation

In this paper we propose and analyze a class of simple Nyström discretizations of the hypersingular integral equation for the Helmholtz problem on domains of the plane with smooth parametrizable boundary. The method depends on a parameter (related to the staggering of two underlying grids) and we show that two choices of this parameter produce convergent methods of order two, while all other stable methods provide methods of order one. Convergence is shown for the density (in uniform norm) and for the potential postprocessing of the solution. Some numerical experiments are given to illustrate the performance of the method.

preprint2012arXiv

Convergence analysis of a high-order Nystrom integral-equation method for surface scattering problems

In this paper we present a convergence analysis for the Nystrom method proposed in [Jour. Comput. Phys. 169 pp. 2921-2934, 2001] for the solution of the combined boundary integral equation formulations of sound-soft acoustic scattering problems in three-dimensional space. This fast and efficient scheme combines FFT techniques and a polar change of variables that cancels out the kernel singularity. We establish the stability of the algorithms in the $L^2$ norm and we derive convergence estimates in both the $L^2$ and $L^\infty$ norms. In particular, our analysis establishes theoretically the previously observed super-algebraic convergence of the method in cases in which the right-hand side is smooth.