Numerical analysis of the spectrum for the highly oscillatory integral equation with weak singularity☆
Introduction
We focus on the spectra of the integral equation where the kernel function is assumed to be even function, such as in which is an oscillatory parameter. Once , the integral equation becomes one with weak singularity. The integral operator has arisen in the singular integrals of the lower dimensional varieties, singular Radon transform [19], [20], and the application of the boundary integral equation of acoustic, electromagnetic wave scattering [3], [4]. In particular, its numerical computation is also proposed as an open problem in Section 6 of Brunner et al. [1]. The integral operator of (1.1) with (1.2) or (1.3) in is compact which inspires the computational way of finite projection operator approximation [14], [21]. The compactness of the operators also indicates that the discrete spectrum of the operator accumulate to the origin.
Highly oscillatory integral arises in many engineering and numerical analysis fields, for instance, acoustic scattering, magnetics, quantum mechanics and so on. In the last decades, much efforts are devoted to the quadrature of the highly oscillatory integral. The asymptotic expansion method is fundamental to the numerical computations but only working for the sufficiently large oscillatory parameter [5]. Filon type method is applied extensively with high accuracy when the moment function can be calculated exactly [9], [12]. Avoiding the computation of moment, Levin method is constructed for the oscillatory integral without the stationary point [15], [16]. Numerical steepest descent method turns the oscillation integral to the non-oscillatory one by changing the integral path for analytic function [11]. The oscillatory integral with Bessel oscillation is explored in [22]. Furthermore, the asymptotic expansion and Filon type methods have been extended to the highly oscillatory singular integral recently [8].
Weak singularity and high oscillations hinder the accurately fast numerical computation for the integral equation. Since the analytic solution is difficult to obtain, numerical computation for (1.1) with (1.2) or (1.3) is a better option. The theoretical analysis of the numerical approximation for a compact operator has been studied in [18]. In [2], several numerical computation methods are presented which are the first step to the numerical computation for the Fox–Li operator. The asymptotic analysis of the spectra of the integral operator with the absolute oscillator out of the singularity is formulated in [1]. In [7], one finite section method with the modified Fourier basis has been presented to analyse the spectra for the Fox–Li integral operator , where and . For the integral equation (1.1) with (1.2), the modified Fourier basis is firstly used to approximate the unknown in which the entries of the matrix are represented in the hypergeometric functions [10]. Here, we are interested in introducing the finite section method to analyse the numerical spectrum (1.1), especially the influence of the singular parameter for the distribution of the spectrum. For the absolute oscillator integral equation with power-law or logarithmic singularities, we construct a unified numerical method which can be implemented easily by considering the affection of two oscillatory factors , and .
In this paper, a finite section method based on the modified Fourier basis is applied to compute the spectrum of the highly oscillatory integral equation with algebraic singularity. The accurate approximation of the non-periodic function facilitates the approximation of the modified Fourier basis for the integral equation [13]. During transforming the problem of the operator spectrum to the algebraic problem, one infinite coefficient matrix is derived whose entries are in the form of bivariate highly oscillatory singular integral. By some algebra computation, the double integral is simplified to be the explicit expression. The asymptotic properties of the entries are proved to construct the truncation strategy. Finally, the trajectories of the spectra are shown for the highly oscillatory singular integral operators.
The paper is organised as follows. In Section 2, for the integral equation (1.1) with (1.2), the modified Fourier basis is applied to obtain the explicit expressions for the resulting matrix. The asymptotic orders are proved for the entries with increasing oscillatory parameters. We construct the truncation scheme based on the asymptotic estimates. Then the numerical results are provided to show how the singular and oscillatory parameters affect the spectra. The similar analysis for the integral (1.1) with (1.3) is explored in Section 3. Section 4 concludes the paper.
Before giving the development of the method in details, we present a basic frame of the method for the integral equation (1.1) with (1.2) or (1.3). Applying an orthogonal system as the trial basis to approximate the unknown eigenfunctions , with the same test function as the trial basis, results in an algebraic system where For the inner integral, it holds that Thus, we have
Furthermore, if the orthogonal system is constituted of the even basis and odd basis , the eigenfunction in (1.1) is expanded as Substituting the expansion into (1.1), we get an infinite-dimensional algebraic eigenvalue problem where the elements of the matrices are and , . The parities of two sets and determine that the integrals in and are zeros. Therefore, the system (1.5) reduces to two uncoupled sub-systems The definitions of (1.6), (1.7) imply that and which means the matrices and are (complex) symmetric. In the next section, we explore the exact expressions and asymptotics of and by introducing the modified Fourier basis.
Section snippets
For the highly oscillatory singular kernel
In this section, we take the modified Fourier basis to discretise the integral equation of (1.1) with the kernel (1.2). We start from the integrals (1.6), (1.7).
For the logarithmic kernel (1.3)
In this section, we focus on the singular integral equation (1.1) with (1.3). Similar to the analysis in Section 2, we commence from the expansion of the special integral where is Euler constant and the function
Conclusions
We introduce one spectra computation method for the highly oscillatory integral equation with algebraic and logarithmic singularity. The advantage of the modified Fourier basis used to discretise the integral equation is that we can derive the explicit expression for the double highly oscillatory singular integral in the form of the special functions. Furthermore, the decay rate of the coefficient matrix can also be obtained by proving the asymptotics of the entries. The truncation strategy is
Acknowledgements
The author would like to thank Arieh Iserles, who gave the valuable content and direction of the paper.
References (22)
- et al.
An O(1) integration scheme for three-dimensional surface scattering problems
J. Comput. Appl. Math.
(2007) An optimal bound for the spectral variation of two matrices
Linear Algebra Appl.
(1985)- et al.
Spectral computation of highly oscillatory integral equations in laser theory
J. Comput. Phys.
(2019) Fast integration of rapidly oscillatory functions
J. Comput. Appl. Math.
(1996)Uniform estimates for oscillatory integral operators
J. Funct. Anal.
(1991)- et al.
The spectral problem for a class of highly oscillatory Fredholm integral operators
IMA J. Numer. Anal.
(2010) - et al.
The computation of the spectra of highly oscillatory Fredholm integral operators
J. Integral Equations Appl.
(2011) - et al.
Numerical-asymptotic boundary integral methods in high-frequency acoustic scattering
Acta Numer.
(2012) - et al.
Computing Highly Oscillatory Integrals
(2018) - et al.
Quadrature methods for highly oscillatory singular integrals
J. Comput. Math.
(2021)
A generalization of Filon-Clenshaw-Curtis quadrature for highly oscillatory integrals
BIT
Cited by (1)
Asymptotics on the Fredholm integral equation with a highly oscillatory and weakly singular kernel
2023, Applied Mathematics and Computation
- ☆
The work is supported by Natural Science Basic Research Plan in Shaanxi Province of China (Grant No. 2020JM-026) and National Natural Science Foundation of China (Grant No. 11871393).