K-e湍流模型

K是紊流脉动动能（J），ε 是紊流脉动动能的耗散率（%）

K越大表明湍流脉动长度和时间尺度越大，ε 越大意味着湍流脉动长度和时间尺度越小，它们是两个量制约着湍流脉动。

但是由于湍流脉动的尺度范围很大，计算的实际问题可能并不会如上所说的那样存在一个确切的正比和反比的关系。在多尺度湍流模式中，湍流由各种尺度的涡动结构组成，大涡携带并传递能量，小涡则将能量耗散为内能。

在入口界面上设置的K和湍动能尺度对计算的结果影响大，

至于k是怎么设定see fluent manual "turbulence modelling"

作一个简单的平板间充分发展的湍流流动，

基于k-e模型。

确定压力梯度有两种方案，一是给定压力梯度，二是对速度采用周期边界条件，压力不管！

k-epsiloin湍流模型参数设置：k－动能能量；epsilon－耗散率；

在运用两方程湍流模型时这个k值是怎么设置的呢epsilon可以这样计算吗Mepsilon＝Cu*k*k/Vt%

这些在软件里有详细介绍。陶的书中有类似的处理，假定了进口的湍流雷诺数。

fluent帮助里说，用给出的公式计算就行。

k－e模型的收敛问题！

应用k－e模型计算圆筒内湍流流动时，网格比较粗的时计算结果能收敛，但是当网格比较密的时候，湍流好散率就只能收敛到10的－2次方，请问大侠有没有解决的办法

用粗网格的结果做初场网格加密不是根本原因，更本的原因是在加密过程中，部分网格质量差注意改进网格质量，应该就会好转.

在求解标准k-e双方程湍流模型时（采用涡粘假设，求湍流粘性系数，然后和N－S方程耦

合求解粘性流场），发现湍动能产生项（雷诺应力和一个速度张量相乘组成的项）出现负

值，请问是不是一种错误现象

如果是错误现象一般怎样避免。另外处理湍动能产生项采

用什么样的差分格式最好。而且因为源项的影响，使得程序总是不稳定，造成k,e值出现负

值，请问有什么办法克服这种现象。

你可以试试这里计算的时候加一个判断，出现负值的时候强制为一个很小的正值。

这可能是因为你采用的数值格式的问题，一般计算程序对k方程都要做一定处理，

以保证k的正定。

比如，强制规定源项与0的关系，以使数值计算稳定。

就ke模型而言。

它是problem dependent.对简单的无弯曲无旋转无...的湍流问题，它能算而且能给出好的结果，但对复杂的流动问题，它就不能使用了。

出现负的ke不仅仅是计算格式的问题，

更重要的是模型问题，没有谁能证明ke模型在任何流动问题中都能保证ke 是正的。

有这么一些办法避免ke出现负值

1。对K＝ln(k)和E=ln(e)求解，问题：壁面ke=0难处理，

2。先用层流计算500步，然后再用ke算

3。各种强制限制办法

4。源项局部线性化

5。算到一定程度，如果k值趋势对了，就干脆不求ke方程

几种湍流模型

解决湍流的模型总计就是那几个方程，Flue nt又从工程和数值的角度进行了整理，下面就是这些湍流模型的详细说明。 FLUENT提供了以下湍流模型： ?Spalart-Allmaras 模型 ?k-e模型 —标准k-e模型 —Ren ormalizatio n-group （RNG^e 模型 —带旋流修正k-e模型 ?k-3模型 —标准k- 3模型 —压力修正k- 3模型雷诺兹压力模型大漩涡模拟模型 几个湍流模型的比较： 从计算的角度看Spalart-Allmaras模型在FLUENT中是最经济的湍流模型，虽然只有一种方程可以解。由于要解额外的方程，标准ke模型比Spalart-Allmaras模型耗费更多的计算机资 源。带旋流修正的k-e模型比标准ke模型稍微多一点。由于控制方程中额外的功能和非线性，RN&七模型比标准k-e模型多消耗10?15%的CPU时间。就像k七模型，k-3模型也是两个方程的模型，所以计算时间相同。 比较一下k◎莫型和k-3模型，RSM模型因为考虑了雷诺压力而需要更多的CPU时间。然而高效的程序大大的节约了CPU时间。RSM模型比k-e模型和k-3模型要多耗费50?60%的CPU 时间，还有15?20%的内存。 除了时间，湍流模型的选择也影响FLUENT勺计算。比如标准k-e模型是专为轻微的扩散 设计的，然而RNGk-e模型是为高张力引起的湍流粘度降低而设计的。这就是RNG莫型的缺点。同样的，RSM模型需要比k-e模型和k-3模型更多的时间因为它要联合雷诺压力和层流。 概念：1?雷诺平均：在雷诺平均中，在瞬态N-S方程中要求的变量已经分解为时均常量和变量。 相似的，像压力和其它的标量 ；（10.2-2） i「 这里??表示一个标量如压力，动能，或粒子浓度。 2. Boussinesq逼近从雷诺压力转化模型：禾U用Bouss in esq假设把雷诺压力和平均速度梯度 联系起来： +茁飞（肚+川亦）也（10 2-O） Boussinesq假设使用在Spalart-Allmaras模型、k-e模型和k- 3模型中。这种逼近方法好处是对计算机的要求不高。在Spalart-Allmaras模型中只有一个额外的方程要解。k-e模型和k-3模型 中又两个方程要解。Bouss inesq假设的不足之处是假设u t是个等方性标量，这是不严格的。

