第8章 偏微分方程数值解.ppt

第8章 偏微分方程数值解;2. 抛物型方程;3．双曲型方程;二、定解问题;双曲型方程初始条件表示弦在两端振动规律为已知：;Poisson方程反应稳定状态的情况，与时间无关，所以不需要提初始条件。边界条件的提法为： 其中? (x, y)为已知边界，s是区域D的边界。 ;本章主要针对几个典型的微分方程介绍常用的差分方法和有限元方法。 这些方法基本思想是： 把一个连续问题离散化 通过各种手法化成有限形式的线性方程组 然后求其解。 ;计算机只能作有限次的加、减、乘、除运算，它既不能求导数，更不能解偏微分方程 如果想在计算机上求得微分方程数值解，它的主要做法是 把偏微分方程中所有的偏导数分别用差商代替 从而得到一个代数方程组——差分方程组 然后对差分方程求解，并以所求的解作为偏微分方程数值解; 因此需要对区域进行剖分，用网格点来代替连续区域，所以差分法亦称“网格法”。; 把整体分割成若干个单元来处理问题的方法在数学上称为“离散化方法” 。 在结点上采用离散化方法（数值微分、数值积分、泰勒展开等）将微分方程的初边值问题化成关于离散变量的相应问题，这个相应问题的解就是方程在点xi上的数值解f(x)，或在点（xi , ti）上的数值解U( xi , ti)。 一般来说，不同的离散化导致不同的方法。;例：取一边长为1的正方形均匀薄板，上下侧面绝热，四周保持恒温，求板内各点的稳定温定分布。;记;得到u (i, k)的近似ui,k，所满足的线性代数方程组： ;简单迭代法;表8.1;0;0;用差分法解偏微分方程需要考虑三个问题： ;8.2 椭圆型方程的差分解法 ;考虑泊松方程第一边值问题：;(一) 矩形网格;（二）五点差分格式;则差分方程可简写成： ;;这四个式子两两相加便有： ;于是可得差分方程的截断误差 ;（三）边值条件的处理 ;（2）线性插值法;则u在这些点上的值有近似关系： ;（3）列不等距差分方程;8.3 抛物型方程的差分解法 ;现以热传导方程为例，介绍抛物型方程的有限差分格式。;8.3.1 矩形网格;8.3.2．古典差分格式;于结点（j, k）处偏导数与差商之间有如下近似的关系：;视为 u (xj, tk) 的近似值。 ;将（6）写成便于计算的格式：;为了提???截断误差的阶，可以利用中心差商：;截断误差阶为o (?2 + h2)，较古典显格式高。;如果利用向后差商 ;8.3.3 六点对称格式;将以上各式代入（3）式得到差分方程： ;整理,得 ;8.3.4 稳定性;考察 Richardson 格式的稳定性。 ; r = ? 时 Richardson 格式的误差传播; r≤ 1/2 时古典显格式的误差传播; 差分格式关于初值稳定的实际含义是：如果其解 在某一层存在误差，则由它引起的以后各层上的误差不 超过原始误差的M倍（M为与? 无关的常数）。;8.4 双曲型方程的差分解法 ;也就是说，在平面 x t上，沿着;在物理上常见的双曲型偏微分方程最简单模型是波动方程 ;8.4.1 矩形网格;a）迎风格式 ;令r= ? / h，得;稳定性的讨论：;当a0时，恒有 ，格式（8.4.7）不稳定 ；;b) Lax-Friedrichs格式;传播因子 ;c）Lax-Wendroff格式;d）古典隐式格式;e）Grank-Nicholson格式;8.4.3 二阶双曲型方程的差分格式;b）隐格式;考试就要到了*******