1. 引言
众所周知,Laplace方程 [1] (又称调和方程)
![](//html.hanspub.org/file/7-2620429x10_hanspub.png)
是最典型的椭圆型偏微分方程之一,它具有广泛的应用背景,例如静电学中的电势及牛顿万有引力理论中的引力势均满足该方程。传统上人们认为有限元法擅长计算复杂区域上椭圆型偏微分方程的适定问题 [2] ,但有限元法计算步骤较为复杂和抽象。与网格方法相比,无网格方法的最大优点是它不依赖于求解区域的网格划分,从而有效避免了网格划分所需要的巨大而繁重的工作量,特别是当数值求解过程中求解区域和边界条件随时间而变动时,无网格方法则具有传统网格方法无法比拟的优势 [3] [4] [5] [6] [7] 。鉴于以上分析,本文考虑二维不规则区域上Laplace方程的无网格方法。
首先利用Laplace算子的基本解作为基函数构造近似解,进而通过配置法利用边界条件得到数值解,数值实验的结果表明了此方法的可行性和精确性。
2. 基本解方法
基本解方法是一种数值求解椭圆型方程边值问题的边界型方法,它用微分算子的基本解的线性组合来逼近问题的解。由于基本解满足偏微分方程,只须考虑边界条件,因此属于边界型方法 [8] 。
考虑边值问题
![](//html.hanspub.org/file/7-2620429x11_hanspub.png)
其中
是微分算子,
是边界条件算子,
是已知函数,
是求解区域
的边界。
假设
是微分算子
的基本解,则基本解方法的步骤可以简述如下:
1) 选取中心点
,其中
是包含求解区域
的封闭曲面(线);
2) 构造
用于近似原问题的解
;
3) 选择边界点
,并令
![](//html.hanspub.org/file/7-2620429x25_hanspub.png)
特别地,当
为线性算子时,有
![](//html.hanspub.org/file/7-2620429x27_hanspub.png)
4) 求解上述方程组得到数值解。
3. Laplace方程的基本解无网格方法
当微分算子
时,考虑二维Laplace方程的Dirichlet边值问题
![](//html.hanspub.org/file/7-2620429x29_hanspub.png)
由文献 [9] 知微分算子
的基本解是
![](//html.hanspub.org/file/7-2620429x31_hanspub.png)
其中
。
此时,为求数值解
,只须求解线性方程组
![](//html.hanspub.org/file/7-2620429x34_hanspub.png)
其中
。
附录中给出了具体的MATLAB求解程序。
4. 数值实验
考虑二维不规则区域上Laplace方程的边值问题
![](//html.hanspub.org/file/7-2620429x36_hanspub.png)
![](//html.hanspub.org/file/7-2620429x37_hanspub.png)
其中
。求解区域的边界
是由一个大椭圆和两个小圆组成。该问题的真解是
。
为衡量数值解的精度,定义最大绝对误差如下
![](//html.hanspub.org/file/7-2620429x41_hanspub.png)
实验过程中,中心点的选择采用如下两种方案
方案I:在外边界(大椭圆)外部选择中心点如图1所示;
方案II:在内边界(两个小圆)的内部和外边界(大椭圆)外部分别选择中心点如图2所示。
为了考察中心点个数对数值解精度的影响,表1和表2分别列出了两种方案下不同中心点个数对应的最大绝对误差,从表中数据可以看出随着中心点个数的增加,数值解精度在不断提高,并且方案I的精度明显优于方案II。
![](//html.hanspub.org/file/7-2620429x42_hanspub.png)
Figure 1. Chosen points, numerical solution and absolute error for scheme I
图1. 标方案I的选点、数值解及误差情况
![](//html.hanspub.org/file/7-2620429x43_hanspub.png)
Figure 2. Chosen points, numerical solution and absolute error for scheme II
图2. 标方案II的选点、数值解及误差情况
![](Images/Table_Tmp.jpg)
Table 1. Comparison of ABE for scheme I
表1. 方案I的最大绝对误差比较
![](Images/Table_Tmp.jpg)
Table 2. Comparison of ABE for scheme II
表2. 方案II的最大绝对误差比较
5. 结论
通过本文结果可以知道,对于不规则区域上椭圆型偏微分方程的边值问题,采用基本解方法求解是方便可行的,同时启示我们构造高精度的无网格方法是值得研究的。
基金项目
大学生创新创业训练计划项目(201610452005),山东省自然科学基金项目 (BS2015DX012)。
附录
for s=1:M
for t=1:M
q(s,t)=sqrt((node1(s,1)-node(t,1))^2+(node1(s,2)-node(t,2))^2);
A(s,t)=log(1/q(s,t));
end
end
alfa=A\f;
for t=1:M
qq=((X-node(t,1)).^2+(Y-node(t,2)).^2).^(1/2);
u=u+alfa(t).*log(1./qq);
end
期刊投稿者将享受如下服务:
1.投稿前咨询服务 (QQ、微信、邮箱皆可)
2.为您匹配最合适的期刊
3.24小时以内解答您的所有疑问
4.友好的在线投稿界面
5.专业的同行评审
6.知网检索
7.全网络覆盖式推广您的研究
投稿请点击:http://www.hanspub.org/Submission.aspx
期刊邮箱:aam@hanspub.org
NOTES
*通讯作者。