煤粉燃烧模拟——湍流破碎模型.

煤粉燃烧——湍流破碎模型（EBU） 简介 该帮助文件主要介绍煤粉燃烧模型的设置和求解，采用湍流破碎模型（EBU）。 EBU燃烧模型，也称涡团破碎模型，假设化学反应的平均速度与化学动力学无关，而只取决于低温的反应物和高温的燃烧产物之间的湍流混合作用。 主要包括： 1）煤粉燃烧模型的建立和求解 2）湍流破碎模型（EBU）的应用 3）选择合适的求解参数 4）计算结果的后处理 问题描述 3D模型的剖面图如图1所示。左侧为两个环形入口，右侧为一个圆形出口。由于模型的对称性，取系统的1/4进行建模。煤粉和携带空气（一次风）从内环进入燃烧室，二次风从外环进入燃烧室，发生燃烧反应，产物从压力出口流出。模型建立和求解 Step1：网格 1、读取mesh文件 2、检查网格：grid check 3、显示网格：display grid a）从列表框中选取所有surfaces b）点击display，并关闭grid display 面板 Step2：模型 1、选择k-ε湍流模型 Define models viscous 2、启动能量方程 Define models energy 3、启动物质输运方程 Define models species transport & reaction

a）选择model列表中的species transport b）选择reaction列表中的volumetric c）选择mixture material 列表中的coal-hv-volatiles-air d）选择turbulence-chemistry interaction（湍流与化学反应的作用）列表中的eddy dissipation（涡流耗散） e）点击ok，关闭species model面板 4、启动Discrete ordinates模型（DO离散坐标系） a）从model列表中选择discrete ordinates b）设置flow iterations per radiation iteration（流动和辐射迭代次数）为1 c）设置angular discretization （角离散化）中的Theta Divisions and Phi Divisions为4 d）设置angular discretization （角离散化）中的Theta Pixels and Phi Pixels 为3 e）点击ok，关闭radiation model 5、启动discrete phase（离散相）模型 a）M ax number of steps（最大步数）40000 b）启动specify length scale（步长），设为0.0025 c）点击ok，关闭discrete phase面板 Step3：injections（喷口） 1、v-1入口截面设9个喷口：define injections a）点击creat，新建喷口 c）9个喷口通用性质，见表 c）点击turbulent dispersion，并启动discrete random walk model

湍流模型的选择依据

