Approximate solution of singular integral equations of the first kind with Cauchy kernel

AbstractIn this work a study of efficient approximate methods for solving the Cauchy type singular integral equations (CSIEs) of the first kind, over a finite interval, is presented. In the solution, Chebyshev polynomials of the first kind, Tn(x), second kind, Un(x), third kind, Vn(x), and fourth kind, Wn(x), corresponding to respective weight functions W(1)(x)=(1−x2)−12, W(2)(x)=(1−x2)12,W(3)(x)=(1+x)12(1−x)−12 and W(4)(x)=(1+x)−12(1−x)12, have been used to obtain systems of linear algebraic equations. These systems are solved numerically. It is shown that for a linear force function the method of approximate solution gives an exact solution, and it cannot be generalized to any polynomial of degree n. Numerical results for other force functions are given to illustrate the efficiency and accuracy of the method

A fast and well-conditioned spectral method for singular integral equations

We develop a spectral method for solving univariate singular integral equations over unions of intervals by utilizing Chebyshev and ultraspherical polynomials to reformulate the equations as almost-banded infinite-dimensional systems. This is accomplished by utilizing low rank approximations for sparse representations of the bivariate kernels. The resulting system can be solved in ${\cal O}(m^2n)$ operations using an adaptive QR factorization, where $m$ is the bandwidth and $n$ is the optimal number of unknowns needed to resolve the true solution. The complexity is reduced to ${\cal O}(m n)$ operations by pre-caching the QR factorization when the same operator is used for multiple right-hand sides. Stability is proved by showing that the resulting linear operator can be diagonally preconditioned to be a compact perturbation of the identity. Applications considered include the Faraday cage, and acoustic scattering for the Helmholtz and gravity Helmholtz equations, including spectrally accurate numerical evaluation of the far- and near-field solution. The Julia software package SingularIntegralEquations.jl implements our method with a convenient, user-friendly interface

An integral equation method for solving neumann problems on simply and multiply connected regions with smooth boundaries

This research presents several new boundary integral equations for the solution of Laplace’s equation with the Neumann boundary condition on both bounded and unbounded multiply connected regions. The integral equations are uniquely solvable Fredholm integral equations of the second kind with the generalized Neumann kernel. The complete discussion of the solvability of the integral equations is also presented. Numerical results obtained show the efficiency of the proposed method when the boundaries of the regions are sufficiently smooth

On the convergence of local expansions of layer potentials

In a recently developed quadrature method (quadrature by expansion or QBX), it was demonstrated that weakly singular or singular layer potentials can be evaluated rapidly and accurately on surface by making use of local expansions about carefully chosen off-surface points. In this paper, we derive estimates for the rate of convergence of these local expansions, providing the analytic foundation for the QBX method. The estimates may also be of mathematical interest, particularly for microlocal or asymptotic analysis in potential theory

Unsteady two dimensional airloads acting on oscillating thin airfoils in subsonic ventilated wind tunnels

The numerical calculation of unsteady two dimensional airloads which act upon thin airfoils in subsonic ventilated wind tunnels was studied. Neglecting certain quadrature errors, Bland's collocation method is rigorously proved to converge to the mathematically exact solution of Bland's integral equation, and a three way equivalence was established between collocation, Galerkin's method and least squares whenever the collocation points are chosen to be the nodes of the quadrature rule used for Galerkin's method. A computer program displayed convergence with respect to the number of pressure basis functions employed, and agreement with known special cases was demonstrated. Results are obtained for the combined effects of wind tunnel wall ventilation and wind tunnel depth to airfoil chord ratio, and for acoustic resonance between the airfoil and wind tunnel walls. A boundary condition is proposed for permeable walls through which mass flow rate is proportional to pressure jump