北航湍流模式大作业.doc

.北 京 航 空 航 天 大 学湍流模式大作业学院名称 航空科学与工程学院 专业方向 飞行器设计 学 号 学生姓名 日 期 2017 年 1月14 日 ;.摘要本文首先介绍了层流与湍流的概念，湍流的几种数值模拟方法及其优缺点。由于时间所限，本文仅对雷诺平均数值模拟(RANS)的封闭模式进行了详细讨论，包括零方程模式，一方程模式，二方程模式和雷诺应力模式。在雷诺方程中建立雷诺应力和平均运动场间的关系，就可以使雷诺方程封闭，此模式称为低阶矩模式。如果以雷诺方程和雷诺应力输运方程联立作为湍流统计平衡的控制方程，这种模式称为二阶矩模式。此外，依据封闭关系式的数学方程模式，湍流模式又可分为代数模式和微分模式。目录摘要11. 介绍31.1 层流和湍流31.2 湍流的几种数值模拟方法及其优缺点32. 湍流的平均化、雷诺粘性应力43. 湍流模型63.1 零方程模型63.2 一方程模型63.3 二方程模型73.4 雷诺应力模式71. 介绍1.1 层流和湍流粘性流体的运动存在着两种完全不同的流动状态：层流状态和湍流状态。雷诺首先于1883年通过做圆管内流动实验观察到层流与湍流现象。当圆管中流动速度较小时，管中的流线之间层次分明，互不掺混，这样的流动称为层流。当流速增大后，流体作复杂、无规律、随机的不定常运动，称为湍流。流动状态与雷诺数、下临界雷诺数和上临界雷诺数有关。当时，流动为层流；当时，流动为不稳定过渡状态；当时，流动为湍流。湍流是在连续介质范畴内流体的不规则运动，它有别于物质分子的不规则运动。具体来说，在极不规则的湍流中，流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度。因此湍流运动产生的质量和能量的输运将远远大于分子热运动产生的宏观输运。1.2 湍流的几种数值模拟方法及其优缺点如果从流动控制方程出发，进行最精细的数值模拟，称为直接数值模拟(DNS, Direct Numerical Simulation)。湍流的直接数值模拟(DNS)不存在封闭性问题，原则上可以求解所有湍流问题。由于计算机资源的限制，迄今为止DNS只能求低雷诺数的简单湍流问题；若只需预测湍流的统计量时，从雷诺平均方程出发的模拟，称为雷诺平均数值模拟(RANS, Reynolds Average Numerical Simulation)。RANS的缺点是缺乏普适性，要使这种方法成为一种简单可靠的方法还需要很多的工作；介于DNS和RANS之间的是大涡模拟(LES, Large Eddy Simulation)。LES的思想是大尺度脉动用数值方法计算，小尺度脉动对大尺度的作用通过模型假设处理，有较宽的适用范围。2. 湍流的平均化、雷诺粘性应力经典的湍流理论认为，湍流是一种完全不规则的随机运动，湍流场中的物理量在时间和空间上呈随机分布，不同的瞬时有不同的值，关注某个瞬时的值是没有意义的。因此，雷诺首创用统计平均方法来描述湍流的随机运动，即对各瞬时量进行平均得到有意义的平均值。从N-S方程出发，利用平均化运算的法则推导平均物理量满足的方程组。只考虑不可压缩流体情形，假设体力可以忽略，此时，N-S方程具有下列形式 (1)运用(1)式中的连续性方程，运动方程可改写为 (2)对方程组(2)中各式两边进行平均化运算，并利用平均化运算法则得到 (3)考虑到方程组(3)的第四式，方程组(3)中的头三个方程可改写成另一种形式，把脉动项移到右边，得到(4)将上式和应力形式的运动方程对比其中是应力张量，则有其中是单位张量。是平均运动变形速度张量，是对称的二阶张量，它的表达式如下 (5)由(5)式得到在湍流运动中除了平均运动的粘性应力外，还多了一项由于脉动所引起的应力，称为雷诺粘性应力。平均流动元所满足的方程(4)称为雷诺方程。上述雷诺方程是不封闭的，方程个数只有4个，而未知函数有10个，即、及六个雷诺粘性应力分量。为了使方程组封闭，必须在雷诺应力即平均速度之间建立补充关系式。3. 湍流模型所谓湍流模型理论，就是根据理论和经验，对雷诺平均运动方程和雷诺应力项建立表达式或方程，然后对雷诺应力方程的某些项提出尽可能合理的模型和假设，使方程组封闭求解的理论。3.1 零方程模型所谓零方程模型，就是在雷诺方程和连续方程以外，不需要另外再加任何方程式来使方程组封闭。即雷诺粘性应力能直接用某些物理量和物理常数表达出来，没有引进高阶微分方程，又称代数模型。(1) Boussinesq涡粘性系数法一般的方法利用Boussinesq假设把雷诺压力和平均速度梯度联系起来： (6)(2) Cebeci-Smith模型该模型通常写成涡粘性系数形式，在内层有 (7)式中是衰减长度系数，当壁面不存在质量引射时为常数。在外层有 (8)式中为边界层位移厚度，是外部自由流速度，为间隙因子。3.2 一方程模型(1) The Spalart-Allmaras模型Spalart-Allmaras模型的变量中是湍流动粘滞率除了近壁区域，方程是： (9)这里是湍流粘度生成的，是被湍流粘度消去，发生在近壁区域。对于解决动力漩涡粘性，S-A模型是相对简单的方程，不必要去计算和剪应力层厚度相关的长度尺度。(2) 简化的二方程模型由于方程难模化和精度低，可以在保留方程的基础上舍弃方程，将表示为 式中，可视为混合长度，必须根据具体流场情况另外给出。3.3 二方程模型(1) 标准模型标准模型是个半经验公式，主要是基于湍流动能和扩散率。方程是个精确方程，方程是个由经验公式导出的方程。湍流动能方程方程为 (10)扩散方程方程为 (11)方程中表示由层流速度梯度而产生的湍流动能，是由浮力产生的湍流动能，由于在可压缩湍流中，过渡的扩散产生的波动。(2) 标准模型标准模型是一种经验模型，是基于湍流能量方程和扩散速率方程。 (12) (13)在方程中，是由层流速度梯度而产生的湍流动能。是由方程产生的，和表明了和的扩散率。和由于扩散产生的湍流。3.4 雷诺应力模式雷诺应力模式简称为RSM，由于雷诺平均运动方程所包含的雷诺应力为一未知项，所以最直接的办法是建立雷诺应力的方程。(1) 雷诺应力方程和湍流动能方程得到的方程步骤为：把N-S方程化成的形式，将瞬时速度表示成平均速度与脉动速度之和，即、，将其代入形式方程，减去平均运动的形式，最终得到雷诺应力方程为 (14)在模式理论中要用到湍动能，在式(14)中令i=j得到 (15)这就是湍动能方程。(2) 经模化后的雷诺应力方程和湍动能方程在雷诺应力方程(14)和湍动能方程(15)中存在未知项，如、等，因此