fluent教程 第三讲,湍流模拟.pptVIP

湍流模拟 什么是湍流? 湍流：非定常，非周期性的三维速度脉动、? 强化物质、动量和能量的输运. 瞬时速度分解为平均速度和脉动速度: Ui(t) ? Ui + ui(t) 压力、温度、组分浓度值具有类似的脉动 平均量输运方程 雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有 ： 其中， 和 分别是平均速度和脉动速度（i=1,2,3） 类似地，对于压力等其它标量，我们也有 ： 把上面的表达式代入瞬时的连续与动量方程，并取平均 （去掉平均速度 上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式 平均量输运方程（续） 湍流模拟的方法 直接数值模拟（DNS）只适合低雷诺数流动。 求解雷诺平均的 Navier-Stokes (RANS) 方程: 其中 (雷诺应力) 时间平均的湍流脉动用实验常数和流场平均速度信息来模拟. 大涡模拟（Large Eddy Simulation，LES）数值求解大涡，小涡用模型. 湍流输运模型 涡旋粘性系数模型 涡旋粘性系数模型是根据湍流应力和分子运动引起的粘性应力相似而提出来的 雷诺应力模型 雷诺应力模型则是从基本方程出发，直接推出雷诺应力的方程，但在雷诺应力的方程中包含有更高阶的相关项，对这些更高阶的相关项再建立相应的模型。 湍流应力 湍流应力 在三维空间中，下标和分别可取为1，2和3，所以湍流应力有9个分量组成，是一个二阶张量，若用矩阵形式表示可写为： 主对角线上的三个分量，和称为湍流正应力，其余的六个分量称为湍流切应力。三个湍流正应力之和是湍流脉动动能的两倍 湍流应力 若以主对角线上的三个分量作为对称轴，则对称的两个切应力分量是相等的，很显然，这是一个对称的二阶张量。很容易可以证明，在各向同性湍流中，湍流正应力的三个分量相等，即 判断湍流 湍流特点 额外应变率 流向曲率 测向分离 加速或减速 有旋 回流 (或分离) 二次流 3D振荡流动 Transpiration (吹风/吸气) 自由湍流 剪切层相互作用 需要作出选择 湍流模型 涡旋粘性系数模型 涡旋粘性系数模型是Boussinesq于1877年最早提出来的。他假定在近似平行的剪切流中，湍流应力张量中的切应力分量和平均速度在横向方向的梯度成正比，其比例系数称为涡旋粘性系数，以后把这个假定再推广到三维的流动。即湍流应力张量和平均流场应变率之间有关系： 仿照分子运动引起的粘性系数，我们可以假定涡旋粘性系数和湍流中含能涡旋的特征长度和速度尺度成正比，即： 而根据确定速度尺度V和长度尺度L方法的不同，又可以细分为各种不同的模型。而其中最简单，V和L的确定方法又是一致的模型应该是双方程模型。 涡旋粘性系数模型 在双方程模型中，假定特征速度V和湍流动能k的平方根成正比，特征长度由湍流动能k和另外一个辅助的量确定。比如在k-e模型中，辅助的量选为湍流动能的耗散率，根据量纲分析得长度尺度 涡旋粘性系数模型 辅助的量是湍流频率 ， 根据量纲分析得长度尺度 Reynolds Stress Terms in RANS-based Models 湍流模型 0方程模型 单方程模型 双方程模型 雷诺应力模型 大涡模拟 单方程模型: Spalart-Allmaras 湍流粘性系数: 求解 的输运方程: The additional variables are functions of the modified turbulent viscosity and velocity gradients. 单方程模型: Spalart-Allmaras 用于计算航空中考虑边界作用的高速流动问题 可以考虑边界层逆压梯度 叶轮机械流动 壁面网格可粗可细 网格细了，可以考虑低雷诺数流动，考虑边界层影响. 网格粗时候，可以得到相对较好的流动结果. 应用举例 无粘流结果 单方程湍流模型结果 双方程模型: 标准 k-? 模型 双方程模型: 标准 k-? 模型 “Baseline model” (Two-equation) 最广泛用于工程流动与换热计算 结果比较透彻，优点、缺点明显 半经验性 k 方程经过严