两类尾迹型剪切流中旋涡位错.pdfVIP

两类尾迹型剪切流中的旋涡位错 凌国灿1 赵红亮2’1 牛家玉3 (1．中国科学院力学研究所LNM国家重点实验室，北京100080 2．解放军理工大学理学院流体力学研究中心，江苏南京211101 3．中国科学院力学研究所LNM国家重点实验室，北京100080) 摘要本文用高精度紧致有限差分与Fourier谱方法，通过对三维不可压N-S方程的直接数 值模拟计算，研究具有展向不同非均匀性的尾迹型流动时空发展的基本特征和转捩性态。所研究 的流动包括来流沿展向呈指数衰减型的局部非均匀流、台阶型非均匀流等两类典型情况。数值结 果表明流动不稳定性的非线性发展及空间非均匀发展导致旋涡脱落在相位、频率和强度上的差异。 在变化率最大的地方出现晨向涡量向流向及垂向的显著转变，形成不同类型的旋涡位错。用流动 关键物理量及涡线详细刻画了旋涡位错的结构特征和动力学过程，描述了旋涡位错中复杂的三维 涡链。利用FFT方法分析了流动中速度的时空变化及旋涡位错对流动转捩的影响。 关键词旋涡位错，尾迹型剪切流 1 前言 在尾迹、混合层等许多实际流动的转捩过程中发现有一系列不同类型大尺度、形态复杂的涡结 构，称之为“旋涡位错”(‘Vortex 性的一个基本特征，是流动转捩及向湍流过渡的一种新的动力学途径和机制。研究旋涡位错及其 产生的空间结构特征，动力学机制对理解流动转捩、流动掺混等具有重要的理论意义和深刻的实际 应用背景。实验显示和测量在对于旋涡位错的形成、特征和动力学机制方面提供了逐渐深化的信 息和认识，但对旋涡位错的结构特征，实验显示难以提供详细信息，在动力学方面还很缺乏深入研 模拟出了尾迹型流动及圆柱绕流尾迹中的旋涡位错结构，对旋涡位错的形成和随机性的旋涡位错 发生和对流场影响作出了重要的数值研究探索，但对于流动的展向涡结构演化成为三维位错结构 的动力学过程和机制的研究还不充分，特别是实际流动旋涡位错中的涡联的复杂性和具体详情还 没有得到清晰的描述。旋涡位错对流场特性及转捩特征的影响也需要进一步探讨。基于这些考虑 2001)工作基础上利用不可压N_s方程的高精度格式对尾迹型来流进行进一步的 我们在文(Ling 直接数值模拟研究，对流动的时间一空阉演化作更充分的数值模拟，即更大的空间计算域，用更细 的分辨作长时间的计算模拟，并提供新的结果。为了避免绕柱流动模拟的大规模计算，也是为了通 过模型流的研究，对旋涡位错有较深入的基本认识，我们研究两类典型的具有展向非均匀的尾迹型 剪切流动，即对于模型来流剖面的亏损强度给以局部指数衰变型扰动，以及对来流剖面速度型的特 征尺度沿展向引入台阶型变化。这两种情况与柱体直径有一个局部小变化或柱径有台阶型变化的 ·186· 流动有一定的相似性。 2数值方法 旋涡位错结构是发生在三维开放型流动中的大尺度流动结构。因此如果用数值方法研究这类 问题就需要用时空模式来模拟。也就是说当大尺度涡结构向下游发展时流向尺度计算域要足够大 和合适的人流、出流边界条件来适应这一过程。为此，这里用紧致差分～谱方法来求解原始变量形 式的不可压缩三维N—S方程组。N—S方程采用三步时间分裂格式 J．一l 矗一∑a五rt ，，一， —}一蚤驴m[：～V)加 兰[A=_堕=一v，户寸1t 】m—q 一， … 盐蛰At二虽：一1Rev卵 Ling(2001)，这里仅给出计算的相关参数。 人流为尾迹型动量亏损剖面为 (2) U(y，z)一1．0一D(z)．(2．0一cosh(b(z)，)2e_‘出)，)2 c(z)组合表示该剖面的特征长度参数。 3 数值结果显示与分析 3．1来流沿展向呈指数衰减非均匀 通过n(z)引入沿展向的局部非均匀，