2012年秋季学期研究生课程考核（读书报告、研究报告）考核科目高等流体力学学生所在院（系）机电工程学院学生所在学科机械制造及自动化学生姓名高强学号12S008123学生类别工学硕士考核结果阅卷人湍流的数值模拟一、流体力学概述流体力学是研究流体的力学运动规律及其应用的学科。

主要研究在各种力的作用下，流体本身的状态，以及流体和固体壁面、流体和流体间、流体与其他运动形态之间的相互作用的力学分支。

除水和空气之外，这里的流体还指作为汽轮机工作介质的水蒸气、润滑油、地下石油、含泥沙的江水、血液、超高压作用下的金属和燃烧后产生成分复杂的气体、高温条件下的等离子体等等。

它的主要基础是牛顿运动定律和质量守恒定律，常常还要用到热力学知识，有时还用到宏观电动力学的基本定律、本构方程和高等数学、物理学、化学的基础知识。

气象、水利的研究，船舶、飞行器、叶轮机械和核电站的设计及其运行，可燃气体或炸药的爆炸，汽车制造，以及天体物理的若干问题等等，都广泛地用到流体力学知识。

许多现代科学技术所关心的问题既受流体力学的指导，同时也促进了它不断地发展。

二、数值计算在流体力学研究中的应用数值计算是研究流体力学的重要方法。

它是针对流体运动的特点，用数学语言将质量守恒、动量守恒、能量守恒等定律表达出来，从而得到连续性方程、动量方程和能量方程。

此外，还要加上某些联系流动参量的关系式(例如状态方程)，或者其他方程。

这些方程合在一起称为流体力学基本方程组。

求出方程组的解后，结合具体流动，解释这些解的物理含义和流动机理。

通常还要将这些理论结果同实验结果进行比较，以确定所得解的准确程度和力学模型的适用范围。

从基本概念到基本方程的一系列定量研究，都涉及到很深的数学问题，所以流体力学的发展是以数学的发展为前提。

反过来，那些经过了实验和工程实践考验过的流体力学理论，又检验和丰富了数学理论，它所提出的一些未解决的难题，也是进行数学研究、发展数学理论的好课题。

按目前数学发展的水平看，有不少题目将是在今后几十年以内难于从纯数学角度完善解决的。

在流体力学理论中，用简化流体物理性质的方法建立特定的流体的理论模型，用减少自变量和减少未知函数等方法来简化数学问题，在一定的范围是成功的，并解决了许多实际问题。

对于一个特定领域，考虑具体的物理性质和运动的具体环境后，抓住主要因素忽略次要因素进行抽象化也同时是简化，建立特定的力学理论模型，便可以克服数学上的困难，进一步深入地研究流体的平衡和运动性质。

20世纪50年代开始，在设计携带人造卫星上天的火箭发动机时，配合实验所做的理论研究，正是依靠一维定常流的引入和简化，才能及时得到指导设计的流体力学结论。

此外。

流体力学中还经常用各种小扰动的简化，使微分方程和边界条件从非线性的变成线性的。

声学是流体力学中采用小扰动方法而取得重大成就的最早学科。

声学中的所谓小扰动，就是指声音在流体中传播时，流体的状态(压力、密度、流体质点速度)同声音未传到时的差别很小。

线性化水波理论、薄机翼理论等虽然由于简化而有些粗略，但都是比较好地采用了小扰动方法的例子。

每种合理的简化都有其力学成果，但也总有其局限性。

例如，忽略了密度的变化就不能讨论声音的传播；忽略了粘性就不能讨论与它有关的阻力和某些其他效应。

掌握合理的简化方法，正确解释简化后得出的规律或结论，全面并充分认识简化模型的适用范围，正确估计它带来的同实际的偏离，正是流体力学理论工作和实验工作的精华。

流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。

20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。

数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了“计算流体力学”。

从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。

数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。

数值计算方法最近发展很快，其重要性与日俱增。

随着计算机CPU 速度的提高, 硬盘及内存容量的持续增大, 以及数值计算方法的不断完善, 计算流体力学( CFD) 已经成为民用飞机气动设计过程中的一个强有力的工具。

然而湍流模型仍然是阻碍人们应用N-S 方程组进行飞机设计的瓶颈之一。

近20 年来, 人们对粘性流动进行了比较深入的研究, 提出了许多湍流模型, 具有代表性且应用较广的模型是: 代数B-L 模型, 半方程J-K 模型, 一方程Spalar-t Allmaras模型, 以及两方程k-E和k-X 模型等。

三、湍流湍流是一种高度复杂的三维非稳态、带旋转的不规则流动。

湍流中流体的各个物理参数，如速度、压力、温度等都随时间与空间发生随机变化。

从物理机理上说，可以把湍流看成是各种不同尺度的涡旋叠合而成的，这些涡的大小及旋转方向分布是随机的。

大尺度的涡旋主要由流动的边界条件所决定，其尺寸可与流场的大小相比拟，是引起低频脉动的原因；小尺度的涡主要是由粘性力所决定的，其尺寸可能只有流场尺度的千分之一，是引发高频脉动的原因。

大尺度的涡破裂后形成小尺度的涡，较小尺度的涡破裂后形成更小的涡。

大尺度的涡从主流获得能量，通过涡间的转化将能量传给小尺度的涡。

最后由于粘性的作用，小尺度的涡不断消失，机械能就转化(即耗散)为流体的热能。

