有限差分法的基本知识.pptVIP

有限差分法的基本知识;一根紧拉着的均匀柔软弦，长为l，两端固定在X轴上O、L两点，当它在平衡位置附近做垂直于OL方向的微小横向振动时，求这根弦上各点的运动规律。;讨论如何将这一物理问题转化为数学上的定解问题。要确定弦的运动方程，需要明确：;简化假设：;横向：;由中值定理：;………一维波动方程;

假设外力在处外力密度为：方向垂直于轴。;弦振动方程中只含有两个自变量：。由于它描写的是弦的振动，因而它又称为一维波动方程。类似可以导出二维波动方程（如膜振动）和三维波动方程，它们的形式分别为：;建立数学物理方程是一个辩证分析的过程。

由于客观事物的复杂性，要求对所研究的对象

能够抓住事物发展的主要因素，摈弃次要因素，

使问题得到适度的简化。;☆均匀杆的纵振动;通过截面x，受到弹性力P(x,t)S的作用

通过截面x+dx受到弹性力P(x+dx,t)S的作用

P(x,t)为单位面积所受的弹性力(应力)，沿x方向为正．;☆静止空气中一维微小压力波的传播;代入;☆静止空气中三维声波方程;8.1.2扩散方程（抛物型方程）;所要研究的物理量：;利用Fourier热力学定律和能量守恒定律来建立热传导方程。;在dt时段内通过微元的两端流入的热量;为此所需的热量为;即：;根据热学中的傅立叶定律;;☆一维浓度扩散方程;8.1.3稳态方程(调和方程);这就是稳态方程，称为泊松方程。;又如在理想势流场中，存在速度势φ(x,y,z)，速度与φ(x,y,z)的关系为：;三类典型的偏微分方程;8.1.3有限差分法的基本知识

1、差分方程

2、截断误差

3、收敛性

4、稳定性;§8.1.3差分方程;有限差分法求解偏微分方程的基本过程是：首先将求解区域划分为差分网格，用有限个网格点代替连续的求解域，将待求解的变量（如密度、速度等）存储在各网格点上，并将偏微分方程中的微分项用相应的差商代替，从而将偏微分方程转化为代数形式的差分方程，得到含有离散点上的有限个未知变量的差分方程组。求出该差分方程组的解，也就得到了网格点上流动变量的数值解。;模型方程;Poisson方程：;;;离??网格点;

;;

;;;;;;;差分方程的建立过程;1.划分网格;;时间导数用一阶向前差商近似代替：;差分方程和其定解条件一起，称为相应微分方程

问题的差分格式。上述初值问题的差分格式可改写为：;显式有限差分模板：;时间推进：;例考虑长度为1的均匀直杆，其表面是绝热的，而且杆截面足够细，可;时间导数用一阶向前差商近似代替：;显式有限差分模板：;;如仍取而为缩短计算时间，时间步长取，则最终的差分方程：;§8.1.4截断误差;;上例中，令表示差分方程的精确解．利用Taylor级数将

上式中邻近节点的解在(i，n)点展开，整理并略去上标后可得

上式就是与差分方程等价的微分方程式。一般地说，任何一个微

分方程的差分方程,其差商都可以用Taylor级数表示，这样都可

以得到一个与差分方程对应的新的微分方程，该微分方程称为差

分方程的修正方程式。;上式中的就是差分方程与微分方程的差别，称之为截断误

差。显然与、成正比，一般情况下，当步长趋向零时,有限差分方程的截断误差是趋向于零的，则称有限差分方程与相应的偏微分方程是相容的。

一个可用的偏微分方程的差分表达式必须是相容的。否则在

、趋近零时，差分方程不能趋于原微分方程，差分方程的解就不能代表微分方程的解，差分求解就失去了意义!

;§8.1.6收敛性;§8.1.7稳定性