可压缩槽道流的标度律和流动结构相关性分析.docVIP

可压缩槽道流的标度律和流动结构相关性分析* 裴杰1，胡凯衡1,2，佘振苏1,3 北京大学工学院，湍流与复杂系统国家重点实验室，北京 100871 中国科学院成都山地灾害与环境研究所, 四川成都, 610041 Department of Mathematics, University of California, Los Angeles, Los Angeles, CA 90095, USA Email: she@pku.edu.cn ,010摘要：本文对直接数值模拟所得到的Ma数为0.8和1.3的槽道湍流流场（Re=3000）进行了统计分析，计算了速度结构函数六个分量的相对标度指数，并将计算结果与流场可视化结构进行对照比较，研究湍流流场的奇异标度特性与流动结构形态的相关性。结果表明，在槽道的中心区湍流发展比较充分，结构性弱，具有很好的ESS标度指数，最接近于各向同性湍流，但仍然不是完全的各向同性湍流；在槽道的近壁区和对数区只有流向速度沿流向和展向的分量（，）的标度行为良好，法向速度分量（）的间歇性最强，结构性也最强。随着Ma数的增加，近壁区的条带结构更有序，间歇性更弱。在过渡区、部分对数区和中心区，间歇性基本不受压缩性影响。在大部分对数区和中心区，压缩性影响具有不确定性。这些分析表明，湍流流场结构形态与奇异标度律的测量值之间存在一定的关系，奇异度和间隙性强弱是随着流动结构的变化而变化。 关键词：可压缩湍流，标度律，间隙性，层次结构 引言 湍流的标度律刻画了湍流不同尺度的脉动结构之间的统计相关特性，是对多尺度复杂系统临界自组织状态的定量刻画。 标度律在K41理论[1,2]出现之后越来越受到人们的重视。该理论认为：当雷诺数充分大时，标度律是普适的。近几十年来的实验[3,4]和数值模拟[5,6]发现许多封闭性流动表现出比均匀各向同性湍流更小的标度指数。从动力学角度来讲，湍流是一种包含着从大尺度到小尺度的能量级串过程的非平衡介质，能量级串引起的多级相关反映在各阶速度结构函数的标度指数具有非线性的依赖关系，这就是著名的反常标度律问题。 湍流的唯象理论模型采用基于能量级串的随机乘子过程来解释反常标度律[7]，但是这个模型忽略了对流动结构，如条带，发夹涡等的分析。湍流的She-Leveque层次理论模型(SL94)[8,9,10,11]引用的层次的概念来描述湍流这个多尺度、多幅度随机场，根据最高激发态这一具体结构从统计和结构上描述了湍流中的有序结构和不规则运动，解释了惯性区的间歇现象。 80年代以来，湍流的直接数值模拟有了很大的发展，但结果大部分是针对不可压缩流动的。而实际应用中的流动有相当一部分是可压缩的，比如航空航天领域涉及到的湍流均是可压缩湍流。可压缩槽道流的物理现象比不可压缩槽道流更复杂，数值计算要求的精度要更高，实验的条件更加严格。在可压缩湍流的研究中，对流场的拟序结构分析是相当重要的一项。 数值模拟 直接数值模拟是使用精确的数值计算方法直接求解控制湍流运动的N-S方程，不附加任何湍流模式。可压缩槽道流作为一种可压缩的壁湍流，包含了声速效应和壁面效应的影响，是一种有代表性的典型流动。 本文中对二种不同Ma数的槽道流开展了计算，计算区域为，网格在流向(x)和展向(z)上划分为的均匀网格，在法向(y)用140个非等距网格以保证近壁处分辨率（详细计算条件见[12,13,14]）。速度场的初始分布为无量纲化的Poiseuille层流解加上一定振幅的均匀随机拢动。 壁面处的边界条件为：无滑移条件和恒温条件。 Ma=0.8和 Ma=1.3的计算结果比较 Morkovin理论认为，湍流场中的压力脉动值与其平均量相比不是很大时，压缩性主要影响平均量[15]， 表1列出的二个流场的初始参数和部分结果参数表明：随着Ma的增大，压缩效应使壁面处的平均密度变大，槽道中心的与壁面的粘性比增大，而壁面摩阻结果相近。 (a)(b) 图1： (a) 平均速度剖面，点线是不可压槽道流的流向平均速度的壁面律和对数律；(b)为平均密度变化曲线。 表1：流场的相关参数比较(Ma=1.5为Coleman等[16]的计算结果) Case Re数 Ma=0.8 3000 18.65 6.21 0.41 5.68 190 1.11 0.0035 1.076 Ma=1.3 3000 20.59 6.86 0.246 7.33 209 1.276 0.0038 1.198 Ma=1.5 3000 19 12 0.1 ___ 222 1.335 0.004 1.252 近壁区、对数区和中心区的结构及标度律分析 条带结构是近壁区的最重要的拟序结构。图2-图3为二个Ma数的可