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

、 、 维普资讯 q，}2、 一 鸦／ 赵永凯 高殿荣 梁启国 尹敏镐 0；7’| 燕山大学 大庆石化 ．巷厂 摘 要 阐述了流体力学有限元分析中应用流函数 一涡量法时典型边界条件的处理方 法．并给出计算实例 关键词 涡 一流函数j去 边界条件 磁 分类号 03571 O 引言 求解不可压缩粘性流体二维流动问题的数值方法有速度 一压力法、涡 一流函数法 和流函数法．其中涡 一流函数法应用较广LI～ ． 在应用涡 一流函数法进行有限元分析时，数值边界条件的处理不但影响解的精度而 且影响解的稳定性 2【】．本文结合工程实际给出了在有限元分析中确定几种典型边界条件 的方法 ． ’ 1 边界条件的处理 如图l所示的沿背部台阶的不可压缩粘性 流体二维流动 ．(1点为角点 ． 在涡 一流函数方程的求解过程中．不仅需 口 要知道边界上的流函数和涡量值．有时还要对 边界上的涡量值不断地进行修正 ．下面分别讨 瀛 论各种边界上的流函数值和涡量值 ．设 Ax和 … 、 A，分别为x方向和 方向的有限元划分网格间 ～ 距．q田，分别为 方向和1’方向的步长指标 ． 1)壁面边界处理 若壁面是不可渗透的．则沿壁面的流函数I}- 图 L 沿背鄙台阶的边界条件 为常数 ．一阶精度的壁涡公式为： +D( ) (1) ( 二阶精度的壁涡公式为： = +D( )2) 2【) 图2 用内点上的 +，nf-，+l表示边界上的涡量值 一二阶壁涡公式在大网格雷诺数和 变网格情况下往往引起数值不稳定和计算不收敛；一阶壁涡公式虽然精度低．但计算稳 1996年9月 1日收到 t本项 目得到浙江大学流体传动及控制国家重点实验室开放基金的资助 ／ 维普资讯 东北重型机械学院学报 997妊 定．易于收敛 ． … 另外一种二阶壁涡公式p为： ： Qf．=—— ——— — 一 +十D((【△r)J L3j【J) “ 2(Ay)‘ 用时间迭代法确定壁涡的方法L。’71为： 根据不可压缩粘性流体二维流动的流函数 一涡量公式 【～101 + ： 一Q