近地面大气湍流平均动能耗散率测量与分析.pdfVIP

维普资讯 第 39卷 第 6 期 力 学 学 报 Vo1．39．NO．6 2007年 l1月 ChineseJournalofTheoreticalandApplied M echanics NOV．、2007 近地面大气湍流平均动能耗散率测量与分析 吴晓庆 。) 聂 群 方 强 (中国科学院安徽光学精密机械研究所大气光学中心，合肥 230031) 摘要 用三维超声风速计测量了合肥地区风速脉动和温度脉动数据，并用二阶和三阶径向风速结构函数以及相 似理论估算了湍流动能耗散率，估算结果基本一致．其与稳定度相关性分析表明．在大气处于中性条件下，湍 流动能耗散率最大，并随 lz／Ll值增大而下降；Kolmogorov耗散尺度 叩与 ￡l／成反比，其与稳定度 z／L的 变化趋势与s随z／L的变化趋势正好相反； 在大气处于中性条件下最小，并随 lz／Ll值的增大而增大， 但在不稳定条件下 增大得更快．热力湍流内尺度不能用湍流动能耗散率 s进行计算，而应引入与温度0梯 度有关的耗散率 ￡． 关键词 大气湍流，平均动能耗散率，稳定度，热力湍流，内尺度 中图分类号： P425．2+1 文献标识码： A 文章编号： 0459—1879(2007)06—0721—06 引 言 确定因素带来的误差．王汉封 【J通过 PIV测量数 据，运用大涡模拟中的亚网格应力模型估算出槽道 有许多外场实验和实验室实验都证实大气近地 内湍流动能耗散率．由于PIV测量的是一个平面上 面层速度场和温度场存在惯性子区，其表现形式是存 的速度分布，该方法 目前只能实现对近似二维的或 在 一5／3功率谱幂率，二阶结构函数与空间两点间距 各向同性流体的湍流动能耗散率的估算．本文比较 离呈2／3指数关系．惯性子区最主要的特征是局地各 了3种湍流平均动能耗散率计算方法．通过三维超 向同性，其物理统计量 (速度结构函数、温度结构函数 声风速计测量的近地面层大气速度和温度脉动信号 等)随尺度 的变化通常是呈幂次律，与大尺度脉动统 得到了湍流动能耗散率，分析了湍流动能耗散率与 计特性、黏性耗散机理和流动环境无关，这一假设是 大气稳定度的相关性，指出热力湍流内尺度不能用 Kolmogorov(K41)湍流理论的基础．但由许多实验证 湍流动能耗散率进行计算． 实，在阶数大于2时，结构函数与空间两点间距离存 在的指数关系与K41理论发生偏离．援引Landau[] 1湍流平均动能耗散率的计算方法 和Lifshitz[2J的说法，偏离的原因之一是由于惯性区 局地湍流平均动能耗散率 e定义为 内存在间歇性，湍流动能耗散率存在起伏．于是出 现了一些对 K41进行修正的模型，如 模型、对数 e 正态模型和其它多重分形模型． Kuznetsov[3_把这 ~=-1j=l(考+Oxi]＼。一 些模型中出现的间歇性称为内间歇性，因为模型只 考虑了耗散率的起伏．湍流动能耗散率的传统测量 ， J分别代表风速分量， (=1．5×10 m2．s )是空 方法是采用热线风速仪，从测量的速度脉动时间序 气分子黏滞系数， ()代表半径为r的空间平均． 列信号中分别得到能谱、空间相关系数或脉动速度 在满足局地各向同性假定下，空间平均简化为 的空间梯度等再计算出湍流动能耗散率．热线风速 线性平均，式 (1)可写为 仪测量需要很高的采样频率，且探针本身结构对流 场的影响以及各热丝间相互影响会带来一定的