层流—湍流转捩的-breakdown-过程的内在机理.pdfVIP

近代空气动力学研讨会论文集 周恒罗纪生王新军 (天津大学力学系南开大学天津大学刘辉应用数学中心，天津30072) 摘要层流一湍流的转捩是一个古老但依然没有解决的问题。研究的成果主要集 中在以下几个方面：1)扰动演化的线性与非线性理论；2)感受性问题，即扰动是如何 产生的；3)转捩的实际预测。转捩的文献很多，这里无法一一列出，但从(AnnualRe viewof FluidMechanics}的综述文章中可追索到该课题的研究历程。尽管对转捩的 研究付出了巨大的努力，但仍然有一个重要的问题还没有涉及到，那就是最终导致转 捩的breakdown的内在机理是什么?传统的看法是，转捩从扰动的放大开始。由于 非线性的作用，当扰动增长时，高次谐波将产生，流动变得越来越复杂，最终导致湍 流。这一过程似乎非常清楚，但却漏掉了一个环节，即转捩的breakdown过程究竟发 生了什么。下面我们将对直接数值模拟的结果进行分析，分析的结果表明，平均流剖 面稳定性特性的改变在breakdown过程中起了关键性的作用。 关键词转捩Breakdown过程直接数值模拟 1已知的结语 首先，我们将列出一些有关湍流和转捩的已知事实，然后给出一些有关槽道流转捩数值模 拟的分析结果，这一结果将揭示转捩的breakdown过程的关键机理。 平均流剖面的变化在转捩中起着关键的作用。一个明显的解释是，层流与湍流具有截然 不同的平均流剖面，但这并没有提供breakdown发生的内在机理。 从流动稳定性理论的观点来看，当雷诺数足够大时，层流的平均流剖面具有不稳定的特征 解。在流动的法向此特征解占据流动的整个区域，因此特征函数具有整体的特性。 与层流相反，湍流平均速度剖面不具有整体的不稳定的特征解。从20世纪60年代发现 相干结构以来，经过30多年的研究，人们已经普遍认识到，对平面槽道流和边界层流等带有固 壁边界的湍流，存在两种不同的流动区域，即壁面区与湍流边界层外区或槽道流的核心区。前 者仅占整体流动区域的-／j,部分，而后者通常占据整个流动区域的90％以上。Tsujimoto与 3的研究结果表明，虽然没有外区的脉动壁面区的湍流将不能维持，但壁面区具有自 Miyake[1 身的动力学特性，如湍流强度、雷诺应力分布等都主要依赖这一动力学特性，而不是依赖外区 脉动的性质。众所周知，壁面区的动力学特征是产生相干结构的根源，而相干结构又从平均流 中吸收能量，并通过所谓的喷射和扫掠与边界层外区相互作用。根据这一事实，周恒”~51与 他的同事们提出了近壁区产生相干结构的动力学模型。 湍流平均速度剖面的一个明显的特征是法向速度梯度在壁面区很大，而在湍流边界层外 区或槽道流核。fi,lX却很小。这两个区域中间为缓冲区，在这里速度梯度从很陡变得平缓。两 个性质截然不同区域导致了不稳定波整体特性的破损。如周恒等所述，本质上，相干结构是由 层流一湍流转捩的“breakdown”过程的内在机理 55 于流动的某种不稳定性产生的，他们用流体不稳定模型描述了相干结构。这类不稳定性具有 较强的局部特性，即描述相干结构的特征解在边界层外区迅速衰减，因而与远离壁面区的扰动 没有联系。因此可将槽道湍流描述为三个截然不同的区域中的扰动相互作用的结果，这三个 区域分别为两个壁面区和一个中心区。湍流能量主要由两个壁面区域中相干结构的生成提 供。中心区充当相对被动的角色，即通过相干结构的喷射从壁面区接受湍流能量，相应地，中 心区也将在壁面区内引起一些扰动，从而产生新的相干结构。性质截然不同的壁面区和边界 层的外区或槽道流的核心区的存在显然是平均流修正的结果。 2槽道流转捩的直接数值模拟 我们用平面槽道流来研究转捩。数值方法为伪谱方法；即在流向和展向采用傅里叶变换， 而在法向采用两点4阶精度的紧致差分格式来解傅里叶展开所产生的方程。 计算起始于具有抛物线型平均速度剖面的层流，然后引入扰动。扰动由三个T—S波组成， 形式如下： (2) 其中z是流向坐标，z是展向坐标，Y是法向坐标捣与岛(i=1，2，3)分别是32和2向的波数。 化条件为maxl“】I=1，“1