新型外燃旋流热风炉内流动与传热过程数值计算

D0I:10.13374/i.issn1001-053x.2004.01.010 第26卷第1期 北京科技大学学报 VoL26 No.1 2004年2月 Jourral of University of Science and Technology Beijing Feb.2004 新型外燃旋流热风炉内流动与 传热过程数值计算 张书臣”郭鸿志》刘向军”程小虎”蔡漳平)赵树治)程克友) 1)北京科技大学机城工程学院，北京1000832)济钢集团公司技术中心，济南2500003)济钢集团公司第一炼铁厂，济南250000 摘要阐明中小型高炉的风温由1000℃左右提高到1250-1300℃的主要措施是将煤气与 空气双预热：提出应用高速烧嘴和新型外燃旋流热风炉来预热助燃空气，烧单一煤气可以实 现高炉1250~1300℃高风温.对炉内流动与传热过程进行的数值计算结果表明，在蓄热球床 表面没有偏流，气流及温度分布均匀，在中试炉上测量的烟气温度分布与计算结果基本吻合， 关键词外燃旋流热风炉：速度场：温度场：数值计算 分类号TF748.41 近10年来，发达国家的高炉风温已达1250℃ 蓄热球床截面烟气气流速度和温度分布均匀. 以上，正向1350℃迈进.目前我国热风炉风温平 4min左右烟气与空气阀门自动切换换向，将原 均仅为1000℃左右，两者相比差距甚大.提高高 热风炉所需空气预热到所需高温(7001000℃)， 炉风温可以大幅度降低焦比，提高喷煤量，提高 同时预热煤气，烧单一煤气，可实现高炉1250- 高炉利用系数.因此，一些大型钢铁企业正积极 1300℃高风温. 采取措施，努力提高高炉风温.济钢技术中心与 北京科技大学合作开发济钢第一炼铁厂高炉热 1数值模型、边界条件及计算方法 风炉1250-~1300℃高风温燃烧技术，将助燃空气 最高能预热到1000℃. 1.1数学模型 热风炉的技术关键是解决偏流问题，使气流 湍流模型采用考虑气体旋转的修正的k一e 分布均匀.传统的外燃式热风炉比内燃式热风炉 双方程模型，采用Boussinesq假设"： 最大的优越性是气流分布均匀.但是传统外燃式 网-=8腰}子uv.n6. 热风炉没有旋流，占地大，造价高，顶燃式热风炉 流动与传热过程控制方程如下： 实质上是新型内燃式热风炉，顶部燃烧产物整流 (1)质量守恒方程. 后进入蓄热体，气流分布均匀，但是由于没有采 0+7-om=0. 用高速烧嘴，热风炉顶部燃烧空间很大，蓄热体 (2)动量方程， 采用传统的耐火材料，因此炉体很大.国内小型 g之+V.om=-VpHV.r-v(径pk+pg, O 石球热风炉，实质是顶燃式热风炉，由于没有采 2 用高速烧嘴，火焰与石球上层接触，常常烧结.热 t=4[(VH(W门，4=utcP 风炉小型化的技术关键是采用高速烧嘴与高温 (3)能量方程， 空气燃烧技术，纯高炉煤气在高速烧嘴内部完全 ohW=7(像h0s 燃烧，高速烧嘴切向安装在新型外燃旋流热风炉 (4)k-σ方程 上部，高温烟气在炉上部空间高速旋转流动，这 -6划=-k.nrT4会-e 就避免了火焰烧结下部小球，同时避免偏流，使 -0a=子e-cpev.4v偿 收稿日期20020708 张书臣男，30岁，硕士

