有限元方法
有限元方法简介
汪老师通过一个简单的常微分方程初边值问题求解的例子引出有限元方法,当真解未知时,通过数值离散方法,求出向量u(u是真解在n+1个点处的值作为分量的向量),然后再利用一些插值方法(比如说分片线性函数插值)解出方程的数值解uh,随后很重要的一步是误差估计。从理论上证明uh是对u的逼近。
简略的写为以下几方面:
- 提出问题P
- 真解u未知
- 通过数值离散方法求出真解u在一些点处的值(近似意义下)
- 利用插值方法求出真解u的数值解uh
- 误差估计(从理论上证明uh是对u的逼近)
常见的数值离散方法有:
- 有限差分法
- 有限元方法
- 有限体积法
等等。
注:这几个“有限”指的都是有限维线性空间
具体的实际问题有:
- 纯扩散问题
比如往一个盛满清水的容器里滴入一滴颜料,颜料的扩散浓度 - 弹性问题
地震时的裂缝穿过弹性材料 - 流体力学-Stokes方程
流体通过阻碍物会产生一系列的旋涡,比如云通过火山口 - 分子束外延薄膜生长
利用特殊的金属枪往半导体材料上铺粒子
老师举的这几个例子都很有意思。
微分方程简介
微分方程分类
讲了三种微分方程的分类
- 根据方程中导数性质
- 常微分方程(ODE)
- 偏微分方程(PDE)
- 根据方程中最高阶导数
- 一阶
- 二阶
- 三阶…
- 根据方程是否线性
- 线性方程:方程关于u和u的导数项都是线性的
- 非线性方程
拟线性方程(关于最高阶导数项是线性的)
完全线性方程(关于最高阶导数项不是线性的)
判断一个方程具体是什么方程,首先应该理解的有自变量、未知函数、已知函数。
二阶线性偏微分方程的分类(以二维为例)
通过考虑二阶项
对二阶线性偏微分方程进行分类,以二维线性空间为例,方程的形式如下:
A
u
x
x
+
B
u
x
y
+
C
u
y
y
+
D
u
x
+
E
u
y
+
F
=
f
Au~xx~+Bu~xy~+Cu~yy~+Du~x~+Eu~y~+F=f
Au xx +Bu xy +Cu yy +Du x +Eu y +F=f
其中
A
,
B
,
C
,
D
,
E
,
F
,
f
A,B,C,D,E,F,f
A,B,C,D,E,F,f均为已知量。
I=AC-(B2)/4
(这里是二阶项的系数矩阵行列式,我不知道怎么输入:D,有木有人帮助一下,)
I>0-----椭圆形方程
I=0-----抛物型方程
I<0-----双曲型方程
若是n维,则看二阶项系数矩阵的特征值情况:
- 若有零特征值,则为抛物型方程
- 若所有特征值同号,则为椭圆型方程
- 若有异号特征值,则为双曲型方程
二阶椭圆型偏微分方程的边值问题
- 第一类边值问题:狄利克雷型边值问题(关于u)
- 第二类边值问题:诺依曼型边值问题(关于u的导数)
- 第三类边值问题:混合型边值问题(关于u和u的导数)