求解热传导方程的高精度隐式差分格式毕业论文.docx

新疆大学毕业论文（设计）XINJIANGUNIVERSITY题目：求解热传导方程的高精度隐式差分格式所属院系：数学与系统科学学院专业：信息与计算科学声明本人郑重声明该毕业论文（设计）是本人在开依沙尔老师指导下独立完成的,本人拥有自主知识产权，没有抄袭、剽窃他人成果，由此造成的知识产权纠纷由本人负责。声明人（签名）：年月曰亚库甫江.买买提同学在指导老师的指导下，按照任务书的内容，独立完成了该毕业论文（设计），指导教师已经详细审阅该毕业论文（设计）。指导教师（签名）：年月日新疆大学毕业论文（设计）任务书班级：信计07-2姓名：亚库甫江.买买提论文（设计）题目：求解热传导方程的高精度隐式差分格式专题：毕业设计论文（设计）来源：教师自拟要求完成的内容：学习和掌握一维热传导方程已有的各种差分格式的基础上，扩散方程对空间变量应用紧致格式离散，对时间变量应用梯形方法，构造热传导方程的精度为+A4）数值格式，讨论格式的稳定性，最后数值例子来验证。发题日期：2012年12月25日完成日期:2012年5月28日实习实训单位：数学学院地点：数学学院论文页数：19页:图纸张数：4指导教师：开依沙尔老师教研室主任院长（系主任）摘要本文首先对热传导方程经典差分格式进行复习和讨论，然后热传导方程对空间变量四阶紧致格式进行离散,时间变量保持不变，把一维热传导方程转化为常微分方程组的初值问题，再利用梯形方法构造热传导方程方程的时间二阶空间四阶精度的一种差分格式，并稳定性进行分析，数值结果与Crank-Nicho1son格式进行比较，数值结果表明，该方法是有效求解热传导方程的数值计算.关键词：热传导方程，高精度紧致格式；梯形方法；两层隐格式；Crank-Nico1son格式ABSTRACTThispaperfirststudyonsomec1assica1finitedifferencefortheheatconductionequation,seconde1yseconde1yweapp1ycompactfinitedifferenceapproximationoffourthorderfordiscretizingspatia1derivativesbut1eavethetimevariab1eContinuous.Thisapproachresu1tsinasystemofODEs,whichcanthenbeusedtrapezodia1formu1aderivedfourthorderinspaceandsecondorderintimeunconditiona11ystab1eimp1icitscheme.thestabi1ityand1oca1truncationerroroftheobtainedmethodareana1ysied.Numerica1experimentsshowsthatthismethodUsefu1,efficientmethodforso1vingdiffusionequationKeywords:Heatconductioneqution;Higher-odercompactscheme;Trapezodia1formu1a;Two-1eve1imp1ictscheme;Crank-Nico1sonscheme引言1预备知识21 .扩散方程的经典差分格式31.1 显式差分格31.1.1 显式的截断误差41.1.2 显式差分格式的稳定性41.2 隐式差分格式51.2.1 隐式差分格式的截断误差51.2.2 隐式差分格式的稳定性61.3 Crank-NicoIson格式61.3.1 Crank-NicoIson差分格式的截断误差71.3.2 Crank-NicoIson差分格式的稳定性82 .高精度格式的构造92.1 梯形方法92.2 本文格式的构造102.3 稳定性分析113 .数值实验13结论17致谢184 考文献19引言热传导方程是一类描述物理量随时间的扩散和衰减规律的抛物型微分方程.自然环境、工程设备及生物机体中的许多物理现象，诸如气体的扩散、液体的渗透、热的传导、以及半导体材料中杂质的扩散等都可用热传导方程来描述.由于物理问题本身的复杂性，其精确解往往不容易求得，因此研究其数值求解方法无疑具有非常重要的理论意义和工程应用价值.求解该问题的数值方法主要有差分法、有限元法、边界元法等，其中有限差分方法数值求解扩散方程的应用广泛的有效地方法之一。目前求解该问题的主要的差分格式有显式格式，隐式格式，Crank-Nico1son格式等值。虽然显式格式计算简单，但是稳定性有所限制，一般隐式格式和Crank-Nic。ISon格式分别为一阶和二阶精度的绝对稳定的隐式格式，还显得误差阶不够高，得到的结果也往往不能令人满意，考虑到这些不足文7中半离散方法构造0(，+/?)格式结果Crank-Nico1son格式进行比较,在文10待定参数法构造精度。卜、/“)的显式格式但是稳定性条件比较苛刻，它文的稳定性条件为旨,本文热传导方程对空间变量应用紧致格式离散，对时间变量应用梯形方法，构造热传导方程的精度为0(/+/)的绝对稳定的隐式差分格式，并讨论了稳定性，数值值结果与经典Crank-NiChO1SOn格式进行比较，数值结果表明，该方法是有效求解扩散方程的数值计算.本文分为三大部分，第一部分简单介绍热传导方程的经典差分格式，第二部分主要介绍热传导方程的高精度格式的构造和稳定性，第三部分给出具体的数值算例，结果与Crank-NiCO1SOn格式，准确值进行比较，最后给出结论。预备知识利用下面的各种数值微分公式得到不同的差分格式一（Xr，乙）=8uTj+h_dx，+（马+|,乙）一（巧|,乙）2h÷（z2）÷（z2）"（x+，乙）一2u（X八乙）+“（xj,乙）截断误差：一般说来，微分方程的解不会精确地满足差分方程。将差分方程中的各个项同时用微分方程的解在相应点的值代入，利用泰勒展开，就会得到一个误差项，这个误差项就是截断误差。相容性：若时间步长r以及空间步长力同时趋于0,截断误差学0,就说差分格式与微分方程是相容的。一个差分格式与一个微分方程相容，则表明当汇,20时，差分算子与微分算子对任一光滑函数的作用是相同的，所以可用相容的差分格式近似相应的微分方程，而截断误差则是对这一近似程度的一个度量。收敛性：考察差分格式在理论上的准确解能否任意逼近微分方程的解。如果当时间步长不以及空间步长力趋于0时，q=N（Xj0,我们称差分格式是收敛的，即时间步长汇以及空间步长力趋于0时，差分格式的解逼近于微分方程的解。稳定性：差分格式的计算是逐层计算的，计算第+1层上的丁时，要用到第层上计算出来的结果。计算；时的舍入误差，必然会影响到片的值，从而就要分析这种误差传播的情况。因此，一个有实用价值的数值方法应该具有能够控制这种误差影响的性能，这就是数值方法的稳定性。精度：如果一个差分格式的截断误差£=。（"+/优），就说差分格式对时间/是4阶精度的，对空间X是P阶精度的。1ax等价定理：给定一个适定的线性初值问题以及与其相容的差分式，则差分格式的稳定性是差分格式收敛性的充分必要条件。定理I(VOnNeI1ma1Ii1条件)微分方程的差分格式稳定的必要条件是当rr0,77对所有kWR有2y(G(r,As)1+Mr，/=1,2,，夕其中G(,Q为增长因子(或增长矩阵)，%(G(,Q)表示G(r,Q的特征值，M为常数。定理2如果差分格式的增长矩阵G«,女)是正规矩阵，则VOnNeUmann条件是差分格式稳定的必要且充分条件。推论2.1当G«,%)为实对称矩阵，酉矩阵，Hermite矩阵时;vonNeumann条件是差分格式稳定的充分必要条件。推论2.2当P=I时，即G(,2)只有一个元素，则VOnNeUmann条件是差分格式稳定的充要条件。定理3如果存在常数K/.使得G(r,Ar)1+K，O<r<r0,则差分格式是稳定的。1.热传导方程的经典差分格式考虑一维热传导方程的初边界问题:u2uCT=°，tOx2M(X,0)=/(x),(")=Q)"S,7)=0(,),a<X<b,t>Oa<x<b0t0我们可以对学用向前差分"(秩=包T÷3，_dtJz用二阶差商x2正得到差分格式为uni+x-uni-ah2一Ui1=-=0(1.1.D1.1.1显式差分格式的截断误差证：(用tay1or展开),.'Wy=u(Xj,t,1)M=u(x9t+)=u(xj,tn)+-,j+2,(0«4«DA33u,xh4r4u+31>+-r1?+?”ZAXZ、，r°“"22U,.uj+="(+)=u(xj,1)+M17+CzV乙CzV(0<921)/A、/、fSu,12。23d%h4rS4Mnw.一=w(x-x)=u(xj,tn)-My÷-2y-÷7r-r1+x2!x3!3x4!x3(0名1)把上述代入差分格式中，得截断误差为：八H23?介2合4h24U%g)=7Jj+5d/-j>÷j+1>÷/”r°2”MI,产°力】+四=f1>1>1+1>=o(+h2).h2,4uhi54mv,11024"W+右1+6(02,1)从上述可知，截断误差为号（XM=。+川），它对空间方向为一阶截断误差而对时间方向为二阶截断误差。112显式差分格式的稳定性证：先把差分格式公式（11D改写为：丁=unj+（uf;+1-2；+；T）其中X=W利用稳定性的Fourier方法，令；=5,并将它代入上式就得到vn+'eikjh=vneikjh+eikj-2÷e-it)vw消去共因子有vw+,=1+eikj-2+e-ikj)vn由此得到增长因子G(r,Ao=1+(eikj-2÷e-ikj)=(1-2)÷eikj÷e-ikj)=(12)÷2coskh=-22(1coskh)khG(r)=1-2(1-cos)H1-4sin1Chkkh=>11-4Asin21=>04sin2222因为k>O,h>O且20,所以必然左边成立，则右边为1rhkh12sin2IKsin21=22<1><-222显然这个格式是相容的。它在2'时稳定的，因为按照1ax定理可知;2它是条件收敛的（收敛条件4'21.2隐式差分格式我们可以噂用向后差分，需用二阶差商，得到差分格式为:unj_u"U-a-h2=0(1.2.1)1.2.1 隐式差分格式的截断误差证：（用tay1or展开）；=”(XM)引Tj(QAOj(“)_吟；+，热片(OWa1)UI4C/fh4f4u不1旅"出(0&1)h44u1zT1病