第 ‘ 卷 第 期 年 月 北 京 科 技 大 学 学 报 珍 。 恤面 血 ， 】 一 新型外燃旋流热风炉 内流动与 传热过程数值计算 张 书 臣 ” 郭鸿志 ” 刘 向 军 ” 程 小 虎 ” 蔡漳平 ” 赵树 治 ” 程克友 ” 北 京科技大学机械工 程 学 院 ， 北 京 济钢集 团公 司 技术 中心 ， 济南 济钢 集 团公 司第一炼铁厂 ， 济南 摘 要 阐明 中小 型 高炉 的风温 由 ℃ 左 右 提高 到 ℃ 的主 要 措施 是 将煤气 与 空气双预热 提 出应用 高速烧 嘴和 新型外燃旋流 热风 炉 来预热助燃空气 ， 烧 单一 煤气 可 以实 现 高炉 一 ℃ 高风 温 对炉 内流动 与传热 过程 进行 的数值计 算结果表 明 ， 在蓄热球床 表面 没 有偏流 ， 气流及温度分布均匀 在 中试炉 上测 量 的烟气温度分布与计算结果基本吻合 关键词 外燃 旋流 热 风 炉 速度场 温度场 数值计 算 分 类号 近 年来 ， 发达 国家 的高炉风温 己达 ℃ 以上 ， 正 向 ℃ 迈 进 目前 我 国热风 炉 风 温 平 均仅 为 ℃ 左 右 ， 两 者 相 比差 距 甚 大 提 高 高 炉 风温 可 以大 幅度 降低 焦 比 ， 提 高喷煤 量 ， 提 高 高炉 利用 系 数 因此 ， 一 些 大 型钢 铁 企 业 正 积 极 采取 措 施 ， 努 力提 高 高炉 风温 济钢 技 术 中心 与 北 京科技 大 学 合 作 开 发济钢 第 一 炼 铁 厂 高 炉 热 风 炉 一 ℃ 高风温 燃烧 技 术 ， 将助燃 空气 最 高能预 热 到 ℃ 热 风 炉 的技 术关键 是解 决偏 流 问题 ， 使气 流 分布 均匀 传统 的外 燃 式热 风 炉 比 内燃 式热 风 炉 最 大 的优越性 是气 流 分布均匀 但 是传 统外 燃 式 热风 炉 没有旋 流 ， 占地 大 ， 造价 高 顶燃 式 热风 炉 实质 上 是新 型 内燃 式热风 炉 ， 顶 部燃烧 产物整流 后 进 入 蓄 热 体 ， 气 流 分布 均 匀 但 是 由于 没 有采 用 高速 烧 嘴 ， 热 风 炉顶 部燃烧 空 间很 大 ， 蓄热 体 采用 传 统 的耐 火材 料 ， 因此 炉体 很 大 国 内小 型 石 球热 风 炉 ， 实质 是 顶 燃 式 热 风 炉 ， 由于 没 有采 用 高速烧 嘴 ， 火焰与石 球上层 接触 ， 常常烧结 热 风 炉 小 型化 的技术 关 键 是 采 用 高速 烧 嘴 与 高温 空气 燃烧技术 ， 纯 高炉煤气 在 高速 烧 嘴 内部完全 燃烧 ， 高速烧 嘴切 向安装在 新 型外燃旋 流 热 风 炉 上 部 ， 高温 烟 气 在 炉 上 部 空 间高速 旋 转流 动 ， 这 就 避 免 了火焰烧 结 下 部 小球 ， 同时避 免偏 流 ， 使 蓄 热 球 床 截 面 烟 气 气 流 速 度 和 温 度 分 布 均 匀 左 右 烟 气 与 空气 阀 门 自动 切 换 换 向 ， 将 原 热 风 炉 所 需 空气 预 热 到所 需 高温 小 ℃ ， 同 时预 热 煤气 ， 烧 单 一 煤气 ， 可 实现 高炉 一 ℃ 高风温 数值模型 、 边 界 条件及 计 算方法 数 学模 型 湍 流 模 型 采 用 考 虑 气 体 旋 转 的 修 正 的裕一￡ 双 方程 模 型 ， 采 用 假 设‘ 、 一 阵 会会 一扣 。 ， 、 ， · 流 动 与传 热 过 程 控 制 方程 如 下 一飞 质 量 守恒 方程 争 ， · 勿的 一 · 动 量 方 程 掣 十 · 。 。 一、 一 令 ， 一 ， · 〔‘甲约 ‘， 的 · ， ， 。 一 ， 祥 能量 方 程 ， · 勿冷，一 甲 · 告 十必 一 方 程 甲 · 勿玖 二 一 甲 · 勿纯 二 一 静 · 。 · 。 令 。 一 。 加 。 。 · 除叫 暗 ，· 收稿 日期 刁 刁 张 书臣 男 ， 岁 ， 硕士 DOI ：10．13374／j ．issn1001－053x．2004．01．010

