【毕业论文】对流占有的特征差分算法.doc

PAGE PAGE 13 目 录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc263170613 前 言 PAGEREF _Toc263170613 \h 1 HYPERLINK \l _Toc263170614 第一章 有限差分法介绍 PAGEREF _Toc263170614 \h 2 HYPERLINK \l _Toc263170615 §1.1 有限差分法基本知识 PAGEREF _Toc263170615 \h 2 HYPERLINK \l _Toc263170616 §1.2 对流占优扩散方程常用的几个差分格式 PAGEREF _Toc263170616 \h 6 HYPERLINK \l _Toc263170617 §1.2.1 中心显示格式 PAGEREF _Toc263170617 \h 6 HYPERLINK \l _Toc263170618 §1.2.1 迎风格式 PAGEREF _Toc263170618 \h 7 HYPERLINK \l _Toc263170619 第二章 特征差分格式 PAGEREF _Toc263170619 \h 9 HYPERLINK \l _Toc263170620 §2.1 对流方程的特征线理论 PAGEREF _Toc263170620 \h 9 HYPERLINK \l _Toc263170621 §2.2 基于线性插值的一种特征差分格式 PAGEREF _Toc263170621 \h 12 HYPERLINK \l _Toc263170622 §2.2.1 的情形 PAGEREF _Toc263170622 \h 13 HYPERLINK \l _Toc263170623 §2.2.2 的情形 PAGEREF _Toc263170623 \h 17 HYPERLINK \l _Toc263170624 §2.2.3 的情形 PAGEREF _Toc263170624 \h 17 HYPERLINK \l _Toc263170625 §2.2.4 的情形 PAGEREF _Toc263170625 \h 21 HYPERLINK \l _Toc263170626 §2.2.5 小结 PAGEREF _Toc263170626 \h 21 HYPERLINK \l _Toc263170627 §2.3 基于双线性插值的改进特征差分格式 PAGEREF _Toc263170627 \h 22 HYPERLINK \l _Toc263170628 §2.3.1 算法构造 PAGEREF _Toc263170628 \h 23 HYPERLINK \l _Toc263170629 §2.3.2 双线性插值 PAGEREF _Toc263170629 \h 25 HYPERLINK \l _Toc263170630 第三章 算例结果 PAGEREF _Toc263170630 \h 27 HYPERLINK \l _Toc263170631 结 论 PAGEREF _Toc263170631 \h 31 HYPERLINK \l _Toc263170632 参考文献 PAGEREF _Toc263170632 \h 32 HYPERLINK \l _Toc263170633 致 谢 PAGEREF _Toc263170633 \h 35 HYPERLINK \l _Toc263170634 附 录 PAGEREF _Toc263170634 \h 36 前 言 随着现代工业的发展和人类活动的加剧，污染物以及悬浮颗粒排入大气、江河、湖泊和海洋中，它们在大气及水体中的变化过程都可以概括为对流扩散问题的各种数学模型。欲刻画这类流体运动的物理现象，就要求解对流扩散方程。差分方法是数值求解的一种经典方法。 常用的算法主要有迎风方法，中心型差分方法以及沿着特征线离散的特征类方法。迎风方法由于在一阶精度的差分格式，所以常出现数值耗散现象。而中心型Crank-Nicolson差分格式关于时间以及空间均有二阶的精度。但是该格式形成的系数矩阵并不一定能满足对角占优条件。 另外，对于比较特殊的对流占优问题，传统的迎风格式或者中心差分格式，会出现数值扩散或数值震荡．因为这时候，方程呈现出双曲型方程的基本特性，因此．需要构造稳定性好，适用于小扩散系数，且基本上能反映双曲型方程的特征性质的数值解法．特征差分方法正是这样一类方法。八十年代Douglas和Russell[1]提出了求解对流占优扩散问题的特征有限差分法．该算法在一定程度上克服了因为对流占优项而引起的数值震荡．其后，这