麦1因变量中及其T、S、S,与a Table 1 The dependent variable and its I.,S,S,,a 因变量 ro S+ Sr 涡量 1:+=:{P,o} dP)亡+Cop⊙时词尺度) px( +P()×Pp+Px(pf 矢性流函数 0 1 -ip8-v3p310 p x书×P是+日 p 数性势函数 0 -Sm 脉动动能 1 B:if/prb G- 0 脉动动能耗散半 1 μeff/p,e C1+Gk /K-Caps/k 0 能量 1 Me ff/pr k 0 -∑TpA Q +C.TSm 气相组份 Y: 1 Hej f /pr s -1s a:∑npkm% 混合物分数 1 H:f i/pr f 0 0 于的脉动均方值 1 Heff/prg G&-C&2pge/k 0 颗粒与气相间的对流传热，ās为相变对气相S组份供献的分数。 中一四法在算法上不需采用交错网格，因而较之原型方程采用SIMPLE系列算法用了交 错网格，在计算三维流场时前者仅用1套网格，后者需用4套网格，且前者因不用交错网格 而减少了伪扩散。 2算 例 梁栋3,4)对国内研制的一种煤粉预燃室一扁平射流煤粉预燃室三维气相冷态场做了数 值模拟。图1示扁平射流煤粉预燃器。图2示燃烧室纵对称面上合速测量值，图3示计算结 果，两图的比例尺不一致，对比两者可以看出，它们在定量上有一定差异，但定性上很相 符。图4示计算所得速度分量V,在3个不同位置横截面上的分布。自图2~图4可分析对稳 图1平扁射流煤粉预燃器〔3〕 图2对称面上实测合速度分布，Y=0〔4们 Fig.1 The flat jct combustor(3) Fig.2 Measured resultant velocity distribution on the symmetrical plane,Y=04) 0?表 因变量价及其 ，、 ，、 ，，与 价 ，， ，， ，，， 。 因 变 量 价 。 厂 。 ， ， ， 涡 量 才 ， ， ， 尸 ， 。 争 尸 万 ， 。 矢性流函 数 了 忿 · 。 才 二 。 ， 了丽着灭正万 ‘ ‘ 知 厂 些 夕 厂 又 书卜 一 共卜 一了夕 夕 叻 一 夕 劝 夕 一 少 了 厂 尸 户梦 厂 一 一 十 今臼 。 。 、 。 粤 尸 数性势函 数 脉动动能 脉动动能 耗散率 能 量 一 子名， ， “ 。 ， ‘ 粼 ， 几 。 一 户。 一 户刃 ， 反 碑左几护 气相组份 混合物分数 的脉动均方值 ，‘ 。 厂 尸 ， 召 。 井 。 尹 ， 一 刃 。 户 、 ‘ ， ， 叮 刃介 ，‘ ，‘ 吞 一 夕。 命 夕一 颗粒氢气想间的对 流传热 ， “ 为相变对气相 组份 供献的分数 。 叻一。 法 在算法上不需采用 交错 网格 ， 因而较之原型方程采用 系列算法用 了交 错 网格 ， 在计算三维流场时前者仅用 套 网 格 ， 后者需用 套网格 ， 且前 者因不用 交错 网 格 而减少 了伪扩散 。 算 例 梁 栋 〔 ” ， ‘ ’ 对 国内研 制的 一种煤粉预燃室一扁平射流煤粉预燃室三维气相 冷态场 做 了数 值模拟 。 图 示扁 平射 流煤粉 预燃 器 。 图 示 燃烧室纵对称面上 合速测量值 ， 图 示计算结 果 ， 两 图的 比 例尺 不一致 ， 对 比 两 者可 以看出 ， 它们在定量上有一定差异 ， 但 定 性 上 很 相 符 。 图 示计算所得速度分量犷 二 在 个 不 同位置 横截面上的分 布 。 自图 一 图 可 分 析对稳 县 一 探之续止止 阵三云月 升 图 平扁 射 流煤 粉预燃 器〔 〕 · 冬 〔 〕 多琴 图 对称面上实侧合速度分布 ， 二 〔 〕 · 认 冬 “ ， 二 〔哆〕