同是，由于边界的作用，扰动及速度梯度的影响，新的涡又不断产生，这就构成了湍流运动。

可见，湍流的一个重要特点是物理量的脉动，非稳态的N-S方程对湍流运动仍是适用的。

湍流是自然界当中普遍存在的一种非常复杂的流动现象, 但人们对湍流机理的认识及其数值模拟方法至今仍处于探索阶段。

包括已故诺贝尔奖获得者Feynman 在内的好几位物理学家认为, 湍流是经典物理学中尚未得到解决的最后一个大难题。

人们用现代的理论和方法系统地研究湍流现象始于19 世纪末,Reynolds 提出统计平均方法是湍流研究的起点。

一个多世纪以来, 尽管在湍流本质认识和实际应用方面, 湍流研究都取得了很大的进步, 但是随着计算流体力学及计算空气动力学方法的不断完善, 计算机性能的不断提高, 湍流的数值模拟方法已成为阻碍人们应用N- S 方程进行水流运动特性分析、管道螺旋流水力输送研究[2] 、飞机设计等的瓶颈之一。

对湍流基础研究的进展, 可以直接促进许多实际工程及科学应用的进步。

目前, 湍流数值模拟的方法有: 直接数值模拟( Direct Numerical Simulation, DNS) 、大涡数值模拟( Large Eddy Simulation, LES)和雷诺平均模拟( Reynolds Averaged Navier-Stokes, RANS)。

1、直接模拟(DNS)直接数值模拟（DNS）特点在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对Navier-Stokes方程直接求解。

这种方法能对湍流流动中最小尺度涡进行求解，要对高度复杂的湍流运动进行直接的数值计算，必须采用很小的时间与空间步长，才能分辨出湍流中详细的空间结构及变化剧烈的时间特性。

基于这个原因，DNS目前仅限于相对低的雷诺数中湍流流动模型。

另外，利用DNS模型对湍流运动进行直接的数值模拟对计算工具有很高的要求，计算机的内存及计算速度要非常的高，目前DNS模型还无法应用于工程数值计算，还不能解决工程实际问题。

2、大涡模拟(LES)大涡模拟（LES）是基于网格尺度封闭模型及对大尺度涡进行直接求解N-S方程，其网格尺度比湍流尺度大，可以模拟湍流发展过程的一些细节，但其计算量仍很大，也仅用于比较简单的剪切流运动及管流。

大涡模拟的基础是：湍流的脉动与混合主要是由大尺度的涡造成的，大尺度涡是高度的非各向同性，而且随流动的情形而异。

大尺度的涡通过相互作用把能量传递给小尺度的涡，而小尺度的涡旋主要起到耗散能量的作用，几乎是各向同性的。

这些对涡旋的认识基础就导致了大涡模拟方法的产生。

Les大涡模拟采用非稳态的N-S方程直接模拟大尺度涡，但不计算小尺度涡，小涡对大涡的影响通过近似的模拟来考虑，这种影响称为亚格子Reynolds应力模型。

大多数亚格子Reynolds模型都是将湍流脉动所造成的影响用一个湍流粘性系数，既粘涡性来描述。

LES对计算机的容量和CPU的要求虽然仍然很高，但是远远低于DNS方法对计算机的要求，因而近年来的研究与应用日趋广泛。

3、Reynolds平均法许多流体力学的研究和数值模拟的结果表明，可用于工程上现实可行的湍流模拟方法仍然是基于求解Reynolds时均方程及关联量输运方程的湍流模拟方法，即湍流的统观模拟方法。

统观模拟方法的基本思想是用低阶关联量和平均流性质来模拟未知的高阶关联项，从而封闭平均方程组或关联项方程组。

虽然这种方法在湍流理论中是最简单的，但是对工程应用而言仍然是相当复杂的。

即便如此，在处理工程上的问题时，统观模拟方法仍然是最有效、最经济而且合理的方法。

在统观模型中，使用时间最长，积累经验最丰富的是混合长度模型和 K-E模型。

其中混合长度模型是最早期和最简单的湍流模型。

该模型是建立在层流粘性和湍流粘性的类比、平均运动与湍流的脉动的概念上的。

该模型的优点是简单直观、无须增加微分方程。

缺点是在模型中忽略了湍流的对流与扩散，对于复杂湍流流动混合长度难以确定。

到目前为止，工程中应用最广泛的是K-E模型。

另外针对K-E模型的不足之处，许多学者通过对K-E模型的修正和发展，开始采用雷诺应力模型（DSM）和代数应力模型（ASM）。

近年来，DSM模型已用来预报燃烧室及炉内的强旋及浮力流动。

很多情况下能够给出优于K-E模型的结果。

但是该模型也有不足之处，首先它对工程预报来说太复杂，其次经验系数太多难以确定，此外，对压力应变项的模拟还有争议。

更主要的是，尽管这一模型考虑了各种应变效应，但是其总精度并不总是高于其它模型，这些缺点导致了DSM模型没有得到广泛的应用。

总之，虽然从本质上讲DSM模型和ASM模型比K-E模型对湍流流场的模拟更加合理，但DSM和ASM中仍然采用精度不高的E方程，模型中常数的通用性还没有得到广泛的验证，边界条件不好给定，计算也比较复杂。

正因为如此，目前用计算解决湍流问题时仍然采用比较成熟的K-E模型。

Reynolds平均法的分类在这类方程中，将非稳态N-S方程对时间作平均，即把湍流运动看成二个流动的叠加，一是时间平均流动，二是瞬时脉动流动。