论文部分内容阅读
层流-湍流的转捩是一个古老但依然没有解决的问题。百余年来,工程实践的需求使得人们对这一问题的研究不断深入。尽管对转捩的研究付出了巨大的努力,但是仍然有一个重要的问题还没有被涉及到,那就是最终导致转捩的“breakdown”的内在机理是什么?传统的看法是,转捩从扰动的放大开始,由于非线性的作用,当扰动增长时,将产生高次谐波,流动变得越来越复杂,最终导致湍流。这一过程似乎非常清楚,但却漏掉了一个关键的环节,即转捩的“breakdown”过程究竟发生了什么。本文对转捩的过程进行了直接数值模拟,对结果的分析表明,平均流剖面稳定性特性的改变在“breakdown”过程中起了关键性的作用。目前,能用于工程中湍流计算的只有湍流模式理论。在大多数的做法中,其雷诺应力最终都归结为通过某种涡粘系数与平均流的变形率联系起来,从而使方程封闭。不同的模式或亚格子应力模式的差别只在涡粘系数的不同算法,由于所用的涡粘系数都是实数,所以无论哪种方法都没有考虑到雷诺应力与速度变形率的相位关系。近年来对湍流涡粘模型的研究表明,目前所采用的湍流涡粘模式中的涡粘系数不能都是实数,这对湍流模式中的涡粘模型提出了新的挑战。为了给出更准确的涡粘系数,本论文采用数值模拟方法,系统研究了具有周期性吹吸边界条件的湍流流场中的雷诺应力与速度变形率的相位特性,并研究了雷诺数、壁面吹吸幅值以及移动吹吸边界条件等因素对这一相位特性的影响。重点对雷诺应力与速度变形率的相位差进行了分析,说明在此类流场中涡粘系数应该考虑相位的特性,并对此条件下的涡粘系数形式进行了讨论。