解决湍流的模型总计就就是那几个方程,Fluent 又从工程与数值的角度进行了整理,下面就就是这些湍流模型的详细说明。 FLUENT 提供了以下湍流模型: ·Spalart-Allmaras 模型 ·k-e 模型 －标准k-e 模型 －Renormalization-group (RNG) k -e 模型 －带旋流修正k -e 模型 ·k-ω模型 －标准k-ω模型 －压力修正k-ω模型 雷诺兹压力模型 大漩涡模拟模型 几个湍流模型的比较: 从计算的角度瞧Spalart-Allmaras 模型在FLUENT 中就是最经济的湍流模型,虽然只有一种方程可以解。由于要解额外的方程,标准k -e 模型比Spalart-Allmaras 模型耗费更多的计算机资源。带旋流修正的k -e 模型比标准k -e 模型稍微多一点。由于控制方程中额外的功能与非线性,RNG k -e 模型比标准k -e 模型多消耗10～15%的CPU 时间。就像k -e 模型,k -ω模型也就是两个方程的模型,所以计算时间相同。 比较一下k -e 模型与k -ω模型,RSM 模型因为考虑了雷诺压力而需要更多的CPU 时间。然而高效的程序大大的节约了CPU 时间。RSM 模型比k -e 模型与k -ω模型要多耗费50～60%的CPU 时间,还有15～20%的内存。 除了时间,湍流模型的选择也影响FLUENT 的计算。比如标准k -e 模型就是专为轻微的扩散设计的,然而RNG k -e 模型就是为高张力引起的湍流粘度降低而设计的。这就就是RNG 模型的缺点。 同样的,RSM 模型需要比k -e 模型与k -ω模型更多的时间因为它要联合雷诺压力与层流。 概念: 1、雷诺平均:在雷诺平均中,在瞬态N-S 方程中要求的变量已经分解位时均常量与变量。 相似的,像压力与其它的标量 )22.10('-+=ΛΛΛi i i φφφ 这里φ表示一个标量如压力,动能,或粒子浓度。 2、 Boussinesq 逼近从雷诺压力转化模型:利用Boussinesq 假设把雷诺压力与平均速度梯度联系起来: Boussinesq 假设使用在Spalart-Allmaras 模型、k -e 模型与k -ω模型中。这种逼近方法好处就是对计算机的要求不高。在Spalart-Allmaras 模型中只有一个额外的方程要解。k -e 模型与k -ω模型中又两个方程要解。Boussinesq 假设的不足之处就是假设u t 就是个等方性标量,这就是不严格的。 1. Spalart-Allmaras 模型(1equ):

Fluent湍流模型选取的准则

Fluent湍流模型选取的准则 湍流模型选取的准则：流体是否可压、建立特殊的可行的问题、精度的要求、计算机的能力、时间的限制。为了选择最好的模型，你需要了解不同条件的适用范围和限制。 FLUENT软件中提供以下湍流模型：1 Spalart-Allmaras 模型；2 k-ε模型； 3 k-ω模型； 4 雷诺应力模型（RSM）； 5 大涡模拟模型（LES）。 1 Spalart-Allmaras 模型 应用范围： Spalart-Allmaras模型是设计用于航空领域的，主要是墙壁束缚 （wall-bounded）流动，而且已经显示出很好的效果。在透平机械中的应用也愈加广泛。 在湍流模型中利用Boussinesq逼近，中心问题是怎样计算漩涡粘度。这个模型被Spalart-Allmaras提出，用来解决因湍流动粘滞率而修改的数量方程。 模型评价： Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，不需要求解当地剪切层厚度的长度尺度；由于没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合；比如平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显等问题。 Spalart-Allmaras模型中的输运变量在近壁处的梯度要比k-ε中的小，这使得该模型对网格粗糙带来数值误差不太敏感。 Spalart-Allmaras模型不能断定它适用于所有的复杂的工程流体。例如不能依靠它去预测均匀衰退，各向同性湍流。 2 k-ε模型 ① 标准的k-ε模型： 最简单的完整湍流模型是两个方程的模型，要解两个变量，速度和长度尺度。在FLUENT中，标准k-ε模型自从被Launder and Spalding提出之后，就变成工程流场计算中主要的工具了。适用范围广、经济、合理的精度。它是个半经验的公式，是从实验现象中总结出来的。 湍动能输运方程是通过精确的方程推导得到，耗散率方程是通过物理推理，数学上模拟相似原型方程得到的。

FLUENT中常用的湍流模型

The Spalart-Allmaras模型 对于解决动力漩涡粘性，Spalart-Allmaras 模型是相对简单的方程。它包含了一组新的方程，在这些方程里不必要去计算和剪应力层厚度相关的长度尺度。Spalart-Allmaras 模型是设计用于航空领域的，主要是墙壁束缚流动，而且已经显示出很好的效果。在透平机械中的应用也愈加广泛。 在原始形式中Spalart-Allmaras 模型对于低雷诺数模型是十分有效的，要求边界层中粘性影响的区域被适当的解决。在FLUENT中，Spalart-Allmaras 模型用在网格划分的不是很好时。这将是最好的选择，当精确的计算在湍流中并不是十分需要时。再有，在模型中近壁的变量梯度比在k-e模型和k-ω模型中的要小的多。这也许可以使模型对于数值的误差变得不敏感。想知道数值误差的具体情况请看5.1.2。 需要注意的是Spalart-Allmaras 模型是一种新出现的模型，现在不能断定它适用于所有的复杂的工程流体。例如，不能依靠它去预测均匀衰退，各向同性湍流。还有要注意的是,单方程的模型经常因为对长度的不敏感而受到批评，例如当流动墙壁束缚变为自由剪切流。 标准k-e模型 最简单的完整湍流模型是两个方程的模型，要解两个变量，速度和长度尺度。在FLUENT中，标准k-e模型自从被Launder and Spalding提出之后，就变成工程流场计算中主要的工具了。适用范围广、经济，有合理的精度，这就是为什么它在工业流场和热交换模拟中有如此广泛的应用了。它是个半经验的公式，是从实验现象中总结出来的。 由于人们已经知道了k-e模型适用的范围，因此人们对它加以改造，出现了RNG k-e模型和带旋流修正k-e 模型。k-ε模型中的K和ε物理意义：k是紊流脉动动能（J），ε是紊流脉动动能的耗散率（%）；k越大表明湍流脉动长度和时间尺度越大，ε越大意味着湍流脉动长度和时间尺度越小，它们是两个量制约着湍流脉动。 RNG k-e模型 RNG k-e模型来源于严格的统计技术。它和标准k-e模型很相似，但是有以下改进： ?RNG模型在e方程中加了一个条件，有效的改善了精度。 ?考虑到了湍流漩涡，提高了在这方面的精度。 ?RNG理论为湍流Prandtl数提供了一个解析公式，然而标准k-e模型使用的是用户提供的常数。 ?然而标准k-e模型是一种高雷诺数的模型，RNG理论提供了一个考虑低雷诺数流动粘性的解析公式。这些公式的效用依靠正确的对待近壁区域 这些特点使得RNG k-e模型比标准k-e模型在更广泛的流动中有更高的可信度和精度。 带旋流修正的k-e模型 带旋流修正的k-e模型是近期才出现的，比起标准k-e模型来有两个主要的不同点。 ?带旋流修正的k-e模型为湍流粘性增加了一个公式。 ?为耗散率增加了新的传输方程，这个方程来源于一个为层流速度波动而作的精确方程。 术语“realizable”，意味着模型要确保在雷诺压力中要有数学约束，湍流的连续性。带旋流修正的k-e模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。带旋流修正的k-e模型和RNG k-e模型都显现出比标准k-e模型在强流线弯曲、漩涡和旋转有更好的表现。由于带旋流修正的k-e模型是新出现的模型，所以现在还没有确凿的证据表明它比RNG k-e模型有更好的表现。但是最初的研究表明带旋流修正的k-e模型在所有k-e模型中流动分离和复杂二次流有很好的作用。带旋流修正的k-e模型的一个不足是在主要计算旋转和静态流动区域时不能提供自然的湍流粘度。这是因为带旋流修正的k-e模型在定义湍流粘度时考虑了平均旋度的影响。这种额外的旋转影响已经在单一旋转参考系中得到证实，而且表现要好于标准k-e模型。由于这些修改，把它应用于多重参考系统中需要注意。 标准k-ω模型 标准k-ω模型是基于Wilcox k-ω模型，它是为考虑低雷诺数、可压缩性和剪切流传播而修改的。Wilcox k-ω模型预测了自由剪切流传播速率，像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可以应用于墙壁束缚流动和自由剪切流动。标准k-e模型的一个变形是SST k-ω模型，它在FLUENT中也是可用的，将在10.2.9中介绍它。 剪切压力传输（SST）k-ω模型

fluent湍流模型

第十章湍流模型 本章主要介绍Fluent所使用的各种湍流模型及使用方法。 各小节的具体内容是： 10．1 简介 10．2 选择湍流模型 10．3 Spalart-Allmaras 模型 10．4 标准、RNG和k-e相关模型 10．5 标准和SST k-ω模型 10．6 雷诺兹压力模型 10．7 大型艾迪仿真模型 10．8 边界层湍流的近壁处理 10．9 湍流仿真模型的网格划分 10．10 湍流模型的问题提出 10．11 湍流模型问题的解决方法 10．12 湍流模型的后处理 10．1 简介 湍流出现在速度变动的地方。这种波动使得流体介质之间相互交换动量、能量和浓度变化，而且引起了数量的波动。由于这种波动是小尺度且是高频率的，所以在实际工程计算中直接模拟的话对计算机的要求会很高。实际上瞬时控制方程可能在时间上、空间上是均匀的，或者可以人为的改变尺度，这样修改后的方程耗费较少的计算机。但是，修改后的方程可能包含有我们所不知的变量，湍流模型需要用已知变量来确定这些变量。 FLUENT 提供了以下湍流模型： ·Spalart-Allmaras 模型 ·k-e 模型 －标准k-e 模型 －Renormalization-group (RNG) k-e模型 －带旋流修正k-e模型 ·k-ω模型 －标准k-ω模型 －压力修正k-ω模型 －雷诺兹压力模型 －大漩涡模拟模型 10．2 选择一个湍流模型 不幸的是没有一个湍流模型对于所有的问题是通用的。选择模型时主要依靠以下几点：流体是否可压、建立特殊的可行的问题、精度的要求、计算机的能力、时间的限制。为了选择最好的模型，你需要了解不同条件的适用范围和限制 这一章的目的是给出在FLUENT中湍流模型的总的情况。我们将讨论单个模型对cpu 和内存的要求。同时陈述一下一种模型对那些特定问题最适用，给出一般的指导方针以便对于你需要的给出湍流模型。 10．2．1 雷诺平均逼近vs LES 在复杂形体的高雷诺数湍流中要求得精确的N-S方程的有关时间的解在近期内不太可能实现。两种可选择的方法用于把N-S方程不直接用于小尺度的模拟：雷诺平均和过滤。

湍流模型

我们知道，描述流体运动（层流）的流体力学基本方程组是封闭的，而描述湍流运动的方程组由于采用了某种平均（时间平均或网格平均等）而不封闭，须对方程组中出现的新未知量采用模型而使其封闭，这就是CF D中的湍流模型。湍流模型的主要作用是将新未知量和平均速度梯度联系起来。目前，工程应用中湍流的数值模拟主要分三大类：直接数值模拟（D NS）；基于雷诺平均N-S方程组（RANS）的模型和大涡模拟（LES）。DNS是直接数值求解N-S方程组，不需要任何湍流模型，是目前最精确的方法。其优点在于可以得出流场内任何物理量（如速度和压力）的时间和空间演变过程，旋涡的运动学和动力学问题等。由于直接求解N-S方程，其应用也受到诸多方面的限制。第一：计算域形状比较简单，边界条件比较单一；第二：计算量大。影响计算量的因素有三个：网格数量、流场的时间积分长度（与计算时间长度有关）和最小旋涡的时间积分长度（与时间步长有关），其中网格数量是重要因素。为了得到湍流问题足够精确的解，要求能够数值求解所有旋涡的运动，因此要求网格的尺度和最小旋涡的尺度相当，即使采用子域技术，其网格规模也是巨大的。为了求解各个尺度旋涡的运动，要求每个方向上网格节点的数量与Re3/4成比例，考虑一个三维问题，网格节点的数量与Re9/4成比例。目前，DNS能够求解Re(104)的范围。 基于RANS的湍流模型采用雷诺平均的概念，将物理量区分为平均量和脉动量，将脉动量对平均量的影响用模型表示出来。目前，基于RANS方程已经发展了许多模型，几乎能对所有雷诺数范围的工程问题求解，并得出一些有用的结果。其缺点在于：第一：不同的模型解决不同类型的问题，

几种湍流模型

解决湍流的模型总计就是那几个方程，Fluent 又从工程和数值的角度进行了整理，下面就是这些湍流模型的详细说明。 FLUENT 提供了以下湍流模型： ·Spalart-Allmaras 模型 ·k-e 模型 －标准k-e 模型 －Renormalization-group (RNG) k -e 模型 －带旋流修正k -e 模型 ·k-ω模型 －标准k-ω模型 －压力修正k-ω模型 雷诺兹压力模型 大漩涡模拟模型 几个湍流模型的比较： 从计算的角度看Spalart-Allmaras 模型在FLUENT 中是最经济的湍流模型，虽然只有一种方程可以解。由于要解额外的方程，标准k -e 模型比Spalart-Allmaras 模型耗费更多的计算机资源。带旋流修正的k -e 模型比标准k -e 模型稍微多一点。由于控制方程中额外的功能和非线性，RNG k -e 模型比标准k -e 模型多消耗10～15%的CPU 时间。就像k -e 模型，k -ω模型也是两个方程的模型，所以计算时间相同。 比较一下k -e 模型和k -ω模型，RSM 模型因为考虑了雷诺压力而需要更多的CPU 时间。然而高效的程序大大的节约了CPU 时间。RSM 模型比k -e 模型和k -ω模型要多耗费50～60%的CPU 时间，还有15～20%的内存。 除了时间，湍流模型的选择也影响FLUENT 的计算。比如标准k -e 模型是专为轻微的扩散设计的，然而RNG k -e 模型是为高张力引起的湍流粘度降低而设计的。这就是RNG 模型的缺点。 同样的，RSM 模型需要比k -e 模型和k -ω模型更多的时间因为它要联合雷诺压力和层流。 概念： 1.雷诺平均：在雷诺平均中，在瞬态N-S 方程中要求的变量已经分解为时均常量和变量。 相似的，像压力和其它的标量 )2 2.10('-+= i i i φφφ 这里φ表示一个标量如压力，动能，或粒子浓度。 2. Boussinesq 逼近从雷诺压力转化模型：利用Boussinesq 假设把雷诺压力和平均速度梯度联系起来： Boussinesq 假设使用在Spalart-Allmaras 模型、k -e 模型和k -ω模型中。这种逼近方法好处是对计算机的要求不高。在Spalart-Allmaras 模型中只有一个额外的方程要解。k -e 模型和k -ω模型中又两个方程要解。Boussinesq 假设的不足之处是假设u t 是个等方性标量，这是不严格的。

湍流模型概述

大多数飞行器都是在高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）涡粘性模式

湍流燃烧模型-PDF

PDF 模型 概率密度函数PDF方法以随机的观点来对待湍流问题，对解决湍流化学反应流的问题具有很强的优势。在湍流燃烧中存在一些非输运量( 如反应速率, 密度, 温度及气相体积分数等) 的湍流封闭问题。尽管这些量没有输运方程, 但它们常常是输运变量的已知函数。平均或者过滤高度非线性的化学反应源项会引起方程的封闭问题。因此，用PDF的方法来解决这些非输运量的湍流封闭问题显然是一个既简单又直接的途径。 PDF方法是一种较为流行的湍流燃烧模型，能够较为精确的模拟任何详细的化学动力学过程, 适用于预混、非预混和部分预混的任何燃烧问题。目前, 确定输运变量脉动概率密度函数的方法有输运方程和简化假定两种, 分别称之为输运方程的PDF和简化的PDF。前者建立输运变量脉动的概率密度输运方程，通过求解该方程来获得输运变量脉动的概率分布。后者假定输运变量脉动的概率密度函数的具体形式, 通过确定其中的一些待定参数来获得输运变量脉动的概率分布。湍流燃烧中, 后者应用最为普遍和广泛。在简化的PDF 中, 输运变量脉动的概率密度函数常常采用双 D 分布、截尾高斯分布和B 函数分布等形式。 PDF在理论上可以精确考虑任意详细的化学反应机理，但是其具体求解时需借助其它的模型和算法，而且计算量相对较大。PDF的方程是由N-S方程推导而来，其中的化学反应源项是封闭的，但压力脉动梯度项以及分子粘性和分子扩散引起的PDF的分子输运项是不封闭的，需要引入模型加以封闭。例如，在速度- 标量-湍流频率PDF中，必须采用小尺度混合模型、随机速度模型和湍流频率模型加以封闭。 模化后的输运方程难以用有限容积、有限差分和有限元等方法来求解，比较可行的一种方法是蒙特卡洛(MonteCarlo)方法,在该方法中输运方程被转化为拉格朗日(Lagrangian)方程,流体由大量遵循Lagrang ian方程的随机粒子的系统来描述, 最后对粒子作统计平均得到流场物理量和各阶统计矩。另有与有限容积法相结合的蒙特卡洛法。 PDF 模型的发展 1969年Lungdren首先推导、计算了速度的联合PDF运输方程，避免了对梯度扩散模型进行模拟，对很简单的流动过程得到了简析解[1]。

湍流模型发展综述

湍流模型发展综述 摘要：在概述了湍流问题的基础上，本文简要介绍了湍流的四种模型，对湍流模型在不同情况下的模拟能力进行了对比，最后简述了湍流模型的发展方向。 关键词：湍流模型；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方程的计算。此两模型已在工程计算中

湍流模型概述

大多数飞行器都是在高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）涡粘性模式

四种湍流模型介绍

由于航发燃烧室中的流动特性极其复杂，要想提高数值计算的预测能力,必须要慎重选择湍流模型。用四种不同的湍流模型对带双径向旋流杯的下游流场进行数值模拟,将计算结果与实验结果作对比，比较各湍流模型的原理和物理基础，优劣，并分析流场速度分布和回流区特性。 涉及的湍流模型： 标准k-ε湍流模型（SＫE) １标准k-ε湍流模型有较高的稳定性,经济性和计算精度，应用广泛,适合高雷诺数湍流，但不适合旋流等各向异性较强的流动。 2简单的湍流模型是两个方程的模型，需要解两个变量，即速度和长度。在fｌｕｅｎt中，标准 k-ε湍流模型自从被Lauｎdeｒand Spalｄing 提出之后，就变成流场计算中的主要工具。其在工业上被普遍应用,其计算收敛性和准确性都非常符合工程计算的要求。 3但其也有某些限制,如ε方程包含不能在壁面计算的项，因此必须使用壁面函数。另外，其预测强分离流,包含大曲率的流动和强压力梯度流动的结果较弱。 它是个半经验的公式,是从实验现象中总结出来的。 动能输运方程是通过精确的方程推导得到,耗散率方程是通过物理推理，数学上模拟相似原型方程得到的。 应用范围：该模型假设流动为完全湍流，分子粘性的影响可以忽略,此标准κ－ε模型只适合完全湍流的流动过程模拟。 可实现的ｋ-ε模型是才出现的,比起标准k-ε模型来有两个主要的不同点：·可实现的k－ε模型为湍流粘性增加了一个公式。 ·为耗散率增加了新的传输方程,这个方程来源于一个为层流速度波动而作的精确方程。 术语“rｅalizablｅ”,意味着模型要确保在雷诺压力中要有数学约束，湍流的连续性。 应用范围： 可实现的k-ε模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。 可实现的k－ε模型和RＮG k－ε模型都显现出比标准k－ε模型在强流线弯曲、漩涡和旋转有更好的表现。由于带旋流修正的k-ε模型是新出现的模型，所以还没有确凿的证据表明它比ＲＮＧｋ－ε模型有更好的表现。但是最初的研究表明可实现的k-ε模型在所有k-ε模型中流动分离和复杂二次流有很好的作用。 该模型适合的流动类型比较广泛,包括有旋均匀剪切流,自由流(射流和混合层）,腔道流动和边界层流动。对以上流动过程模拟结果都比标准k-ε模型的结果好，特别是可再现k－ε模型对圆口射流和平板射流模拟中,能给出较好的射流扩张。

湍流模型

湍流模型概述 湍流是一种复杂的非稳态三维流动,通常把瑞流定义为具有随机性、扩散性、高雷诺数、三维祸量脉动性、耗散性及连续性特征的复杂流动。虽然瑞流具有多种特性,但瑞流不是流体本身具有的某些特征而是流体流动的特征,仍是一种连续流动,仍然同层流一样满足流动的基本方程。从数学的观点看,瑞流是N-S方程的 通解,求解端流与求解层流无本质区别,目前己具有足以求解瑞流问题的有关方程式。端流还可以看作是由多种大尺度祸流和小尺度祸流组成的特殊流动。大尺度的祸流主要由流动的边界条件和流动区域的几何形状所决定,是引起流场中低频 脉动的主要原因;小尺度的祸流主要是點性力所决定,是引起流场中高频脉动的主要原因。瑞流的物理量的脉动特点就是由于流体内各种不同尺度祸流的随机运动造成。 用数值方法直接计算瑞流单元运动规律时,计算网格尺寸要小于瑞流单元 尺度,并在瑞流单元尺度内计算N-S方程的通解。但是在实际工程中具有重要意 义的不是端流的精细结构,而是瑞流对于时间的平均(时均)效应。因此,雷诺首先提出了将N-S方程对某一时间比例尺取平均,得到时均N-S方程。虽然瑞流的N-S 方程经过时均化处理后方程式的形式可以保持不变,但是出现了脉动应力项(雷诺应力),因此需要提出相应的端流模型(一个或一组数学方程)使时均方程得到封闭。这种方法按雷诺应力方程模型化方法的不同可分为两类:一类是直接就雷诺应力 建立模型化方程的雷诺应力方程模型;另一类是在雷诺应力与局部时均速度梯度 成比例的Boussinesq假设下引入的瑞流黏度系数模型。另一种瑞流数值计算方法是亚网格尺度模拟,即大祸模拟(LES),由N-S方程出发直接模拟大尺度祸流,小尺度祸流的影响可以通过近似模型来考虑。但是由于大祸模拟计算量仍很大,也只能 模拟一些简单的情况。 工程上通常需要深入了解的是温度场、时均速度场、瑞流脉动时均特性等, 并不需要了解瑞流产生和发展的详细过程。因此,利用雷诺提出的时均值的概念 来研究瑞流运动的方法是一种有效的简化,从N-S方程导出瑞流平均运动方程和 雷诺方程,还导出了连续性方程和能量方程等基本方程。雷诺平均法将瑞流物理 量代入不可压缩瞬态连续性方程、动量方程得到端流平均运动的连续性方程和动量方程。但是在雷诺时均方程组中除了瞬态连续性方程和动量方程外还有一项是

湍流燃烧及其数值模拟

湍流燃烧及其数值模拟研究 1. 湍流燃烧 1.1湍流燃烧基本概念 当流动雷诺数数较小时，由于流体粘性的作用，流体呈层流流态。当流动的特征雷诺数超过相应的临界值，流动从层流转捩到湍流。湍流燃烧是指湍流流动中可燃气的燃烧，在能源、动力、航空和航天等工程领域,经常遇到的实际燃烧过程几乎全部都是湍流燃烧过程。湍流燃烧实质是湍流，化学反应和传热传质等过程相耦合的结果。湍流对燃烧的影响与湍流强度和湍流涡旋尺度有关。小尺度湍流通过湍流扩散使火焰区内的输运效应增加，从而使化学反应速率增加。但气流脉动不会火焰面产生皱褶，只能把火焰变成波纹状。大尺度湍流对火焰内部结构没有影响，但使火焰阵面出现皱褶，增加其燃烧面积，造成火焰表现传播速度增加。当湍流强度及湍流尺度均较大时，火焰前沿不再连续而分裂成四分五裂。 燃烧对湍流的影响主要表现在燃烧释放的热流流团膨胀，影响气体的密度和运动速度，从而影响当地的涡旋，湍流强度和湍流结构。 1.2湍流燃烧分类 湍流燃烧按其燃料和氧化剂的初始混合状态可以分类为：湍流非预混燃烧、预混燃烧和部分预混燃烧。在湍流非预混燃烧燃料和氧化剂事先是分离的，燃料和氧化剂一边混合一边燃烧，燃烧速率主要受湍流混合过程控制，而在湍流预混燃烧中，燃料和氧化剂在进入核心燃烧区以前已经充分混合，化学反应的速率由火焰前缘从炽热的燃烧区向冷态无反应区的传播所控制。上面两种燃烧方式是湍流燃烧的两个极限情形，很多情况下两种燃烧模式是并存的，称为部分预混燃烧。部分预混燃烧可出现在下列情形中叫：(1)在一个完全以非预混燃烧为配置的燃烧装置发牛了局部熄火；(2)当预混火焰前缘穿过非均匀的混气时；(3)射流非预混火焰发生抬举，其根部是一个典型的部分预混火焰。这三种部分预混燃烧情形涉及了经常受到关注的燃烧研究话题如局部熄火、火焰稳定等，它们对研究湍流燃烧过程的机理有很大意义。 在湍流燃烧中,湍流流动过程和化学反应过程有强烈的相互关联和相互影响.湍流通过强化混合而影响着时平均化学反应速率,同时化学反应放热过程又影响着湍流,如何定量地来描述和确定这种相互作用是湍流燃烧研究的一个重要内容. 湍流是非常复杂的，它包括湍流问题，湍流与燃烧的相互作用，流动参数与化学动力参数之间的耦合机理等问题。因此湍流燃烧是工程科学中最复杂的领域之一。 湍流燃烧的研究已进行多年，研究的方法有试验研究，理论分析和数值模拟等。计算流体力学和计算机技术的发展，数值模拟由于它的廉价性和可操作性在国际上受到越来越多的重视，得到了广泛的应用。 2.湍流燃烧数值模拟 2.1湍流燃烧数值模拟简介 湍流燃烧数值模拟(Numerical Simulation of Turbulent Combustion)是指应用计算机为工

几种湍流模型

解决湍流的模型总计就是那几个方 程， 就是这些湍流模型的详细说明。FLUENT 提供了以下湍流模型： ·Spalart-Allmaras 模型Fluen t 又从工程和数值的角度进行了整理，下面 ·k- e 模型 －标准k-e模型 －Renormalization-group(RNG)k-e 模型－带旋流修正k-e模型 ·k-ω模型 －标准k-ω模型 －压力修正k-ω模型 雷诺兹压力模型 大漩涡模拟模型 几个湍流模型的比较： 从计算的角度看Spalart-Allmaras模型在FLUENT中是最经济的湍流模型，虽然只有一 种方程可以解。由于要解额外的方程，标准k-e模型比Spalart-Allmaras模型耗费更多的计 算机资源。带旋流修正的k-e模型比标准k-e模型稍微多一点。由于控制方程中额外的功能和 非线性，RNGk-e模型比标准k-e模型多消耗10～15%的CPU时间。就像k-e模型，k-ω模型 也是两个方程 的模型，所以计算时间相同。 比较一下k-e模型和k-ω模型，RSM模型因为考虑了雷诺压力而需要更多的CPU时间。然 而高效的程序大大的节约了CPU时间。RSM模型比k-e模型和k-ω模型要多耗费50～60%的CPU 时间，还有15～20%的内存。 除了时间，湍流模型的选择也影响FLUENT的计算。比如标准 k-e模型是专为轻微的扩散 设计的，然而RNGk-e模型是为高张力引起的湍流粘度降低而设计的。这就是RNG模型的缺点。 同样的，RSM模型需要比k-e模型和k-ω模型更多的时间因为它要联合雷诺压力和层流。 概念： 1.雷诺平均：在雷诺平均中，在瞬态N-S方程中要求的变量已经分解为时均常量和变量。 相似的，像压力和其它的标量 ' (10.22) i i i 这里表示一个标量如压力，动能，或粒子浓度。 2.Boussinesq逼近从雷诺压力转化模型：利用Boussinesq假设把雷诺压力和平均速度梯度 联系起来： Boussinesq假设使用在Spalart-Allmaras模型、k-e模型和k-ω模型中。这种逼近方法好处是对 计算机的要求不高。在Spalart-Allmaras模型中只有一个额外的方程要解。k-e模型和k-ω模型中又两个方程要解。 Boussinesq假设的不足之处是假设ut是个等方性标量，这是不严格的。

fluent中燃烧模型分类

FLUENT燃烧模型 化学反应 模拟方法 方法描述计算反应的选择 有限速率模型需要求解组分质量分数的 输运方程，化学反应机理 由用户自己定义。反应速 率在组分输运方程中作为 源项，并由阿累尼乌斯公 式计算。应用范围最广泛。 应用：模拟化学组分混合、 输运和反应的问题；壁面 或粒子表面反应问题 层流有限 速率模型 使用Arrhenius公式计算化学源项，忽略湍流脉动的影响。对于化学动力学控制的燃烧（如层流燃 烧），或化学反应相对缓慢的湍流燃烧是准确的。但对一般湍流火焰中Arrhenius化学动力学的高度 非线性一般不精确；对于化学反应相对缓慢、湍流脉动较小的燃烧（如超音速火焰）可能可以接受。 漩涡破碎模型 Eddy Dissipation 大部分燃料快速燃烧，整体反应速率由湍流混合控制。复杂且常是未知的化学反应动力学速率可以 完全的被忽略掉。化学反应速率由大尺度涡混合时间尺度k/ε控制。只要k/ε（湍流）出现，燃烧 即可进行，不需要点火源来启动燃烧。（缺点：未能考虑分子输运和化学动力学因素的影响） 适用条件：高雷诺数湍流预混燃烧过程。 EBU-Arrehenius 模型 EDC模型 假定化学反应都发生在小涡中（精细涡），反应时间由小涡生存时间和化学反应本身需要的时间共 同控制。EDC模型能够在湍流反应中考虑详细的化学反应机理。但是他们的数值积分计算开销很大。 使用条件：只有在快速化学反应假定无效的情况下才能使用这一模型（如快速熄灭火焰中缓慢的 CO烧尽、选择性非催化还原中的NO转化问题）。 非预混燃烧模型不求解每个组分的质量分数输运方程，求解混合分数输运方程和一个或两个守恒标量的方程，然后从预测的混合分数公布推导出每一个组分的浓度。通过概率密度函数或PDF来考虑湍流的影响。 应用：主要用于模拟湍流扩散火焰的反应系统。这个系统要求接近化学平衡，氧化物和燃料以两个或者三个进口进入计算域。 预混燃烧模型主要用于单一、完全预先混合好的燃烧系统。反应物和燃烧产物被火焰前沿分开。求解出反应发展变量来预测前沿的位置。湍流的影响通过湍流火焰速度计算。 部分预混燃烧模型描述非预混燃烧完全预混燃烧相结合的系统。结合混合分数方程和反应物发展变量来分别确定组分浓度和火焰前沿位置。适用于计算域内具有变化等值比率的预混火焰情况。通过求解混合分数方程和反应过程参数来确定火焰峰面的位置。 PDF输运方程模型结合CHEMKIN可以考虑详细的化学反应机理，高度的非线性化学反应项是精确模拟，无须封闭模型，可以合理的模拟湍流和详细化学反应动力学之间的相互作用，是模拟湍流燃烧的精确模拟方法。但计算量特别大。 优点：可以计算中间组分；考虑分裂影响；考虑湍流-化学反应之间的作用；无需求解组分输运方程 缺点：系统须满足（靠近）局部平衡；不能用于可压缩或非湍流流动；不能用于预混燃烧。