控制体有限元方法中的一种高阶插值格式.pdfVIP

航空动力学报990223 航空动力学报 JOURNAL OF AEROSPACE POWER 1999年　第14卷　第2期　Vol.14　No.2　1999 控制体有限元方法中的一种高阶插值格式 ** 王旭 　谷传纲 　　【摘要】　针对控制体有限元方法对流动方向不平行于坐标轴方向的回流问题数值模拟 采用各种迎风格式都容易产生假扩散的问题，本文给出了一种由插值单元及其迎风相邻单元 所确定的高阶插值格式。通过对Re=1000时三种不同形状的二维四边形空穴流模拟，表明该格 式与局部斜迎风格式的适当组合可有效地减小控制体有限元方法对回流问题数值模拟的假扩 散。 　　主题词：　有限元　插值法　流动 　　分类号：　O357 　　控制体有限元方法(CVFEM)是用于流动与传热问题数值模拟的有效方法，便于引入迎风 ［1］ 格式及采用速—压同位网格的原始变量法，适用于各种复杂几何域 。但当流动方向与坐 标轴方向不一致且网格Reynolds数较大时，各种迎风格式仍会产生严重的假扩散。本文给出 控制体有限元方法在二维三角形单元网格域上的一种高阶插值格式，并通过对Re=1000时，三 种几何形状上二维空穴流的模拟验证该方法对回流问题的适用性。 1　控制方程及其控制体有限元方法 　　设稳态的粘性不可压流体流动的控制方程为：动量方程 ． （ρv u-μ u）＝S - p i i i ／ x i (1) 连续性方程 ． （ρv）＝0 (2) 其中：ρ为流体密度，μ为粘性或扩散系数，S 为源项。 i file:///E|/qk/hkdlxb/hkdl99/hkdl9902/990223.htm（第 1／6 页）2010-3-22 22:44:44 航空动力学报990223 图1　插值单元上的子控制面及其相邻单元 　　连接三角形单元网格域上每一单元的形心点与边中点，对每一结点P构成一个控制体V 。在V 上对 P P 方程(1)和(2)积分得： (3) (4) 其中δS 为V 界面，n为δS上微元ds的单位外法向量。 P P P 　　对(3)式中的扩散项变量采用多项式插值，对流项变量采用局部斜迎风插值及高阶插值。 　　局部斜迎风格式：如图1所示，定义质量流率为： (5) 其中：δSr 表示子控制面 的面积(长度)，n为按图1规定的单位法向量。 r