TJU计算流体力学第2章.doc

第2章 计算流体力学的有限差分方法 2．1差分格式 用差分来替代微分、用差商来替代导数是最容易被接受的数学近似概念，随着采用符号的变化,它表示的数学意义也发生了变化，它们之间的关系为或或，下面我们来探求它的数学原理。 设存在连续函数,将它用图线来表示 图2-1 拟序函数序列示意图 作为轴的自变量具有原点、方向和长度单位的特征，依的方向和长度单位可以依次选取,,，并拟序以从小到大排列，那么对应的序列，必然存在,,的函数值。同样将函数值依照的序列作拟序处理得到,,或，那么与的区别在什么地方呢? 表达的仅仅是一个数值，而表达的是函数关系及对应的数值，如,当,,,时 ，而表示的是没有包含函数表达式。一种情况是首先要知道函数形式才能求得函数值，另一种情况是只需知道函数值。当点取得足够致密时，序列的顺序连接可以再现的变化趋势。如果满足某种精度要求，可以认为可以完全反映的变化情况，其中表示的序列集合，由此可以看出表示的是数值变化，而表示的是函数变化， 表示的是数值变化趋势，而表示的是函数变化趋势。 数值计算的基本思想就是要从离散的序列中，通过对数值计算成果的分析找出原函数的变化趋势，数值变化趋势是否反映原函数的变化趋势，取决于离散点选取的致密程度；计算数值结果的可靠性，取决于微分方程的离散、计算方法的选用和计算精度控制。可以从分布点中任意选取相邻的几个点研究离散点与原函数之间的关系，建立离散点与原函数间的变化特征。 根据泰勒公式，在点处连续可微（如图2-2），在的邻域内。 图2-2 则在邻域上可以展开为 （2-1） 其中，称为展开点，称为表征点。 当 取邻近的时，上式为 （2-2） 当与充分接近时，的幂次越高，其在（2-2）式右端求和中所占的比重越小，直至可以忽略不计，对这样的小量之和表示为或，这里。当时，用表示，当时，用表示。式（2-2）可以整理为 （2-3） 忽略项，可以看出式（2-3）的近似程度， 近似程度受一次方以上各项的影响，对这样的近似式，称（2-3）为具有的一阶精度，近似式可以表示为（对单变量函数不区分偏导数或全导数） ==或== （2-4） 为表示方便引入符号表示沿正方向的空间步长，即。当时取为，用表示时间步长。 a级数组合差分法 用表示一维空间函数，在，和，处用泰勒级数展开为点的级数式（如图2-3），以数字乘座标为下标表示导数，数值或座标个数表示导数阶数，用加数字表示位置，得 图2-3差分网格示意图 在时间坐标处，空间区间上，以为展开点，为表征点的空间差分为（）的下标表示对的导数，前的系数表示导数次数。 （2-5） 在时间坐标处，空间区间上的空间差分为以为展开点，为表征点的空间差分为 （2-6） 在时间坐标处，空间区间上的空间差分为以为展开点，为表征点的空间差分为 （2-7） 在时间坐标处，空间区间上的空间差分为以为展开点，为表征点的空间差分为 （2-8） 类似的，在空间坐标处，时间区间上的时间差分为以为展开点，为表征点的时间差分为 一阶偏导数可以由以上四式组合而成，由（2-5）得 为在时间坐标处，空间区间上具有一阶精度的一阶偏导向前差分。 由（2-6）得 为在时间坐标处，空间区间上具有一阶精度的一阶偏导向后差分。 由（2-5）减（2-6）得 为在时间坐标处，空间区间上具有二阶精度的一阶偏导中心差分。 由4×（2-5）减（2-7）得 为在时间坐标处，空间区间上具有二阶精度的一阶偏导前差。 由（2-8）减4×（2-6）得 为在时间坐标处，空间区间上具有二阶精度的一阶偏导后差。 二阶偏导数也可以由（2-5）至（2-8）四式组合而成，（2-5）加（2-6）得 为在时间坐标处，空间区间上具有二阶精度的二阶偏导中心差分。 将的一阶向前差分代入式（2-7）一阶导数中得 为在时间坐标处，空间区间上具有一阶精度的二阶偏导向前差分。 空间四点后差为 空间五点后差为 空间四点偏后差为 空间五点偏后差为 从以上的差分格式可以看出，不同的差分构造格式具有不同的精度阶数，也就是先得到构造的差分形式，而后才得到了它的精度。差分精度的阶数取决与略去式中的自变量差值的最低阶次。同价偏导选择点数越多，差分精度越高。分母中的步长方次决定偏导的阶次。余函数中的步长方次决定差分的精度。 b待定系数差分法 根据以为中心展开的泰勒级数 （2-9） 得 可以看出，要保证一阶精度需要有两个结点的物理量组合，在区间上点处对一阶导数可以表示为 （2-10） 展开等号右侧第一项代入（2-10）得 （2-11） 比较（2-11）