Kansa方法

✍ dations ◷ 2024-12-23 11:15:14 #数值微分方程

20世纪90年代,E. J. Kansa将用于散乱数据处理和函数近似的径向基函数用于处理偏微分方程,并提出一种强格式的配点方法。Kansa所提出的径向基函数配点方法是真正的无网格方法,具有易于编程、数学形式简单、方便掌握等优点。该方法提出后不久,被学术界称之为Kansa方法(Kansa method)。

由于径向基函数是采用无需考虑维数的一维欧几里德距离作为变量,Kansa方法适用于高维的和形状复杂的问题。Kansa方法是一种区域型方法,不仅在边界离散使其满足边界条件,同时内部配点需要满足控制方程。

此外,还有一类以径向基函数为核函数的边界型径向基函数配点方法(Boundary-type RBF collocation method),如基本解法(英语:Method of fundamental solution)、边界节点法、奇异边界法、边界粒子法、和正则化无网格法(Regularized meshless method)等。这类方法选取的基函数(也被称为核函数),通常选取控制方程的基本解或通解,因而满足控制方程。因此只需要在边界离散满足边界条件即可。

Kansa方法选取的径向基函数不需要满足控制方程,因此选取基函数有更大的自由空间。多元二次曲面(Multiquadric, MQ)函数是Kansa方法最常用的径向基函数,如果选择了恰当的形参数可以获得谱收敛的精度。

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

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

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

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

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

其中,

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

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

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

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

对于MQ函数,其插值误差取决于自身的形状参数,如何选取恰当的形状参数和关于MQ径向基函数的一些数学理论可以参见以下文献 。

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

相关

  • 角弓反张角弓反张是指项背强直,身体向后反折如弓状的临床表现。急性症状可由番木鳖碱毒素、破伤风、啡噻肼中毒等引起。 慢性或复发性发作可因癫痫、僵硬人症候群、歇斯底里等引起。
  • 见内文。鳐科(学名Rajidae)是鳐形目的典型科,本科鱼身体扁平,胸鳍宽大,腹鳍有一个明显的缺刻,几乎将腹鳍分为两叶,两个背鳍非常小,紧靠尾部,尾细小,通常无尾鳍,如有也没有鳍条支持。鳐
  • 化疗脑化疗脑,在学界又称化疗后认知异常(Post-chemotherapy cognitive impairment,简称PCCI)或化疗相关认知异常(Chemotherapy-Related Cognitive Impairments,简称CRCIs)等,所指的是因为
  • 藻华水华(Water bloom)或藻华(Algal bloom),通常为学术所称“水体富营养化”而造成,是发生在淡水中,由水体中氮磷含量过高导致藻类,细菌或浮游生物突然性过度增殖的一种自然现象,同时也
  • 便当盒饭是指使用盒子盛装饭菜或面条等各种食物以便于携带的盒装餐食,主要流行于亚洲以稻米作为主食的地区。“便当”一词最早源于南宋时期的俗语,意思是“便利的东西、方便、顺利
  • 鲱属鲱属(学名:Clupea)是生活在北太平洋和北大西洋(包括波罗的海)温带浅水中的一个多油的鱼属。这个属有两个种:大西洋鲱和太平洋鲱。这两个种又可以各分数个亚种。鲱鱼是成群的鱼,它们
  • 牙买加君主牙买加君主,牙买加的君主称号。牙买加的国家元首,由英国君主兼任。1962年8月6日一国牙买加独立后,王位设立,作为英联邦内的独立君主国。牙买加君主履行宪法职务,政治中立。由于英
  • 参差调谐参差调谐(staggered tuning)是用于多级调谐放大器的一种技术,其中每一级调谐到的频率有微小的差异。相比同步调谐(每一级都以相同频率进行调谐),参差调谐可以产生更宽的频带,但是其
  • 玛丽昂·内斯特尔玛丽昂·内斯特尔(Marion Nestle), 哲学博士, 是一名在纽约大学获得波莱特·戈达德荣誉奖的营养学教授。她在1988年到2003年期间在纽约大学营养学担任主席。她在加州大学伯
  • 毕波·伊斯特万毕波·伊斯特万(匈牙利语:Bibó István,1911年8月7日-1979年5月10日),出生于布达佩斯,匈牙利律师,政治理论家。有评论认为他是20世纪欧洲最深刻的政治思想家之一。在1956年匈牙利事