环形域Laplace方程的边界积分方程快速傅里叶变换求解与数值模拟研究

《Partial Differential Equations in Applied Mathematics》：Boundary element method for Laplace equation in a ring domain

针对环形域上Laplace方程的混合边值问题，传统数值方法存在效率瓶颈。本研究提出一种基于边界积分方程（BIE）与标准Galerkin方法的快速求解方案。通过引入一周期B-样条离散化，将边界积分算子转化为具有循环块结构的稠密线性系统，并利用快速傅里叶变换（FFT）实现矩阵-向量乘法的加速。结果表明，该方法在保持高精度的同时，显著降低了计算复杂度（O(N ln N)量级），为工程与物理中的多连通区域问题提供了高效、稳定的数值工具。

　　

在科学与工程领域，许多物理现象，如热传导、电磁场分布和流体流动，都可通过偏微分方程描述。其中，Laplace方程作为一类经典的椭圆型方程，在稳态场问题中扮演着核心角色。当求解区域呈环形或多连通结构时，问题变得尤为复杂——边界不连通、内外场耦合等因素使得传统数值方法（如有限元法）面临计算效率低、存储需求大的挑战。尤其在处理混合边界条件（如部分区域已知势函数、部分区域已知通量）时，如何高效、精确地求解成为数值计算中的难点。

为此，研究人员在《Partial Differential Equations in Applied Mathematics》上发表了题为“A fast Fourier transform solution of the Laplace equation in a ring domain with mixed boundary conditions”的论文，针对二维环形域上的Laplace方程，提出了一种基于边界积分方程（Boundary Integral Equation, BIE）与快速傅里叶变换（Fast Fourier Transform, FFT）结合的创新数值方法。该方法不仅克服了传统方法的瓶颈，还通过巧妙的离散化技术，实现了计算复杂度从O(N²)到O(N ln N)的跨越，为多连通区域问题的求解提供了新思路。

研究中采用的关键技术方法主要包括：1. 边界积分方程（BIE）与Steklov-Poincaré算子理论，将原偏微分方程转化为边界上的积分方程，降低问题维度；2. 标准Galerkin离散化，利用一周期B-样条（B-splines）作为基函数，对边界积分算子进行离散，得到具有循环块结构的稠密线性系统；3. 快速傅里叶变换（FFT）算法，利用离散Fourier矩阵作为预条件子，加速矩阵-向量乘法与线性系统求解，显著提升计算效率；4. 特征系统分析，通过参数化边界推导边界积分算子的特征值与特征函数，为离散矩阵的特征值估计提供理论依据。

研究结果

2.1. 问题建模与边界积分方程建立

针对由两个同心圆边界Γ₀、Γ₁构成的环形域Ω，建立了Laplace方程的Dirichlet边值问题与混合边值问题的数学模型。通过Green表示公式与迹算子，推导出矩阵形式的边界积分方程组Vt=Kg，其中V为单层位势算子（single-layer potential operator）、K为双层位势算子（double-layer potential operator），t为未知的法向导数，g为已知边界数据。

2.2. 边界积分算子的性质与特征系统

证明单层算子V是自伴、H^-1/2(Γ)椭圆且有界的，确保变分问题解的唯一存在性。通过引入角度参数化，得到边界积分算子的显式积分核，并利用Fourier级数展开，解析推导出其特征值与特征函数。例如，对于同心圆边界，V_Γ?Γ?的特征值为λ_n=R_?/(2|n|) (n≠0)，为后续离散矩阵特征值分析奠定基础。

3.1. Galerkin离散化与矩阵结构

采用一周期B-样条基函数对变分形式进行Galerkin离散，得到形如Ax=b的线性系统。矩阵A由四个子块构成，分别对应不同边界间的相互作用。由于圆周参数化的周期性与B-样条的平移不变性，离散后的矩阵具有循环块结构（block-circulant structure），使得FFT可应用于矩阵-向量乘法。

3.2. 快速算法设计与计算复杂度

利用FFT实现离散卷积运算，将矩阵-向量乘法复杂度从O(N²)降至O(N ln N)。通过离散Fourier矩阵对单层算子矩阵进行预条件，改善系数矩阵的条件数，结合共轭梯度法（conjugate gradient method）求解线性系统，在保持精度的同时大幅减少计算时间。

4. 数值算例与误差分析

以Dirichlet边值问题为例，选取解析解u(r,θ)=r²cos(2θ)进行测试。在不同网格尺寸下比较了共轭梯度法与直接法的计算时间。结果显示，当N=1024时，FFT加速的共轭梯度法求解时间仅为直接法的1/20。通过计算L²误差与收敛阶，验证了该方法具有最优收敛阶O(h^θ+1)，其中θ为B-样条阶数。

5. 混合边值问题的求解与应用

将方法推广至混合边值问题（Γ₀上为Dirichlet条件、Γ₁上为Neumann条件）。通过Steklov-Poincaré算子将问题重新表述为边界积分方程，并采用相同的Galerkin-FFT框架求解。数值实验表明，该方法能稳定处理混合边界条件，且计算效率与纯Dirichlet问题相当。

结论与讨论

本研究发展了一种基于边界积分方程与快速傅里叶变换的高效数值方法，用于求解环形域上Laplace方程的Dirichlet与混合边值问题。通过理论分析证明边界积分算子的椭圆性与有界性，确保数值格式的稳定性；利用B-样条离散与FFT加速，实现计算复杂度的显著降低。数值实验验证了方法在计算精度、效率与适用性方面的优势。

该方法的重要意义体现在：1. 为多连通区域问题提供了维数约简（dimension reduction）的高效求解框架，克服了区域网格生成的困难；2. 通过FFT利用循环结构，大幅降低存储与计算成本，使大规模问题求解成为可能；3. 可扩展至其他椭圆型方程（如Helmholtz方程、Stokes方程）及更复杂边界形状，为工程计算（如电磁散射、微流体模拟）与物理建模（如势场分析）提供通用工具。未来工作可进一步研究非同心圆边界、三维环形域及非线性问题的扩展。