圆柱尾迹中的转捩波及其和卡门涡相互作用的特征研究.pdfVIP

圆柱尾迹中的转捩波及其与卡门涡 相互作用的特征研究 杨毅 凌国灿 (中科院力学所非线性力学国家重点实验室(LNM)，北京100080) 摘 要 本文采用对二维圆柱绕泷N—S方程直接数值模拟的方法详细计算了霄诺数 结出了转捩渡频率^及其能量随Re变化的规律，I；l及^能量在空闻上的变化规律；本文通 过计算发现转捩波频率^与卡门涡脱落频率，的比值，f腻随尺e的麦化规律在雷诺数 幂。本文还对上述计算结果以及流场涡结构特征提出了一些解释。 美键词 非线性相互作用，二雏圆柱尾连流动，转捩坡 1 前言 圆柱尾迹随雷诺数增加，分离剪切层不稳定性发展将产生频率高于s。数频率的转捩波 (Transition Waves)，形成周期性产生的小涡．剪切层转捩方面的研究工作最早是由Roshko (1 954)开展的，其后在实验研究上最早由B|oor(1964)探测到剪切不稳定渡的存在。在转 捩发生的区域中，Bloor通过数一个基波波长内的转捩波波数得到转捩波频率与基频的比 etal(1987)的实验 值，这个比例与R一“成正比．即，，／，JOCRe“2，这一结果得到了Kourta 研究工作的支持。而Wei 稳定波频率与RP“”成正比。Prasad 稳定性几乎完全是二维的，而且是间歇性地显示的．他用自己的实验结果结合了前人的结果 eta1(1992)风洞实验 通过数据拟合给出的幂次关系是：^，，，=0．0235RP’”。J．Sheridan 则发现在低雷诺数(1600，4000)下的实验结果与RP“5的幂次关系符合得很好，而在高雷诺 数下的实验结果则与0．5的幂次关系有较大的偏差。在理论数值计算方面，Braza(1990)在 尺Po 算结果接近WeiSmith的0．87的幂次关系。 从以上各人的实验以及理论计算结果可以看出圆柱尾迹转捩波频率与交替脱落涡颧率 比随雷诺数变化规律从Bloor的0．5次幂到Wei 0．67次幂各不相同。反映出人们对这两类不稳定性相互作用的机制与特征还未得到一致的 认识和结果，也直接表明人们对于较大雷诺数流动中圆柱尾迹从层流向湍流的转捩的发生、 发展的过程还没有得到清晰的认识，因而对此问题进行进一步的研究具有非常重要的理论 30 和实际意义(Ling，G．C．8dMonkewitz，P．1999)． 2．1数值方法 本文数值求解的控制方程是对二维原始变量形式的Navier--Stokes方程，采用圆柱半 径R、无穷远来流速度U。分别作为长度、速度尺度进行无量纲化，无量纲变量定义为： z=丢．，=丢一警，一一磊 “=矗，”一矗m=学 无量纲化后的二维原始变量形式的Navier--Stokes方程为： 口·V=0 詈+(V·V)“一一差+恚V2“ (2) 害+(y·V)”一凳+毛V2。 (3) 压力Po[sson方程为： (4) V：P警+毛Vn[警+2哿+等] 其中： D=塞+爵By 在本文数值计算中，坐标系固联在圆柱上，初始速度场【k=1。在圆柱表面速度采用无 滑移、无穿透边界条件。来流边界的速度取U。=1．时间项的离散采用二阶精度的分裂步法 求解。对流项和牯性项采用二阶精度的中心差格式离散。本文计算在柱坐标下进行． 采用本文所用的计算程序，计算了雷诺效Re=i000一li000范围内的多个雷诺数下 的流场，对每个雷诺数下计算所得的速度分量“，”及压力分量p分别作快速傅里叶变换 作功率谱计算，样本