偏微分方程的数值解(差分法).pptVIP

解: (1) 作平行于坐标轴的的两族直线, 间隔: (4) 差分方程为: (3) 满足稳定性条件: (2) 根据边界条件与初始条件写出边界上各点的函数值 (5) 解: 计算实例： * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * (1)分割时空网格： 若时间导数用一阶向前差商近似代替，即 空间导数用一阶中心差商近似代替，即 则在 点的对流方程就可近似地写作 (2)差商近似代替 (1)用时间向前差商代替时间导数时的误差为 用空间中心差商代替空间导数时的误差为 (2)对流方程与对应的差分方程之间的误差为： 这也可由Taylor展开得到。因为： 2、截断误差 与时间相关的物理问题,必须给定初始条件,形成一个完整的初值问题。对流方程的初值问题为 差分方程必须有初始条件： 3.初值条件 4.收敛性与稳定性 一个差分格式能否在实际中使用，最终要看能否任意地逼近微分方程的解。这样对于每一个差分格式，人们便从两个方面加以考虑：一是引入收敛性的概念，考察差分格式在理论上的准确解能否任意逼近微分方程的解；二是引入稳定性的概念，考察差分格式在实际计算中的近似解能否任意逼近差分方程的解。 (1)收敛性 (2) 稳定性 1.离散求解域； 2.用时间向前差分和空间中间差分格式代替控 制方程的对应项； 3.差分化; 4.选择适当的计算方法求解线性代数方程组； 5.将求解结果用云图、等值线、动画等方式展 示出来，以供实际应用参考。 二 有限差分法求解应用问题的一般步骤 §4 差分法求解Laplace方程 一.Laplace方程的差分表示 二. Laplace方程在流体力学中的物理意义 不可压流体、稳定、无旋、无源流体，速度势函数满足Laplace方程 三. Laplace方程的差分解法 由于区域不规则,不能用前面方法解 Laplace方程的差分解法 Laplace方程的差分解法 Laplace方程的差分解法 Laplace方程的差分解法 Laplace方程的差分解法 Laplace方程的差分解法 例: 用差分方法求下列定解问题的近似解 1. 作平行于坐标轴、间隔为1的两族直线, 根据边界条件写出边界上各节点的函数值 解: U值图 x y 2.由Laplace差分方程可得: U值图 x y 3.用消去法或迭代法解方程组 同步迭代法:①给定上式中每个未知数Uij以一个值作为其零次近似值,如: 边界点上值的平均值 ②将平均值代入方程的右边即可得到解的一次近似值. ③ 比较第k次迭代值与第k+1次的迭代值, 若小数点后三位数字相同,就可取Uij(k+1)作为原定解问题的近似值 (2) 异步迭代法: 计算U2,1(k+1)时就可以利用已算的U1,1(k+1)的,同理,计算U3,1(k+1)时就可以利用已算的U2,1(k+1) §5 差分法求解热传导方程 考虑下列定解问题： 初始边界已知,且: 解: 1.在xot平面上用二族平行于坐标轴的直线来分割D 将[0,l]等份为N份, [0,t]等份为J份 x t 解: 2. 求出矩形区域上左右边界、底边的节点上的函数值 解: 3. 将热传导方程转化为差分方程 解: 4. 解热传导方程转化为的差分方程 由上式可以解出 式中: 解: 5. 解的稳定性讨论 稳定的条件: 可以证明: 当 时,差分方程是稳定的. 习题2: 写出下列波动方程的差分格式: 解: (1) 作平行于坐标轴的的两族直线, 间隔: (3) 差分方程为: (2) 设: (4) 边界条件为: 习题3: 设区域Ω是边长为1,中心在原点的正方形,用差分解法(取步长h=0.1)求Laplace方程的解的第一次近似值Ui,j: (取零次近似值为 ,采用同步迭代法) 解:(1)作平行于坐标轴的间隔为0.1 的两族直线,根据边界条件写出边界上各点的函数值,域内各节点(ih, jh)的函数值记为: (2)由Laplace的差分方程可知, 域内每个节点的函数值为其上、下、左、右四个相邻节点的函数值的平均值，且已知其零次近似值为0。 (3)用同步迭代法解得: 习题4: 热传导方程第一边值问题： 求: 题中取: 第七章 偏微分方程数值解(差分法) Finite difference method (FDM) §7.1 数值天气预报概述 §7.2 差分原理及逼近