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