湍流模型理论.docVIP

湍流模型理论

湍流模型理论

§3.1引言

自然界中的实际流动绝大部分是三维的湍流流动，如河流，血液流动等。湍流是流体粘性运动最复杂的形式，湍流流动的核心特征是其在物理上近乎于无穷多的尺度和数学上强烈的非线性，这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。回顾计算流体力学的发展，特别是活跃的80年代，不仅提出和发展了一大批高精度、高分辨率的计算格式，从主控方程看相当成功地解决了Euler方程的数值模拟，可以说Euler方程数值模拟方法的精度已接近于它有效使用范围的极限；同时还发展了一大批有效的网格生成技术及相应的软件，具体实现了工程计算所需要的复杂外形的计算网格；且随着计算机的发展，无论从计算时间还是从计算费用考虑，Euler方程都已能适用于各种实践所需。在此基础上，80年代还进行了求解可压缩雷诺平均方程及其三维定态粘流流动的模拟。90年代又开始一个非定常粘流流场模拟的新局面，这里所说的粘流流场具有高雷诺数、非定常、不稳定、剧烈分离流动的特点，显然需要继续探求更高精度的计算方法和更实用可靠的网格生成技术。但更为重要的关键性的决策将是，研究湍流机理，建立相应的模式，并进行适当的模拟仍是解决湍流问题的重要途径。

要反映湍流流场的真实情况，目前数值模拟主要有三种方法：1.平均N-S方程的求解，2.大涡模拟（LES），3.直接数值模拟（DNS）。但是由于叶轮机械内部结构的复杂性以及目前计算机运算速度较慢，大涡模拟和直接数值模拟还很少用于叶轮机械内部湍流场的计算，更多的是通过求解平均N-S方程来进行数值模拟。因为平均N-S方程的不封闭性，人们引入了湍流模型来封闭方程组，所以模拟结果的好坏很大程度上取决于湍流模型的准确度。自70年代以来，湍流模型的研究发展迅速，建立了一系列的零方程、一方程、两方程模型和二阶矩模型，已经能够十分成功的模拟边界层和剪切层流动。但是，对于复杂的工业流动，比如航空发动机中的压气机动静叶相互干扰问题，大曲率绕流，激波与边界层相互干扰，流动分离，高速旋转以及其他一些原因，常常会改变湍流的结构，使那些能够预测简单流动的湍流模型失效，所以完善现有湍流模型和寻找新的湍流模型在实际工作中显得尤为重要。

§3.2湍流模型概述

§3.2.1湍流模型的引入

湍流模式理论或简称湍流模型，就是以雷诺平均运动方程与脉动运动方程为基础，依靠理论与经验的结合，引进一系列模型假设，而建立起的一组描写湍流平均量的封闭方程组。湍流运动物理上近乎无穷多尺度漩涡流动和数学上的强烈非线性，使得理论实验和数值模拟都很难解决湍流问题。虽然N-S方程能够准确地描述湍流运动地细节，但求解这样一个复杂的方程会花费大量的精力和时间。实际上往往采用平均N-S方程来描述工程和物理学问题中遇到的湍流运动。当我们对三维非定常随机不规则的有旋湍流流动的N-S方程平均后，得到相应的平均方程，此时平均方程中增加了六个未知的雷诺应力项，从而形成了湍流基本方程的不封闭问题。根据湍流运动规律以寻找附加条件和关系式从而使方程封闭就促使了几年来各种湍流模型的发展，而且在平均过程中失去了很多流动的细节信息，为了找回这些失去的流动信息，也必须引入湍流模型。目前虽然许多湍流模型已经取得了某些预报能力，但至今还没有得到一个有效的统一的湍流模型。同样，在叶轮机械内流研究中，如何找到一种更合适更准确的湍流模型也有待于进一步研究。

§3.2.2湍流模型的发展历程

模型理论的思想可追溯到100多年前，为了求解雷诺应力使方程封闭，早期的处理方法是模仿粘性流体应力张量与变形率张量关联表达式，直接将脉动特征速度与平均运动场中速度联系起来。十九世纪后期，Boussinesq提出用涡粘性系数的方法来模拟湍流流动，通过涡粘度将雷诺应力和平均流场联系起来，涡粘系数的数值用实验方法确定。到二次世界大战前，发展了一系列的所谓半经验理论，其中包括得到广泛应用的普朗特混合长度理论，以及G.I泰勒涡量传递理论和Karman相似理论。他们的基本思想都是建立在对雷诺应力的模型假设上，使雷诺平均运动方程组得以封闭。1940年，我国流体力学专家周培源教授在世界上首次推出了一般湍流的雷诺应力输运微分方程；1951年在西德的Rotta又发展了周培源先生的工作，提出了完整的雷诺应力模型。他们的工作现在被认为是以二阶封闭模型为主的现代湍流模型理论的最早奠基工作。但因为当时计算机水平的落后，方程组实际求解还不可能。70年代后期，由于计算机技术的飞速发展，周培源等人的理论重新获得了生命力，湍流模型的研究得到迅速发展。建立的一系列的两方程模型和二阶矩模型，已经能十分成功地模拟边界层和剪切层流动，但是对于复杂的工业流动，比如大曲率绕流，旋转流动，透平叶栅动静叶互相干扰等，这些因素对湍流的影响还不清楚，这些复杂流