偏微分方程的有限差分法.pptx

单击此处编辑母版标题样式,Harbin Institute of Technology Yangkun kyang,#,/75,计,算,物,理,学,单击此处编辑母版文本样式,第二级,Harbin Institute of Technology Yangkun kyang,1,/75,第四章 偏微分方程的有限差分法,4.1,有限差分法原理,4.2 热传导方程的差分解法,4.3,波动方程的,差分解法,Harbin Institute of Technology Yangkun kyang,2,/75,4.1,有限差分法原理,抛物线形,双曲型,椭圆形,不可逆过程,可逆过程,平衡过程,热传导方程,波动方程,位势方程,物理学中许多物理规律都用偏微分方程描述，偏微分方程主要分为以下三类：,上式中,a,c,f,以及未知函数,u,为定义在求解区域上的实（复）函数,Harbin Institute of Technology Yangkun kyang,3,/75,4.1,有限差分法原理,有限差分,解法,差分近似代替微分，差商近似代替微商,这样就把求解区域内连续分布函数离散化成求网络节点上的分立函数值，从而把所需求解的,微分方程,变为一组相应的,差分方程,，进一步可以求解离散节点上的函数值。,数学基础,泰勒（,Taylor,）展开,Harbin Institute of Technology Yangkun kyang,4,/75,4.1,有限差分法原理,差商公式的构造,利用泰勒级数展开定义差商,Harbin Institute of Technology Yangkun kyang,5,/75,4.1,有限差分法原理,误差为,O(h),差商公式：,一阶向前差商：,一阶向后差商：,Harbin Institute of Technology Yangkun kyang,6,/75,4.1,有限差分法原理,二阶,向前差商：,式,(2)-,式,(1)X2,误差为,O(h,2,),差商公式：,Harbin Institute of Technology Yangkun kyang,7,/75,4.1,有限差分法原理,二阶,向后差商,:,式,(2)-,式,(1)X2,Harbin Institute of Technology Yangkun kyang,8,/75,4.1,有限差分法原理,一阶向前差商：,Harbin Institute of Technology Yangkun kyang,9,/75,4.1,有限差分法原理,一阶向后差商：,Harbin Institute of Technology Yangkun kyang,10,/75,4.1,有限差分法原理,一阶中心差商：,Harbin Institute of Technology Yangkun kyang,11,/75,4.1,有限差分法原理,二阶中心差商：,Harbin Institute of Technology Yangkun kyang,12,/75,4.1,有限差分法原理,差分格式的收敛性和稳定性,收敛性：,稳定性：,当步长,h,0,时，差分方程的解趋向于微分方程的解。,误差在运算过程中不会失控，即累计误差不会无限增加。,Harbin Institute of Technology Yangkun kyang,13,/75,4.1,有限差分法原理,从数学上讲，没有限制的微分方程会有无穷多个解，不能构成一个定解问题。,从物理上讲，描述物理问题的微分方程仅适用于描述在一个连续体或物理场的内部发生的物理过程，仅靠这些微分方程不足以确定物理过程的,具体,特征。,因此，要想解决实际的物理问题，必须知道一个连续体或物理场的初始状态和边界受到的外界影响。,Harbin Institute of Technology Yangkun kyang,14,/75,4.1,有限差分法原理,初始条件：,与时间相联系,边界条件：,边界,受到外界的影响,偏微分方程的定解条件,常见的物理问题可以归结为三大类边界条件,Harbin Institute of Technology Yangkun kyang,15,/75,4.1,有限差分法原理,2,第二类边界条件,（诺依曼,Neumann,）,1,第一类边界条件,（狄利克雷,Dirichlet,）,热传导,问题：,边界,上,温度分布已知,热传导,问题：,通过,边界,单位面积上的热流量已知,n,表示,的外法线,q,0,定义在,上的已知函数,Harbin Institute of Technology Yangkun kyang,16,/75,4.1,有限差分法原理,由热力学傅立叶定律得：,热流量：,单位面积上的,热流量：,K:,热传导系数,单位时间内通过给定截面的热量，正比例于垂直于该界面方向上的温度变化率和截面面积，而热量传递的方向则与温度升高的方向相反。,Harbin Institute of Technology Yangkun kyang,17,/75,4.1,有限差分法原理,3,第三类边界条件,（洛平,Robin,）,热传导,问题：,边界表面,与外界之间的热量交换已知,a,0,b,0,.c,0,定义在,上的已知函数,外界温度为,u,0,，热交换规律遵循,热传导实验定律,：,单位时间内，从边界单位面积传递给周围的热流量正比于边界表面和外界的温度差。,Harbin Institute of Technology Yangkun kyang,18,/75,4.1,有限差分法原理,对于实际物理问题，边界条件往往是很复杂的，可能是一种或不同边界区域几种边界条件的组合，甚至不能用这三类边界条件描述。,：热交换系数,u,：边界温度,单位面积上的,热流量：,Harbin Institute of Technology Yangkun kyang,19,/75,4.2,热传导方程的差分解法,物理学中对热传导、热辐射以及气体扩散现象的描述，常可以归结为同一类型的抛物线型方程，通常采用二阶偏微分方程描述，这类方程统称为热传导方程。,Harbin Institute of Technology Yangkun kyang,20,/75,4.2,热传导方程的差分解法,一维各向同性、均匀介质，且无热源的热传导方程：,为了求解u,(x,t),，还必须利用边界条件和初始条件。,定解条件,：,边界条件,和,初始条件,。,定解问题,：解存在、唯一并且连续依赖初始条件。,4.2.1一维热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,21,/75,4.2,热传导方程的差分解法,对于一维热传导问题,（第一类边界条件）,数值解就是在求解区域,中某些,离散点（,x,i,t,i,）,上求出,u(,x,i,t,i,),足够近似的解。,Harbin Institute of Technology Yangkun kyang,22,/75,4.2,热传导方程的差分解法,1,把求解区域离散化（确定离散点）,T,l,Harbin Institute of Technology Yangkun kyang,23,/75,4.2,热传导方程的差分解法,2,推导差分递推公式,在节点（,x,i,t,k,）上,二阶,向前差商,O(h,2,),Harbin Institute of Technology Yangkun kyang,24,/75,4.2,热传导方程的差分解法,同样，在节点（,x,i,t,k,）上,一,阶,向前差商,O(h),Harbin Institute of Technology Yangkun kyang,25,/75,4.2,热传导方程的差分解法,一维热传导方程可以近似为,令,O(h),Harbin Institute of Technology Yangkun kyang,26,/75,4.2,热传导方程的差分解法,初始条件,边界条件,Harbin Institute of Technology Yangkun kyang,27,/75,4.2,热传导方程的差分解法,一维热传导方程显示差分递推公式为,Harbin Institute of Technology Yangkun kyang,28,/75,4.2,热传导方程的差分解法,显示差分递推公式的稳定性：,对于一维热传导方程，差分格式为稳定差分格式的充分条件是,即,Harbin Institute of Technology Yangkun kyang,29,/75,4.2,热传导方程的差分解法,为了提高数值解的精度，必须减小,相应就要变小，这必然增加计算量。这就是显示差分格式的缺点，但它的优点是计算简单。,Harbin Institute of Technology Yangkun kyang,30,/75,4.2,热传导方程的差分解法,差分格式计算步骤：,给定,计算,计算初始值：,计算边界值：,用差分格式计算,Harbin Institute of Technology Yangkun kyang,31,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,32,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,33,/75,4.2,热传导方程的差分解法,X=diag(v,k),若,v,为,n,个元素向量，返回一个阶数为,n+abs(k),的方阵,X,，将,v,作为,X,的第,k,个对角元，,k=0,代表主对角元，,k0,表示在主对角元之上，,k0,表示在主对角元以下。,v=ones(1,5);,X1=diag(v),X2=diag(v,1),X3=diag(v,-1),Harbin Institute of Technology Yangkun kyang,34,/75,4.2,热传导方程的差分解法,（,1-2,*,）*,diag(ones(1,N-1)+,*（,diag(ones(1,N-2),1)+diag(ones(1,N-2),-1),）,Harbin Institute of Technology Yangkun kyang,35,/75,4.2,热传导方程的差分解法,例,4.2.1,求热传导方程混合问题,Harbin Institute of Technology Yangkun kyang,36,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,37,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,38,/75,9.2,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,39,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,40,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,41,/75,4.2,热传导方程的差分解法,二维各向同性、均匀介质，且无热源的热传导方程：,初始条件：,边界条件：？,4.2.2,二,维热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,42,/75,4.2,热传导方程的差分解法,同一维类似，把求解区域离散化,Harbin Institute of Technology Yangkun kyang,43,/75,4.2,热传导方程的差分解法,在节点（,x,i,y,j,t,k,）上,Harbin Institute of Technology Yangkun kyang,44,/75,4.2,热传导方程的差分解法,在节点（,x,i,y,j,t,k,）上,Harbin Institute of Technology Yangkun kyang,45,/75,4.2,热传导方程的差分解法,二维热传导方程可以近似为,令,Harbin Institute of Technology Yangkun kyang,46,/75,4.2,热传导方程的差分解法,差分递推公式为,Harbin Institute of Technology Yangkun kyang,47,/75,4.2,热传导方程的差分解法,边界条件：,一、右图中阴影部分为绝热壁，此绝热壁可以用第二类边界条件描述，即热流量为零。,第二类边界条件：,通过,边界表面单位面积上的热流量已知,Harbin Institute of Technology Yangkun kyang,48,/75,4.2,热传导方程的差分解法,差分近似为,Harbin Institute of Technology Yangkun kyang,49,/75,4.2,热传导方程的差分解法,递推公式为：,Harbin Institute of Technology Yangkun kyang,50,/75,4.2,热传导方程的差分解法,二、,i=0,边界，,M,1,jM,2,区域为与高温恒温热源相连接的口，温度可取归一化值,1,。,j=0,和,j=M,边界与低温恒温热源相连，温度始终为,0,。,Harbin Institute of Technology Yangkun kyang,51,/75,4.2,热传导方程的差分解法,二维热传导方程显示差分递推公式为,Harbin Institute of Technology Yangkun kyang,52,/75,4.2,热传导方程的差分解法,稳定差分格式的充分条件是,即,Harbin Institute of Technology Yangkun kyang,53,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,54,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,55,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,56,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,57,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,58,/75,4.2,热传导方程的差分解法,view(az,el),az:,方位角,el:,仰角,view(0,90),view(-37.5,30,）,Harbin Institute of Technology Yangkun kyang,59,/75,4.2,热传导方程的差分解法,Harbin Institute of Technology Yangkun kyang,60,/75,4.3,波动方程的差分解法,一维均匀弦线的自由振动方程：,波动方程的差分解法也利用构造网格节点的方法,Harbin Institute of Technology Yangkun kyang,61,/75,4.3,波动方程的差分解法,用二阶中心差分近似方法得：,Harbin Institute of Technology Yangkun kyang,62,/75,4.3,波动方程的差分解法,代入波动方程得：,令：,Harbin Institute of Technology Yangkun kyang,63,/75,4.3,波动方程的差分解法,边界条件,初始条件,Harbin Institute of Technology Yangkun kyang,64,/75,4.3,波动方程的差分解法,a,初始条件：,对于初始时刻速度，也须用差分格式给出：,Harbin Institute of Technology Yangkun kyang,65,/75,4.3,波动方程的差分解法,a1,向前差分：,误差：,Harbin Institute of Technology Yangkun kyang,66,/75,4.3,波动方程的差分解法,a2,中心差分：,由,得（,k=0,）,+,Harbin Institute of Technology Yangkun kyang,67,/75,4.3,波动方程的差分解法,误差：,整理得：,Harbin Institute of Technology Yangkun kyang,68,/75,4.3,波动方程的差分解法,b,边界条件：,一维波动方程的差分格式有如下两种,Harbin Institute of Technology Yangkun kyang,69,/75,4.3,波动方程的差分解法,第一种：,误差：,Harbin Institute of Technology Yangkun kyang,70,/75,4.3,波动方程的差分解法,第二种：,误差：,Harbin Institute of Technology Yangkun kyang,71,/75,4.3,波动方程的差分解法,第二种差分格式精度要高于第一种，是经常采用的方法。,理论上可以证明，两种差分格式稳定条件是：,Harbin Institute of Technology Yangkun kyang,72,/75,4.3,波动方程的差分解法,波动方程差分格式的计算步骤如下：,;,计算,;,计算,;,计算初值和边值,;,给定,计算,Harbin Institute of Technology Yangkun kyang,73,/75,4.3,波动方程的差分解法,例：计算下列一维波动方程,Harbin Institute of Technology Yangkun kyang,74,/75,4.3,波动方程的差分解法,Harbin Institute of Technology Yangkun kyang,75,/75,4.3,波动方程的差分解法,Harbin Institute of Technology Yangkun kyang,76,/75,4.3,波动方程的差分解法,Harbin Institute of Technology Yangkun kyang,77,/75,4.3,波动方程的差分解法,Harbin Institute of Technology Yangkun kyang,78,/75,4.3,波动方程的差分解法,Harbin Institute of Technology Yangkun kyang,79,/75,4.3,波动方程的差分解法,Harbin Institute of Technology Yangkun kyang,80,/75,上机,4,编程计算下列一维扩散方程的解,要求：,1,推导差分递推公式,2,图形显示计算结果，并与解析解比较。,a x-t-u,三维曲线,b t=0 0.5 1,时刻,x-u,曲线,c x=0.5,处,t-u,曲线,