流体力学有限元中的边界条件处理解析.doc

 流体力学有限元分析中的边界条件处理 梁启国 尹敏镐 赵永凯 高殿荣 燕山大学 大庆石化总厂 摘 要 阐述了流体力学有限元分析中应用流函数 - 涡量法时典型边界条件的处理方 法 ,并给出计算实例 . 关键词 有限元 涡 - 流函数法 边界条件 分类号 O357 . 1 引言 求解不可压缩粘性流体二维流动问题的数值方法有速度 - 压力法 、涡 - 流函数法和 流函数法 ,其中涡 - 流函数法应用较广1 ～4 . 在应用涡 - 流函数法进行有限元分析时 ,数值边界条件的处理不但影响解的精度而 且影响解的稳定性2 . 本文结合工程实际给出了在有限元分析中确定几种典型边界条件 的方法. 0 边界条件的处理 如图 1 所示的沿背部台阶的不可压缩粘性流 体二维流动. C 点为角点. 在涡 - 流函数方程的求解过程中 , 不仅需要 知道边界上的流函数和涡量值 ,有时还要对边界 上的涡量值不断地进行修正. 下面分别讨论各种 边界上的流函数值和涡量值. 设Δ x 和Δy 分别为 x 方向和 y 方向的有限元划分网格间距 , i 和 j 分别 为 x 方向和 y 方向的步长指标. 1) 壁面边界处理 若壁面是不可渗透的 ,则沿壁面的流函数 Ψ 为常数. 一阶精度的壁涡公式为 : 1 图 1 沿背部台阶的边界条件 2 (Ψ i , j - Ψ i , j +1 ) Ωi , j c c + O (Δy) ( 1) = (Δy ) 2 c 二阶精度的壁涡公式为 : 3 (Ψ i , j - Ψ i , j +1 ) Ωi , j + 1 + O ( (Δy ) 2 ) Ωi , j c c c ( 2) = (Δy ) 2 2 c 图 2 用内点上的 Ψ i , j +1 , Ωi , j +1 表示边界上的涡量值 . 二阶壁涡公 c c 式在大网格雷诺数和变网格情况下往往引起数值不稳定和计算不收敛 ; 一阶壁涡公式虽 然精度低 ,但计算稳定 、易于收敛. 另外一种二阶壁涡公式5 为 : 8Ψ i , j +1 - 7Ψ i , j - Ψ i , j +2 + O ( (Δy ) 2 ) Ωi , j c c c ( 3) = 2 (Δy ) 2 c 用时间迭代法确定壁涡的方法6 ,7 为 : 根据不可压缩粘性流体二维流动的流函数 - 涡量公式8 ～10 92Ψ 92Ψ = - Ω + 9x 2 9y2 92Ω 92Ω 9Ω 9Ω 9Ω 1 + + u + v = 9x 2 9y 2 9t 9x 9y Re 对上述方程进行有限元离散 ,得总体矩阵方程 KΨ = MΩ - q MΩQ+ ( K + A )Ω = s ( 4) ( 5) 式中 ,ΩQ= 9Ω/ 9t , K , M , A , q 和 s 为系数矩阵. 在求解式 (4) 和 (5) 时 ,为确定壁面涡量 , 首先由式 (4) 和 t = t 0 时刻的涡量初值来计算 t = t 0 +Δt 时刻的流函数 ,有 Ψ t 0 +Δt = K - 1 ( MΩt 0 - q) ( 6) 然后再用式 (4) 和新的流函数值 Ψ t 0 +Δt 求得壁涡 Ωt 0 +Δt = M - 1 ( KΨ t 0 +Δt + q) ( 7) 最后 ,从式 (5) 出发 ,并利用新的壁涡值可求 t 0 +Δt 时刻其内部各点的涡量值 . 以后每增 加一个时间步长都按上述方法确定壁涡 ,所以也称壁面涡的时间迭代法. 边界 B 6 的处理和 B 1 相同. 2) 对称边界处理 由于流动的对称性 ,图 1 中的边界 B 5 为对称边界 ;而对称线必为流线 ,因此在对称边 界上的流函数 ΨB 为常数. 5 对称边界是滑移边界 ,显然 9 9u v vB = 0 = 0 = 0 9x 9y 5 B B 5 5 由此可求对称边界上的涡量 9 9u v ΩB = = 0 - 9x 9y 5 B 5 3) 上边界处理 上边界 B 3 一般为以下两种情形 :其一