不可压缩流动问题的无网格局部PetrovGalerkin方法的数值模拟.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：093492 不可压缩流动问题的无网格局部 Petrov-Galerkin 方法的数值模拟 1,2 1 1 1 吴学红 戴艳俊 李增耀 陶文铨 1 西安交通大学能源与动力工程学院多相流国家重点实验室 2 郑州轻工业学院机电工程学院 Tel：029 E-mail： wqtao@mail.xjtu.edu.cn 摘要:本文将MLPG 方法应用于计算二维不可压缩流动问题，提出用SUPG 方法克服非线性对流项所产生 的不稳定性，并用混合公式去克服或减轻产生棋盘型的压力场。通过对顶盖驱动流、突扩后台阶流动 两例子的计算表明：本文提出的对流项处理方法有较强的健壮性，即使在较高的Re 情况下也能给出稳 定的收敛解；并且与文献值相比，本文构建的方法都能给出较好的计算结果。 关键词:无网格方法，迎风格式，不可压缩流, 数值模拟 1 引言 无网格方法是从计算固体力学发展起来的一类新的数值方法，最近十几年一些研 究者尝试把无网格方法引入到计算流体力学和传热学领域中。目前无网格方法在计算流 体力学和传热学中的应用还不成熟与完善，还无法与有限容积法或有限元方法相比。 从有限容积法可知[1,2] ，对于不可压缩流动问题，由于非线性对流项和压力梯度项 的存在，易产生振荡速度场和压力场，因此必须研究对流项与压力项的专门处理方法来 克服这两个困难。对于消除非线性对流项引起计算结果的振荡问题，以往有限元方法与 有限容积法基本上是采用迎风思想的格式来加以克服或减轻，如：有限元方法中的 SUPG[3]方法；有限容积法中的一阶迎风格式、二阶迎风、SGSD （Stability guaranteed [4] second-order difference ） 、QUICK格式等。为了克服压力项离散造成棋盘式压力分布 的困难，有限容积方法发展了一系列的专门算法，如SIMPLE系列[2]算法、CLEAR算法 [5,6] [7,8] 、IDEAL算法 等。有限元方法则用所谓的混合公式法去解决压力项所产生的振荡。 混合公式被应用去解决由1971年Babuška[9]和1974年Brezzi[10]建立起来的压力稳定性条 件，后来称这种压力稳定性条件为B-B条件。B-B条件是使解存在唯一性的保证。 由于无网格方法主要是从有限元方法发展起来，有很多相似的地方，一些有限元 方法中的计算方法，可以很方便地应用于无网格计算方法之中，所以本文在构建无网格 方法的离散计算式时用SUPG方法去满足对流的稳定性，同时应用混合公式去满足B-B 条件。在本文应用该种无网格方法计算两个有基准解的算例，去验证该种方法的稳定性 与计算精度。 下面主要对目前文献中的常用于无网格方法计算流体问题中的迎风格式进行阐 述；对于压力项处理的混合公式方法可见文献[11] 。 2 对流项的处理方法 目前应用于计算流体动力学中的无网格方法主要有配点型无网格方法（如FPM 、 SPH 方法）、RKPM 、MLPG 方法、EFG 方法等。下面介绍一下目前在无网格方法中常 见的各类迎风格式。 2.1 有限增量微积分方法 有限增量微积分方法由Oñate 等[12]于1998 年提出。它的基本思想是：如对流扩散 方程，将各项用Taylor 级数展开，将瞬态项和源项展开至一阶，对流项和扩散项展到二 资助项目：国家自然科学基金重点项目（No.）的资助。 阶（理论上FIC 方法能够建立更高的阶数），然后代入对流扩散方程。Oñate 及其合作者 [13,14]把FIC 方法应用到对流扩散方程中，并用最小二乘配点方法对方程进行离散，发展 了有限点法（FPM ）。从FIC 的特点可以看出，这种稳定性的方法相似于有限差分方法， 它不满足一致性条件，特征尺寸也不容易选取，而