学位论文 > 优秀研究生学位论文题录展示

解二维Laplace方程Dirichlet问题直接边界积分方程的Galerkin边界元法

作　者: 董海云
导　师: 祝家麟
学　校: 重庆大学
专　业: 计算数学
关键词: Laplace 方程 Dirichlet 问题直接边界积分方程第一类Fredholm 积分方程 Galerkin 边界元法
分类号: O241.83
类　型: 硕士论文
年　份: 2005年
下　载: 285次
引　用: 3次
阅　读: 论文下载

内容摘要

Laplace 方程是最典型,最简单但应用广泛的椭圆型偏微分方程。用边界元法解边值问题,由不同的边界归化方法可以得到不同的边界积分方程,数值求解边界积分方程也有好几种方法。本文考虑用Green 公式和基本解推导得出直接边界积分方程来求解二维Laplace 方程的Dirichlet 问题,该直接边界积分方程是第一类Fredholm 积分方程。对二维问题,一般的带对数积分核第一类Fredholm 积分方程并不总是唯一可解的,特别是对外边值问题,解在无穷远处的形态有很大的影响。人们在用直接边界元方法进行计算时,并不刻意去考虑积分方程的可解性,但可解性的问题是不能回避的,这涉及到原问题的解与边界积分方程的解的等价性问题。事实上,对内边值问题,第一类Fredholm 直接边界积分方程的可解性条件是自然得到满足的,本文对此做了验证。对外边值问题,考虑到二维Dirichlet问题的解应当在无穷远处有界,故解的边界积分表达式要做修正,对积分方程的解要有约束,这样去解边界积分方程得出的解才等同于原问题的解。一般来说,直接边界积分方程可以很方便的用配点法求解,还未见有实际用Galerkin边界元来解的报道。本文采用Galerkin边界元方法求解直接边界积分方程,是为了验证这两种方法的效率和精度,且Galerkin 法易于进行收敛性分析。Galerkin边界元方法是把积分方程转化为等价的边界变分方程,经用边界元离散后,通过求解线性代数方程组和计算解的离散的积分表达式求得原问题的数值解,该方法需要在边界上计算重积分。本文推出了第一重积分的解析计算公式,对外层积分则采用高斯数值积分。对外边值问题,第一类Fredholm 积分方程的解要附加在边界上积分为零的条件,本文采用Lagrange 乘子放松这个约束,求解扩展的变分方程时,可同时得出解在无穷远的值。本文采用常单元和线性元这两种离散方式,分别用Fortran90 编写了计算程序,对误差与边界元的数量的关系做了数值实验。数值算例表明本文的计算方法是可行的、有效的,用Galerkin 方法求解边界积分方程比用配点法精度高。

全文目录

中文摘要  4-5
英文摘要  5-8
1 绪论  8-15
  1.1 研究现状评述  8-13
    1.1.1 边界元法的历史回顾和发展现状  8-11
    1.1.2 用边界元法解调和方程的几种主要边界积分方程和解的表达式  11-13
  1.2 本文主要研究内容和研究方法  13-15
2 Laplace 方程直接边界积分方程的 Galerkin 边界元解法  15-28
  2.1 Laplace 方程Dirichlet 问题的直接边界积分方程  15-16
  2.2 关于第一类Fredholm 积分方程的可解性问题  16
  2.3 直接边界积分方程的可解性问题  16-19
    2.3.1 关于内边值问题可解性的讨论  16-18
    2.3.2 关于外边值问题可解性的讨论  18-19
  2.4 变分方程及其边界离散  19-25
    2.4.1 直接边界积分方程的变分方程  19-20
    2.4.2 用常单元离散变分方程  20-22
    2.4.3 用线性元离散变分方程  22-25
  2.5 矩阵元素的处理和计算  25-28
3 程序实现  28-56
  3.1 程序框架简介  28-33
  3.2 源程序  33-56
4 算例和数值实验  56-67
5 结论  67-68
致谢  68-69
参考文献  69-72
附录:作者在攻读硕士学位期间发表的论文目录  72-73
独创性声明  73
学位论文版权使用授权书  73

解二维Laplace方程Dirichlet问题直接边界积分方程的Galerkin边界元法

内容摘要

全文目录

相似论文