Efficient quadrature for highly oscillatory integrals involving critical points

This paper based on the Levin collocation method and Levin-type method together with composite two-point Gauss–Legendre quadrature presents efficient quadrature for integral transformations of highly oscillatory functions with critical points. The effectiveness and accuracy of the quadrature are tested.     

Marching schemes for inverse acoustic scattering problems

Summary For the numerical solution of inverse Helmholtz problems the boundary value problem for a Helmholtz equation with spatially variable wave number has to be solved repeatedly. For large wave numbers this is a challenge. In the paper we reformulate the inverse problem as an initial value problem, and describe a marching scheme for the numerical computation that needs only n² log n operations on an n × n grid. We derive stability and error estimates for the marching scheme. We show that the marching solution is close to the low-pass filtered true solution. We present numerical examples that demonstrate the efficacy of the marching scheme.     

On the evaluation of Cauchy principal value integrals of oscillatory functions

The problem of the numerical evaluation of Cauchy principal value integrals of oscillatory functions , where −1<τ<1, has been discussed. Based on analytic continuation, if f is analytic in a sufficiently large complex region G containing [−1, 1], the integrals can be transformed into the problems of integrating two integrals on [0,+∞) with the integrand that does not oscillate, and that decays exponentially fast, which can be efficiently computed by using the Gauss-Laguerre quadrature rule. The validity of the method has been demonstrated in the provision of two numerical experiments and their results.     

Controllability method for acoustic scattering with spectral elements

We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improvements in efficiency due to the higher order spectral elements. For a given accuracy, the controllability technique with spectral element method requires fewer computational operations than with conventional finite element method. In addition, by using higher order polynomial basis the influence of the pollution effect is reduced.     

Wavelet moment method for the Cauchy problem for the Helmholtz equation

The paper is concerned with the problem of reconstruction of acoustic or electromagnetic field from inexact data given on an open part of the boundary of a given domain. A regularization concept is presented for the moment problem that is equivalent to a Cauchy problem for the Helmholtz equation. A method of regularization by projection with application of the Meyer wavelet subspaces is introduced and analyzed. The derived formula, describing the projection level in terms of the error bound of the inexact Cauchy data, allows us to prove the convergence and stability of the method.     

Numerical evaluation for Cauchy type singular integrals on the interval

The singular integral (SI) with the Cauchy kernel is considered. New quadrature formulas (QFs) based on the modification of discrete vortex method to approximate SI are constructed. Convergence of QFs and error bounds are shown in the classes of functions H^α([−1,1]) and C¹([−1,1]). Numerical examples are shown to validate the QFs constructed.     

Complex Gaussian quadrature of oscillatory integrals

We construct and analyze Gauss-type quadrature rules with complex- valued nodes and weights to approximate oscillatory integrals with stationary points of high order. The method is based on substituting the original interval of integration by a set of contours in the complex plane, corresponding to the paths of steepest descent. Each of these line integrals shows an exponentially decaying behaviour, suitable for the application of Gaussian rules with non-standard weight functions. The results differ from those in previous research in the sense that the constructed rules are asymptotically optimal, i.e., among all known methods for oscillatory integrals they deliver the highest possible asymptotic order of convergence, relative to the required number of evaluations of the integrand.     

On rates of acceleration of extrapolation methods for oscillatory infinite integrals

The purpose of this work is to complement and expand our knowledge of the convergence theory of some extrapolation methods for the accurate computation of oscillatory infinite integrals. Specifically, we analyze in detail the convergence properties of theW- and -transformations of the author as they are applied to three integrals, all with totally different behavior at infinity. The results of the analysis suggest different convergence and acceleration of convergence behavior for the above mentioned transformations on the different integrals, and they improve considerably those that can be obtained from the existing convergence theories.     

Integration matrix based on arbitrary grids with a preconditioner for pseudospectral method

A matrix representation of integration for arbitrary grids is introduced. Suitable results are then obtained to be used along with differentiation matrix preconditioner to implement Pseudospectral method on integro-differential equations using arbitrary grids. Numerical examples are given to clarify the efficiency of the new method.     

A collocation method for high-frequency scattering by convex polygons

We consider the problem of scattering of a time-harmonic acoustic incident plane wave by a sound soft convex polygon. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the computational cost required to achieve a prescribed level of accuracy grows linearly with respect to the frequency of the incident wave. Recently Chandler–Wilde and Langdon proposed a novel Galerkin boundary element method for this problem for which, by incorporating the products of plane wave basis functions with piecewise polynomials supported on a graded mesh into the approximation space, they were able to demonstrate that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency. Here we propose a related collocation method, using the same approximation space, for which we demonstrate via numerical experiments a convergence rate identical to that achieved with the Galerkin scheme, but with a substantially reduced computational cost.     

