求解对流扩散方程用的几种欧拉_拉格朗日模型之比较.docxVIP

求解对流扩散方程用的几种欧拉一 拉格 朗 日 模 型 之 比 较 〔葡萄牙 〕 A M 巴普蒂斯泰 〔美国〕 E E 阿达姆斯等 摘 要 本文就 对流扩 散方程的 各种欧拉一拉格朗日解法所使用的插值格式的精确度和稳定性从理论分析和数值试 验上 进行了比较 己 l  言 欧拉一拉 格朗 日方法 (简作 E L 人D 为传输方程 (对流扩散方程 ) 提供了经济 而准确的解法 〔 〕 其 中有些 方法 (E L M / C ) 直 接把浓度 处理成主要因变量〔 ” ) 而 另 外一些则追综质点进而和浓度相联系在各种E LM/C方法中多数把方程分裂成对流方程和 扩散方 程 并采用 各自合 适的技术 求解 但事实上整个传输方程是可以 同时求解的〔3 〕 当按 照分裂法处理时大多数情况是将传输方程从时间域进行离散其实也可以直接将微分方程分裂( 〕 典型的 E L M / C 方法用向后特征方法解对流方程 (或在可能的情况下 解 整 个输运方程) 其 中包括一个追综步和一个求特征线根部浓度的插值步 已经提出了许多种插值格式〔78)及木文各具有不同的稳定比阻尼和耗散特性最终所有的ELM方法都与某种给定的离散单元技术联系起来这些技术或者贯穿于求扩散方程的数值解过程 或者 仅仅通过对流场的 定义 大多属于 有限差分法 〔 〕或有限元法〔” ” ) 由于 种离散技术均有各自的优点就希望对ELM有一个系统的比较本文即致力于这样的比较工作侧重于E LM/C方法解纯对流方程时可以使用的关于浓度的各种插值格式 二 问题的定义和对插值格式简要评述 考察下面形式 的对流方程 黔〔二器+ 歹 卫口C 任- 、= O (i = 1 2 3 ) (1) 口X 诬 J 上式意味着 C 梦 = C 梦 P 是由下式 定义的特征线的根部: d X (2 ) 以 t = 一 谊 其起点为位于时刻—的节点 夕 (见图 1 ) 假定时 刻 。 的所有各节 点的 浓 度 都 已知 则 梦可以通过 点附近适当的节点 (独立于位置 声 的插 值求的 C P ) ( a ) 一维 ( b ) 二维 图1 衬流 方程解格式杀 意图 最简单的插值格式是 线性格式 采用 (在一维情况下) 两个最靠近 P 的节点 然而 这 种格式会引起过大的数值阻尼就象在常规逆风格式里出现的一样因此早期的ELM‘)采用了P 点周围的三节点 (一维问题) 二次拉格朗日多项式 相同的格式 即三节点拉格朗日插值 (简作 3 P一L l ) 曾被本魁等〔 ) 和 巴 普蒂斯泰等〔 〕用于 有限元 (以 及文献〔1 。〕提出的不受柯朗稳定性约束限制的格式) 3 P L l 格式比传统的欧拉方法精确 但仍然会 出 现某些数值阻尼和分散增高拉格朗日插值多项式的阶次似乎可能改善精度可是应用四节点格式 (4 P L D 的 结果导致 显著的数值分散和不稳定状 态 关 于这一点 下文还要论述 霍雷和普列士曼〔〕使用了另一种基于埃尔米特多项式的精确四阶格式使用两点埃尔米特插值的 ZPH l包含最少的插值点数 (一维 问题是 2 点) 但要所需变量除浓度本身外还需要浓度由于对流和扩散而引起的空间导数在一维问题里做到这点并不难但对二维 以上的间题则较为复杂 霍 雷和小松〔 〕及小松等〔 〕试图在保持 ZP H l 精度的俞提下不把浓度导 数作为 独立变 他们 提出了八节点和六节点的准埃尔米特插值格式 (分别以S P Q H 6 P Q H 表示 ) 在这些格式中通过由在 K K 一/1 2 和 K 一 1 三点的浓度及其导 数定义的 (一维 问 题 见 l a ) 两个三次多项式的平均得到 C 梦 = C 彗 至于导 数值则通过顺序 四个节点上的 五个 (对SP Q H ) 或三个 (对 6 P Q H ) 三次拉格朗日多项式所得值并加权平均定 出 最佳” 权重由数值试验和使数学公式 中由于泰勒级数展开引起的人工扩散项变得最 小而得到的 按上述简单方法用三阶或二阶拉格朗日多项式估算K和K一1点处的浓度导数并直 接代入 Z P H l 算式 中去后 我们增补了一些另外的格式 4 P H l S P H l 和 7 P H l 由于这 些格式尚未发表我们在表1中说明了其形成原理 面我们把重点放在一维问题上因为对正交网格来说扩展到二维或三维问题从概念上说是容易的只要求把一维插值公式适当地相乘〔58〕不幸的是对非