6.流动数学模型的改进.pptVIP

§4 近壁区数学模型使用问题探讨  1、近壁区流动特点： 1）壁面区：沿壁面法线不同距离上的流动分为： ①粘性底层：粘性力在动量，热量及质量交换中起主导作用， 湍流切应力可忽略，u呈线性分布； ②过渡层：粘性力可忽略，湍流切应力起主导作用，u 呈对数分布； 2）核心区：完全湍流区 2、近壁区数学模型引入原因： 标准的Κ-ε模型，只是针对充分发展的湍流才是有效的。但对于固体壁面附近的粘性底层内的流动和换热问题，特别是在粘性底层内，湍流应力几乎不起作用。故不能用以上标准的模型，对于这个特殊区域，可以采用低Re数模型或壁面函数法。 3、解决近壁区途径 1）壁面函数法：不对粘性影响明显区域（粘性底层及过渡层）求解，而用一组半经验公式（壁面函数）将壁面上的物理量与湍流核心区内相应物理量联系起来的方法。 2）采用低Re数Κ-ε模型：求解粘性影响较明显区域，要求在壁面区域划分较细的网格。 4、壁面函数法的基本思想： 1）对于湍流核心区流动使用Κ-ε模型求解，近壁区不求解均流场或湍流量的偏微分方程。这是与低Re数Κ-ε模型不同之处； 2）在粘性底层内不布置任何节点，把与壁面相邻的第一个节点布置在旺盛湍流区域内（对数律层）。近壁影响全部集中在第一个节点所在控制体内；此时壁面上的切应力与热流密度仍然按第一个内节点与壁面上的速度和温度的一阶差分来计算； 3）关键是使用半经验公式将壁面上的物理量与湍流核心区内相应物理量联系起来，合理选择方程在边界节点处的有效扩散系数Г以及Κ ， ε的边界条件，以使计算所得的切应力与热流密度能与实际情形基本相符。这种方法能节省内存和计算时间，在工程湍流计算中应用较广。 4.1 壁面函数法单层模型： 特点： 1）近壁的一个网格深入到旺盛湍流区； 2）湍流边界层流动要有一些特征： ①壁面附近速度按对数律分布。 ②从壁面到第一个节点间的切应力均匀。 ③在此区域内脉动动能的产生和耗散相平衡。 1．假设在所计算问题的壁面附近粘性底层以外的地区，无量纲速度与温度服从对数分布律： 其中 （B=5.0~5.5）； σL ， σT ---分子 Pr 数及湍流 Pr 数； E-------壁面光滑系数 9.0 A---为常数，对于光滑管取为26； K----冯.卡门系数 0.4~0.42， --------壁面摩擦速度； （1） （2） （即壁面附近湍动能生成率与耗散率相抵） 下式中的三个无量纲量定义式中，引入了湍流参数k ，同时 又保留了壁面切应力τw 及热流密度 qw，。正是后面的两个量 是工程计算中的主要求解对象。 2．在划分网格时，把第一个内节点 P 布置到对数分布律成立 的范围内，既配置到旺盛湍流区域。 3．第一个内节点上的与壁面区域的当量粘性系数ηt 及当量 导热系数λt 按下式确定： （3） （4） （5） 将y+,T+表达式代入（2）中令T=Tp,再联合由（4）可得第一 个内节点上与壁面相平行的温度满足对数分布律为： （6） 将y+,u+表达式代入（1）中令u=up,再联合（3）可得第一个内节 点上与壁面相平行的流速满足对数分布律为： 再由（3）（5）可得 第一个内节点P上 （7） 其中η为分子粘性系数 再由（4）（6）可得第一个内节点P上 其中λ为分子导热系数 说明：壁面函数法主要是要确定壁面上温度的当量导热系数λt 及流速的当量粘性系数ηt，从而确定出壁面切应力τw 和热流密度qw。 （8） 4．对第一个内节点P上Κp和 εp的确定方法的选择。 1） Κp值按Κ方程计算，边界条件取为 2） 令Гw=0 即 aw=0(p控制容积中，k的产生与耗散都较向壁面的扩散要大得多。） εp值不需求解离散方程，采用代数方程 （9） 5、标准 Κ-ε模型/壁面函数法数值计算中的注意事项 应满足 数值计算实践表明： ， 的值对传热特性的影响较大， 只有在一个合适的范围，计算结果与实验数据才能吻合的好。 1）第一个内节点与壁面间的无量纲距离 2）非线性耦合关系的处理 由于各个变量间强烈的非线性耦合关系，模拟过程中应当 采用亚松弛迭代方法，以利于非线性问题迭代的收敛。 Κ，ε及ηt 的松弛因子在开始时可取0.4-0.5左右。 3） Κ，ε方程中源项的处理方法 4）动量方程源项的离散（采用一系列插值） 先把Bous