VoL.26 No.1 张书臣等：新型外燃旋流热风炉内流动与传热过程数值计算 35 (cx-VVtcpe). 两个高速烧嘴出口在截面1上，截面1距蓄热球 其中，c1=1.44,c2=1.92,c=1.0,=1.0,,=1.3. 床上表面2.1m,截面2距蓄热床表面1.2m,截面 1.2边界条件母 3距蓄热球床上表面0.15m,截面4为通过圆心与 在小型高效旋流热风炉高速烧嘴出口，按具 高速烧嘴出口的纵剖面，高炉煤气的组分（质量 体工况给定4，y,w,T,取k=0.01(w+),0.5l, 分数)：C02,8.4%;02,0.1%;C0,22.5%;CH,0.1%; 特征尺寸I取为高速烧嘴出口的当量直径. H2,1.6%:N2,67.3%. 在固体壁面上，速度取无滑移条件，== 高速烧嘴入口 w0,在邻近固体壁面区，采用壁函数： 截面1： Y*=11.63 截面4 Tw= Bpctkiv Y11.63 截面2 高速烧嘴入口 In(EY) 其中，Y”=pc"ky,q=,k=Y2 截面3 P(vi-pr)" p=g受-1-（经”. 图1物理模型 其中，脚标p为邻壁节点，，为近壁节点速度，y,为 Fig.1 Physical model 节点到壁面距离.式中冯·卡门常数-0.4187， 2.2 结果与分析 壁面粗糙度E取9379，t,9w分别为壁面切应力及 (1)速度分布.图2~4分别为横截面1,2,3的速 热流通量. 度分布.从图2可以看出，截面1的速度分布很不 13计算方法 均匀，最大速度在两个高速烧嘴出口处附近，为 流动与传热过程控制方程的差分方程为： 119.55m/s,最小速度在圆心附近，为12.978m/s. C中p=Csp:+Cm中m十Cs中stCw中w+Cr中e+Cs中atD. 从图3可见，在截面2的圆周上切向速度已趋于 流场计算采用SIMPLER算法，柱坐标系.网 均匀，最大速度为22.91m/s.在圆的R2R的区 格划分采用18×18×18均匀网格，对流项的离散 域，切向速度都比较大，从0-R2区域，其切向速 化采用上风格式，扩散项的离散化采用中心差分 度与截面1的情况类似.从图4可见，在截面3的 格式，源项线性化，代数方程的求解采用亚松驰 圆周方向上，切向速度已趋向于均匀，但在圆心 Gauss--Seidel迭代法， 附近区域，速度很高，达63.42m/s.图5为过轴线 2计算结果与分析 与两高速烧嘴出口纵剖面轴向速度分布，从图中 可见，在图的下部，即靠近蓄热球上表面的R3~R 2.1物理模型 的区域，速度分布比较均匀，在0~R3的区域，速 图1为新型外燃旋流热风炉上部空间的物理 度较大，其中在轴线附近速度很大，可达107.7 模型，亦即数值计算的区域，内径2.2m,高2.4m. ms.图6为矢量速度分布，即总体速度分布图. G3-%w-1 。1 图2计算截面1速度分布 图3计算截面2速度分布 图4计算截面3速度分布 Fig.2 Velocity distribution on Fig.3 Velocity distribution on Fig.4 Velocity distribution on section 1 section 2 section 3

