径向湍流旋转射流的计算研究.pdfVIP

中固工程热撕理擘管 热机气哥鼎力学 第 十 ‘ 丰 量 壕号：012∞1 径向湍流旋转射流的计算研究 吕盘明 (中国科技大学热科学和能源工程系．230026，pmlu@ustc．edu．cn) 一．引言 径向湍流旋转射流在工程上有很多应用，比如喷雾，干燥，加热，炉膛上煤，上 料，液体与固体颓粒的溶化等等，都可因为射流旋转的帮助而进行得更加有效。实验研 究有文献[1】，计算研究有文献[2]，用的湍流模型是k—F EvM，但计算所得的射流扩 张速率与实验相比小了很多，计算是0．I，而实验是0．139。文献[2]怀疑这是由于标准 的占一方程并不完全适合于有旋转的流动．作者用能加入湍流结构的边界条件的湍流模 型．}-gLAEⅦ已经计算过各种平面的和轴对称的自由射流，发现加入合理的上游边 条后，计算出的射流扩张速率和实验符台得都很好[3，4】。至于通道内的旋转流动也用 七一占I．ADSM和露一sIAEVM计算过[5】。因此作者相信七一sLAEVM对改进径向湍流 旋转射流的预铡也应当有效。 二，数值求解 为求解径向湍流旋转射流，我们用时间平均的N—s方程组(用茁卡儿张量表示)： (p矿』WJ=0 (1) (pUUjI，=一P。也∥。I，州，沪硐， (2) · 并用能加^湍流结构的边界条件的二方程湍流模型表示式中的Reynolds应力“，打， 即让 一pu．“J==一，FJ。ji=+∥，【u。l，+u，Ii】 (3) 其中 “=q球2／6 (4) p筹一v．譬雅】+矿一班 (5) p告再【等V小PK(C,。P-c∥) 【6) + 本文由国家自然科学基金项目资助．N 积 2+c与r 计盘 ☆。 一加’… P叫饼c争2+c争cy2H詈+事2+c争2巾嘉禹r】2)∽ 方程(3)中这个Y“是称作湍流局部结构张量Y的协变分量，它由湍流结构的上游边 界条件及同壁边界条件等决定，计算)0的具体方法详见参考文献[5，6]，它的基本 做法都是给定上游射流中心的y妒和射流边沿的儿，，然后由(线性)内插值决定其它 地方的yF·并假定yF沿平均流线不再改变。其它通用的模型常数见表一。 表一： 湍流模型常数 f K q c￡- ％ ak 瓯 1 0．42 0．09 1．44 1．92 1．O 1．3 以上方程(1)(2)(5)(6)可以表示成以下统一的形式，然后联立求解． (8) 砉(∥7脚)+砉(∥’啪=砉【，，‘知+舌【，’‘争∥ 其中r的上标j取值1(因是圆柱坐标系)．交换系数r．和源项s。的表达式见表二。 计算区问：可有两种取法：1．可取如图一1由ABcDGH为顶点的相连成一片的两个四 边形区域分两次分别计算。 2．也可取如图一1由BcGH为顶点的一个四边形区域一次计 算。这后一种方法本文作者试过，但是还未成熟，