紧致差分格式求解二维非稳态不可压N-S方程.pdfVIP

卷 期 燕山大学学报 年 月 文章编号：1007-791X (2008) 03-0278-05 紧致差分格式求解二维非稳态不可压N-S 方程 王涛 ，杨丽芳 （ 兰州交通大学数理与软件工程学院，甘肃 兰州 兰州交通大学机电工程学院，甘肃 兰州 ） 摘 要：采用泰勒展式系数匹配的方法构造基于非等距网格的紧致差分格式并得出了它的截断误差。紧致差分 格式能够很好的模拟不同时刻流场的变化情况，网格系统的选择对精度的影响很大，基于非等距网格的紧致差 分方法是一种比经典差分方法精度更高的求解非稳态纳维斯托克斯方程的有效算法。 关键词 紧致有限差分方法； 方程；非等距网格；截断误差 中图分类号：O357.1 文献标识码：A 引言 （如图 ），模型方程为 由于实际流动问题的复杂性和计算技术的高 速发展，人们对流场的计算提出了更高的要求，不 满足于一阶、二阶精度的解，希望得到更准确的流 场解，紧致差分格式是一类高精度的格式，能在更 宽的谐波范围内逼近真实解的谐波成份，更好的反 映流动的小尺度结构 。 对一类 格式 图 非等距网格系统 进行了总结，提出了紧致型差分格式 。早期 方程数值求解 用谱方法 。和谱方法相比，紧致 一阶导数的紧致差分格式 有限差分方法是有限差分方法的一种，有限差分法 具有格式简单，计算量小，对边界条件处理灵活等 设 ′是 的差分近似，内节点的紧致差分格 优点。和经典的差分格式相比，紧致格式能够在相 式可写成如下形式 同的基架点上构造出更高精度的格式，更适用于求 解非稳态的纳维斯托克斯方程。紧致有限差分法是 ′ ′ ′ 用网格点上函数值的线性组合表示网格点上导数 值的线性组合，通过求解方程组获得每个网格点上 根据泰勒展式系数匹配的方法 可得如下关 的导数值 。付德薰和马延文将迎风机制引入紧 系式 致型差分格式，构造了迎风紧致格式，该格式具有 很高的精度外还能对高频非物理震荡进行衰减，抑 制计算中的混淆误差 ，他们还将迎风紧致格式应 用于可压平面混合层到湍流的直接数值模拟，取得 了很好的效 。 构造基于非等距网格的紧致差分格式 设网格为 ，网格间距为 其中 收稿日期： 作者简介：王 涛（ ），男，山东枣庄人。硕士。主要研究方向为偏微分方程数值解。 期 王 涛 等 紧致差分格式求解二维非稳态不可压 方程 空间精度 求解一阶偏导数内节点和边界节点紧致差分 满足边界节点的紧致差分格式可写如下形式 格式的截断误差分别为 ′ ′ 根据泰勒展式系数匹配的方法可得 求解二阶偏导数内节点和边界节点的紧致差 分格式的截断误差分别为 二阶导数的紧致差分格式 设 ′′是 的差分近似，求解二阶偏导数的内 节点紧致差分格式可写成如下形式 ′′ ′′ ′′ 由以上可知，求解一阶偏导数内节点的紧致差 分格式为四阶精度，二阶偏导数