可压缩各向同性湍流湍动能衰减规律.pdfVIP

李虎，刘旭亮等：可压缩各向同性湍流的湍动能衰减规律 可压缩各向同性湍流的湍动能衰减规律 李虎，刘旭亮，张树海 (中国空气动力研究与发展中心，空气动力学国家重点实验室，四川绵阳，621000) 摘要：本文采用五阶WENO格式直接模拟了初始Taylor微尺度雷诺数为72～153，初始湍流马赫数为0．3～ 1．5的可压缩各向同性湍流，目的是研究可压缩均匀各向同性湍流的湍动能衰减规律。研究表明，可压缩均 匀各向同性湍流湍动能存在一个幂指数形式的衰减，相应的衰减指数n=1．4～1．6，接近不可压缩均匀各向 同性湍流统计平衡态前的弛豫过程，也就是湍动能初期衰减指数的理论预测值1．43：这表明压缩性对可压 缩各向同性湍流的湍动能衰减指数的影响不大。 关键词：可压缩湍流，各向同性，DNS，WENO，幂律衰减 0引言 自相似性解得到的衰减指数n=l。Clark和Zhou 等【_7】的研究也支持这一结论，并指出完全自相似 从Batchelortl】和Saffman[21等人的经典工作 我们知道，目前存在两种标准类型的均匀各向同 解对应的初始能谱形式为E(扣O)=七。George[8】 性湍流问题，具体的能谱形式为能谱形式为 声称找到了适用于全范围Reynolds数的完全自 相似解，指出相应的湍动能衰减指数是依赖于初 以后_÷0)—G矿，其中0-=2是Saffman型湍流，o--4 是Batchelor型湍流。然而，究竟湍流属于Saffman始条件的。George的理论是基于对能谱动力学方 型还是Batchelor型要取决于初始条件，也就是湍 流初场的生成机制。 Leonardp]也发展了一种新的关于二阶和三阶速 在均匀各向同性湍流的衰变过程中，湍动能 度相关函数的完全自相似假设，该假设适用于中 的衰减将呈现为负指数的幂函数形式，即 低Rel数，并且满足Kannan．Howarth方程． OC “2 t～，／／为衰减指数。根据Batchelor的理论如果假设相关函数是部分自相似的，也就是 p巧】，这一衰变过程可划分为三个阶段：初期衰减， 相关距离只有在一定范围内时相关函数才是自相 过渡期和后期衰减。在各向同性湍流衰减的初期， 似的，可以得到Karman．Howarth方程的部分自 湍流还处在发展之中，此时惯性力占支配地位； 相似性解，而这种情况下得到的衰减指数并不为 在衰减的后期，在全波数范围内粘性力将超越惯 1。放宽完全自相似假设中长度尺度必须为Taylor 性力起主导作用；在过渡期，湍流运动机理逐渐 演变，能量的衰减由初期形式过渡到后期形式。 对应于各向同性湍流的初期和后期衰减，湍动能 的衰减指数存在两种不同条件下的极限值，分别 为Rex-÷∞和Rex_O。 量，推导出相应的衰减指数为n=6／5。 基于不可压缩流动的Karman．Howarth方程，在均匀各向同性湍流衰减的后期，湍动能的 结合二阶和三阶速度相关函数的完全或者部分自 幂律衰减指数是能谱形态的函数，由初始能谱 相似假设，在Rex_0的限制条件下，可以演绎出 以幼的卜÷0的低波数成分决定。也就是说，各向 第十五届全国计算流体力学会议论文集 同性湍流后期的衰减规律是依赖于初始条件的 ¨1|。因此在通过DNS研究均匀各向同性湍流的 的格栅湍流实验得到的涡量均方根的幂律衰减指 后期衰减规律时，需要区分初始时刻的湍流类型。 数为1．5+0．2，相应的湍动能衰减指数为2。Skrbek Batchelor和TownsendN指出Batchelor型均 和Stalpll91的液氦实验结果