T型管三通流体动力学分析题目T 型管三通流体动力学分析
小组成员:
学院、专业班:
时间： _____ 指导教师:
目录
摘要 (2)
关键字 (2)
前言 (2)
正文 (2)
一、建立模型 (3)
1、.................................................. 绘制模型立体图
3
2、........................................................ 划分网格
4
3、........................................................ 输入参数
4
4、............................................................ 计算
6
二、分析 (7)
1、............................................................ 压强
7
2、............................................................ 速度
8
3、............................................................ 温度
10
三、总结 (11)
摘要
为了加深对工程流体力学基本概念和基本理论的理解，本组依照指导进行了此次实验。为明确冷水、热水在管内如何混合，混合后的运动状态和特性，我们选取了三通管为模型。将三通管注入冷水和热水，汇合后通入大气中。通过假设将其简化后研究其内部流体运动状态的变化。结合压强、管道的属性、水的速度温度及能量损失等问题，应用软件模型计算得到了三通管道内部的流场分布，并对三通管内道内流动的特性进行分析，得出了三通管道紊流流动的计算结果。
关键字
三通管、混合过程、运动状态、能量损失
、八、亠
刖言
工程流体力学是研究流体受力及其运动规律的一门学科，侧重于应用流体力学的基本原理、理论与方法研究解决实际问题。它以流体为研究对象，是研究流体平衡和运动规律的科学。流体力学在水利、航空、电力、机械、冶金、化学、石油、土木等工业技术中有广泛的应用。对口于本专业的机械工业中的润滑、冷却、液压传动、气力输送以及液压和气动控制问题的解决，都必须应用流体力学的理论。因此它是我们理解掌握现代化工程勘测、设计、运行与管理的知识基础，也是我们继续深造及将来从事研究工作的重要工具。
为深入学习流体力学，培养建模能力和分析实例的能力，培养理论联系实际、实事求是、严格认真的科学态度，本组成员积极配合开展了此次实验。
正文
本组研究的流体类型为水，研究围绕三通管内的冷热水混合进行。为节省实验研究的时间和经费，我们采用数值计算方法来研究该问题。通过软件模拟对其进行定性分析。主要研究黏性流体在等速有温差的
条件下产生的局部损失和沿程损失及其动量变化。以期在实验中更加具象的了解连续性方程、伯努利方程、动量方程和达西-巴赫公式的内涵和应用。为展示研究过程及结果，特在此进行系统陈列。
在进行软件模拟之前首先对流体进行以下假设：
1、质量守恒
2、动量守恒
3、连续介质假设
4、不可压缩流体
一、建立模型
模型简介：本装置模型为"T”型三通管，X方向直管管径50mm，Z方向管径30mm。
流动过程简述：向直径为50mm直管的一端管口通入80 C流速为5m/s的热水，直径30mm管口通入10C流速为2m/s的冷水，水流在管中混合后由50mm直管的另一端管口流入大气?（本次实验不对流出部分液体做研究）。
1、绘制模型立体图
图1.1用Geometry绘制的立体图
2、划分网格
Detaik of ^Inflation* - Inflation
图1.2 in flation参数设置
图1.3生成网格
3、输入参数
对面进行命名，两个输入口为inlet_hot （热水管），inlet_cool （冷水管）。出口为outlet。参数如下: 冷水入口端：直径30mm 温度10C 流速2m/s
热水入口端：直径50mm 温度80C 流速5m/s
出口端：直径50mm 压强1atm
湍流强度：l=0.16X Re A（-1/8）
图1.4 冷水入口设置
图1.5热水入口设置
图1.7计算过程
Outiine Boundary: outlet
□
Sources
Option
图1.6 出口设置
4、计算
数学方法：K-&湍流模型与能量方程 （迭代次数：最小值 1最大值100） Details of outlet n Default Domain in Flow Analysis 1
Boundary Details Basic Setting a Plot Options □
Pressure Averag ing Average Over Whole Outlet
大多数飞行器都是在高Re数下飞行，表面的流态是湍流。为了准确地确定湍流流态下的摩阻、热流，湍流成为一个重要而困难的研究课题。 （一）DNS 目前处理湍流数值计算问题有三种方法，第一种方法即所谓直接数值模拟方法（DNS方法），直接求解湍流运动的N-S方程，得到湍流的瞬时流场，即各种尺度的随机运动，可以获得湍流的全部信息。随着现代计算机的发展和先进的数值方法的研究，DNS方法已经成为解决湍流的一种实际的方法。但由于计算机条件的约束，目前只能限于一些低Re数的简单流动，不能用于工程应用。目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数(Re～200)和非常简单的流动外形，如平板边界层、完全发展的槽道流，以及后台阶流动等。用直接数值模拟方法处理工程中的复杂流动问题，即使是当前最先进的计算机也还差三个量级。 （二）LES 另一种方法称做大涡模拟方法（LES方法）。这是一种折衷的方法，即对湍流脉动部分直接地模拟，将N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。小涡对大涡的影响会出现在大涡方程中，再通过建立模型（亚格子尺度模型）来模拟小涡的影响。由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件，难以找出普遍的湍流模型来描述具有不同的边界特征的大涡结构，宜做直接模拟。相反地，小尺度涡对边界条件不存在直接依赖关系，而且一般具有各向同性性质。所以亚格子模型具有更大的普适性，比较容易构造，这是它比雷诺平均方法要优越的地方。自从1970年Deardorff第一次给出具有工程意义的LES计算以来，LES方法已经成为计算湍流的最强有力的工具之一，应用的方向也在逐步扩展，但是仍然受计算机条件等的限制，使之成为解决大量工程问题的成熟方法仍有很长的路要走。 （三）RANS 目前能够用于工程计算的方法就是模式理论。所谓湍流模式理论，就是依据湍流的理论知识、实验数据或直接数值模拟结果，对Reynolds应力做出各种假设，即假设各种经验的和半经验的本构关系，从而使湍流的平均Reynolds方程封闭。随着计算流体力学的发展，湍流模式理论也有了很大的进步，有了非常丰硕的成果。从对模式处理的出发点不同，可以将湍流模式理论分类成两大类：一类称为二阶矩封闭模式，另一类称涡粘性封闭模式。 （1）雷诺应力模式 所谓二阶矩封闭模式，是从Reynolds应力满足的方程出发，将方程右端未知的项（生成项，扩散项，耗散项等）用平均流动的物理量和湍流的特征尺度表示出来。典型的平均流动的变量是平均速度和平均温度的空间导数。这种模式理论，由于保留了Reynolds应力所满足的方程，如果模拟的好，可以较好地反映Reynolds应力随空间和时间的变化规律，因而可以较好地反映湍流运动规律。因此，二阶矩模式是一种较高级的模式，但是，由于保留了Reynolds应力的方程，加上平均运动的方程整个方程组总计15个方程，是一个庞大的方程组，应用这样一个庞大的方程组来解决实际工程问题，计算量很大，这就极大地限制了二阶矩模式在工程问题中的应用。 （2）涡粘性模式
1、数值的耗散与频散： 在数值解中出现的振幅衰减波长加宽的现象叫数值耗散，与高阶偶次空间偏导数有关；在数值解中出现解得主波后有一系列频及传播速度不等的尾波的现象叫数值频散，与高阶奇次偏导数有关。 2、湍流模型理论：湍流模式理论或简称湍流模型，就是以雷诺平均运动方程与脉动运动方程为基础，依靠理论与经验的结合，引进一系列模型假设，而建立起得一组 描写湍流平均量的封闭方程组。 3、修正的偏微分方程：与差分方程相等价的微分方程称之为修正的微分方程。 4、自适应网格:为了计算具有高雷诺数的流场，必须将流场内的网格加密，但是实际计算中并不需要对全流场的网格所有部分同样加密，只需在某些部分，如物面附近、尾流区等得网格加密即可。因此需要事先估计一些变化较快的区域，但这种估计又是是正确的。有时则不正确。特别是不定常流动，流动过程本身就是变化的，所以需要不断的调整网格的位置和疏密，这样就产生了自适应网格。 5、CFL 条件：定义t C x μ ?=? ，不等式1C ≤ 称为CFL 条件，此条件一般应用于双曲线偏微分方程的显式格式。物理意义：即在时间步长内，波的位移应小于空间步长。 数学意义：差分方程解的依赖区域包含微分方程解得依赖区域。 1、简答CFD 方法求解流动问题的基本步骤 答：①确定流动模型；②计算区域离散化；③用离散节点变量代替场；④将控制方程中偏导数进行离散，得到线性方程组；⑤边界条件和初值条件离散化；⑥离散的线性方程组求解，得到离散值；⑦计算结果数据处理。 2、简述离散偏微分方程的三个原则及LAX 定理 三原则；相容性、稳定性、收敛性。 LAX 定理：对于一个选定的线性偏微分方程的初值问题，对应的差分方法是相容的，则差分方程解得收敛性和稳定性事等价的或者说稳定性是收敛性的充要条件。 3、简述差分格构造的基本规律，并应用规律方程 0t x μμ λ??+=?? 利用网格点（） （）（）构造方程的差分格式，并验证其离散格式的精度等级。 答：构造的基本规律 ：①为保证均匀流场，差分的分子各项系数之和为零 ②分母向量级与微分的阶数一致 ③构造差分级指明针对哪点构造 ④差分格式的精度 由网格点（）（）（）规律方程（）构造得 1 11 1 0n n n n j j j j x x μμμμλ +++---+=?? 令112j j x k k x μμμ-+= ? 用泰勒公 式展开的23 126j j x xx x x x x μμμμμ-??=-?+- 所以12101k k k +=?? -=? 得12 11k k =-??=? 所以1j j x x μμμ--+=? 所以具有一阶精度 4、简要概括流动的数值计算对网格的基本要求 答：①计算域边界上的网格节点都应在边界上 ②物理域上的特点与计算域上的节点要求一一对应 ③网格应尽量尺寸匀称，相邻网格长度比应小于2 ④物理域网格夹角不宜太小（≥45°） ⑤流动参数梯度大的地方网格要加密，否则稀疏。 5、简述人工压缩方法（时间相关法）的基本思想 答：用非定常流动方程来求解定常流动问题，用其稳态求解定常流动的解，将不可压缩的粘性流动的连续方程，添加到可压缩项。则与动量方程构成定常粘性流动时间相关方程，可把非定常流动的稳态解作为非定常流动的解。
湍流模型发展综述 摘要：在概述了湍流问题的基础上，本文简要介绍了湍流的四种模型，对湍流模型在不同情况下的模拟能力进行了对比，最后简述了湍流模型的发展方向。 关键词：湍流模型；Navier-Stokes方程组；J-K模型 Abstract：On the basis of introducing the problems of turbulence, this paper briefly analyzed four kinds of turbulence models and compared their ability of simulation in different situations. At last, the paper expounded the development direction of the turbulence model. Key words：Turbulence model; Navier-Stokes equations; J-K model 一、引言 湍流又称紊流,是自然界中常见的一种很不规则的流动现象。当粘性阻尼无法消除惯性的影响时,自然界中的绝大部分流动都是湍流。 湍流运动的实验研究表明，虽然湍流结构十分复杂，但它仍然遵循连续介质的一般动力学规律，湍流流动的各物理量的瞬时值也应该服从一般的N-S方程。对粘性流体服从的N-S方程进行时均化，就可以得到雷诺平均方程。与定常的N-S方程相比，不同之处是在该式右边多了九项与脉动量有关的项，这脉动量的乘积的平均值与密度的乘积是湍流流动中的一种应力，称为湍流应力或雷诺应力。其中，法向雷诺应力和切向雷诺应力各有三个。 湍流问题就是在给定的边界条件下解雷诺方程。由于雷诺平均方程中未知数个数远多于方程个数而出现了方程不封闭的问题，这就需要依据各种半经验理论提出相应的补充方程式，即各种湍流模型。一般按照所用湍流量偏微分方程的物理含义或者数量进行区分，分别称为梅罗尔—赫林方法和雷诺方法。而后者又将湍流模型分成四类。（1）零方程模型；（2）一方程模型；（3）二方程模型；（4）应力方程模型。下面就对这些模型进行简单的描述。 二、湍流模型简介 1、零方程模型 最初的湍流模型只考虑了一阶湍流计算统计量的动力学微分方程，即平均方程，没有引进高阶统计量的微分方程，因而称之为一阶封闭模式或零方程模型。零方程模型又称为代数模型，代数模型又可以分成以下几种模型：（1）Cebeci —Smith 模型，（2）Baldwin—Lomax 模型，（3）Johnson—King 模型。 其中，B-L与C-S模型的不同之处在于外层湍流粘性系数取法不同。后者适用于湍流边界层，而前者则可用于 N-S方程的计算。此两模型已在工程计算中
汽车外部气体流动模拟 振动和噪声控制研究所 1.模型概述 在汽车外部建立一个较大的长方体几何空间，长度约为30m，宽度和高度约为5m，在空间内部挖出汽车形状的空腔，汽车尺寸参照本田CRV为4550mm*1820mm*1685mm。由于汽车向前开进，气体从车头流向车尾，因此将汽车前方空间设为气体入口，后方空间设为气体出口，模拟气体在车外的流动。另外为了节省计算成本将整个模型按1:100的比例缩小，考虑到模型和流体均是对称的，因此仅画出几何模型的一半区域，建立对称面以考虑生成包含理想气体的流体域。在Catia中建立的模型如图1.1所示。 图1.1几何模型 2.利用ICEM CFD进行网格划分 a)导入有Catia生成的stp格式的模型； b)模型修复，删除多余的点、线、面，允许公差设为0.1； c)生成体，由于本模型仅为流体区域，因此将全部区域划分为一个体，选取方法可以 使用整体模型选取； d)为了后面的设置边界方便，因此将具有相同特性的面设为一个part,共设置了in， out，FreeWalls，Symmetry和Body； e)网格划分，设置Max element=2，共划分了1333817个单元，有225390个节点； f)网格输出，设置求解器为ANSYS CFX，输出cfx5文件。 3.利用ANSYS CFX求解 a)生成域，物质选定Air Ideal Gas，参考压强设为1atm，浮力选项为无浮力模型，
域运动选项为静止，网格变形为无；流体模型设定中的热量传输设定为Isothermal，流体温度设定为288k，湍流模型设定为Shear Stress Transport模型，壁面函数 选择Automatic。 b)入口边界设定，类型为Inlet，位置选定在in，质量与栋梁选定Normal Speed，设 定为15m/s，湍流模型设定类型为Intensity and Length Scale=0.05，Eddy Len.Scale=0.1m。 c)出口边界设定，边界类型为Outlet，位置选out。质量与动量选项为Static Pressure，相对压强为0pa。 d)壁面边界设定，边界类型为Wall，位置选在FreeWalls。壁面边界详细信息中指定 WallInfluence On Flow为Free Slip。 e)对称边界设定，边界类型为Symmetry，位置选在Symmetry。 f)汽车外壁面设定，边界类型为Wall，位置设在Body，壁面详细信息选项中指定Wall Influence On Flow为No Slip，即汽车壁面为无滑移壁面。 g)初始条件设定，初始速度分量设为U方向为15m/s，其他两个方向的速度为零。 h)求解设置，残差类型选为RMS，残差目标设定为1e-5，当求解达到此目标时，求解 自动终止。求解之前的模型如图3.1所示。 图3.1求解之前的模型 4.结果后处理 从图4.1中可以看出计算收敛。
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。 从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解析解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
我所知道的计算流体力学(CFD)大牛们 (1) Jameson的故事 Jameson是当今CFD届的超级大牛。偶的超级偶像哦。Jameson是个英国人，出生在军人世家。从小随老爹驻守印度。于是长大了也抗起枪到海外保卫日不落帝国，军衔是Second Lieutenant。无奈“日不落”已落，皇家陆军已经不需要他了。大概有什么立功表现把，退役后就直接进了剑桥大学。在那里拿到博士学位。辗转间从英国来到了美国，从工厂又到了学校。成了Princeton的教授。在那里提出了著名的中心差分格式和有限体积法。就是在这里，发表了他那篇著名的中心差分离散的有限体积法。中心差分格式，大家都知道，是二阶，但是稳定范围特别小，Pe不能超过2，于是就得加人工粘性（一听这名字，数学家就倔嘴巴，不科学嘛），这是大学生都知道的事，怎么加就是学问了。Jameson用二阶项做背景粘性，用四阶项抑制激波振荡（也亏他想得出来），配合他提出的有限体积法，获得了极大的成功，很快风靡世界，工程界几乎无一例外在使用他的方法，原因很简单，他的方法乐百氏，而且又有相当精度。从此大行于市，座上了P大的航空系系主任，也确立了CFD界第一大牛人的地位。 Jameson发文章有个特点，喜欢发在小会议上或者烂杂志上，反正是SCI检索不到地方。包括后来关于非结构网格，多重网格等等经典的开创性文章，都是这样。（如果按照清华的唯SCI论的评判标准，我估计在清华最多只能给他评一个副教授当当。）牛牛的人总是遭人忌妒，哪里都这样。看着Jameson的有限体积方法这么受欢迎，有些人就红眼了。于是说，有限体积方法不错，可惜只适合于定常问题计算，非定常计算就不怎么样嘛。Jameson那里能容忍别人对他的得意之做胡说。于是，灵机一动，想出了一个双时间尺度的方法，引进一个非物理时间，把非定常问题变成了一个定常问题计算，还真好使，又风靡世界，从此天下太平。97年，Jameson年龄到了，就从P大退休了，结果又被聘请到Standford大学当Thomas V. Jones Professor搞起了湍流来。前不久偶导师见他回来，对欧们边摇头边说，“几年不见，老得快不行了”，言下之意，我们如果想多活几年，不要去搞什么湍流。 (2) Steven A. Orszag的故事 Steven A. Orszag是一个天才级别的人物啦。在直接数值模拟，谱方法，湍流模型等等许多方面都有开创性的贡献。天才嘛，总是有缺陷的，不是生活不能自理，就是不懂得处理人际关系。前者还好办，只是lp不舒服，后者嘛，让同事和同行不舒服，可麻烦就大了。 不幸的是，Orszag属于后者。对于他的恃才傲物，有人早就恨得牙根痒痒，报复的机会终于来了。 三十年前，湍流模型的先驱们，是通过数值试验，再连懵带猜的确定下了双方程湍流模型的参数。20年前，Orszag突发奇想，能否用RNG（重整化群理论）从理论上推导这些参数呢？RNG理论在相变上取得了很大的成功，发明者也在81年获得了Nobel奖。牛人就是牛人很快居然真从理论上推出了这些参数。这下湍流模型界可炸开了锅，这岂不是要砸掉很多人的饭碗？这不等于说那些老家伙几十年前的工作一钱不值么？这帮大学霸可不是省油的灯。环顾地球之大，Orszag居然找不到一本杂志愿意接受他这篇文章。Orszag这个郁闷呀，这个生
湍流模型 目前计算流体力学常用的湍流的数值模拟方法主要有以下三种： ?直接模拟（direct numerical simulation, DNS） 直接数值模拟（DNS）特点在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对Navier-Stokes方程直接求解。这种方法能对湍流流动中最小尺度涡进行求解，要对高度复杂的湍流运动进行直接的数值计算，必须采用很小的时间与空间步长，才能分辨出湍流中详细的空间结构及变化剧烈的时间特性。基于这个原因，DNS目前仅限于相对低的雷诺数中湍流流动模型。另外，利用DNS模型对湍流运动进行直接的数值模拟对计算工具有很高的要求，计算机的内存及计算速度要非常的高，目前DNS模型还无法应用于工程数值计算，还不能解决工程实际问题。 ?大涡模拟(large eddy simulation, LES) 大涡模拟（LES）是基于网格尺度封闭模型及对大尺度涡进行直接求解N-S方程，其网格尺度比湍流尺度大，可以模拟湍流发展过程的一些细节，但其计算量仍很大，也仅用于比较简单的剪切流运动及管流。大涡模拟的基础是：湍流的脉动与混合主要是由大尺度的涡造成的，大尺度涡是高度的非各向同性，而且随流动的情形而异。大尺度的涡通过相互作用把能量传递给小尺度的涡，而小尺度的涡旋主要起到耗散能量的作用，几乎是各向同性的。这些对涡旋的认识基础就导致了大涡模拟方法的产生。Les大涡模拟采用非稳态的N-S方程直接模拟大尺度涡，但不计算小尺度涡，小涡对大涡的影响通过近似的模拟来考虑，这种影响称为亚格子Reynolds应力模型。大多数亚格子Reynolds模型都是将湍流脉动所造成的影响用一个湍流粘性系数，既粘涡性来描述。LES对计算机的容量和CPU的要求虽然仍然很高，但是远远低于DNS方法对计算机的要求，因而近年来的研究与应用日趋广泛。 ?应用Reynolds时均方程(Reynolds-averaging equations)的模拟方法 许多流体力学的研究和数值模拟的结果表明，可用于工程上现实可行的湍流模拟方法仍然是基于求解Reynolds时均方程及关联量输运方程的湍流模拟方法，即湍流的统观模拟方法。统观模拟方法的基本思想是用低阶关联量和平均流性质来模拟未知的高阶关联项，从而封闭平均方程组或关联项方程组。虽然这种方法在湍流理论中是最简单的，但是对工程应用而言仍然是相当复杂的。即便如此，在处理工程上的问题时，统观模拟方法仍然是最有效、最经济而且合理的方法。在统观模型中，使用时间最长，积累经验最丰富的是混合长度模型和K-E模型。其中混合长度模型是最早期和最简单的湍流模型。该模型是建立在层流粘性和湍流粘性的类比、平均运动与湍流的脉动的概念上的。该模型的优点是简单直观、无须增加微分方程。缺点是在模型中忽略了湍流的对流与扩散，对于复杂湍流流动混合长度难以确定。到目前为止，工程中应用最广泛的是k-ε模型。另外针对k-ε模型的不足之处，许多学者通过对K-E模型的修正和发展，开始采用雷诺应力模型（DSM）和代数应力模型（ASM）。近年来，DSM模型已用来预报燃烧室及炉内的强旋及浮力流动。很多情况下能够给出优于k-ε模型的结果。但是该模型也有不足之处，首先它对工程预报来说太复杂，其次经验系数太多难以确定，此外，对压力应变项的模拟还有争议。更主要的是，尽管这一模型考虑了各种应变效应，但是其总精度并不总是高于其它模型，这些缺点导致了DSM模型没有得到广泛的应用。总之，虽然从本质上讲DSM模型和ASM模型比k-ε模型对湍流流场的模拟更加合理，但DSM和ASM中仍然采用精度不高的E方程，模型中常数的通用性还没有得到广泛的验证，边界条件不好给定，计算也比较复杂。正因为如此，目前用计算解决湍流问题时仍然采用比较成熟的K-E模型。 需要注意的是： 1、大涡模拟有自己的亚格子封闭模型，这和k-ε模型完全是两回事。LES的亚格子模型表现的是过滤掉的小涡对大涡的影响（这种影响是相互的）。而Reynolds时均方程的k-ε是建
湍流模型理论 §３．1 引言 自然界中的实际流动绝大部分是三维的湍流流动，如河流,血液流动等。湍流是流体粘性运动最复杂的形式,湍流流动的核心特征是其在物理上近乎于无穷多的尺度和数学上强烈的非线性，这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。回顾计算流体力学的发展，特别是活跃的８０年代,不仅提出和发展了一大批高精度、高分辨率的计算格式,从主控方程看相当成功地解决了Euler方程的数值模拟，可以说Euler方程数值模拟方法的精度已接近于它有效使用范围的极限；同时还发展了一大批有效的网格生成技术及相应的软件，具体实现了工程计算所需要的复杂外形的计算网格；且随着计算机的发展,无论从计算时间还是从计算费用考虑,Ｅｕlｅｒ方程都已能适用于各种实践所需。在此基础上,8０年代还进行了求解可压缩雷诺平均方程及其三维定态粘流流动的模拟。９0年代又开始一个非定常粘流流场模拟的新局面，这里所说的粘流流场具有高雷诺数、非定常、不稳定、剧烈分离流动的特点，显然需要继续探求更高精度的计算方法和更实用可靠的网格生成技术。但更为重要的关键性的决策将是，研究湍流机理,建立相应的模式，并进行适当的模拟仍是解决湍流问题的重要途径。 要反映湍流流场的真实情况，目前数值模拟主要有三种方法：１．平均Ｎ－S方程的求解，２.大涡模拟（LＥS），3.直接数值模拟(DNS）。但是由于叶轮机械内部结构的复杂性以及目前计算机运算速度较慢，大涡模拟和直接数值模拟还很少用于叶轮机械内部湍流场的计算,更多的是通过求解平均N-S方程来进行数值模拟。因为平均N-S方程的不封闭性,人们引入了湍流模型来封闭方程组，所以模拟结果的好坏很大程度上取决于湍流模型的准确度。自7０年代以来，湍流模型的研究发展迅速，建立了一系列的零方程、一方程、两方程模型和二阶矩模型,已经能够十分成功的模拟边界层和剪切层流动。但是,对于复杂的工业流动,比如航空发动机中的压气机动静叶相互干扰问题,大曲率绕流,激波与边界层相互干扰,流动分离，高速旋转以及其他一些原因,常常会改变湍流的结构,使那些能够预测简单流动的湍流模型失效,所以完善现有湍流模型和寻找新的湍流模型在实际工作中显得尤为重要。 §3.2 湍流模型概述 §3.２.1 湍流模型的引入