关于对流扩散方程的一些花样的探讨.pdfVIP

摘 要 本文讨论了以下一维对流扩散初边值问题数值求解的几个差分格式 0掰L0￡了1 鬻+6器=a貉+m朋 u(。，0)=，(z) 0oL u(o，t)=9l(f) 0fr U(L，f)=卯(￡) 0￡丁 先是基于中心型差分格式的探讨．第一章介绍了最简单的显式格式以及隐式 格式，并且在此基础上推导得常用的Crallk．Nicols咖型差分格式(六点隐格式，简记 为lygs)．第二章旨在克服中心型差分格式的缺点之一：系数矩阵并不总能满足对角 占优．在Crank．Njcolson型差分格式的基础上，提出了两种适应于变系数问题的改 进格式(罚格式以及迎风格式，分别简记为f西，y龟s)以保证其所形成的系数矩阱：对 角元占优，可以更加有效的求解．对改进格式，分别讨论了截断误差以及精度．并且 对系数给出一定的限制条件，以方便用能量估计法来证明差分解按照离散如范数 收敛于Sobolev空间内的精确解． 一般来说，应用中心型差分格式进行求解，结果是令人满意的．但是对对流占优 的对流扩散方程，如果利用中心型差分格式(以及其他一般的差分格式)进行求解干丁 时会出现数值震荡．而特征有限差分方法则可以有效地克服这一缺点第三章对基 于线性插值的特征差分格式进行了改进．采用了文【19}的利用最近点插值的思想， 参照文f1 51对边界点的处理办法，将原来的基于线性插值的特征差分方法予以改进。 使得算法更加简便，并且从理论上证明了格式的收敛性：差分解按照最大模范数收 敛于Sobolev空间内的精确解． 在理论分析之后给出了具体算例，证实了六点隐格式(1ygs)的精度最高，迎风杆 式(yf98)次之，罚格式(fgS)在精度上最低．对于一般的问题，迎风格式比罚格式更好． 因为前者比后者的精度更加高，适用范围也更加广．但是对于对流占优的陡峭前R} 问题，利用迎风格式或六点隐格式求解都发生了数值震荡，而罚格式则可以避免产 生数值震荡，虽然有一点数值扩散的迹象，仍不失为一个适用于对流占优问题的算 法．至于第三章所提出的基于线性插值豹改进特征差分格式，由于其构造简单，而 且能有效消除数值震荡现象，适用于求解对流占优的对流扩散方程． 关键词：Crar止Nicol80n型差分格式；数值震荡；线性插值；特征差分格式 ABSTRACT jnl、一 This discusses80medifrerencescllemesfororledimensionalconwrl paper djffusion asfollow8： initial_boundaryprobkms otr 鬻+6鬻一n器+脚，t)ozL u(霉，O)=，(o) OS茗≤￡ U(0，t)=91(￡) O≤tS丁 c，(L，f)=出(￡) 0冬t≤r We central withthe dj跽renceschemes． withthe begin Combiningsiml)1e：=i1 8clleme and usu出Crank-Nicolsondj行erence expIicit implicitscheme，the schelne(sjx． isdenotedto for deductedin pointimplici