Worst case scenario analysis for elliptic problems with uncertainty

This work studies linear elliptic problems under uncertainty. The major emphasis is on the deterministic treatment of such uncertainty. In particular, this work uses the Worst Scenario approach for the characterization of uncertainty on functional outputs (quantities of physical interest). Assuming that the input data belong to a given functional set, eventually infinitely dimensional, this work proposes numerical methods to approximate the corresponding uncertainty intervals for the quantities of interest. Numerical experiments illustrate the performance of the proposed methodology.     

Pseudospectral method for quadrilaterals

In this paper, we investigate the pseudospectral method on quadrilaterals. Some results on Legendre-Gauss-type interpolation are established, which play important roles in the pseudospectral method for partial differential equations defined on quadrilaterals. As examples of applications, we propose pseudospectral methods for two model problems and prove their spectral accuracy in space. Numerical results demonstrate the efficiency of the suggested algorithms. The approximation results and techniques developed in this paper are also applicable to other problems defined on quadrilaterals.     

Feynman path integrals as infinite-dimensional oscillatory integrals: Some new developments

We discuss a basic mathematical approach to Feynman path integrals as infinite-dimensional oscillatory integrals. We present new results on asymptotics of such integrals which exploit recently developed approximation techniques via finite dimensional oscillatory integrals. Applications are also given, namely to the study of the trace of the time evolution operator in quantum mechanics and to the interpretation of Gutzwiller's trace formula as a leading term in an asymptotic expansion around classical periodic orbits.The second named author is an Alexander von Humboldt Stiftung fellow.     

Nonlinear elliptic equations with measure data

In this paper we prove the existence of solutions of nonlinear equations of the type-div(a(x, u, Du)+H(x, u, Du)=f, wherea andH are Caratheodory functions andf is a bounded Radon measure. We remark that the operator can be not coercive. We give also some regularity results.     

The superconvergence of the Newton-Cotes rule for Cauchy principal value integrals

We consider the general (composite) Newton-Cotes method for the computation of Cauchy principal value integrals and focus on its pointwise superconvergence phenomenon, which means that the rate of convergence of the Newton-Cotes quadrature rule is higher than what is globally possible when the singular point coincides with some a priori known point. The necessary and sufficient conditions satisfied by the superconvergence point are given. Moreover, the superconvergence estimate is obtained and the properties of the superconvergence points are investigated. Finally, some numerical examples are provided to validate the theoretical results.     

Existence of solutions to quasi-linear elliptic problems with generalized Robin boundary conditions

We characterize the L^q-solvability of a class of quasi-linear elliptic equations involving the p-Laplace operator with generalized nonlinear Robin type boundary conditions on bad domains. Some uniqueness results are also given.     

Application of Cauchy integrals and singular integral equations in scattered data problems

A new approach is suggested for constructing “smooth” surfaces which contain discontinuities in functional values or derivatives at prescribed locations. The approach is based on solving singular integral equations with Cauchy-type kernels. It is applicable in the interpolation or approximation of scattered data. We investigate, in particular, several two-dimensional biharmonic and triharmonic problems which have smooth solutions except on given line segments, across which different types of discontinuities occur. We also discuss some issues concerning the application of the approach in practical problems.     

Convergence rates of iterated Tikhonov regularized solutions of nonlinear III — posed problems

Summary In this paper we investigate iterated Tikhonov regularization for the solution of nonlinear ill-posed problems. In the case of linear ill-posed problems it is well-known that (under appropriate assumptions) then-th iterated regularized solutions can converge likeO(²² /(2n+1)), where denotes the noise level of the data perturbation. We give conditions that guarantee this convergence rate also for nonlinear ill-posed problems, and motivate these conditions by the mapping degree. The results are derived by a comparison of the iterated regularized solutions of the nonlinear problem with the iterated regularized solutions of its linearization. Numerical examples are presented.Supported by the Austrian Fonds zur Förderung der wissenschaftlichen Forschung,project P-7869 PHY, and by the Christian Doppler Society     

Asymptotic expansions of Mellin convolution integrals: An oscillatory case

In a recent paper [J.L. López, Asymptotic expansions of Mellin convolution integrals, SIAM Rev. 50 (2) (2008) 275-293], we have presented a new, very general and simple method for deriving asymptotic expansions of for small x. It contains Watson’s Lemma and other classical methods, Mellin transform techniques, McClure and Wong’s distributional approach and the method of analytic continuation used in this approach as particular cases. In this paper we generalize that idea to the case of oscillatory kernels, that is, to integrals of the form , with c∈R, and we give a method as simple as the one given in the above cited reference for the case c=0. We show that McClure and Wong’s distributional approach for oscillatory kernels and the summability method for oscillatory integrals are particular cases of this method. Some examples are given as illustration.