A new spectral method approach for singular integral equations

本研究计划的核心是引入一种新的、快速的、光谱精确的算法来求解复杂一维边界上的一般奇异积分方程，该算法允许在两个空间维度上表示椭圆型偏微分方程的解。奇异积分方程在电磁学、地震成像、断裂力学、流体力学和波束物理的声散射中有着丰富的历史。******该计划的结果将在一个用Julia编程语言编写的开源包中实现。边界密度端点奇点的理论确定允许直接求解器在不使用hp自适应细化的情况下获得光谱精确的全局解。通过成功地进一步开发一类新的直接求解器，该软件包将以稳定和及时的方式求解广泛的奇异积分方程。******最近引入的直接求解器将与基于递归块对角化的分层求解器结合使用。这将特别利用由椭圆型偏微分方程的强制奇异积分算子引起的层次非对角低秩结构。分层求解器涉及一个独立于强迫项的预计算阶段。一旦这种预计算分解了算子，求解许多强迫项的复杂性就会降低，因此只需要一小部分时间。******本程序还将考虑定义在一类重要边界上的奇异积分方程：那些是从单位区间和圆的多项式映射的方程。我们将通过谱映射定理对带奇异积分算子进行大量的分析。求解具有混合边界条件或多重连通轮廓的奇异积分方程会导致结点处具有复杂代数奇异结构的分段定义解。这些困难将通过设计完全捕捉在连接处产生的复杂单一结构的碱基来解决。******新的谱法将应用于求解Stokes流问题、双调和方程和断裂力学的应力应变计算。新的谱方法与稳定和高阶时间步进格式的结合将允许探索Helmholtz方程的时域积分公式和RayleighTaylor不稳定性的模拟。它还将允许在重要应用中进行实验和潜在的新现象发现，如纳米尺度的光学元，逆散射问题的解决，以及深水内波的本杰明诺方程的模拟。

This research programme is centred around the introduction of a new, fast, and spectrally accurate algorithm for solving general singular integral equations on complicated one-dimensional boundaries, which allows for a representation of the solution of elliptic partial differential equations in two spatial dimensions. Singular integral equations have a rich history in acoustic scattering for electromagnetics and seismic imaging, fracture mechanics, fluid dynamics, and beam physics.******Results of the programme will be implemented in an open source package written in the Julia programming language. Theoretical determination of the endpoint singularities of the boundary densities allows for the direct solver to obtain spectrally accurate global solutions without the use of h-p adaptive refinement. By successfully furthering the development of a new class of direct solvers, the software package will solve a wide range of singular integral equations in a stable and timely manner.******The recently introduced direct solver will be combined with a hierarchical solver based on recursive block diagonalization via Schur complements. This will specifically exploit the hierarchically off-diagonal low-rank structure arising from coercive singular integral operators of elliptic partial differential equations. The hierarchical solver involves a pre-computation phase independent of the forcing term. Once this pre-computation factorizes the operator, the solution to many forcing terms has a lower complexity and therefore takes a fraction of the time.******This programme will also consider singular integral equations defined on an important class of boundaries: those that are polynomial maps from the unit interval and circle. A considerable analysis will be performed to again obtain banded singular integral operators via the spectral mapping theorem. Solving singular integral equations with either mixed boundary conditions or multiply connected contours leads to piecewise-defined solutions with complicated algebraic singular structure at the junctions. These difficulties will be approached by designing bases that fully capture this complicated singular structure arising at the junctions.******The new spectral method will be applied to solve problems of Stokes flow, the biharmonic equation and stress and strain computations for fracture mechanics. Combination of the new spectral method with stable and high-order time-stepping schemes will allow for the exploration of time-domain integral formulations of the Helmholtz equation and the simulation of RayleighTaylor instability. It will also allow experimentation and potential discovery of new phenomena in important applications such as optical metacages at the nanoscale, the solution of inverse scattering problems, and simulation of the BenjaminOno equation for internal waves in deep water.