Source author record

V. Dominguez

V. 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

7works
1topics
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

7 published item(s)

preprint2016arXiv

High order Nyström methods for transmission problems for Helmholtz equation

We present superalgebraic compatible Nyström discretizations for the four Helmholtz boundary operators of Calderón's calculus on smooth closed curves in 2D. These discretizations are based on appropriate splitting of the kernels combined with very accurate product-quadrature rules for the different singularities that such kernels present. A Fourier based analysis shows that the four discrete operators converge to the continuous ones in appropriate Sobolev norms. This proves that Nyström discretizations of many popular integral equation formulations for Helmholtz equations are stable and convergent. The convergence is actually superalgebraic for smooth solutions.

preprint2016arXiv

Sobolev estimates for constructive uniform-grid FFT interpolatory approximations of spherical functions

The fast Fourier transform (FFT) based matrix-free ansatz interpolatory approximations of periodic functions are fundamental for efficient realization in several applications.In this work we design, analyze, and implement similar constructive interpolatory approximations of spherical functions, using samples of the unknown functions at the poles and at the uniform spherical-polar grid locations. The spherical matrix-free interpolation operator range space consists of a selective subspace of two dimensional trigonometric polynomials which are rich enough to contain all spherical polynomials of degree less than $N$. The spherical interpolatory approximation is efficiently constructed by applying the FFT techniques with only ${\mathcal{O}}(N^2 \log N)$ complexity. We describe the construction details using the FFT operators and provide complete convergence analysis of the interpolatory approximation in the Sobolev space framework. We prove that the rate of spectrally accurate convergence of the interpolatory approximations in Sobolev norms (of order zero and one) are similar (up to a log term) to that of the best approximation in the finite dimensional ansatz space. Efficient interpolatory quadratures on the sphere are important for several applications including radiation transport and wave propagation computer models. We use our matrix-free interpolatory approximations to construct robust FFT-based quadrature rules for a wide class of non-, mildly-, and strongly-oscillatory integrands on the sphere. We provide numerical experiments to demonstrate fast evaluation of the algorithm and various theoretical results presented in the article.

preprint2014arXiv

High-order Nyström discretizations for the solution of integral equation formulations of two-dimensional Helmholtz transmission problems

We present and analyze fully discrete Nyström methods for the solution of three classes of well conditioned boundary integral equations for the solution of two dimensional scattering problems by homogeneous dielectric scatterers. Specifically, we perform the stability analysis of Nyström discretizations of (1) the classical second kind integral equations for transmission problems [KressRoach, 1978], (2) the single integral equation formulations [Kleinman & Martin,1988], and (3) recently introduced Generalized Combined Source Integral Equations [Boubendir,Bruno, Levadoux and. Turc,2013]. The Nyström method that we use for the discretization of the various integral equations under consideration are based on global trigonometric approximations, splitting of the kernels of integral operators into singular and smooth components, and explicit quadratures of products of singular parts (logarithms) and trigonometric polynomials. The discretization of the integral equations (2) and (3) above requires special care as these formulations feature compositions of boundary integral operators that are pseudodifferential operators of positive and negative orders respectively. We deal with these compositions through Calderon's calculus and we establish the convergence of fully discrete Nyström methods in appropriate Sobolev spaces which implies pointwise convergence of the discrete solutions. In the case of analytic boundaries, we establish superalgebraic convergence of the method.

preprint2013arXiv

A Nystrom flavored Calderón Calculus of order three for two dimensional waves

In this paper we present and test a full discretization of all elements of the Calderón Calculus (layer potentials and integral operators) for the Helmholtz equation in smooth closed curves in the plane. The resulting integral equations provide approximations of order three for all variables involved. Test are shown for a wide array of direct, indirect and combined field integral equation at fixed frequency and for a Convolution Quadrature based approximation in the time domain.

preprint2013arXiv

On the discontinuous Galerkin method for solving boundary value problems for the Helmholtz equation: A priori and a posteriori error analyses

We apply the local discontinuous Galerkin (LDG for short) method to solve a mixed boundary value problems for the Helmholtz equation in bounded polygonal domain in 2D. Under some assumptions on regularity of the solution of an adjoint problem, we prove that: (a) the corresponding indefinite discrete scheme is well posed; (b) there is convergence with the expected convergence rates as long as the meshsize h is small enough. We give precise information on how small h has to be in terms of the size of the wavenumber and its distance to the set of eigenvalues for the same boundary value problem for the Laplacian. We also present an a posteriori error estimator showing both the reliability and efficiency of the estimator complemented with detailed information on the dependence of the constants on the wavenumber. We finish presenting extensive numerical experiments which illustrate the theoretical results proven in this paper and suggest that stability and convergence may occur under less restrictive assumptions than those taken in the present work.

preprint2012arXiv

Filon-Clenshaw-Curtis rules for highly-oscillatory integrals with algebraic singularities and stationary points

In this paper we propose and analyse composite Filon-Clenshaw-Curtis quadrature rules for integrals of the form $I_{k}^{[a,b]}(f,g) := \int_a^b f(x) \exp(\mathrm{i}kg(x)) \rd x $, where $k \geq 0$, $f$ may have integrable singularities and $g$ may have stationary points. Our composite rule is defined on a mesh with $M$ subintervals and requires $MN+1$ evaluations of $f$. It satisfies an error estimate of the form $C_N k^{-r} M^{-N-1 + r}$, where $r$ is determined by the strength of any singularity in $f$ and the order of any stationary points in $g$ and $C_N$ is a constant which is independent of $k$ and $M$, but depends on $N$. The regularity requirements on $f$ and $g$ are explicit in the error estimates. For fixed $k$, the rate of convergence of the rule as $M \rightarrow \infty$ is the same as would be obtained if $f$ was smooth. Moreover, the quadrature error decays at least as fast as $k \rightarrow \infty$ as does the original integral $I_{k}^{[a,b]}(f,g)$. For the case of nonlinear oscillators $g$, the algorithm requires the evaluation of $g^{-1}$ at non-stationary points. Numerical results demonstrate the sharpness of the theory. An application to the implementation of boundary integral methods for the high-frequency Helmholtz equation is given.

preprint2012arXiv

Interpolation and cubature approximations and analysis for a class of wideband integrals on the sphere

We propose, analyze, and implement interpolatory approximations and Filon-type cubature for efficient and accurate evaluation of a class of wideband generalized Fourier integrals on the sphere. The analysis includes derivation of (i) optimal order Sobolev norm error estimates for an explicit discrete Fourier transform type interpolatory approximation of spherical functions; and (ii) a wavenumber explicit error estimate of the order $\mathcal{O}(κ^{-\ell} N^{-r_\ell})$, for $\ell = 0, 1, 2$, where $κ$ is the wavenumber, $N$ is the number of interpolation/cubature points on the sphere and $r_\ell$ depends on the smoothness of the integrand. Consequently, the cubature is robust for wideband (from very low to very high) frequencies and very efficient for highly-oscillatory integrals because the quality of the high-order approximation (with respect to quadrature points) is further improved as the wavenumber increases. This property is a marked advantage compared to standard cubature that require at least ten points per wavelength per dimension and methods for which asymptotic convergence is known only with respect to the wavenumber subject to stable of computation of quadrature weights. Numerical results in this article demonstrate the optimal order accuracy of the interpolatory approximations and the wideband cubature.