数值分析第八节.pdfVIP

数值分析 Numerical Analysis 7 李 元庆 东南大学数学系 2012 年 8 月29 日 李元庆 (东南大学数学系) 数值分析 Numerical Analysis 7 2012 年 8 月29 日 1 / 68 偏微分方程 （PDEs） 偏微分方程 （PDEs）：含有多元未知函数偏导数的方程。 描述自然界连续变化的重要的PDEs： 1 Maxwell’s equations （电磁场领域）； 2 Navier-Stokes equations （流体力学领域）； 3 Linear elasticity equations （弹性力学）； 4 ¨ Schrodinger’s equations （量子力学）； 5 Einstein’s equation （广义相对论）； 6 等等。。。 李元庆 (东南大学数学系) 数值分析 Numerical Analysis 7 2012 年 8 月29 日 2 / 68 偏微分方程的求解 相 比常微分方程，偏微分的求解更为困难 。 原因有很多，比如说： 求解问题的区域更加不规则 初边值条件更为复杂等。 除少数几种特殊情况，求出偏微分方程的精确解是很困难 的。 即使是数值求解，PDEs的求解依然是当今的热门领域 。 李元庆 (东南大学数学系) 数值分析 Numerical Analysis 7 2012 年 8 月29 日 3 / 68 偏微分方程数值求解的热门方法 有限元方法； 边界元方法； 有限差分方法； 基于积分方程的快速多极算法； 区域分解算法； 等等。。 李元庆 (东南大学数学系) 数值分析 Numerical Analysis 7 2012 年 8 月29 日 4 / 68 三类常见的偏微分方程 二阶常系数偏微分方程。 一般形式： auxx + buxy + cuyy + dux + euy + fu + g = 0 根据 b2 4ac 符号的不同，二阶常系数偏微分方程常常分为三类 。 李元庆 (东南大学数学系) 数值分析 Numerical Analysis 7 2012 年 8 月29 日 5 / 68 三类常见的偏微分方程：椭 圆型 椭圆型方程：如 2 维Poisson方程 2 u + 2 u = f (x y ) (x y ) Ω R2 x 2 y 2 u(x y ) = (x y ) Ω 此时，a = c = 1 b = 0 b2 4ac = 4 0 。 简单地说，椭 圆型方程描述的是已经达到了稳定状态的系统 （或者说是均衡系统），因而同时间无关。 典型代表：位势方程 （描述电场、引力场和流场等物理现象）。 李元庆 (东南大学数学系) 数值分析 Numerical Analysis 7 2012 年 8 月29 日 6 / 68 三类常见的偏微分方程：抛物型 抛物型方程：如 u 2 u = f (x t) x (ab) t (0 T) t x 2 u(x 0) = (x ) x [ab] u(a t) = (t) u(b t) = (t) t [0 T] 此时，