张 书 臣等 新型外 燃旋流 热风炉 内流 动 与传热过程 数值 计算 含 · 价 ‘ 沪。 其 中 ， ， 伪 ， 伪 ， 二 ， 氏 · 边 界 条件 ’卜 在 小型 高效旋流 热 风 炉 高速烧 嘴 出 口 ， 按 具 体工 况 给 定 ， ， ， ， 取 汁 ， 于 护叨 ， 特 征尺 寸 取 为高速烧 嘴 出 口 的 当量 直 径 在 固 体 壁 面 上 ， 速 度 取 无 滑 移 条 件 ， “ 赃 ， 在邻 近 固体壁 面 区 ， 采用壁 函数 两 个 高速 烧 嘴 出 口 在截 面 上 ， 截 面 距 蓄 热球 床 上 表 面 ， 截 面 距 蓄 热床 表 面 ， 截 面 距 蓄热球床 上表 面 ， 截面 为通过 圆心 与 高速 烧 嘴 出 口 的纵 剖 面 高炉 煤气 的组 分 质 量 分 数 认 ， ， ， 凡 ， 玫 ， 凡 ， 、 高速烧嘴大口 一 截面 口 爪 、 高速烧嘴入 口 截厂面︷ 二 一 隘 截面 其 中 ， 一 ， 。 ， 一 税 “ 东鹊 万 一 截面 介一 呼 一 今 一 口 口 其 中 ， 脚标尸为邻 壁 节 点 ， 咋 为近壁 节 点速度 ，外 为 节 点到壁 面距 离 式 中冯 · 卡 门常数介 ， 壁 面 粗糙度 取 ， ‘ ， 如 分 别 为壁 面 切 应 力及 热 流 通 量 计 算方 法 流 动 与 传 热 过 程控 制 方程 的差 分 方程 为‘ ， 沪 价沙 神砂 劝 刊叻砂 劝砂几叻 习七 流 场 计 算采 用 算 法 ， 柱 坐 标 系 网 格 划 分采 用 均 匀 网格 ， 对 流 项 的离散 化采 用 上 风 格 式 ， 扩 散项 的离散化采用 中心 差 分 格 式 ， 源 项 线 性 化 ， 代 数 方 程 的求 解采用 亚 松 驰 一 迭 代法 图 物 理模型 啥 妙 计 算结 果 与分 析 物理 模 型 图 为新型 外燃 旋流热 风 炉上 部 空 间 的物理 模 型 ， 亦 即数 值 计 算 的 区 域 ， 内径 ， 高 结 果 与分 析 速度 分布 图 分 别 为横截 面 ，， 的速 度分 布 从 图 可 以看 出 ， 截面 的速 度分 布很 不 均 匀 ， 最 大速 度在 两 个 高速烧 嘴 出 口 处 附近 ， 为 耐 ， 最 小速 度 在 圆心 附近 ， 为 而 从 图 可 见 ， 在 截 面 的 圆周 上 切 向速度 己趋 于 均 匀 ， 最 大速 度 为 公 在 圆 的招 的 区 域 ， 切 向速 度 都 比较 大 ， 从 区 域 ， 其 切 向速 度 与截面 的情 况类似 从 图 可 见 ， 在 截 面 的 圆周 方 向上 ， 切 向速 度 己趋 向于 均 匀 ， 但 在 圆心 附近 区域 ， 速 度 很 高 ， 达 岁 图 为过 轴线 与两 高速烧 嘴 出 口 纵 剖面轴 向速度分布 ， 从 图 中 可 见 ， 在 图 的下 部 ， 即靠 近 蓄 热 球 上表 面 的招 的 区 域 ， 速 度 分 布 比较 均 匀 ， 在 的 区 域 ， 速 度 较 大 ， 其 中在 轴 线 附近速 度 很 大 ， 可 达 州 图 为矢 量速 度 分 布 ， 即 总 体速 度 分 布 图 一二三一 洪蒸然﹃雏竹︸ 图 计 算截面 速度 分 布 图 计 算截面 速 度 分 布 · 图 计 算截面 速 度分 布 幻改

·36· 北京科技大学学报 2004年第1期 与测量结果的最大误差不超过5%.这是因为本 课题是采用的是未加修正的k一ε双方程模型.如 果想减小计算结果的误差，应使用修正的k一ε双 方程模型. 图5计算截面4速度分布 Fig.5 Velocity distribution on section 4 图8计算戴面4温度分布 Fig.8 Temperature distribution on section 4 860 840 。一剥量值 ◆一计算值 820 800 图6总体速度分布 -1.22-0.88-0.4400.440.881.22 Fig.6 General velocity distribution in the furnace 半径/m 图9测量与计算的截面3与截面4交界线上的温度分布 可以看出，这种流线图谱是比较理想的，气流分 Fig.9 Measured and calculated temperature distribution on 布比较均匀，避免了偏流 the cross line of section 3 and secton 4 (2)温度分布 图7为截面1的温度分布，可以看出，这个 3 结论 截面上的温度分布并不均匀，火焰的温度达到 1140℃，而轴线附近的温度仅792℃. (1)阐明占我国铁产量23的中小高炉的风 图8为截面4的温度分布.可以看出，在该图 温由1000℃左右提高到1250~1300℃主要措施 下部的温度分布比较均匀，在870900℃左右.图 是将煤气与空气双预热，同时提高原热风炉耐火 9为中试炉上用热电偶测量的与计算的截面3与 材料质量.提出应用高速烧嘴和高温空气燃烧技 截面4交界线上的温度分布.由图9可看出，计算 术研制新型外燃旋流热风炉来预热高温助燃 空气，烧单一高炉煤气，可以实现高炉1250- 1300℃高风温. (2)数值计算了新型外燃旋流热风炉内的速 度场与温度场，计算结果表明，在蓄热球床表面， 烟气没有偏流，烟气气流分布比较均匀，烟气温 度分布均匀.在中试炉上测量的温度分布与计算 基本吻合.数值计算的结果为新型外燃旋流热风 炉的设计与研制提供了参考依据， 参考文献 图7计算截面1温度分布 Fig.7 Temperature distribution on section 1 】郭鸿志，张欣欣，刘向军.传输过程数值模拟M.北

