基于湍流自旋理论的二维方腔旋转流动的数值模拟.pdfVIP

第 30卷第5期 华 北 水 利 水 电 学 院 学 报 Vo1．30No．5 2009年 lO月 JournalofNoahChinaInstituteofWaterConservancyandHydroelectricPower 0et．2009 文章编号 ：1002—5634(2009)05—0004一O4 基于湍流 自旋理论的二维方腔旋转流动的数值模拟 董 荭，刘 云 (华北水利水 电学院，河南 郑州450011) 摘 要 ：基于湍流 自旋理论对旋转方腔 内不可压流体的流动进行 了二维数值模拟 ，用有限容积法在交叉网格 上对控制方程进行离散 ，采用 SIMPLEC算法进行计算 ，对流项 的离散方面采用 了延迟修正 的 QUICK格式 ，最 后采用二维的 TDMA算法进行迭代求解 ，并对边界条件处理采用 了Block块技术．同时与经典理论 即Navier— Stokes(N．S)方程 的模拟情况做 了比较．得 出下列结论 ：N—s方程表 明方腔旋转时内部流体只会作刚体旋转 ；而 自旋模型表明该类流动具有不平凡性和不稳定性 ：即流体 的内聚力不足 以使流体跟随固体边界作刚体运动 ， 旋转边界必然在流体中诱导出小尺度耗散涡旋 ，这些小尺度涡旋在固壁附近组织出正负相 间的拓扑缺陷 (涡 旋场强)，这一流体定向奇性会与速度场耦合 ，并容易形成复杂的大尺度 的涡旋运动． 关键词 ：湍流 自旋理论 ；方腔 ；旋转流动 ；数值模拟 中图分类号 ：0351．2 文献标识码 ：A 湍流流动是工程技术领域与 自然界中常见的流 _2 叫 v2V (AA + ) 一 一 动现象 ¨。 ．传统的湍流统计理论从细观层次来研 究小涡的局部级 串、相似结构和普适规律等．湍流运 (2AaV2+VzA) (1) 动的 自旋理论 则从宏观层 次来研 究小涡的宏观 +2toV。一 v V2+ (AA +AB) ：一 2P+ 效应 、大涡结构以及它们与速度场 的耦合动力学过 程．自旋理论将湍流的描述建立在统计平均意义上 ， (2A Vl+ 1A) (2) 即用统计平均的场量和如粘性系数等宏观物性参数 2．自旋场方程为 描述湍流，而且认为在流体单元内部分布着大量随 +2∞：AV 一AOdA (3) 时间涨落的耗散小涡 ，由于这些耗散小涡的统计平 3．弯扭场方程为 均效应 ，使得流体单元具有除速度 以及速度梯度带 VllV2一 1Vl+ Al=AO(1 一O,A1)一 来 的旋率以外的独立变化的涡旋量． ／l：(vh9。A：一 A) (4) 二维粘性方腔流为不可压缩流体 ，边界为正方 V12 — 2Vl+ A2=A (2 一O,A2)+ 形物面边界 ，在动坐标系上研究其流动现象 ．笔 A,9( a。A2一 02A。) (5) 者用 Navier．Stokes方程 和 自旋理论对二维坐标下 由上述方程可看出，在旋转标架下得出的控制方 的不可压缩方腔旋转 流动进行计算，考察涡旋场 程比在固定标架下得 出的控制方程只多 出了科氏加 如何产生、定常流动时涡旋场分布及对速度场的影