求解对流扩散方程的紧致修正方法.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093466 求解对流扩散方程的紧致修正方法 张昆，杨茉，康张扬 （上海理工大学动力工程学院，上海，200093） (Tel: 021 E-mail:zhangkun@126.com) 摘要：提出了求解对流扩散方程的紧致修正方法，该方法是在低阶离散格式的源项中，引入紧致修正 项，从而构造高阶紧致修正格式进行求解。发展了求解一阶和二阶导数的基于非等距网格的紧致格式， 先求得方程的低精度数值解，代入紧致格式后得出紧致修正项。采用紧致修正方法对典型的对流扩散 方程进行计算，结果表明，紧致修正方法虽然与二阶经典差分方法建立在相同的结点数上，但其具有 四阶精度，所以计算误差较小。在相同的结点数下，非等距网格可以得出比等距网格精度更高的解， 紧致修正方法比经典差分方法对网格的依赖性更强。 关键词：紧致格式；紧致修正方法；对流扩散方程； 非等距网格。 0 前言 流体流动现象可采用对流扩散方程来描述，求解该方程的经典差分格式由于精度较 低，不能满足实际计算和模拟的需求，而高精度的谱方法对边界条件的限制又十分苛刻。 所以紧致差分格式的研究受到大量研究者的关注。求解结点导数值的紧致差分格式是由 Lele 于 1991 年提出[1]。该格式是一种高精度，基架点少的隐式差分格式。采用紧致格 式求解结点的导数值时需要求解关于结点导数和结点值所构成的代数方程组，所以紧致 格式不能像经典格式那样直接应用于对流扩散方程中导数项的离散。目前求解对流扩散 方程的紧致（紧凑）方法多是通过摄动原理建立，即对现有的低阶格式进行截断误差分 析，然后修正原低阶格式的系数，消除截断误差中的低阶项，从而提高截断误差阶而得 到的。陈国谦、Zhang 分别建立了二维和三维对流扩散方程的均分网格上的四阶紧致格 式[2-3]，非均分网格上的高阶紧致差分格式由田振夫建立[4]。这些应用摄动原理所构造 的紧致格式虽然简单明了，但因计算复杂，实际应用此思想来提高截断误差阶并不容易 [5]。同时，以上二维对流扩散方程的四阶紧致格式是建立在9 个结点上，而三维是建立 在 27 个结点上，所以摄动原理所得出的紧致格式与经典格式具有相同阶数的截断误差 时，实质上是建立在同样数目的结点上。 为使紧致格式能够灵活的应用于对流扩散方程的数值求解，又能建立在较少的网格 结点上，提出了一种新的紧致格式的修正方法，简称紧致修正方法。该方法对于那些原 来采用低阶格式的离散方程，在源项中加入紧致修正项即可纳入高阶紧致格式，求解过 程中，需要先计算低精度的数值解，利用已得到的数值解，通过已建立的有关导数的紧 基金项目：国家自然科学基金资助项目（No ），上海市教委发展基金项目资（NO.J50501 ）， 上海理工大学创新基金 致格式求解导数值，然后将其代入修正项后，再计算高精度的数值解。紧致修正方法与 延时修正方法[6,7]中高阶离散格式的引入思想相同，所以可增加迭代求解代数方程组的 稳定性，但与延时修正方法的求解过程不同，这样做既可避免延时修正方法的反复迭代 过程，又可以保证数值解的高精度。本文所建立的求解二维对流扩散方程的紧致格式是 建立在5 个结点上，三维建立在7 个结点上，所以该格式的基架点数明显少于以往的紧 致格式。有利于增加迭代求解代数方程组的收敛性。 1 数值方法 1.1 基于非等距网格的紧致格式 以往的紧致格式多数建立在等距网格上，应用于非等距网格时需经过坐标变换。以 下建立的紧致格式可以直接应用于非等距网格的数值求解中。对于函数φ ，建立如图一 所示的计算网格{x , i 1,2,..., N −1, N } ，网格间距h x −x 。 i i i i−1 图一：一维网格系统示意图 在i 点处，具