北 京 科 技 大 学 学 报 年 第 期 冷旧 ， 仲到 ， 之 功口习 甘乙翻 与测 量 结 果 的最 大 误 差 不 超 过 这 是 因为本 课题 是采用 的是 未加 修正 的 一 双 方程模型 如 果想 减 小计算 结果 的误 差 ， 应 使用 修正 的 一 双 方 程模 型 且口目曰剐 幽翔带姗溯川则沼彻，皿 湘枷洲﹄﹃翻一满艘 一三一 土 、 一 、勘 ， 几 图 计 算截面 速 度分 布 咭 伙 ，上 二 图 计 算截面 温 度分布 面 了…朴 侧、明咖 图 ‘ 总 体 速度分 布 落 】 曲州 可 以看 出 ， 这 种流 线 图谱 是 比较 理 想 的 ， 气 流 分 布 比较均匀 ， 避免 了偏流 温 度 分布 图 为截 面 的温度 分布 ， 可 以看 出 ， 这 个 截面 上 的温 度 分布 并 不 均 匀 ， 火 焰 的温 度 达 到 ℃ ， 而轴 线 附近 的温 度 仅 ℃ 图 为截面 的温度 分布 可 以看 出 ， 在 该 图 下部 的温度分布 比较均 匀 ， 在 ℃ 左 右 图 为 中试炉 上用 热 电偶 测 量 的与计 算 的截 面 与 截面 交界线 上 的温 度分布 由图 可看 出 ， 计 算 一 ， 刁 刁 科 科 半径 图 ，测 与计算的截面 与截面 交界线 上 的温 度分布 落 比 肠 邝 邝 曲 翻 灯 皿 即 ‘亡一 图 计 算截面 温 度分布 曲杭 。 结论 阐 明 占我 国铁产 量 的 中小 高炉 的风 温 由 ℃ 左 右 提 高到 ℃ 主 要 措 施 是将煤气 与空气 双预 热 ， 同时提 高原热风 炉耐火 材料质量 提 出应用 高速烧 嘴和 高温 空气燃烧技 术 研 制 新 型 外 燃 旋 流 热 风 炉 来 预 热 高 温 助 燃 空气 ， 烧 单 一 高 炉 煤气 ， 可 以实现 高炉 卜 ℃ 高风 温 数 值 计 算 了新 型 外 燃 旋 流 热 风 炉 内的速 度场 与温度 场 计算 结果表 明 ， 在 蓄热球床表面 ， 烟气 没 有偏 流 ， 烟 气 气 流 分布 比较均 匀 ， 烟 气 温 度 分布均匀 在 中试 炉 上测量 的温度分布 与计 算 基 本 吻合 数值 计算 的结果 为新型外燃 旋流热风 炉 的设计 与研 制提供 了参考 依 据 参 考 文 献 郭鸿志 ， 张欣欣 ，刘 向军 传输过 程数值模拟四」北 兹圈目赢曰 翩﹃翻︸州川姗脚初柳朋枷卿

