湍流预溷火焰模型.ppt

设τ等于τ-的时间分数为α，则τ等于τ+的时间分数为(1 - α) ，根据平均值的定义 (3-62) (3-63) 因此 (3-64) (3-65) 或者按下述关系求解 α 、τ+和τ- 与 同τmax、 同τmin的相对大小有关 对于 和 对于 对于 根据α 、τ+和τ-，求出 和 求温度的均值、瞬时值和脉动均方根值 化学反应只是把滞止焓中的化学能转化成其中的内能和动能，在一定条件下滞止焓本身的脉动可以忽略 (3-66) (3-67) (3-68) (3-69) 湍流预混燃烧的速率的求解过程 →湍流反应度脉动模型 → T ? 式(3-68) → S 式(3-57) →拉切滑模型的燃烧速率表达式(3-54) 或(3-55) →引用式(3-48) 和(3-49) ，便完成了对 方程源项的模拟 。 2 控制方程的统一形式 座标系中控制方程的统一形式 实际求解的微分方程(3-71)  湍流模型常数 其他湍流模型常数值同表3-2。 3 求解程序 在矩形截面管道中心部分有一个直径为2.5毫米的圆柱形火焰稳定器（其轴线垂直于主流方向），该矩形截面在垂直于火焰稳定器轴线方向上的宽度是7.8厘米，来流是丙烷和空气的预混物。 1) 边界条件 原则：根据试验状态确定计算的边界条件。 为了把回流区排除在积分区域之外，积分区域自距离火焰稳定器（其直径为d）2.5d的下游开始。 在积分区域的开始截面上 速度、滞止焓和燃料质量分数的分布，可直接参照Howe等人实验数据给定。 k、ε和gτ，没有对应的实验数据，根据经验确定。 k、ε和gτ k取为局部时间平均速度百分之二的平方值（即假设湍流强度为2%）； 燃料质量分数脉动的均方值（计算gτ时需要）取为局部燃料质量分数均值的1%； ε的初值： (3-72) lm是局部的混合长度(Nikurades公式)。 积分区域的边界 积分区域的内边界 管道的对称平面 对称性，函数梯度为0的形式(所有待解函数的边界条件 ) 积分区域的外边界 管道壁 没有质量穿透，近似考虑成等温， 速度和三个湍流量（k、ε和gτ）在固体壁面上的值均认为等于0。 2) 壁面函数 前述湍流模型适用于完全湍流的区域。 在固体壁面附近，流体的输运特性变化很大，当十分靠近壁面时，会出现层流粘性作用大于湍流输运作用的情况。 方法： 壁面函数（wall-function method） 比较经济，且对于利用经验公式比较方便。 低雷诺数模型（low-Reynolds-number modelling method）。 壁面函数 壁面函数是一组特殊的关系式，它关联流体在近壁区域和壁面上的状态。 与边界条件的区别： 不是为了确定在壁面上的状态。 作用：体现在修正近壁区域各变量的有效交换系数或者给出变量的数值。 可以用分析或者经验的办法导出来。 壁面函数的确定 速度、滞止焓和燃料质量分数的壁面函数 由解控制近壁区域的Couette流动的常微分方程得到 k、ε和gτ的壁面函数 给出适于近壁区域的k、ε和gτ的控制方程相当困难 采用经验加假设的方法来确定 k、ε和gτ的壁面函数 认为： 在近壁区域，以下考虑近似符合实际 平行于壁面的速度分量的分布遵守对数分布律； 湍流量k、ε和gτ分别处于各自的局部平衡状态，即它们各自的产生率和耗散率相等 联系k和ε的湍流长度尺度（length scale）与到壁面的垂直距离成正比。 (3-73) (3-74) (3-75) τw表示壁面剪应力， K是实验得到的常数0.4， 下标p表示紧靠壁面的节点。 k、ε和gτ的壁面函数 4 计算结果分析 1) 火焰扩展角 平面管道中火焰稳定器后面的火焰面沿轴向的分布，它是各横截面上是最大温度脉动值所处位置的轨迹 根据这个轨迹，可以方便地计算出火焰扩展角火焰扩展角的绝对值大约为2度。与实验事实相符。 各种不同来流状态下火焰扩展角 火焰面的位置 虽然火焰扩展角相同，但火焰面的位置却有所不同。 来流的温度越低，则火焰面越接近轴线；大的来流速度使得火焰面靠近轴线。 2) 温度脉动和湍流动能的分布随轴向距离的变化 曲线1和2分别取自其与火焰稳定器的距离为90d和180d的截面。 每个曲线都有一个峰值；这个峰值的位置随着与火焰稳定器的距离的增加而逐渐靠近管道壁面； 温度脉动在管道中心部分比较平稳，在火焰表面达到极大值。 不同轴向位置的湍流动能的分布 剖面的形状各不相同；k的值开