两点第三边值问题的非均匀网格四阶有限体积方法.pdf

thefourth-orderfinite volumeschemeonnon-uniformmesh． value ofthird Keywords：two-pointboundaryproblem volumemethodwith kind；finite nonuniform mesh；Richardson extrapolation order technique；convergenceanalysis；fourth accuracy III 目录 第一章引言…………………………………………………………………．1 第二章两点第三边值问题非均匀网格二阶有限体积方法的外推……………………．3 §2．1二阶有限体积格式的构造及收敛性分析…………………………………．3 §2．2外推辅助方程构造及收敛性分析………………………………………．10 §2．3数值算例……………………………………………………………17 第三章两点第三边值问题非均匀网格四阶精度有限体积方法……………………．．19 §3．1四阶非均匀网格有限体积格式的建立…………………………………．．19 §3．2收敛性分析…………………………………………………………21 §3．3导数恢复……………………………………………………………31 §3．4数值算例……………………………………………………………36 参考文献……………………………………………………………………．40 致谢…………………………………………………………………………43 IV 第一章引言 偏微分方程数值解法在数值分析中占有重要地位，很多科学问题的数值计算需要通过求 解偏微分方程来实现．有限差分，有限元和有限体积是三种常用的方法．有限体积方法It-5]， 亦被称作广义差分法(GDMs)【5一，从微分方程的积分守恒形式出发，通过选取有限元子空 间作为试探函数空间，分片常数空间作为检验函数空间来离散方程，进而导出计算格式．该 方法既最大限度地保持了差分方法的简单性，又兼有有限元法的精确性，且能保持物理量的 局部守恒性，已被广泛应用于计算流体力学中． 在过去的几十年里，在运用高精度紧差分格式求解偏微分方程问题上有了很大的进步， 条件的处理大概被分为三类，第一类为用通过在定义域两边引入虚点的办法使用中心差分离 于这两种方法，要得到高精度必要的先决条件为方程在边界处成立．我们知道这一假设并不 是对所有方程都适用．为了避免这种假设，第三种方法为在靠近实际点离散边界条件[19-22】． 关于格式的收敛性分析，文献【23]使用能量方法研究格式的收敛性且证明结果显示文献[19】 三阶精度而文献【14】和【21]中的格式有四阶精度．这告诉我们有关紧差分格式的理论分析 还可以进一步完善． 目前已有的高精度紧致格式大部分都是在均匀网格上提出的110-23]，然而在实际问题计 算中，经常遇到待求物理量如浓度，温度等急剧变化或空间分布不均匀，大梯度问题或者局 部奇性问题，使用均匀网格格式计算精度受到很大影响．比较合理的做法是在大梯度或边界 层区域内多布置些计算节点而在小梯度或物理量变化比较平缓的区域内少布置一些计算节 点，这样既可以兼顾算法的稳定性和计算结果的精度性又可以节省计算量．因此，发展非均 匀网格上的高精度紧致格式就具有重要的理论意义和实际应用价值． 本文针对两点第三边值问题使用非均匀网格提出两种方法使其精度达到四阶．第二章针 对两点边值问题非均匀网格二阶有限体积方法进行外推，第三章直接构造四阶紧有限体积格 式，并给出了节点导数值的高