旋转库特流中的有限元法.docVIP

2002 年 11 月JOURNAL OF SICHUAN UNIVERSITY ( ENGINEE 2002 年 11 月 JOURNAL OF SICHUAN UNIVERSITY ( ENGINEERING SCIENCE EDITION) Nov. 2002 文章编号 :100923087 (2002) 0620071204 旋转库特流中的有限元法 柳1 ,陈文梅2 ,褚良银2 ,汤潍蔚3 杨 (1 . 四川大学 生物医学工程系 ,四川 成都 610065 ;2 . 四川大学 化工学院 ,四川 成都 610065 ;3 . 四川省 纺织工业研究所 ,四川 成都 610072) 摘 要 : 讨论了旋转管式膜器环隙间定常流动的有限元解法 。首先给出了定常流动的解析解 ;其次 ,根据流体力学 方程建立了元素方程组 ,并将元素方程组汇总成总体方程组 ,归纳了旋转库特流定常流动的有限元算法 。最后 ,通 过给定的参数计算了旋转库特流切向速度分布的数值解 ,并绘制了相应的切向速度分布曲线 。 关键词 :旋转库特流 ;有限元法 ;元素 中图分类号 : TQ051. 84 文献标识码 :A The Finite Element Method of Rotary Couette2Poiseuille Flo w YAN G L iu1 , CHEN Wen2mei2 , CHU L iang2yin2 , TAN G Wei2wei3 (1 . Dept . of Biomedical Eng. ,Sichuan Univ. ,Chengdu 610065 ,China ; 2 . School of Chem. Eng. , Sichuan Univ. , Chengdu 610065 , China ; 3 . Sichuan Textile Industry Inst . , Chengdu 610072 ,China) Abstract : The finite element method of rotary Couette2Poiseuille flow in annular gap of rotary tubular membrane separator was discussed in this paper . At first , the analytical solution of the steady flow was obtained. Then the whole equations were obtained by constructing the element equations of the fluid equation , and the algorithm was given out at the same time . Finally , the numerical solutions of the azimuthal velocity were calculated and the corresponding distribution curves were draw out by assuming some parameters of the flow. Key words :rotary Couette2Poiseuille flow ; finite element method ; element 有限元法起始于 20 世纪 50 年代 ,70 年代开始 在流体力学中应用 ,至今已较为完善 。国外已用有 限元法来解亚音速 、超音速和跨音速的欧拉方程 、全 速度位方程和 N - S 方程 ,并广泛应用到航海 、航空 和航天领域 。有限元法的特点是 : 网格的划分可以 是任意的 ,且网格的划分可以比较大 ,它弥补了有限 元差分法的不足之处 。本文讨论了旋转管式膜器环 隙间定常流动的有限元解法 ,并编制了算法程序 ,得 出旋转库特流切向速度分布的数值解 。因为旋转管 式膜分离器内的流场可以近似表为两同轴圆筒环隙 间的流场 ,所以该算法可为进一步研究旋转管式膜 器的流场作准备 。 旋转库特流的解析解 先给出同轴旋转圆筒环隙间的流体在定常流动 下的速度表达式 ,简述如下 :两半径分别为 R1 和 R2 ( R2 R1) 的无限长旋转同轴圆筒环隙间充满粘性 1 不可压缩流体 , 内筒和外筒的转速分别为 ω1 和 ω2 。 采用柱坐标系 ( r ,θ, z ) , 设 Vθ 表示流场中任一点处 的切向速度 , 则定常流动时流体运动的 N - 表示如下 : S 方程 V 2 d p =ρ θ 收稿日期 :2001212210 基金项目 :