高雷诺数时旋转圆柱绕流的ＳＵＰＧ有限元解法.pdfVIP

高雷诺数时旋转圆柱绕流 的SUPG有限元解法 夏南 (上海大学，上海市应用数学和力学研究所，上海200072) 摘要 本文用SUPG方法求解ir旋转圆柱非定常绕流问题．研究一个无限长圆柱在不可压缩粘性流 体中垂直于其轴突然起动，同时绕轴突然旋转的情况．此流动产生有趣的非定常分离现象，在紧靠圆柱 表面的边界层发展成一环流区，后破碎形成尾迹．分离现象从运动壁面出现该问题存在三种基本的涡脱 落方式。依赖于参数胎和圆柱表面旋转速度相比于平动速度的参数函这样的a≠0的情况比萨0复杂- 在理论上：有重要意义．本文进{亍了Re：1500情况下的数值计算，给出了涡脱落的图像，获得了满意的结果． 关键词有限元，不可压牯性流，非定常流 1前言 对于N—S方程．应用有限元法通常是从加权残数法出发，目前使用最广泛的是 缩N—S方程时，往往只能计算到雷诺数几百左右，很难再算上去．自1976年以来，有限 提出了各种迎风有限元理论，对权函数进行修正，采用不对称的权函数。增大节点上游方 向的权函数，而减小下游方向的权函数。Chrisie的讨论限于一维对流扩散方程，后来 嗣格也能得到比较精确的解，但若将它们用：F多维问题或有浓问题时，解常在与流向相垂 年提出了SUPG(Stream]jng—Upwind 动方向起作用的摄动项．取得了迎风法所具备的长处，同时又避免了迎风法存在的虚假扩 散问题。本文用SUPG方法求解了旋转圆柱非定常绕流问题。研究一个无限长圆柱在不可 压缩粘性流体中垂直于其轴突然起动，同时绕轴突然旋转的情况。此类问题有两个基本无 量纲参数，一个是雷诺数，另一个是代表圆柱表面旋转速度与平动速度之比的参数口，定 义萨a孵7比．近年来，绕圆柱的非定常粘性流动更多的引起了人们的关注，因为由此流动 产生有趣的非定常分离现象，在紧靠圆柱表面的边界层发展成一环流区，后破碎形成尾 迹，分离现象从运动壁面出现．这样的删的情况比a=O复杂，在理论上有重要意义，这 的情况，得出一个结论：该问题存在三种基本的涡脱落方式，依赖于参数Re和口。本文进 行了R萨1500，萨l情况下的数值计算，获得了满意的结果。 ·326· 2基本方程 对圆柱问题．RP=型，D为直径。方程组为 卜 }1 V晓 塑反 影瓦 塑彦 2～以 (I) ff 塑a吻 塑砂咄 (2) 。鎏霎鳘兰鼋条件：远场条件，为均匀来流条件，有狄里赫利和诺依曼两种提法．本文中 ‘’一 取诺依曼形式。以s分别表示法向和切向。 警睁一_ (3) 固壁条件，yk=0，反映壁面无穿透。 门) 涡量边界条件：远场条件，取∞lL=O，计算中发现，涡量在远场衰减很快，外边界条 件取害：o和国。o差别甚微。 固壁条件：采用壁涡公式．对静止固壁： 嘶壁涡公式：‰2寿(‰l毗) (6) 二阶壁涡公式朋v (7) 2寿(‰嘿+J)一r-C2w+型 对旋转圆柱问题：壁面速度圪=(半)。。 国w2砑3({f，。q△n一‰。)一等 (8) 圆柱问题采用的网格为同心圆周线和径向射线构成的四边形网格，周向等分，径向南外逐 }坚 渐变宽，表达式为：，j=g“，女=0,1…2…， 3有限元方法 对涡流函数解组-Ⅳ-s方程以矿=w十p为权函数，泊松方程仍以矿为权函数。 f￡(害+矿