对论文“任意精度的三点紧致格式及其在CFD中的应用”的评论.pdfVIP

对论文“任意精度的三点紧致显格式 1 及其在 CFD 中的应用”的评论 张红娜 1* 王艺 2 宇波2 魏进家 3 李凤臣 1 1. 哈尔滨工业大学能源科学与工程学院，黑龙江 哈尔滨 150001 2. 中国石油大学石油天然气工程学院，城市油气输配技术北京市重点实验室，北京 102249 3. 西安交通大学能源与动力工程学院，陕西 西安 710049 *通信作者：E-mail:lala12587@ 摘 要：本文证明了林建国 2007 年发表在应用数学和力学上“任意精度的三点紧致显格式及 其在 CFD 中的应用”[1] 中提出的紧致显格式与传统的差分格式实质相同，是传统差分格式的 另一表达形式，并不具有紧致格式的优点。尽管如此，但这种表达形式更紧凑，推导获得高 精度的差分表达式相对于传统的泰勒展开求待定系数的方法也更加简单。 关键词：显式紧致 传统差分 中心格式 迎风格式 中图分类号：TK01 1．引言 隐式紧致格式在湍流的数值计算中得到了广泛的应用，为了避免隐式紧致格式联立代数 方程组求解，林建国等人[1]利用逐阶迭代的思想，提出了一种关于三点的任意精度“显式紧 致”格式，以提高计算效率。从形式上看，这种显式格式的表达简单，易于编程。然而在展 开表达式中的差商算子后，发现这种格式与传统的差分格式一致，只是表达方式不同，并且 在顶盖驱动流算例中，分别采用这种格式与同精度的传统差分格式编程计算，得到的结果完 全一致。因此得出结论，这种“显式紧致”格式其实质是传统的差分高精度格式的另一种表达 形式。尽管与传统差分格式实质相同，但这种格式更容易得到高精度的表达式，利用传统方 法得到高精度表达形式的推导过程非常繁琐，而采用这种方法，提高精度时只需增添表达式 中的某一项，易于扩展。此外，采用高精度差分格式时，这种“显式紧致”形式更适合于编程。 本文中，首先给出文献[1]中得到的“显式紧致”格式，接下来通过理论推导证明这种格 式与传统的差分格式实质相同，并且在顶盖驱动流算例中，通过比较说明采用这种格式与传 统差分格式得到的计算结果一致。 2 ．均分网格下空间导数三点“显式紧致”格式 下面给出文献[1]中得到的二阶导数和一阶导数的表达式。首先定义均分网格上算子和差 商： k −1 k −1 f −f k δ f −δ f i+1 i−1 x i+1 x i−1 δf = ， δ f = ， x i x i 2 2 k k −1 k −1 + + + + δ f =f −f ， δ f =δ f −δ f x i i+1 i x i x i+1 x i k 1 k 1 k −1 k −1 k −1 δ f =f +f −2f ，δ f =δ f =δ f +δ f −2δ f xx i i +1 i −1 i xx i xx i xx i +1 xx i −1 xx i 其中，f i 为i 点的函数值，k ≥2 ，Δx 为空间网格步长 。两竖线“| |”表示算子，上 标表示用该算子进行运算的次数。利用逐阶迭代以及忽略高阶小量的方法，最终可以得到二 1本文受到教育部科学技术研究重点