两个非同心旋转圆柱间粘性流动的基本法流及其数值解.pdfVIP

第26卷第4期 工 程 数 学 学 报 v01．26No．4 2009年08B OF MATHEMATICSAug．2009 CHINESEJOURNALENGINEERING 文章编号：1005-3085(2009)o垂0637．11 两个非同心旋转圆柱间粘性流动的基本流及其数值解术 张引娣1，一，李开泰2 (1一长安大学理学院，西安710064；2-西安交通大学理学院，西安710049) 具，得到两个非同心旋转圆柱之间粘性流体流动的基本流所满足的方程。这一基本流由两部分组 成，一部分是以角速度u旋转的同心圆柱所产生的基本流，而另一部分是带有高阶小参数的扰动 流动。把扰动项用三角函数展开得到一满足非齐次边界条件的常微分方程组。在此基础上，我们 得到了基本流的数值求解方法，并给出数值算例。 关键词：修正双极坐标系；张量分析；渐近分析方法；Navier-Stokes方程 中图分类号：0242．21 文献标识码：A 分类号：AMS(2000)76D05；65N99 1 引言 轴承润滑理论中，利用两个非同心旋转圆柱之间流体流动产生的升力来达到工程设计的目 的。非同心旋转圆柱之间流体的流动，当内外圆柱半径之差与圆柱半径相比很小时，是一个 等[9-13】。在轴承润滑理论中，李开泰应用变分不等式等原理，建立了Reynolds方程【8J。 理论，在修正双极坐标系下，把偏心率作为一个小参数，运用张量形式下的摄动方法，可推导 出两个非同心旋转圆柱之间粘性流动的基本流所满足的方程，它是两个同心旋转圆柱之间粘性 程组，对得到的常微分方程组我们必须先进行边界条件的处理，最后化为齐次边界条件。在此 基础上，我们将对基本流所满足的方程数值求解，并且将给出数值计算例子。 2 两个非同心旋转圆柱之间的基本流 两个非同心旋转圆柱间流体的流动，当内外圆柱半径之差与圆柱半径之比很小时，是一种 薄区域上的流动。这个区域在半径方向上的尺寸与周围方向长度之比很小。在这种情形下，人 们自然提出一个问题，非同心圆柱之间流动是否可以视为两个同心旋转圆柱之间流动的一个扰 动，为讨论这个问题，采用一种“修正双极坐标系”并利用张量分析，可以得到一个正面的回 答。并且可以有效地得到基本流以及二次流。 收稿日期：2008-05-09．作者简介：张引娣(1962年5月生)，女，副教授．研究方向；偏微分方程数值方法． ’基金项目：国家自然科学基金10771167)． 万方数据 垡塑————————————————j巳—』堕堂堂塑 苎!!堂 。…。。”…… 、。设竺圆柱的半径分别为n，6，且设6一n《口，偏心距记为e，间隙比率记为6，偏心率记 为E，那么 ’ 6=1b-一a，e：E6，o≤E1．o 一 ～1，(1) ；，修正翠磬卷尊。曼外嘤柱为固定，内圆柱旋转角速度为u，外圆柱的旋转轴为圆柱坐标 。一” 系(r，0，z)的z轴。那么双极坐标系(p，lp，名)，它与(r，0，名)有如下关系 跏，加掣，专钏坩，e：砂． (2) 这是垡箩煮实擎，0=o，妒=o是同一条射线。当E_o时，∞，妒)与∽曰)重合。 ～……… ，。，翟：圾坐标系(p，垆，z)下，两个非同心圆变成同心圆，其内面硅酌半荏焉1，外圆柱的半径 为厣，即