环空内非牛顿流体依时性Couette流动的谱法研究.pdfVIP

2010年 11月 四川师范大学学报 (自然科学版) NOV．．2010 第 33卷 第 6期 JoumalofSichuanNorma1．University(NaturalScience) Vo1．33．No．6 环空内非牛顿流体依时性 Couette流动的谱法研究 付 强 (西南民族大学理科学报编辑部，四川 成都610o41) 摘要：采用谱方法研究了环形套管内OldroydB流体依时性Couette流动的速度场分布．该问题可归结 为高阶微方程的初值边值 问题的求解 问题．采用构造的Chebysehv多项式为基的谱方法，成功地得到该问 题的解析解答．其结果与差分数值解一致． 关键词：OldroydB流体；Couette流动；谱方法 中图分类号：0373 文献标志码：A 文章编号：1001—8395(2010)06—0808—03 doi：10．3969／j．issn．1001—8395．2010．06．018 环形套管形成环空内的流体流动是化学工业 Q (Y)=2Ky一(m+1)K， 和石油工业 中的重要流动现象，环空 内的Couette Q2(Y)=8K2y 一8(，n+1) )，+ 流动就是其中一种典型流动，因此，对其研究有着 2(m +1) 一1， f1) 重要的理论意义和应用价值．现有资料 中，已有不 Q3(Y)=32K3y 一48K3(m+1)y2+ 少的文献对该问题进行研究 ¨ ．文献 [1]研究了 [24(m+1) 一6K]y一 两相流的旋转流动，文献 [2]研究了同心球问的 4(m +1)K3+3(m +1)K， Couette流动，韩式方 用符号运算与流体模型相 结合研究了上随体 Maxwell流体的旋转流动，尹析 其中，m= ，r为内圆柱半径，r为外圆柱半径，K 明。。。采用解析法研究了非牛顿流体的旋转流动． 1 本文将采用谱方法研究OldroydB流体在环空内的 一 m 一1‘ Couette流动． 2 Couette流动的谱方法解 2．1 基本方程 采用 OldroydB流体模型的本构 方程 ，设流体是不可压缩的，采用柱坐标系统(，0， =)，流场具有轴对称形式 u，=0， u日=u(r，t)， ：=0． (2) OldroydB模型本构方程为 图1环形套管管内Couette流动示意图 S+A1S=T1o(D+A2D)， (3) Fig．1 A schematicillustrationofCouetteflow inrnigsleeve 其中，S为偏应力张量，D为变形速度张量，A。为松