高Re 数下NS 方程有限元数值解法研究 - 河海期刊.PDFVIP

第 卷第 期 河 海 大 学 学 报 !% ’ 8G1 +!% ,G + ’ 年 月 !$ # =?@,)A B CC)D ?,D8E@.D7F H0I !$ 高 数下 方程有限元数值解法研究 ! !# $ ! $ ! 郭新贵 ，倪福生 ，洪锡军 ，胡沛成 （ 上海交通大学塑性成形工程系，上海 ； 河海大学机电工程学院，江苏 常州 ） $+ !’ ! + !$’!! 摘要：采用原参数形式 方程直接解法，研究高 数下 方程数值解法 波阵技术可节省内 ,-. ! ,-. + 存，提高计算效率；混合插值函数可避免求解过程中的压力振荡；时间推进解法和简化迎风有限元 技术可防止解的非物理振荡；低雷诺数的收敛解作为高雷诺数的初场，定常解作为非定常解初场， 可提高解稳定性，加快收敛速度+ 最后以驱动方腔为计算实例，对解法进行了验证，得到的结果与经 典结果吻合较好，为高 数下 方程的数值解提供了参考 ! ,-. + 关键词： 方程；高 数； 有限元；时间推进法 ,-. ! /0123456 中图分类号： 文献标识码： 文章编号： （ ） 78$9’ + $ ) $ $%: !$ ’ ;# 9 ! ! ! 方程的数值求解一直是流体力学的首要难题 特别是在高 数条件下，方程中粘性项作用减弱，椭 ,-. + ! 圆方程退化为双曲方程，难以得到稳定解 本文采用原参数形式的 方程的直接解法，对高 数下 方 + ,-. ! ,-. 程数值解法进行了研究 在直接解法中，散度条件作为压力方程给出，动量方程与压力方程同步求解，避免了 + ［］ 动量方程与压力方程迭代求解，从而使计算程式得到简化 $ + 最后采用/0123456 有限元迭代解法和时间推进 解法，对高 数下驱动方腔内剪切流动进行数值模拟，得到了较好的结果，为高 数下 方程的数值解