高等计算流体力学讲义.doc

《 高等计算流体力学讲义.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、§10紧致格式一、紧致格式的定义及构造1.紧致格式的定义及待定系数的确定方法许多物理现象中包含有大范围的时间和空间尺度，湍流就是一个常见的例子，对于这些多尺度的物理现象的数值模拟要求适当地体现出所有相关的尺度，特别是小尺度的演化过程。这些要求促进了谱方法的研究，但是谱方法只能局限于简单计算区域和简单边界。为了克服谱方法的局限性，人们试图采用有限差分方法和谱元法，并在这方面取得了一定的进展。对于显式的高精度、高分辨率的有限差分法，必然需要较多的计算模板点，这样将带来边界点处理的困难。以均匀网格上的有限差

2、分问题为例，在求某一点上物理量的导数（为简单起见，仅考虑一阶导数的情况）时，一般是利用该点附近的有限个网格点上的函数值构造导数的差分近似，即，(1)其中代表对处导数的差分近似。一般来说，差分算子L是线性的。此时，如果，(2)则。(3)也就是说，要构造阶精度的差分近似，至少需要个模板点。目前，湍流的直接数值模拟、大涡模拟、计算气动声学领域对数值计算的精度提出了很高的要求。如果用常规的方法构造高阶格式，需要的模板点数也相应增加，这给边界附近数值方法的实施和边界处理带来了困难。紧致格式是一种可以用较少模板点

3、逼近导数到高阶精度的差分格式，它除了具有边界处理相对简单的优点外，更重要的是具有类似谱方法的高分辨率；这一点对于计算具有复杂多尺度性质的流场（如湍流和声场）具有重要的意义。所谓紧致格式，是一种隐式的数值微分计算格式。即(4)一般情况下，Q也是线性算子。于是，i点上物理量的导数值不再能单独求出，而需要同网格线上的其它点（在附加的边界条件下）一起求解。下面讨论一般的线性紧致格式的构造过程。线性紧致格式可以写为：。(5)如何确定式的系数，使得逼近到指定的阶数呢？我们把左右两侧在处作Taylor展开，有(6)

4、对式求导，有(7)(8)……把式乘以与式合并，有(9)把式乘以与式合并，有(10)依次类推，可以消去式中的高阶导数项。由式，要使逼近到一阶精度，需要满足：要使要使逼近到二阶精度，则应有要使要使逼近到三阶精度，则应有可见，如果要构造阶精度的格式，只需式两端前个对应项的系数相等即可。在式中有k+l+m+n+2个待定参数，按照上面的分析，该紧致格式最高可达k+l+m+n+1阶精度。当要求的精度低于k+l+m+n+1时，中可以存在无穷多种满足精度要求的系数。特别的，当要求的精度为k+l+m+n阶时，可以有一个

5、系数可以任意给定；当要求的精度为k+l+m+n-1阶时，可以有两个系数可以任意给定。其余可依此类推。在确定系数后，实际求解的过程中，式的右端是已知的。左端的未知量构成了一个有k+l+1个对角元素的方程组，需要结合边界条件联立求解。一般三对角线方程可以用高效的追赶法求解，所以常取k=1,l=1。2．中心型紧致格式如果在式中且，则称相应的紧致格式为中心型紧致格式。例如左侧三点、右侧5点的中心型紧致格式为(11)式中为网格间距，方程组的系数矩阵为三对角阵。将上式中的各项进行Taylor展开，并令各阶项（）的

6、系数相等，得到满足第一式时有二阶精度，满足前两式时有四阶精度，三式都满足时有六阶精度。四阶精度的格式满足，其中可以任取。其截断误差为如果，则格式退化为四阶中心差分格式，即截断误差为如果，那么第三式满足，格式具有六阶精度，此时截断误差为。从上面的例子可以看出，在相同模板点下，常规的差分方法精度最高为四阶，而紧致格式的精度可以达到六阶。2．迎风型紧致格式迎风型紧致格式是一种偏心（左右不对称）紧致格式。考虑方程，时空间导数用下式计算(12)而时计算方法为：。(13)这种类型的格式即为迎风格式。以的情况为例分

7、析其精度，利用Taylor展开并令各阶项（）的系数相等，得到满足前两式具有两阶精度，满足所有三式时具有三阶精度。由前两式得到截断误差为如果，则格式退化为二阶迎风格式，即截断误差为如果，那么第三式满足，格式具有三阶精度，此时截断误差为时，容易知道的关系与前述情况相同。二、紧致格式误差的Fourier分析1．分辨率分析设是周期为L的周期性函数，把[0，L]上划分为N等分，网格间距h=L/N，对作Fourier级数展开：(14)定义无量纲化波数与坐标分别为和，则解可以表示为：(15)对求一次导数后，得到：(

8、16)假定用紧致差分得到的导数可以写为：(17)于是可以用来衡量差分后的导数值逼近真实导数值的程度。如果规定某个值，则总可以找到一个区间，使得位于这个区间内的所有波数，逼近真实导数的误差都小于。所以代表了格式在条件下的分辨能力的有效范围，越大，格式的分辨率越高。一般称为分辨效率。下面推导的公式。为了简便起见，以三点中心型紧致格式为例，即：(18)将方程当中各项作Fourier展开，只考虑一个频率成分，可得：其中：。把上面的关系代入式，化简得到：对于我们讨