Kansa方法

Kansa方法，也称为改进MQ方法或MQ配点法，源于著名的MQ插值。该方法的有效性和高效性已经得到了许多问题的验证。此外，由于部分问题的基本解和通解不存在，如变系数问题和非线性问题，Kansa方法比边界型径向基函数配点方法拥有更加广阔的应用范畴。

在d维物理区域 内考虑以下边值问题，

${\displaystyle Lu(X)}$  =${\displaystyle f(X)}$  ,${\displaystyle X\in \Omega ,}$  (1)

${\displaystyle u(X)}$  =${\displaystyle g(X)}$  ,${\displaystyle X\in \partial \Omega _{D},}$ (2)

${\displaystyle {\frac {\partial u(X)}{\partial n}}}$ =${\displaystyle h(X)}$  , ${\displaystyle X\in \partial \Omega _{N},}$  (3)

其中L代表微分算子，d为问题的维数，${\displaystyle \partial \Omega _{D}}$  , ${\displaystyle \partial \Omega _{N}}$ 分别代表狄利克雷边界和诺伊曼边界且${\displaystyle \partial \Omega _{D}\cup \partial \Omega _{N}=\partial \Omega }$ 。 Kansa方法通过径向基函数的线性组合来逼近待求的函数，即：

${\displaystyle {{u(X)}^{*}}}$ =${\displaystyle \sum \limits _{i=1}^{N}\alpha _{i}\phi \left(r_{i}\right)}$ , (4)

其中${\displaystyle {{\alpha }_{i}}}$ 为待求参数，${\displaystyle \phi \left(r_{i}\right)}$ 代表径向基函数，如MQ函数。 为确保所求函数的唯一性，在上式右端添加一组多项式：

${\displaystyle {{u(X)}^{*}}}$ =${\displaystyle \sum \limits _{i=1}^{N}\alpha _{i}\phi \left(r_{i}\right)}$ +${\displaystyle \sum \limits _{k=1}^{M}\alpha _{k+N}\gamma _{k}\left(X\right)}$ , (5)

其中 ${\displaystyle {{\gamma }_{k}}(X)}$  为多项式。径向基插值形式(4)和(5)都经常应用在计算中。(4)式的形式简单易于掌握且在大多数情况下都能获取较好的计算结果，因此在工程领域使用广泛；(5)式的形式严谨且理论基础坚实，所以数学工作者偏于采用后者。 将(4)式或(5)式代入方程组(1)-(3)可得以下线性方程组：

${\displaystyle A\alpha =b}$ , (6)

其中，

${\displaystyle \mathbf {A} =\left({\begin{matrix}L({\phi })&L({\gamma })\\{\phi }&{\gamma }\\{\frac {\partial {\phi }}{\partial n}}&{\frac {\partial {\gamma }}{\partial n}}\\{\gamma }&0\\\end{matrix}}\right)}$ , ${\displaystyle \mathbf {b} =\left({\begin{matrix}f\\g\\h\\0\\\end{matrix}}\right)}$ ,${\displaystyle {\phi }}$ =${\displaystyle \phi \left(x_{i},x_{j}\right)}$ ,${\displaystyle {\gamma }}$ =${\displaystyle \gamma _{k}\left(X_{i}\right)}$ . (7)

通过求解以上线性方程组，可求解待定参数${\displaystyle {{\alpha }_{i}}}$ ，根据(4)式或(5)式即可得到待求函数。

偏微分方程的数值求解通常采用有限差分法，有限单元法或边界单元法。有限差分法通常需要规则的网格系统，难以处理不规则区域问题。比之有限差分法，有限单元法能适于处理更复杂的形状，但网格的划分及其再划分在计算时依旧不可避免。边界单元法在处理一些工程问题效果显著，比如反问题、无限域问题和薄壁结构问题。然而，边界单元法受限于控制方程的基本解难以确定，使其应用范畴受到约束。

近来几十年，由于标准有限单元法和边界单元法在处理高维、移动边界和复杂边界等问题需要耗费庞大的计算成本，无网格或无单元方法受到极多关注。Kansa方法^[1][2]是一种真正的无网格方法，不需要划分网格和单元而是通过径向基函数（如MQ函数）在配置的节点处满足相关条件即可。

虽然经过诸多学者的研究，但依旧缺乏对Kansa方法严谨的数学证明^[3]。另外，混合边界会破坏插值矩阵的对称性。文献^[4][5]提及的对称埃尔米特径向基函数插值方案(Hermite RBF collocation scheme)其可解性具有可靠的数值分析。其中，Kansa方法和对称埃尔米特方法都存在一个共同的问题，即相邻边界节点的数值解精度比内部节点低1-2个数量级。边界偏微分方程配点(The PDE collocation on the boundary, PDECB))方案可以消除这一缺陷。然而，这一方案缺乏数学上的理论支持且需要在边界附近的区域内或区域外设置一系列节点，因此在处理复杂区域或多连通问题时非常复杂。随后提出的一种相似处理方式^[6]，在相同的边界节点同时满足控制方程和边界条件，而其缺陷在于产生的插值矩阵是不对称的且方法本身同样缺乏明确的理论基础。通过使用第二格林公式，改进的Kansa方法^[7][8]可以弥补以上缺陷。

对于MQ函数，其插值误差取决于自身的形状参数，如何选取恰当的形状参数和关于MQ径向基函数的一些数学理论可以参见以下文献 ^{[9][10][11][12]}。

Kansa方法广泛应用于计算科学。^[1]中Kansa方法用于求解椭圆型、双曲型和抛物型三类偏微分方程。Kansa近来也应用于求解各类常微分和偏微分方程，包括两相和三相混合模型的组织工程问题^[13][14]，冲击波下的一维非线性Burger方程^[15]，潮汐和海流模拟中的浅水方程 ^[16]，热传导方程^[17]，自由边界问题^[18]，分数阶扩散方程^[19]。

• 径向基函数
• 基本解法（英语：Method of fundamental solution）
• 边界节点法
• 边界粒子法
• 奇异边界法

• Modified Kansa method