流体力学Petrov-Galerkin有限元法研究进展.PDFVIP

第15卷第4期 计 算 力 学 学 报 V o l. 15 N o. 4 1998年11月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS N ovem ber 1998 流体力学 有限元法研究进展 Petrov Galerk in 王建军　　陆明万　　张　雄 (清华大学工程力学系, 北京, 100084) ( 摘　要　详细评述了在流体力学分析中 方法 主要是 方法和 Petrov Galerk in SU PG ) GL S 方法 的研究进展, 并提出了今后的发展方向和研究重点。 关键词　流体分析; 法; 有限元分析 Petrov Galerk in 分类号　 241 7 O 1　引　　言 近20年来, 利用 格式直接从描述粘性不可压流体的 方程出发, 建 Galerk in N avier Stokes 立离散化的数值求解格式进行有限元分析, 已经受到人们越来越多的关注。但是, 在用试函数 与权函数相同的经典 Galerk in 法求解时, 往往会引起寄生振荡(Spu riou s O scillation s) 导致求 解失败。因此, 从80年代初期以来, 以解决数值振荡为目标, 围绕增添求解过程的人工耗散, 以 试函数与权函数采用不同形式的 原理为基础, 提出了许多以增加权函数摄动 Petrov Galerk in 项为主要策略的求解方法。 本文从说明经典 方法的数值振荡原因入手, 评述了各种 方法的 Galerk in Petrov Galerk in 基本概念、主要特点、发展历程以及它们之间的关系, 最后提出了这类方法今后应研究的重点 和发展方向。 2　数值振荡的原因 经典有限元法, 又称为 有限元法, 其试函数与权函数是相同的函数。众 Bubnov Galerk in 所周知, 用这种方法求解流体问题会产生数值振荡。 一般认为, 产生这种数值振荡有两个原因。一是对流项的存在, 该项的数值离散, 解除了相 邻节点间原有的耦合关系, 引起数值求解劣质的散射特性, 从而使求解中局部振荡不断扩散、 增大, 使数值计算失败。 产生数值振荡的另一原因, 是离散过程中速度场和压力场有限元插值函数的不当组合, 使 [ 1, 2 ] ( ) 得相应的协调条件 又称为 条件 得不到满足, 从而引起压力场的数值振 Babu ska B rezzi 荡。解决这种振荡的传统方法是速度场和压力场采用不等阶插值。一般认为, 速度场插值函数 至少要比压力场的高一阶, 才能保证不产生这种振荡。另一种方法是罚函数法[ 3, 4 ] , 这种方法利 用罚