Vol.26 No.1 张书臣等：新型外燃旋流热风炉内流动与传热过程数值计算 37· 京：冶金工业出版社，1998.256 9郭鸿志，姬超月，李有章，连铸中间包熔池流场矢 2帕坦卡.传热与流体流动的数值计算（张政译）[M) 性流函数-祸量法三维数值模拟[).化工冶金，1997， 北京：科学出版社，1980.188 18(5):68 3陶文铨.数值传热学M.西安：西安交通大学出版 10朱洁，郭鸿志，傅杰.DC-EAES-LF熔池电弧电渣传 社，995.255 热过程三维数值模拟[刀.北京科技大学学报，1997， 4郭鸿志.矢性流函数-涡量法充型过程数值模拟[门. 193):254 北京科技大学学报，1999,21(1)23 11冯妍卉.反向凝固器内带有固液相变的湍流流动与 5郭鸿志.75 t DCEAF熔池传输现象的仿真仞.北京 传热传质的数值模拟D].北京：北京科技大学，2000 科技大学学报，1998,20(2：126 14 6 Guo H Z.3D-Numerical Modelling of Electromagnetic 12 Guo HZ.Math modelling of electromagnetic flow Field Velocity and Temperature Fields in DC Electric Arc Elec- and heat transfer in DC electric arc furnace bath [J].J Univ tric Slag Ladle Fumaces[M].Beijing:Science Press,1997. Sci Technol Beijing,1995,2(2):133 165 3郭鸿志.DC-EAF熔池电磁搅拌与传热的数值模拟 7郭鸿志，朱洁，傅杰。直流电弧电渣钢包炉内电磁场 [J.北京科技大学学报，1997,19(3：284 速度场与温度场的三维数值模拟)化工治金，1997， 14刘向军.四角切向燃烧煤粉锅炉内燃烧过程的数值 18(2):162 模拟D].北京：清华大学，1997.12 8李晨曦.压力下真空密封铸造及矢性流函数-涡量 15李晨曦，郭鸿志.涡流旋风粉煤燃烧器的温度场与 法充型过程数值模拟D].北京：北京科技大学，1998. 速度场[J.北京科技大学学报，1995,19(5)：156 44 Numerical Calculation of Flow and Heat Transfer in an Outer-Combustion Swirl- Flowing Hot Blast Stove ZHANG Shuchen,GUO Hongzhi,LIU Xiangjun,CHENG Xiaohu,CAI Zhangping,ZHAO Shuzhi, CHENG Keyou 1)Mechanical Engineering School,University of Science and Technology Beijing,Beijing 100083,China 2)Technology Center,Jinan Steel Group,Jinan 250000,China 3)First Iron-Smelting Factory,Jinan Steel Group,Jinan 250000,China ABSTRACT It is clarified that the important method to improve the blast temperature from 1000C to 1250~1 300C is to preheat both auxiliary air and coal gas.The blast furnace temperature can be reached to 1250~ 1 300'C by burning single blast stove gas if using high speed burners and a type of external-combustion swirl-flo- wing hot blast stove to preheat auxiliary combustion air.The calculating results,similar to those determinating ex- perimental data in a middle experiment hot blast stove,proved that there is no eccentric flow in the face of the store energy balls and the flow and temperature fields even. KEY WORDS external-combustion swirl-flowing hot blast furnace;flow field;temperature field;numerical cal- culation

匕 张书 臣 等 新型外 燃旋流 热风炉 内流 动 与传热 过 程数值计算 京 冶 金 工 业 出版社 ， 帕坦 卡 传热 与流 体流动 的数值计算 张政 译 』 北京 科学 出版社 ， 陶文 锉 数值传 热 学 西 安 西 安交通 大学 出版 社 ， 郭鸿 志 矢性 流 函数一 涡 量 法 充型 过 程 数值模拟 明 北 京科技 大学 学报 ， ， 郭鸿志 熔池传 输现 象 的仿真闭 北 京 科技大学学报 ， ， ” 了七 ， 郭鸿 志 ， 朱 洁 ， 傅杰 直 流 电弧 电渣钢 包 炉 内电磁 场 速度场与温度 场 的三维数值模拟 化 工冶金 ， ， 李晨 曦 压 力下 真 空 密 封铸造 及矢 性 流 函数一 涡 量 法 充型 过程数值模拟 北 京 北 京科技 大 学 ， 郭鸿 志 ， 姬超 月 ， 李 有章 连铸 中间包熔池 流 场 矢 性流 函 数一 涡 量法 三 维数值模拟 闭 化工 冶金 ， ， 朱洁 ， 郭鸿志 ， 傅杰 一 熔池 电弧 电渣 传 热 过程 三 维数值模拟 阴 北 京科 技 大 学学报 ， ， 冯 妍 卉 反 向凝 固器 内带 有 固液相 变 的湍 流流动与 传热传 质 的数值模拟 田 北京 北 京科技大学 ， ， ， 郭鸿 志 一 熔 池 电磁 搅拌 与 传 热 的数值模拟 田 北京科技 大 学学报 ， ， 刘 向军 四角切 向燃烧煤粉锅 炉 内燃烧过程 的数值 模拟 田 北 京 清华大 学 ， 李晨 曦 ， 郭鸿 志 涡 流旋 风粉煤燃烧器 的温度场 与 速度 场 北 京科技大学学报 ， ， 一 月只 ，，， 〔 ，，， 柳盯 ，，， 月五刃 如 ，，， 月了 爬护动犷 ， 乙阮咬 产，， 〔 万王刃 砂 ，， 阴 ， 盯 ， ， 毛戈助。 ， ， 的 的 一 ， ， ， 而 犷 一 一 一 一 ， 以 ， 一 一