边界拟合坐标涡量一流函数绕物体流动计算.PDFVIP

第25卷 第6期 重庆 建筑大 学学报 V01．25 No．6 2003年12月 Journal of Chongqing Jianzhu University Dec．2003 文章编号：1006—7329(2003}06—0058—04 边界拟合坐标涡量一流函数绕物体流动计算 甘 孜 (重庆交通学院 河海建筑工程系，重庆 400074) 摘要：列出了用边界拟合坐标表示的涡量一流函数绕流计算偏微分方程组，及其交替隐式 有限差分和显式有限差分方程组。用自动调控边界计算了矩形区域单园柱绕流运动，具 有无穷域绕流特征。绘制了速度矢量分布图，和流函数，涡量，两个分速度及压强的等值 线分布图。 关键词：涡量一流函数；边界拟合坐标；有限差分；绕流计算；无穷域；自动调控边界 中图分类号：0351．2 文献标识码：A Around Body Flow Computation with Vortex—Flow—stream Function in Boun dary Fitted Coordinates GAN Zi (Depamnent of River and Ocean Construction Engine I】g，C ng Institute of Communications，Cho~qing 400074，P．R．China) Abstract：IIl this paper．partial differential equations of vortex flow—stream fimction in boundary fitted~oor- dinates and related ADI equations were listed ．The flow around a cylinder in quadrilateral area by auto control boundary was computed．The diagram of distribution of speed alTOW，constant lines of flow—stream function， vortex，two speed parts and pressure were studied． Keywords：vo~ex—flow—strealTl function；boundary fitted coordinates；finite difference；flow computation； irrfinite region；automatic boundary 目前绕流偏微分方程计算普遍存在无穷域流动边界有限问题。一般提出计算矩形域两侧边界条件 为光滑固壁流线即第一类边界条件(见参考文献[23第114页)，与无穷域绕流问题有一些差异。为 缩小差异接近无穷域绕流，一般方法是扩大边界。这导致计算量过大。在无穷远来流绕矩形区域 内单园柱体或物体群的流动计算的力学模型中，本文采取涡量一流函数方法，偿试使用主控方程自 动调控边界流动，企图实现无穷域绕流流动条件，即涡量耗损扩散性透过外边界直达无穷远域。计 算结果显示外边界流动具有无穷域绕流流动特征。 偏微分方程自动生成自适应边界拟合坐标，物理域网线分布合理，弯曲弧线边界网格点均匀分 布，计算域边界网格点排列整齐，有利于多物体群绕流各式边界处理及涡量一流函数计算稳定。有 关边界拟合坐标网系的生成见参考文献[5]。 拟合坐标涡量方程(1)式是对流扩散方程，流函数方程(2)式是伯松方程，压强方程(3)式是伯 松方程。外边界大部和内部域的流函数计算用(10)一(11)式。内部涡量计算用方程(6)一(7)式。 全部边界和内部压强计算用(13)一(14)式．以上各计算式组成交替隐式差分方程组，是无条件稳定 的。外边界及内部速度计算按差分方程(17)式。压强项不出现在(1)式中，根据速度即可计算压 收稿日期：2O03—08—28 作者简介：甘孜(1960一)，男，重庆人，硕士，讲师，从事水力学教学及研究工作。 第6期