湍流理论和湍流模型(博士课程)要点详解.ppt

S-A模型的数值求解 湍流模型求解与N-S方程组的求解可采用“松耦合”的方式，即在同一次时间推进中它们的求解是相对独立的。S-A湍流模型方程的对流项和耗散项采用中心格式进行有限体积离散。 这样，在第 i 个网格单元上有 （左端第一项） （左端第二项） （右端第一项） （右端第二项） （右端第三项） （右端第四项） 经过上述空间离散后，湍流模型方程可以写为下面的半离散形式 其中， 是湍流模型方程的残值项。 采用隐式欧拉方法进行时间离散求解，表达式如下， 将方程右端的残值项进行一阶泰勒展开，略去高阶项后得到 式中的偏导数 可以通过数值求解雅可比矩阵的方法求得，即 其中 为一个很小的正数，一般取 经过整理后，方程变为 这个方程采用高斯-赛德尔迭代进行求解。 边界条件方面，固壁面边界上 ；自由流处 。 初始条件可将 取定为自由流的值。 （三）二方程模式(标准 k-ε 模式) k- ε 模式是在涡粘模式的基础上发展起来的，它和代数模式的主要区别是在于k- ε模式的涡粘系数包含部分历史效应。具体来说，它把涡粘系数和湍动能及湍动能耗散联系在一起。 湍动能耗散率的输运方程为 用量纲分析，涡粘系数可以写成 在湍流模式中，湍动能耗散是最难构造准确模型的。 通常采用的ε模式是依据类比方法，基本思想是：湍动能耗散的生成、扩散以及消耗等项与湍动能方程中的对应项（生成、扩散和耗散）有类似的机制和公式。 k-ε方程的封闭方程为 参考文献：Abid, R., “Evaluation of Two-Equation Turbulence Models for Predicting Transitional Flows,” International Journal of Engineering Science, Vol. 31, pp. 831-840, 1993. （四）二方程模式 (k-ω 模型) 参考文献：Wilcox, D., Turbulence Modeling for CFD, DCW Industries, Inc., La Canada, California, 1993. 湍流粘性系数为 式中 相关常数为 （五）二方程模式 (Menter SST 模型) Menter 使用一个混合函数将和模型优缺点进行综合，克服了各自的缺点，得到了一种更鲁棒的湍流模型，通常称为Menter SST模型。物面附近，混合函数趋于0，模型趋于方程；远离物面则趋于1，模型趋于方程。相对而言，模型提高了强压力梯度的预测精度，增强了压力引起的边界层分离的计算能力。 Menter SST模型的湍流涡粘系数为 Menter SST模型方程可以写成 参考文献：Menter F.R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications. AIAA Journal, 1994, 32: 1598-1605. 其中基本常数为 其它常数的计算方法是 式中 （六）雷诺应力模式 雷诺应力模式简称为RSM（Reynolds Stress Model）。 由于雷诺平均运动方程所包含的雷诺应力为一未知数，所以最直接的办法是建立雷诺应力 的方程。 得到 的步骤： （1）把 N-S 方程化成 的形式 （2）将瞬时速度表示成平均速度与脉动速度之和 （3）代入 形式的方程，再减去平均运动的形式 ，即可得 雷诺应力 的方程为： ， 式中已标明了相应的物理意义。 课程考查试题 （1）由张量表达形式的N-S方程出发，推导雷诺平均运动方程。 （2）从流体微元物理运动的角度，简要推导湍流流动中附加 切应力-雷诺应力的表达形式。 （3）简要推导基于普朗特混合长度理论建立的湍流粘性系数 表达式。 （4）简述二维湍流边界层的结构。 （5）简要对比分析各种湍流模型。 Thank you！ 简单剪切湍流 近代湍流研究的重大进展之一：发现剪切湍流中存在拟序结构。 湍流脉动并非完全不规则的随机过程，而是在不规则的脉动中包含可辨认的有序大尺度运动。 剪切湍流： 壁湍流，即壁面附近的剪切湍流，例如槽道、圆管、边界层湍流； 自由剪切湍流，例如射流、混合层、远场尾流。 简单剪切湍流：平均剪切流动是平行流动或准平行流动。 研究雷诺数很高的壁湍流， H是直槽宽度之半, 或圆管半径, 或边界层的平均名义厚度, Um是平均特征速度。 设定坐标系：x为流动方向，y为垂直壁面方向，z为平均流动的展向 进一步假定直槽沿展向无限长，流向单位长度上的平均压降是常数，则 平均运动是定常的单向平行直线运动