基于边界积分方程的Galerkin无网格方法

论文摘要

无网格方法在近年来得到广泛关注,其基本特点是场函数建立在独立的节点上,节点之间无需网格联接。边界积分方程能使所考虑问题的维数降低一维,是求解线性问题和外部问题的一种有效工具。基于边界积分方程的无网格方法是无网格方法的一个重要分支。本文首先回顾了无网格方法的发展历史和研究现状,综述了无网格方法数学理论的研究进展,介绍了无网格方法的基本原理,总结了无网格方法的特点、优越性以及目前无网格方法的难点和存在的问题。然后在大量前人工作的基础上,提出了一种新的基于边界积分方程的Galerkin无网格方法——Galerkin边界点法,并成功地将其应用于位势问题、弹性力学问题和流体力学问题的求解。在Galerkin边界点法中,首先将边值问题归结为边界积分方程的弱形式或变分公式,然后利用移动最小二乘近似构造变分公式中的试探函数和检验函数,进而得到近似的有限维逼近空间。Galerkin边界点法利用了移动最小二乘近似的无网格思想和边界积分方程的降维特性,因此它的输入数据只是求解域边界上的离散分布的点。由于引入了变分公式,Galerkin边界点法能保持变分问题的对称性和正定性,该性质使得Galerkin边界点法是耦合有限元方法或者其它已经建立的区域型无网格方法(如无单元Galerkin法)的理想方法,这种耦合技术非常适合求解无限域问题。另外,虽然用移动最小二乘近似构造的形函数不具有插值特性,但是通过把边界函数与检验函数相乘并在边界上积分,Galerkin边界点法中的边界条件仍能容易地精确满足。本论文针对Galerkin边界点法进行了理论分析和数值应用,具体研究工作如下:首先研究了任意边界上的移动最小二乘近似,给出了移动最小二乘近似的性质。当节点和权函数满足一定条件时,证明了移动最小二乘近似的近似函数在Sobolev空间中的最优误差估计。误差结果表明,移动最小二乘近似的逼近误差与节点间距密切相关。其次给出了Galerkin边界点法求解作为拟微分算子方程的边界积分方程的算法。边界积分方程首先被转化为相应的变分形式,然后用移动最小二乘近似的形函数构造解空间。为了计算积分,我们在边界上构造了背景网格。基于移动最小二乘近似的误差公式和拟微分算子理论,推导了用Galerkin边界点法求解边界积分方程的解的误差公式。从误差分析的过程中可以看出,Galerkin边界点法的误差主要来自两个方面:一是用背景网格上的积分去近似边界积分;二是用移动最小二乘近似去逼近边界变量。我们还考虑了边界积分方程的未知量需要满足一定约束条件的情形,此时我们采用Lagrange乘子放松这个约束,并给出了相应的数值实施过程和误差估计。第三,用Galerkin边界点法求解了Laplace问题、双调和问题、线弹性问题和Stokes问题。把边值问题归化为等价的第一类边界积分方程,再用Galerkin边界点法进行求解。我们给出了数值求解过程,并就一般情形详细地进行了误差分析,得到了最佳渐进误差估计。当积分用的背景网格和边界重合时,我们进一步得到了近似解的能量模估计。数值算例表明了这种方法的有效性,并且数值结果和理论分析吻合。

论文目录

• 摘要
• ABSTRACT
• 1 绪论
• 1.1 引言
• 1.2 无网格方法的研究历史和现状
• 1.2.1 基于全局Galerkin 弱式的无网格方法
• 1.2.2 基于Petrov-Galerkin 弱式的无网格方法
• 1.2.3 基于配点强式的无网格方法
• 1.2.4 基于边界积分方程的无网格方法
• 1.3 无网格方法数学理论研究进展
• 1.4 无网格方法的基本原理
• 1.4.1 无网格方法的基本近似
• 1.4.2 无网格方法的积分方案
• 1.4.3 边界条件的处理
• 1.5 无网格方法评述
• 1.6 本文主要工作
• 2 移动最小二乘（MLS）近似
• 2.1 引言
• 2.2 预备知识
• 2.3 MLS 近似的基本原理
• 2.4 MLS 近似的误差估计
• 2.5 本章小结
• 3 边界积分方程的无网格近似和误差分析
• 3.1 引言
• 3.2 拟微分算子方程
• 3.3 GALERKIN 过程
• 3.4 MLS 近似
• 3.5 逼近空间的建立和近似解的构造.
• 3.5.1 积分背景网格的构造
• 3.5.2 建立逼近空间
• 3.5.3 问题的近似解
• 3.6 误差分析
• 3.7 背景网格与边界重合时的误差分析
• 3.8 约束条件的处理
• 3.8.1 问题的近似解
• 3.8.2 误差分析
• 3.8.3 背景网格与边界重合
• 3.9 本章小结
• 4 位势问题的GALERKIN 边界点法
• 4.1 引言
• 4.2 LAPLACE 方程DIRICHLET 问题
• 4.2.1 解的积分表示及变分公式
• 4.2.2 近似解及数值实施
• 4.2.3 误差分析
• 4.2.4 背景网格与边界重合
• 4.2.5 数值算例
• 4.3 LAPLACE 方程NEUMANN 问题
• 4.3.1 解的积分表示及变分公式
• 4.3.2 近似解及数值实施
• 4.3.3 误差分析
• 4.3.4 背景网格与边界重合
• 4.3.5 数值算例
• 4.4 双调和问题
• 4.4.1 解的积分表示及变分公式
• 4.4.2 近似解及数值实施
• 4.4.3 误差分析
• 4.4.4 数值算例
• 4.5 本章小结
• 5 弹性问题的GALERKIN 边界点法
• 5.1 引言
• 5.2 解的积分表示及变分公式
• 5.3 近似解及数值实施
• 5.4 误差分析
• 5.5 背景网格与边界重合
• 5.6 数值算例
• 5.7 本章小结
• 6 STOKES 问题的GALERKIN 边界点法
• 6.1 引言
• 6.2 解的积分表示及变分公式
• 6.3 近似解及数值实施
• 6.4 误差分析
• 6.5 背景网格与边界重合
• 6.6 数值算例
• 6.7 本章小结
• 7 总结与展望
• 7.1 总结
• 7.2 展望
• 致谢
• 参考文献
• 附录