北大考研-工学院研究生导师简介-郑一

清华大学经济管理学院硕士生导师简介-郑捷郑捷经济系助理教授办公室：伟伦楼561凯程教育是五道口金融学院和清华经管考研黄埔军校,在2014年,凯程学员考入五道口金融学院28人,清华经管11人,五道口状元武xy出自凯程, 在2013年,凯程学员考入五道口金融学院29人,清华经管5人,状元李少h出在凯程, 在凯程网站有很多凯程学员成功经验视频,大家随时可以去查看. 2016年五道口金融学院和清华经管考研保录班开始报名!个人简介研究成果研究项目郑捷，清华大学经济管理学院经济系助理教授。

2003年获得清华大学经济学学士学位，同年获北京市高校优秀毕业生、清华大学优秀毕业生称号。

2005年获得清华大学经济学硕士学位。

2007年获得美国华盛顿大学经济学硕士学位，同年获华盛顿大学经济系最佳助教奖、最佳第三年论文奖。

2011年获得华盛顿大学经济学博士学位。

研究领域包括微观经济理论、产业组织、数理经济学、金融经济学、实验经济学。

当前研究主要为理性泡沫、健忘与福利、自我控制、网络效应、软件产业发展以及开放源代码现象的经济学研究。

曾在世界计量经济学会亚洲年会(AMES2011)、世界计量经济学会北美夏季年会(NASMES2011)、世界计量经济学会第10届全球大会(ESWC2010)、国际经济学会第16届全球大会(IEA2011)、英国皇家经济学会年会（RES2011）、第23届博弈论国际会议（ICGT2012）、美国中西部宏观经济学会议（MMM2010)、美国中西部经济理论与贸易会议（MWIEGSpring2009,Fall2010）、中国留美经济学会中国年会（CES2010,2011,2013）、中国博弈论与实验经济学学会首届会议(CGTEEA2010)、中国博弈论与应用国际会议（CICGTA2012）等国际会议上宣讲学术论文。

2011年与钟笑寒教授、连暐虹教授共同建立清华大学经济管理学院经济科学与政策实验室（ESPEL），并担任实验室管理员至今。

北京理工大学计算机学院导师简介及研究方向计算机技术方向丁刚毅，北京大学理学学士，北京理工大学工学博士。

北京理工大学后备中青年学术带头人。

95年至96年在美国俄亥俄大学做高级访问学者。

2002年3月任北京理工大学软件学院副院长，负责科研、实验室、对外合作等工作。

指导硕士生、博士生十余名，发表论文40余篇。

承担科研项目十几项，主要涉及武器系统仿真设计，高性能计算，软件测试等，十五期间负责的科研、开发项目总经费800余万元。

曾获部级二等奖二项，校级一级科研及论文奖三项。

学术科研专长为计算机仿真、软件工程、机电工程分析、武器效应仿真评估、分布式并行计算等，被聘为总装备部仿真技术专业组成员、中国计算机仿真协会副理事长、《计算机仿真》杂志编委。

地址：北京市海淀区中关村南大街5号北京理工大学软件学院502室邮编：100081电话：010-******** 68913182 68912741传真：************EMAIL：***********.cn胡昌振胡昌振，男，1967年9月生，教授.博导，研究方向：信息安全。

联系方式：68914977-23李凤霞，女，教授，北京市教学名师。

研究方向包括虚拟现实技术、计算机仿真技术、运动模拟技术等，曾获部级二等奖1项，部级三等奖5项，先后主持和参与国家级项目十余项。

目前正承担的项目包括国防科工委基础研究项目1项，总装重点研究项目4项。

联系方式：************;135****1901;***********************.cn李小平教授性别：男学历：博士职务：北京理工大学远程教育学院副院长电话：68915623 手机：130****9599办公地点逸夫楼继续教育学院407室刘玉树教授、博导程成副教授主要学习经历: 美国UIUC大学访问学者。

中科院软件所获得博士学位，专业：计算机应用技术，研究方向：虚拟环境人机交互技术;山东大学计算机学院获得硕士学位，专业：计算机应用技术，研究方向：图形学和辅助设计;南开大学计算机与系统科学系获得学士学位，专业：控制论。

北大郑也夫教授作品目录（全）(2013-10-19 00:28:07)转载▼一、专著2013，《吾国教育病理》，北京：中信出版社。

2013，《语镜子》，（待出）2011，《半开放社会》，香港：文化中国出版有限公司。

2011，《沙葬》，上海 : 复旦大学出版社。

2009，《神似祖先》，北京 : 中国青年出版社。

2008，《与本科生谈：论文与治学》，济南：山东人民出版社。

2007，《后物欲时代的来临》，上海 : 上海人民出版社。

2007，《抵抗通吃》，济南 : 山东人民出版社。

2004，《知识分子研究》，北京 : 中国青年出版社。

2004，《阅读生物学札记》，北京 : 中国青年出版社。

2004，《被动吸烟者说》，北京 : 中国青年出版社。

2003，《中国社会中的信任》（与彭泗清等合著），北京 : 中国城市出版社。

2002，《城市社会学》，北京 : 中国城市出版社。

2001，《信任论》，北京 : 中国广播电视出版社。

1998，《忘却的纪念》，北京 : 中央编译出版社。

1998，《世界杯断想》，北京 : 今日中国出版社。

1997，《游戏人生》，海口 : 海南出版社。

1995，《走出囚徒困境》，北京 : 光明日报出版社。

1995，《代价论》，北京 : 三联书店。

1993，《礼语咒词官腔黑话——语言社会学丛谈》，北京 : 光明日报出版社。

1993，《命运守护神: 九集电视片文本》，北京 : 中国广播电视出版社。

1991，《中国足球的出路》，北京 : 知识出版社。

二、编著2013年，《科场现形记》，（待出）2013，《双城记：京沪众生素描》（与陈映芳合编），上海：上海交通大学出版社。

2012，《北大清华人大社会学硕士论文选编》（与沈原、潘绥铭合编），北京 : 中国社会科学出版社。

2011，《北大清华人大社会学硕士论文选编》（与沈原、潘绥铭合编），北京 : 中国青年出版社。

2010，《北大清华人大社会学硕士论文选编》（与沈原、潘绥铭合编），北京 : 中国城市出版社。

郑义—中国商业组织力第一导师
中国组织力研究专家
企业组织Z型管理模式创始人
北京大学民营经济研究院学部专家
企业营销组织管理专家
中国员工心理健康工程特约专家
中华讲师网品牌讲师
前沿讲座特邀讲师
中国讲师网著名讲师
历任：
天津交易所财经委员
华易期货总经理
加拿大投资资询有限公司总经理
天泽衡器有限责任公司投资人
影响力教育训练集团济南公司总经理
影响力教育训练集团执行顾问（2006年进入影响力集团）
资历背景：
有20年的市场营销、企业管理经验，任山东省经信委，山东省品牌发展研究会专家及专家委员会委员。

山东大学管理学院特聘教授，济南大学教授、研究生导师。

擅长企业团队组织管理教练，营销主管经理人的组织打造及组织力提升训练，辅导企业组织建设，组织力提升，实现组织管理团队三倍速成长。

大型项目和案例：
1、担任十期影响力集团PMBA导师，成功培养山东能源、立白、海信、鲁能、宏济堂制药、济南青年汽车、保龄宝生物、中国“姜老大”、中国邮政、益方中联等700余家企业中高管理干部和职业经理人。

2、参加中国八届培训论坛，并连续4届成为“人力资源战略与开发国际论坛”主讲嘉宾。

案例1山东青岛金岭电器设备出口欧美企业，连续五年不盈利，郑老师执行顾问进入后三个月扭亏为盈实现利润730万
案例2济南PMBA学员企业顶好商贸请郑老师为其全国服装经销商组织培训，实现现场订单销售增长3倍；
案例3另一家学员企业香港轩美诗当年实现培训成果展现经销商销售倍增业绩。

案例4顾问企业山东能源集团2012年实现营业额1600亿，2013年8个月（顾问期间）突破2000亿，具世界500强460位。

第２０卷第２期材　料　与　冶　金　学　报Ｖｏｌ ２０Ｎｏ ２　收稿日期：２０２１ ０１ ０８．　基金项目：国家重点研发计划项目（２０１７ＹＦＢ０３０４０００）．　作者简介：郑占一（１９９５—），男，硕士研究生，Ｅ ｍａｉｌ：ｚｈｅｎｇｚｙ＠ｓｔｕｍａｉｌ ｎｅｕ ｅｄｕ ｃｎ．　通讯作者：齐凤升（１９８０─），男，副教授，Ｅ ｍａｉｌ：ｑｉｆｓ＠ｍａｉｌ ｎｅｕ ｅｄｕ ｃｎ．２０２１年６月ＪｏｕｒｎａｌｏｆＭａｔｅｒｉａｌｓａｎｄＭｅｔａｌｌｕｒｇｙＪｕｎｅ２０２１ｄｏｉ：１０ １４１８６／ｊ ｃｎｋｉ １６７１－６６２０ ２０２１ ０２ ００２转底炉内冶金粉尘还原过程数值模拟郑占一，齐凤升，刘中秋，李宝宽（东北大学冶金学院，沈阳１１０８１９）摘　要：基于计算流体力学方法并以收缩核模型为基础建立了转底炉内燃烧、烟气流动、气体与冶金粉尘球团传热传质及冶金粉尘球团化学反应的全耦合数学模型，计算了中径３６ｍ的转底炉内流场、温度场及冶金粉尘球团内铁氧化物的还原反应，重点分析了球团内部各种铁氧化物浓度及球团的金属化率．采用文献中球团在高温硅钼炉内进行的还原实验验证了模型的可靠性．结果表明，在本文工况下，经过一个工作周期（２５ｍｉｎ），炉膛内烟气流速随流动方向逐渐增大，转底炉中径处球团温度为１４１６ ７Ｋ，铁的浓度由３４７７ ５０ｍｏｌ／ｍ３增长至９７１９ ９４ｍｏｌ／ｍ３，冶金粉尘球团的金属化率最高可达９０ ８５％，平均金属化率为８１ ４２％．关键词：转底炉；冶金粉尘球团；收缩核模型；直接还原；金属化率中图分类号：ＴＦ０６２ 文献标识码：Ａ 文章编号：１６７１ ６６２０（２０２１）０２ ００８５ ０７ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔｒｅｄｕｃｔｉｏｎｐｒｏｃｅｓｓｉｎｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅＺｈｅｎｇＺｈａｎｙｉ，ＱｉＦｅｎｇｓｈｅｎｇ，ＬｉｕＺｈｏｎｇｑｉｕ，ＬｉＢａｏｋｕａｎ（ＳｃｈｏｏｌｏｆＭｅｔａｌｌｕｒｇｙ，ＮｏｒｔｈｅａｓｔｅｒｎＵｎｉｖｅｒｓｉｔｙ，Ｓｈｅｎｙａｎｇ１１０８１９，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ａｆｕｌｌｙｃｏｕｐｌｅｄｍａｔｈｅｍａｔｉｃａｌｍｏｄｅｌｏｆｃｏｍｂｕｓｔｉｏｎ，ｇａｓｆｌｏｗ，ｈｅａｔａｎｄｍａｓｓｔｒａｎｓｆｅｒ，ｃｈｅｍｉｃａｌｒｅａｃｔｉｏｎｉｎｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔｐｅｌｌｅｔｓｗａｓｅｓｔａｂｌｉｓｈｅｄｂａｓｅｄｏｎｃｏｍｐｕｔａｔｉｏｎａｌｆｌｕｉｄｄｙｎａｍｉｃｓｍｅｔｈｏｄａｎｄｓｈｒｉｎｋｉｎｇｃｏｒｅｍｏｄｅｌ．Ｔｈｅｆｌｏｗａｎｄｔｅｍｐｅｒａｔｕｒｅｆｉｅｌｄ，ｒｅｄｕｃｔｉｏｎｒｅａｃｔｉｏｎｏｆｉｒｏｎｏｘｉｄｅｏｆａ３６ｍｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅｗｅｒｅｃａｌｃｕｌａｔｅｄｂｙｔｈｉｓｍａｔｈｅｍａｔｉｃａｌｍｏｄｅ．Ｔｈｅｍｏｌａｒｉｔｙｏｆｉｒｏｎｏｘｉｄｅｓａｎｄｔｈｅｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｒａｔｅｏｆｃｏｍｐｏｓｉｔｅｐｅｌｌｅｔｓｗｅｒｅａｎａｌｙｚｅｄ．Ｔｈｅｒｅｓｕｌｔｓｏｆｒｅｄｕｃｔｉｏｎｅｘｐｅｒｉｍｅｎｔｉｎｔｈｅｌｉｔｅｒａｔｕｒｅｗａｓｃａｒｒｉｅｄｏｕｔｔｏｖｅｒｉｆｙｔｈｅｒｅｌｉａｂｉｌｉｔｙｏｆｔｈｅｍｏｄｅｌ．Ｄｕｒｉｎｇｏｎｅｗｏｒｋｉｎｇｃｙｃｌｅ（２５ｍｉｎ），ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｔｈａｔＴｈｅｖｅｌｏｃｉｔｙｏｆｇａｓｉｎｔｈｅｆｕｒｎａｃｅｉｎｃｒｅａｓｅｄｇｒａｄｕａｌｌｙｗｉｔｈｔｈｅｆｌｏｗｄｉｒｅｃｔｉｏｎ．Ｔｈｅｔｅｍｐｅｒａｔｕｒｅｏｆｔｈｅｐｅｌｌｅｔｓａｔｔｈｅｍｉｄｄｌｅｄｉａｍｅｔｅｒｏｆｔｈｅｒｏｔａｒｙｈｅｒａｔｈｆｕｒｎａｃｅｗａｓ１４１６ ７Ｋ，ａｎｄｔｈｅｉｒｏｎｍｏｌａｒｉｔｙｉｎｃｒｅａｓｅｓｆｒｏｍ３４７７ ５０ｍｏｌ／ｍ３ｔｏ９７１９ ９４ｍｏｌ／ｍ３，ｔｈｅｈｉｇｈｅｓｔｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｒａｔｅｏｆｔｈｅｐｅｌｌｅｔｓｗａｓ９０ ８５％，ａｎｄｔｈｅａｖｅｒａｇｅｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｒａｔｅｗａｓ８１ ４２％．Ｋｅｙｗｏｒｄｓ：ｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ；ｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔｐｅｌｌｅｔｓ；ｓｈｒｉｎｋｉｎｇｃｏｒｅｍｏｄｅｌ；ｄｉｒｅｃｔｒｅｄｕｃｔｉｏｎ；ｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｒａｔｅ 钢铁行业是我国经济的支柱性产业，其生产过程会产生大量的冶金粉尘，产生量约为粗钢产量的８％～１２％［１－２］．２０２０年我国钢铁行业粗钢产量为１０ ６５亿ｔ，冶金粉尘产量至少为８５１８万ｔ．钢铁企业冶金粉尘的含铁量（质量分数）一般在３０％～７０％［３－４］，还含有ＺｎＯ，Ｐｂ，ＫＣｌ，ＮａＣｌ等成分．转底炉十余年来从加热炉转变为冶炼设备，既可用于铁精矿的煤基直接还原，又可处理钢铁企业的冶金粉尘［５］，逐渐成为处理冶金粉尘的主要设备．转底炉还原冶金粉尘的工作过程涉及炉底球团直接还原、炉内烟气流动、传热传质、煤气燃烧等复杂过程，因此对转底炉工作过程的研究十分困难．一些学者进行了转底炉数学模型的研究，主要是转底炉热平衡计算和炉内状态模拟［６－７］及对转底炉的加热制度和加热设备的模拟计算［８－９］．刘颖等［１０－１１］以球团为研究对象，建立了转底炉还原冶金粉尘球团过程一维非稳态数学模型，研究了影响球团金属化率的主要因素，按重要程度排序依次为：炉膛温度＞球团直径＞反应时间＞碳氧比．Ｗｕ等［１２－１３］建立了转底炉直接还原过程的集成模型，将转底炉的三维ＣＦＤ模型与球团内部直接还原的一维模型进行迭代，描述金属氧化物的还原过程．Ｄａｓｇｕｐｔａ等［１４］在转底炉还原球团矿的数学模型中将单球团模型扩展为多层球团模型，给出了时间－温度和时间－温度－化学吸热等值线，以及多床层系统产生的净热流和一氧化碳产生量．这些对球团的研究模型能够反映球团内部组分的化学反应状况及浓度变化，但缺少球团化学反应与转底炉内部过程的耦合计算，不能反映球团在转底炉各个位置的状态．本文采用数值模拟方法建立了转底炉内燃烧与冶金粉尘球团中铁氧化物还原的全耦合数学模型，分析了冶金粉尘球团在随炉底转动过程中的温度变化，以及金属氧化物浓度、金属化率等参数．该数学模型解决了冶金粉尘球团运动与炉膛加热的传热传质问题，以及转底炉中的冶金粉尘球团中铁氧化物的还原问题，为转底炉工业应用提供理论指导．１　数学模型１ １　几何模型根据实际尺寸建立转底炉几何模型，如图１所示，转底炉中径为３６ｍ，炉宽５ ２７ｍ，炉高１ ６１５ｍ．烧嘴布置在距炉底０ ８０７５ｍ处，内侧布置烧嘴２６个，外侧布置烧嘴３８个，各区域角度及出口、入口如图１（ａ）所示．对计算区域进行网格划分，考虑计算量、计算速度和时间成本，经网格无关性验证，确定网格数量为１５０万个，炉膛上方燃烧区域为非结构网格，炉底料层区域为结构化网格，如图１（ｂ）所示．１ ２　控制方程１ ２ １　基本控制方程在转底炉工作中伴随着燃烧、传热传质及化学反应等过程，这些物理化学变化在转底炉工作过程中相互作用．转底炉内部烟气流动、传热传质及化学反应过程满足质量、动量及能量守恒．各个过程的守恒方程如下：连续性方程：ρｔ＋ ·（ρ珒ν）＝０（１）动量方程：ｔ（ρ珒ν）＋ ·（ρ珒ν　珒ν）＝!"#烧嘴还原一区还原四区还原二区还原二区均热区预热区布料排料区烟气出口物料出口物料入口烧嘴还原三区$%&' $%&'$%&''(&($)&(*(&( *+&'排烟区,-.图１　转底炉几何模型及网格划分Ｆｉｇ １　ＧｅｏｍｅｔｒｉｃｍｏｄｅｌａｎｄｇｒｉｄｓｏｆＲＨＦ（ａ）—几何模型；（ｂ）—网格划分－ ｐ＋·μ 珒ν＋ 珒ν()Ｔ－２３ ·珒ν[]Ｉ（２）能量方程：ｔρ()Ｅ＋ ·珒νρＥ＋()[]ｐ＝ ·ｋｅｆｆ Ｔ－∑ｉｈｉ珒Ｊ()ｉ（３）式（１）～（３）中，ρ为密度，ｋｇ／ｍ３；ｔ为时间，ｓ；珒ν为速度矢量，ｍ／ｓ；ｐ为压力，Ｐａ；μ为黏度，Ｐａ·ｓ；Ｉ为单位张量；Ｅ为总能量，Ｊ／ｋｇ；Ｔ为温度，Ｋ；ｋｅｆｆ为有效传热系数，Ｗ／（ｍ·Ｋ）；ｈｉ为显焓，Ｊ／ｋｇ；Ｊｉ为扩散通量，ｋｇ／（ｍ３·ｓ）．１ ２ ２　湍流模型对于转底炉内烟气的湍流流动，采用标准ｋ－ε模型．湍动能ｋ和耗散率ε的控制方程为：ｔρ()ｋ＋ ·ρｋ珒()ν＝·μ＋μｔσ()ｋ[]ｋ＋６８材料与冶金学报 第２０卷Ｇｋ＋Ｇｂ－ρε－ＹＭ（４）ｔ（ρε）＋ ·（ρε珒ν）＝·μ＋μｔσ()ε[]ｋ＋Ｃ１εεｋＧｋ－Ｃ２ερε２ｋ（５）式（４）～（５）中，ｋ为湍动能，ｍ２／ｓ２；μｔ为湍流黏度，Ｐａ·ｓ；ε表示湍动能耗散率；Ｇｋ表示速度梯度产生的湍动能，Ｊ／（ｍ３·ｓ）；Ｇｂ表示浮力产生的湍动能，Ｊ／（ｍ３·ｓ）；ＹＭ表示波动和扩张对总耗散率的影响；模型常数分别为Ｃ１ε＝１ ４４，Ｃ２ε＝１ ９２，σｋ＝１ ０，σε＝１ ３．１ ２ ３　燃烧模型组分输运模型是通过求解混合物中各个组分的对流、扩散和反应确定的守恒方程，可以描述化学物质的混合和传输过程．本文使用组分输运模型模拟各组分的质量分数：ｔ（ρＹｉ）＋ ·（ρ珒νＹｉ）＝－ ·珒Ｊｉ＋Ｒｉ（６）燃烧模型采用基于涡耗散模型的湍流－化学相互作用模型，反应产物的净生成率由式（７）和式（８）计算结果的最小值表示：Ｒｉ，ｒ＝ｖ′ｉ，ｒｍｉＡρεｋｍｉｎｗＲｖ′Ｒ，ｒｍ()Ｒ（７）Ｒｉ，ｒ＝ｖ′ｉ，ｒｍｉＡＢρεｋ∑ＰｗＰ∑Ｎｊｖ″ｊ，ｒｍｊ（８）式（６）～（８）中，Ｙｉ为组分ｉ的质量分数；ｗＰ为生成物组分的质量分数；Ｒｉ为化学反应源项，ｋｇ／（ｍ３·ｓ）；Ａ和Ｂ为经验系数，Ａ＝４，Ｂ＝０ ５；ｖ′ｉ，ｒ为反应物的化学计量数；ｖ″ｊ，ｒ为生成物的化学计量数；ｍｉ为反应物ｉ的分子质量；ｍｊ为生成物ｊ的分子质量；ｗＲ为任一反应物的质量分数；ｍＲ为任一反应物的分子质量．１ ２ ４　辐射模型离散坐标辐射模型求解范围涵盖整个光学深度，有较高的精确度，且适用于滑移网格的计算，可表示为：·（Ｉ（珒ｒ，珒ｓ）珒ｓ）＋（ａ＋σｓ）Ｉ（珒ｒ，珒ｓ）＝ａｎ２σＴ４π＋σｓ４π∫４π０Ｉ（珒ｒ，珒ｓ′）Φ（珒ｓ·珒ｓ′）ｄΩ′（９）式（９）中，珒ｒ为位置向量；珒ｓ为方向向量；珒ｓ′为散射方向矢量；ａ为吸收系数；ｎ为折射率；σｓ为散射系数；σ为Ｓｔｅｆａｎ Ｂｏｌｔｚｍａｎｎ常数，５ ６７×１０－８Ｗ／（ｍ２·Ｋ４）；Ｉ为辐射强度，Ｗ／ｓｒ；Φ为相函数；Ω′为立体角，ｓｒ．１ ２ ５　多孔介质模型转底炉炉底为一层冶金粉尘球团，料层空间被流体与固体混合物占据，并随着炉底转动．本研究将料层假设为一层多孔介质，用以描述流体在料层区域流动时产生的压降，同时将多孔介质区域的温度作为冶金粉尘球团反应前沿面温度．通过源项的方式实现化学反应过程中球团与炉内烟气的传热传质．多孔介质区域控制方程如下：动量方程：（γρｆ珒ν）ｔ＋ ·（γρｆ珒ν珒ν）＝－γ ｐ＋ ·（γτ）－μα珒ν＋Ｃ２１２ρ｜珒ν｜珒()ν（１０）能量方程：ｔγρｆＥｆ＋（１－γ）ρｓＥ[]ｓ＋ ·珒ν（ρｆＥｆ＋ｐ[]）＝ ·ｋｅｆｆ Ｔ－（∑ｉｈｉＪｉ）＋（τ·珒ν[]）＋∑ｊＲｊ·ΔＨｊ（１１）式（１０）～（１１）中，γ为多孔介质的孔隙率；ρｆ和ρｓ分别为流体和固体的密度，ｋｇ／ｍ３；Ｅｆ和Ｅｓ为流体与固体的能量，Ｊ／ｋｇ；Ｒｊ为化学反应速率，ｍｏｌ／（ｍ３·ｓ）；ΔＨｊ为化学反应焓变，Ｊ／ｍｏｌ．１ ２ ６　收缩核模型对于冶金粉尘球团内部铁金属氧化物的还原，真正的还原剂为固体碳．固体碳直接还原铁氧化物可以看作铁氧化物的一氧化碳间接还原反应和碳气化反应的加和，铁氧化物的还原遵循Ｆｅ２Ｏ３→Ｆｅ３Ｏ４→ＦｅＯ→Ｆｅ的逐级还原规律．本文模型中考虑的化学反应如下：碳的气化反应：Ｃ＋ＣＯ２＝２ＣＯ铁氧化物的还原：３Ｆｅ２Ｏ３＋ＣＯ＝２Ｆｅ３Ｏ４＋ＣＯ２Ｆｅ３Ｏ４＋４ＣＯ＝３Ｆｅ＋４ＣＯ２，Ｔ＜８４３ＫＦｅ３Ｏ４＋ＣＯ＝３ＦｅＯ＋ＣＯ２，Ｔ＞８４３ＫＦｅＯ＋ＣＯ＝Ｆｅ＋ＣＯ２在生产中，将冶金粉尘球团布置在转底炉炉底，球团在炉底转动过程中接受烟气与炉壁的辐射热量，温度升高，然后在热力学条件允许时发生一系列的化学反应．本模型中使用收缩核模型描述球团内部进行的铁金属氧化物的还原反应，以７８第２期 郑占一等：转底炉内冶金粉尘还原过程数值模拟气固相反应动力学模型计算球团化学反应速率［１５］，通过自定义标量输运方程的形式与炉膛内的控制方程进行耦合求解，计算球团中各组分收缩核半径．控制方程为：ｍｊ ｔ＝－ｋ·π·ｄ２ｊ·Ｍｊ·Ｒｓ，ｊ可简化为：ｒｊｔ＝－ｋ·Ｍｊρｊ·Ｒｓ，ｊ（１２）式中，ｒｊ为各组分收缩核半径，ｍ；Ｍｊ为各组分摩尔质量，ｋｇ／ｍｏｌ；ρｊ为各组分密度，ｋｇ／ｍ３；Ｒｓ，ｊ为各化学反应界面反应速率，ｍｏｌ／（ｍ２·ｓ）．碳的气化反应速率为：Ｒｓ，Ｃ＝ｋＣｅ－ＥＣＲＴρＣ（ｐＣＯ２－ｐｅｑＣＯ２）（１３）铁氧化物还原反应速率为：Ｒｓ，ＦｅｘＯｙ＝ｋＦｅｘＯｙｅ－ＥＦｅｘＯｙＲＴρＦｅｘＯｙ（ｐＣＯ－ｐＦｅｘＯｙ，ｅｑＣＯ）（１４）式（１３）～（１４）中：Ｒｓ，Ｃ为碳气化反应速率，ｍｏｌ／（ｍ２·ｓ）；ｋＣ为碳气化反应指前因子，ｍｏｌ／（ｍ·ｋｇ·ｓ·Ｐａ）；ρＣ为碳的质量浓度，ｋｇ／ｍ３；ｐＣＯ２为反应体系中ＣＯ２分压，Ｐａ；ｐｅｑＣＯ２为碳气化反应达到平衡时ＣＯ２分压，Ｐａ；Ｒｓ，ＦｅｘＯｙ为各铁氧化物还原反应速率，ｍｏｌ／（ｍ２·ｓ）；ｋＦｅｘＯｙ为铁氧化物还原反应指前因子，ｍｏｌ／（ｍ·ｋｇ·ｓ·Ｐａ）；ρＦｅｘＯｙ为铁氧化物的质量浓度，ｋｇ／ｍ３；ｐＣＯ为反应体系中ＣＯ分压，Ｐａ；ｐＦｅｘＯｙ，ｅｑＣＯ为碳气化反应达到平衡时ＣＯ分压，Ｐａ；ＥＣ为碳气化反应表观活化能，Ｊ／ｍｏｌ；ＥＦｅｘＯｙ为铁氧化物还原反应表观活化能，Ｊ／ｍｏｌ；Ｒ为理想气体常数，８ ３１４Ｊ／（ｍｏｌ·Ｋ）．１ ３　边界条件燃料和助燃气体入口为烧嘴出口，形状分别为圆形和与该圆同心的圆环．入口类型为速度入口，入口速度由气体流量折算．表１为流量２７０００ｍ３／ｈ、预热温度５２３Ｋ下的燃料成分．助燃空气流量为９５００ｍ３／ｈ，富氧用氧气流量为３０００ｍ３／ｈ，预热温度为７７３Ｋ．烟气出口类型为压力出口，转底炉各壁面为恒定温度，炉顶为１２０℃，炉墙及炉底为９０℃．多种冶金粉尘与黏结剂通过配比后混合，通过造球机制作成冶金粉尘球团．球团的主要成分如表２所示，本模型中将冶金粉尘球团假设为半径８ｍｍ的圆球团，球团进入转底炉前的温度为３１０Ｋ．通过滑移网格方法实现冶金粉尘球团随炉底在炉内的转动，转动速度为０ ００３５ｒａｄ／ｓ．表１燃料成分（体积分数）Ｔａｂｌｅ１　Ｆｕｅｌｃｏｍｐｏｓｉｔｉｏｎ（ｖｏｌｕｍｅｆｒａｃｔｉｏｎ）％　ＣＯＯ２ＣＯ２Ｈ２Ｎ２２３ ３５００ ８２９２９ １７５１ ４９７４５ １４９表２　冶金粉尘球团主要成分（质量分数）Ｔａｂｌｅ２　Ｃｈｅｍｉｃａｌｃｏｍｐｏｓｉｔｉｏｎｏｆｔｈｅｐｅｌｌｅｔ（ｍａｓｓｆｒａｃｔｉｏｎ）％　ＴＦｅＦｅ２Ｏ３ＦｅＯＭＦｅＣａＯＭｇＯＣＺｎ其他４２ ８０２７ １２１２ ７３１３ ９１９ ０１１ ８０１２ ９５１ ９１２０ ５７２　模型验证采用文献［１２］中的实验数据对模型进行验证，将４，１０，１６，２２ｍｉｎ时球团金属化率的模拟值与实验值进行对照．金属化率为转底炉还原冶金粉尘球团的一个主要的技术指标，其计算公式为：η＝ＭＦｅＴＦｅ×１００％（１５）表３为转底炉工作不同时间后球团金属化率的模拟和测量结果．由于初始阶段的球团成分不同，转底炉工作１０ｍｉｎ内球团金属化率模拟和测量结果相差较大．１０ｍｉｎ后模拟的转底炉状态接近实验状态，可用于模型验证．最终预测误差在７ ８２％以内（一般误差在１０％以内被认为准确性较好），验证了数学模型的可靠性．表３　球团金属化率Ｔａｂｌｅ３　Ｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｒａｔｅｏｆｔｈｅｐｅｌｌｅｔｓｔ／ｍｉｎ实验值模拟值相对误差％４７ ２３２ １—１０３７ ９３５ ７－６ １６１６６３ ４５８ ８－７ ８２２２７７ ６８１ ５４ ７８８８材料与冶金学报 第２０卷３　结果与讨论３ １　转底炉内流场和温度场分布特征图２为转底炉运行一个周期（２５ｍｉｎ）后，炉膛烧嘴处（距炉底０ ８０７５ｍ）速度场矢量图．结果表明，燃气与助燃气体以恒定速度经烧嘴喷入炉膛内部，炉膛内烟气逆时针流向物料入口，烟气流速随流动方向逐渐增大，最后从烟气出口流出．图３为炉底冶金粉尘球团的温度分布图．温度为３１０Ｋ的冶金粉尘球团随着炉底的转动进入转底炉内．冶金粉尘球团顺时针运动接受烟气与炉壁的辐射热量，温度升高．在转底炉中径处球团升温最快，靠近转底炉内侧及外侧墙壁的球团则升温较缓．冶金粉尘球团在出口处被加热至１４１６ ７Ｋ，在出口处靠近内侧及外侧墙壁的球团温度则介于１２５０～１３００Ｋ之间．３ ２　铁氧化物的还原图４为转底炉运行一个周期后冶金粉尘球团内部各种铁氧化物的收缩核半径云图，从中可以看出当冶金粉尘球团达到临界反应温度后，铁氧化物以Ｆｅ２Ｏ３→Ｆｅ３Ｏ４→ＦｅＯ的顺序逐级进行还原反应．图４（ａ）表明Ｆｅ２Ｏ３的收缩核半径在满足Ｆｅ２Ｏ３还原反应的条件后迅速减小，这是由于Ｆｅ２Ｏ３的还原反应所需的热力学和动力学条件较为简单，因此在转底炉中径处球团中Ｆｅ２Ｏ３的收缩核半径在４５０ｓ内减小至０．达到临界温度８４３Ｋ后，Ｆｅ３Ｏ４与ＦｅＯ均参与反应，收缩核半径开始减小．由于ＣＯ还原ＦｅＯ需要较高的热力学及动力学条件，从图４（ｂ）和（ｃ）中可以看出，Ｆｅ３Ｏ４的收缩核半径减小较快，ＦｅＯ的收缩核半径在还原过程中减小得较为缓慢，两者均未完全反应．!"#$%!&#%'!!#'($#$)*+%$'#'&"#$,&#&%!#''-./0图２　转底炉烧嘴处速度矢量图Ｆｉｇ ２　ＶｅｌｏｃｉｔｙｖｅｃｔｏｒｄｉａｇｒａｍａｔｂｕｒｎｅｒｏｆＲＨＦ!"#"$%&'"((&%$)*+),-$.""&/-))""..-).".0*-1)/10-1.,$(-()&,1-(.&.)-0)$0$-*/)1&-*11$,-..(,/-.1(./-//图３　转底炉内冶金粉尘球团温度分布ｇ ３　ＴｅｍｐｅｒａｔｕｒｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｐｅｌｌｅｔｓｉｎＲＨＦ!"#""$""#""%&"#""'("#""')"#""*%"#""*+"#""&'"#""&""#"")&"#"",("#"",)"#""+%"#""++"#"""'"-./-01-21图４　转底炉内铁氧化物收缩核半径Ｆｉｇ ４　ＲａｄｉｕｓｏｆｉｒｏｎｏｘｉｄｅｕｎｒｅａｃｔｅｄｃｏｒｅｉｎＲＨＦ（ａ）—Ｆｅ２Ｏ３；（ｂ）—Ｆｅ３Ｏ４；（ｃ）—ＦｅＯ ９８第２期 郑占一等：转底炉内冶金粉尘还原过程数值模拟 图５为冶金粉尘球团在转底炉工作一个周期后，球团内铁氧化物及铁的浓度云图．由图中可以看出，反应发生后Ｆｅ２Ｏ３的浓度减小，Ｆｅ的浓度在还原区域不断增大，而Ｆｅ３Ｏ４和ＦｅＯ在还原区域由于存在相互转化，所以浓度先上升后下降．同时，冶金粉尘球团温度分布的不均匀导致在转底炉径向上的铁氧化物的浓度分布不均匀．图６为转底炉运行一个周期后中径处的铁氧化物浓度变化曲线图．从图中可以看出，炉底在运动至距转底炉入口约６５°时，Ｆｅ２Ｏ３开始反应，运行至距入口１５０°时Ｆｅ２Ｏ３的反应基本完成，其浓度由２３７３ｍｏｌ／ｍ３减少到０．Ｆｅ３Ｏ４和ＦｅＯ的浓度是一个先升高后降低的过程，炉底运动至距入口约１４５°时Ｆｅ３Ｏ４的浓度升高至最大，为１０９２ １８ｍｏｌ／ｍ３．ＦｅＯ在距入口约１００°时开始富集增多，在距入口１８０°时浓度达到最大，为３９８１ ８０ｍｏｌ／ｍ３．转底炉中ＦｅＯ的浓度最大为４６１２ ９５ｍｏｌ／ｍ３，出现在靠近转底炉的侧壁处，这是由于在侧壁处的温度较低，ＦｅＯ的反应速率较小，造成ＦｅＯ的富集时间较长，富集量较大．!"#!$#%&#%'#()*)+,,((,)+-,(,).+,,/01.+-,/12-+,,/-(-+-,/)-1+,,//01+-,/,/*+,,0.*+-,1*0+,,-,0+-,))2+,,/12+-,,34563).1/(+2-.).2+.(.,0-+02)0((+)-)--0+0()(2-+(2),)/+*/(*10+()(-+.+12((./+/1/2**+1)/*/.+/,/.-,+-2//0*+,.2()+-,34563)2*()+022(**+*(00)/+--0)0-+)0*2)2+(/*.2)+,.*,.1+0*11,,+121/-.+-(-*,0+)--(1(+/0.0/1+,/.)12+0.)2()+1*).**+-,34563)/,2(+/0/,/.+/12)1+/-0-0+/.*0,+/)*,(+//)(.+/,-.1+,2.10+,0)2,+,1)/(+,-().+,./-1+,)*0+,/,34563)图５　转底炉内铁及铁氧化物的摩尔浓度Ｆｉｇ ５　ＭｏｌａｒｉｔｙｏｆｉｒｏｎａｎｄｉｒｏｎｏｘｉｄｅｉｎＲＨＦ（ａ）—Ｆｅ２Ｏ３；（ｂ）—Ｆｅ３Ｏ４；（ｃ）—ＦｅＯ；（ｄ）—Ｆｅ 图７展示了转底炉内冶金粉尘球团在炉底中径处铁的浓度及球团的金属化率．结果表明，铁氧化物在进入还原区域、经过逐级反应后，浓度不断增大，在反应后期浓度增长放缓．这是由于在反应过程中收缩核半径不断减小，反应界面的面积不断减小，使反应放缓．在转底炉工作一个周期后，转底炉中径处的冶金粉尘球团中铁的浓度由３４７７ ５０ｍｏｌ／ｍ３增长至９７１９ ９４ｍｏｌ／ｍ３；同时，在转底炉的中径处球团金属化率由３２ ５０％增大至９０ ８５％．图８展示了转底炉工作一个周期后，出口处的冶金粉尘球团的金属化率．结果表明，在出口中心处球团的金属化率最大，为９０ ８５％；中０９材料与冶金学报 第２０卷心两侧的球团金属化率逐渐减小，出口处的冶金粉尘球团平均金属化率为８１ ４２％．!"#$% !"%$& !"$'(')''#('#''*(+%''角度,- .&(++&+++%(++%+++#(++#+++*(++*+++(+++ ,-/01 /2%.图６　转底炉中径处铁氧化物浓度变化曲线Ｆｉｇ ６　ＭｏｌａｒｉｔｙｏｆｉｒｏｎｏｘｉｄｅｉｎｔｈｅｍｉｄｄｌｅｄｉａｍｅｔｅｒｏｆＲＨＦ!"#$% #&'(角度!" )*+,++-*+-++,*+'++浓度球团金属化率,++.+/+0+1+*+2+'+球团金属化率!3,++++.+++/+++1+++*+++2+++'+++图７　转底炉中径处铁的浓度及球团金属化率Ｆｉｇ ７　ＭｏｌａｒｉｔｙｏｆｉｒｏｎａｎｄｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｉｎｔｈｅｍｉｄｄｌｅｄｉａｍｅｔｅｒｏｆＲＨＦ球团金属化率!"#$#%#&'&%$距离!(&'')*)'+*+',*,'-*图８　转底炉出口处球团金属化率Ｆｉｇ ８　ＴｈｅｉｒｏｎｍｅｔａｌｌｉｚａｔｉｏｎｏｆｔｈｅｐｅｌｌｅｔｓａｔｔｈｅｏｕｔｌｅｔｏｆＲＨＦ４　结　论（１）在本文工况下的转底炉中径处，Ｆｅ２Ｏ３的浓度在４５０ｓ内由２３７３ｍｏｌ／ｍ３降低至０，Ｆｅ３Ｏ４和ＦｅＯ的浓度则是先升高后降低，Ｆｅ３Ｏ４的浓度在达到最大值１０９２ １８ｍｏｌ／ｍ３后也迅速减小，ＦｅＯ经历了１０４７ｓ的反应后浓度为９２３ ５０ｍｏｌ／ｍ３．　（２）球团直径对球团还原的影响分为两个阶段，在７５０～１２５０ｓ的反应阶段含碳球团直径较大，金属化率升高得较快；在１２５０ｓ之后，球团金属化率随着球团直径的减小而升高得缓慢．（３）冶金粉尘球团在转底炉内经过一个周期（２５ｍｉｎ）的工作过程后，金属化率最高达９０ ８５％，转底炉出口处的平均金属化率为８１ ４２％．参考文献：［１］佘雪峰，薛庆国，王静松，等．钢铁厂含锌粉尘综合利用及相关处理工艺比较［Ｊ］．炼铁，２０１０，２９（４）：５６－６２．（ＳｈｅＸｕｅｆｅｎｇ，ＸｕｅＱｉｎｇｇｕｏ，ＷａｎｇＪｉｎｇｓｏｎｇ，ｅｔａｌ．Ｃｏｍｐｒｅｈｅｎｓｉｖｅｕｔｉｌｉｚａｔｉｏｎａｎｄｒｅｌａｔｉｖｅｔｒｅａｔｍｅｎｔｏｆｚｉｎｃ ｃｏｎｔａｉｎｉｎｇｄｕｓｔｉｎｉｒｏｎａｎｄｓｔｅｅｌｗｏｒｋｓ［Ｊ］．Ｉｒｏｎｍａｋｉｎｇ，２０１０，２９（４）：５６－６２．）［２］陈砚雄，冯万静．钢铁企业粉尘的综合处理与利用［Ｊ］．烧结球团，２００５，３０（５）：４２－４６．（ＣｈｅｎＹａｎｘｉｏｎｇ，ＦｅｎｇＷａｎｊｉｎｇ．Ｏｎｔｈｅｃｅｎｔｒａｌｉｚｅｄｔｒｅａｔｍｅｎｔａｎｄｃｏｍｐｒｅｈｅｎｓｉｖｅｕｔｉｌｉｚａｔｉｏｎｏｆｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔ［Ｊ］．ＳｉｎｔｅｒｉｎｇａｎｄＰｅｌｌｅｔｉｚｉｎｇ，２００５，３０（５）：４２－４６．）［３］ＣａｎｔａｒｉｎｏＭＶ，ＦｉｌｈｏＣＤＣ，ＭａｎｓｕｒＭＢ．Ｓｅｌｅｃｔｉｖｅｒｅｍｏｖａｌｏｆｚｉｎｃｆｒｏｍｂａｓｉｃｏｘｙｇｅｎｆｕｒｎａｃｅｓｌｕｄｇｅｓ［Ｊ］．Ｈｙｄｒｏｍｅｔａｌｌｕｒｇｙ，２０１２，１１１／１１２：１２４－１２８．［４］ＳｅｎｋＤ，ＧｕｄｅｎａｕＨＷ，ＧｅｉｍｅｒＳ，ｅｔａｌ．Ｄｕｓｔｉｎｊｅｃｔｉｏｎｉｎｉｒｏｎａｎｄｓｔｅｅｌｍｅｔａｌｌｕｒｇｙ［Ｊ］．ＩＳＩＪＩｎｔｅｒｎａｔｉｏｎａｌ，２００６，４６（１２）：１７４５－１７５１．［５］熊华文，戴彦德．转底炉直接还原技术对钢铁行业资源综合利用的意义及发展前景分析［Ｊ］．中国能源，２０１２，３４（２）：５－７，１３．（ＸｉｏｎｇＨｕａｗｅｎ，ＤａｉＹａｎｄｅ．Ｓｉｇｎｉｆｉｃａｎｃｅｆｏｒｒｅｓｏｕｒｃｅｃｏｍｐｒｅｈｅｎｓｉｖｅｕｔｉｌｉｚａｔｉｏｎｏｆｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅｄｉｒｅｃｔｒｅｄｕｃｔｉｏｎｉｎｓｔｅｅｌｉｎｄｕｓｔｒｙａｎｄａｎａｌｙｓｉｓｏｆｉｔｓｄｅｖｅｌｏｐｉｎｇｐｒｏｓｐｅｃｔｓ［Ｊ］．ＥｎｅｒｇｙｏｆＣｈｉｎａ，２０１２，３４（２）：５－７，１３．）［６］徐萌．转底炉煤基热风熔融炼铁工艺的基础性研究［Ｄ］．北京：北京科技大学，２００６．（ＸｕＭｅｎｇ．Ｆｕｎｄａｍｅｎｔａｌｒｅｓｅａｒｃｈｏｎｃｏａｌｈｏｔ－ａｉｒｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅｐｒｏｃｅｓｓ［Ｄ］．Ｂｅｉｊｉｎｇ：ＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＢｅｉｊｉｎｇ，２００６．）［７］高金涛，周春芳，朱荣，等．转底炉分区域供热研究［Ｊ］．北京科技大学学报，２０１４，３６（Ｓ１）：１１０－１１６．（ＧａｏＪｉｎｔａｏ，ＺｈｏｕＣｈｕｎｆａｎｇ，ＺｈｕＲｏｎｇ，ｅｔａｌ．Ｒｅｓｅａｒｃｈｏｎｔｈｅｈｅａｔｓｕｐｐｌｙｏｆｄｉｆｆｅｒｅｎｔｓｅｃｔｉｏｎｓｉｎａｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ［Ｊ］．ＪｏｕｒｎａｌｏｆＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＢｅｉｊｉｎｇ，２０１４，３６（Ｓ１）：１１０－１１６．）（下转第９６页）１９第２期 郑占一等：转底炉内冶金粉尘还原过程数值模拟二者呈线性关系．（２）较大的气化剂流速会影响流体和颗粒间的换热时间，恶化换热效果．对于气固顺流式移动床应当合理地控制气化剂流速．（３）出口气体中ＣＯ的质量分数随着气化剂流速的增加而降低，减小焦炭直径有助于加快气化反应的速率，出口气体中ＣＯ的质量分数随着焦炭直径的减小而增大．参考文献：［１］毛艳丽，曲余玲，王涿．高炉熔渣处理及显热回收工艺的研究进展［Ｊ］．上海金属，２０１３，３５（３）：４５－５０．（ＭａｏＹａｎｌｉ，ＱｕＹｕｌｉｎｇ，ＷａｎｇＺｈｕｏ．Ｒｅｖｉｅｗｏｆｂｌａｓｔｆｕｒｎａｃｅｍｏｌｔｅｎｓｌａｇｔｒｅａｔｍｅｎｔａｎｄｓｅｎｓｉｂｌｅｈｅａｔｒｅｃｏｖｅｒｙｔｅｃｈｎｏｌｏｇｉｅｓ［Ｊ］．ＳｈａｎｇｈａｉＭｅｔａｌｓ，２０１３，３５（３）：４５－５０．）［２］ＬｉＰ，ＱｉｎＱ，ＹｕＱＢ，ｅｔａｌ．Ｆｅａｓｉｂｉｌｉｔｙｓｔｕｄｙｆｏｒｔｈｅｓｙｓｔｅｍｏｆｃｏａｌｇａｓｉｆｉｃａｔｉｏｎｂｙｍｏｌｔｅｎｂｌａｓｔｆｕｒｎａｃｅｓｌａｇ［Ｊ］．ＡｄｖａｎｃｅｄＭａｔｅｒｉａｌｓＲｅｓｅａｒｃｈ，２０１０，９７／１０１：２３４７－２３５１．［３］杨世亮．流化床内稠密气固两相流动机理的ＣＦＤ－ＤＥＭ耦合研究［Ｄ］．杭州：浙江大学，２０１４．（ＹａｎｇＳｈｉｌｉａｎｇ．ＣＦＤ ＤＥＭｃｏｕｐｌｉｎｇｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｄｅｎｓｅｔｗｏ ｐｈａｓｅｆｌｏｗｍｅｃｈａｎｉｓｍｓｉｎｆｌｕｉｄｉｚｅｄｂｅｄｓ［Ｄ］．Ｈａｎｇｚｈｏｕ：ＺｈｅｊｉａｎｇＵｎｉｖｅｒｓｉｔｙ，２０１４．）［４］ＹａｎＬＢ，ＣａｏＹ，ＺｈｏｕＨ，ｅｔａｌ．Ｉｎｖｅｓｔｉｇａｔｉｏｎｏｎｂｉｏｍａｓｓｓｔｅａｍｇａｓｉｆｉｃａｔｉｏｎｉｎａｄｕａｌｆｌｕｉｄｉｚｅｄｂｅｄｒｅａｃｔｏｒｗｉｔｈｔｈｅｇｒａｎｕｌａｒｋｉｎｅｔｉｃｔｈｅｏｒｙ［Ｊ］．ＢｉｏｒｅｓｏｕｒｃｅＴｅｃｈｎｏｌｏｇｙ，２０１８，２６９：３８４－３９２．［５］ＬｉｕＤＹ，ＣｈｅｎＸＰ，ＺｈｏｕＷ，ｅｔａｌ．ＳｉｍｕｌａｔｉｏｎｏｆｃｈａｒａｎｄｐｒｏｐａｎｅｃｏｍｂｕｓｔｉｏｎｉｎａｆｌｕｉｄｉｚｅｄｂｅｄｂｙｅｘｔｅｎｄｉｎｇＤＥＭ ＣＦＤａｐｐｒｏａｃｈ［Ｊ］．ＰｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅＣｏｍｂｕｓｔｉｏｎＩｎｓｔｉｔｕｔｅ，２０１１，３３（２）：２７０１－２７０８．［６］ＫｕＸＫ，ＬｉＴ，Ｌ ｖ ｓＴ．ＣＦＤ ＤＥＭｓｉｍｕｌａｔｉｏｎｏｆｂｉｏｍａｓｓｇａｓｉｆｉｃａｔｉｏｎｗｉｔｈｓｔｅａｍｉｎａｆｌｕｉｄｉｚｅｄｂｅｄｒｅａｃｔｏｒ［Ｊ］．ＣｈｅｍｉｃａｌＥｎｇｉｎｅｅｒｉｎｇＳｃｉｅｎｃｅ，２０１５，１２２：２７０－２８３．［７］ＬｉｕＭＬ，ＣｈｅｎＭ，ＬｉＴＪ，ｅｔａｌ．ＣＦＤ－ＤＥＭ－ＣＶＤｍｕｌｔｉ ｐｈｙｓｉｃａｌｆｉｅｌｄｃｏｕｐｌｉｎｇｍｏｄｅｌｆｏｒｓｉｍｕｌａｔｉｎｇｐａｒｔｉｃｌｅｃｏａｔｉｎｇｐｒｏｃｅｓｓｉｎｓｐｏｕｔｂｅｄ［Ｊ］．Ｐａｒｔｉｃｕｏｌｏｇｙ，２０１９，４２：６７－７８．［８］ＬｉｕＤＹ，ＢｕＣＳ，ＣｈｅｎＸＰ．ＤｅｖｅｌｏｐｍｅｎｔａｎｄｔｅｓｔｏｆＣＦＤ ＤＥＭｍｏｄｅｌｆｏｒｃｏｍｐｌｅｘｇｅｏｍｅｔｒｙ：ａｃｏｕｐｌｉｎｇａｌｇｏｒｉｔｈｍｆｏｒＦｌｕｅｎｔａｎｄＤＥＭ［Ｊ］．Ｃｏｍｐｕｔｅｒｓ＆ＣｈｅｍｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ，２０１３，５８：２６０－２６８．［９］ＣｒｏｗｅＣＴ，ＳｏｍｍｅｒｆｅｌｄＭ，ＴｓｕｊｉＹ．Ｍｕｌｔｉｐｈａｓｅｆｌｏｗｓｗｉｔｈｄｒｏｐｌｅｔｓａｎｄｐａｒｔｉｃｌｅｓ［Ｍ］．Ｆｌｏｒｉｄａ：ＣＲＣＰｒｅｓｓ，１９９８．［１０］ＧｉｄａｓｐｏｗＤ．Ｍｕｌｔｉｐｈａｓｅｆｌｏｗａｎｄｆｌｕｉｄｉｚａｔｉｏｎ：ｃｏｎｔｉｎｕｕｍａｎｄｋｉｎｅｔｉｃｔｈｅｏｒｙｄｅｓｃｒｉｐｔｉｏｎ［Ｊ］．ＪｏｕｒｎａｌｏｆＮｏｎ ＮｅｗｔｏｎｉａｎＦｌｕｉｄＭｅｃｈａｎｉｃｓ，１９９４，５５（２）：２０７－２０８．［１１］王帅．流化床内稠密气固两相反应流的欧拉 拉格朗日数值模拟研究［Ｄ］．杭州：浙江大学，２０１９．（ＷａｎｇＳｈｕａｉ．Ｅｕｌｅｒｉａｎ Ｌａｇｒａｎｇｉａｎｓｉｍｕｌａｔｉｏｎｏｆｄｅｎｓｅｒｅａｃｔｉｖｅｇａｓ ｓｏｌｉｄｆｌｏｗｓｉｎｆｌｕｉｄｉｚｅｄｂｅｄｓ［Ｄ］．Ｈａｎｇｚｈｏｕ：ＺｈｅｊｉａｎｇＵｎｉｖｅｒｓｉｔｙ，２０１９．）（上接第９１页）［８］ＬａｎｄｆａｈｒｅｒＭ，ＳｃｈｌｕｃｋｎｅｒＣ，ＰｒｉｅｌｅｒＲ，ｅｔａｌ．ＤｅｖｅｌｏｐｍｅｎｔａｎｄａｐｐｌｉｃａｔｉｏｎｏｆａｎｕｍｅｒｉｃａｌｌｙｅｆｆｉｃｉｅｎｔｍｏｄｅｌｄｅｓｃｒｉｂｉｎｇａｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅｕｓｉｎｇＣＦＤ［Ｊ］．Ｅｎｅｒｇｙ，２０１９，１８０：７９－８９．［９］赵凯，宫晓然，胡长庆，等．转底炉用蓄热式烧嘴的模拟［Ｊ］．材料与冶金学报，２０１５，１４（２）：１２１－１２５．（ＺｈａｏＫａｉ，ＧｏｎｇＸｉａｏｒａｎ，ＨｕＣｈａｎｇｑｉｎｇ，ｅｔａｌ．Ａｓｉｍｕｌａｔｉｏｎｆｏｒｎｏｚｚｌｅｏｆｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ［Ｊ］．ＪｏｕｒｎａｌｏｆＭａｔｅｒｉａｌｓａｎｄＭｅｔａｌｌｕｒｇｙ，２０１５，１４（２）：１２１－１２５．）［１０］刘颖．转底炉内冶金粉尘含碳球团直接还原过程数学模型研究［Ｄ］．北京：北京科技大学，２０１５．（ＬｉｕＹｉｎｇ．Ｍａｔｈｅｍａｔｉｃａｌｍｏｄｅｌｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｄｉｒｅｃｔｒｅｄｕｃｔｉｏｎｏｆｃａｒｂｏｎ ｃｏｎｔａｉｎｉｎｇｐｅｌｌｅｔｓｍａｄｅｏｆｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔｉｎａｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ［Ｄ］．Ｂｅｉｊｉｎｇ：ＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＢｅｉｊｉｎｇ，２０１５．）［１１］ＬｉｕＹ，ＳｕＦ，ＷｅｎＺ，ｅｔａｌ．ＣＦＤｍｏｄｅｌｉｎｇｏｆｆｌｏｗ，ｔｅｍｐｅｒａｔｕｒｅ，ａｎｄｃｏｎｃｅｎｔｒａｔｉｏｎｆｉｅｌｄｓｉｎａｐｉｌｏｔ ｓｃａｌｅｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ［Ｊ］．ＭｅｔａｌｌｕｒｇｉｃａｌａｎｄＭａｔｅｒｉａｌｓＴｒａｎｓａｃｔｉｏｎｓＢ，２０１４，４５（１）：２５１－２６１．［１２］ＷｕＹＬ，ＪｉａｎｇＺＹ，ＺｈａｎｇＸＸ，ｅｔａｌ．Ｍｏｄｅｌｉｎｇｏｆｔｈｅｒｍｏｃｈｅｍｉｃａｌｂｅｈａｖｉｏｒｉｎａｎｉｎｄｕｓｔｒｉａｌ ｓｃａｌｅｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅｆｏｒｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔｒｅｃｙｃｌｉｎｇ［Ｊ］．ＭｅｔａｌｌｕｒｇｉｃａｌａｎｄＭａｔｅｒｉａｌｓＴｒａｎｓａｃｔｉｏｎｓＢ，２０１７，４８（５）：２４０３－２４１８．［１３］ＷｕＹＬ，ＪｉａｎｇＺＹ，ＺｈａｎｇＸＸ，ｅｔａｌ．Ｐｒｏｃｅｓｓｏｐｔｉｍｉｚａｔｉｏｎｏｆｍｅｔａｌｌｕｒｇｉｃａｌｄｕｓｔｒｅｃｙｃｌｉｎｇｂｙｄｉｒｅｃｔｒｅｄｕｃｔｉｏｎｉｎｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ［Ｊ］．ＰｏｗｄｅｒＴｅｃｈｎｏｌｏｇｙ，２０１８，３２６：１０１－１１３．［１４］ＤａｓｇｕｐｔａＳ，ＳａｌｅｅｍＳ，ＳｒｉｒａｎｇａｍＰ，ｅｔａｌ．Ａｃｏｍｐｕｔａｔｉｏｎａｌｓｔｕｄｙｏｎｔｈｅｒｅｄｕｃｔｉｏｎｂｅｈａｖｉｏｒｏｆｉｒｏｎｏｒｅ／ｃａｒｂｏｎｃｏｍｐｏｓｉｔｅｐｅｌｌｅｔｓｉｎｂｏｔｈｓｉｎｇｌｅａｎｄｍｕｌｔｉ ｌａｙｅｒｂｅｄｒｏｔａｒｙｈｅａｒｔｈｆｕｒｎａｃｅ［Ｊ］．ＭｅｔａｌｌｕｒｇｉｃａｌａｎｄＭａｔｅｒｉａｌｓＴｒａｎｓａｃｔｉｏｎｓＢ，２０２０，５１（２）：８１８－８２６．［１５］华一新．冶金过程动力学导论［Ｍ］．北京：冶金工业出版社，２００４：１６２－１６５．（ＨｕａＹｉｘｉｎ．Ｉｎｔｒｏｄｕｃｔｉｏｎｔｏｋｉｎｅｔｉｃｓｏｆｍｅｔａｌｌｕｒｇｙｐｒｏｃｅｓｓ［Ｍ］．Ｂｅｉｊｉｎｇ：ＭｅｔａｌｌｕｒｇｉｃａｌＩｎｄｕｓｔｒｙＰｒｅｓｓ，２００４：１６２－１６５．）６９材料与冶金学报 第２０卷。

782018年4月总第287期ISSN1672-1438CN11-4994/T 工程教育专业认证体系下的自动化专业课程建设研究胡哲源 张焕龙 姜利英 郑安平郑州轻工业学院电气信息工程学院 河南郑州 450002摘 要：在全国深入开展工程教育专业认证的背景下，自动化专业存在适应工程教育专业认证的需求和相关课程建设问题。

以单片机与嵌入式系统课程为例，通过研究工程教育专业认证的需求特点，结合郑州轻工业学院自动化专业人才培养的现状，从培养目标、培养方案、教学内容、教学方法等方面探讨面向工程教育专业认证体系的课程建设思路。

关键词：工程教育认证；自动化；单片机；课程建设作者简介：胡哲源，工学硕士，助教；通讯作者：张焕龙，工学博士，讲师；姜利英，工学博士，教授；郑安平，工学硕士，教授。

基金项目：郑州轻工业学院第四批青年教师教学改革与研究项目“工程教育专业认证体系下的自动化专业建设问题研究”(编号：0035-12005034000)。

2016年，我国正式成为《华盛顿协议》第18个成员。

《华盛顿协议》是国际工程教育互认协议中最具影响力的协议之一，加入其中对我国工程教育领域具有重要意义。

工程教育专业认证以学生为中心，强调成果导向的理念，通过质量监控和反馈机制持续促进教育质量的提高[1]。

自动化专业是工科中传统的宽口径专业，培养的人才具有适应多种工业及科研领域的能力。

随着智能制造、工业4.0等一系列工业新形势的诞生，传统的教学培养方式已经不能满足当代产业人才的需求。

高校需要在教育模式上积极改进和创新，以培养理论扎实、富于创新、实践力强的新型复合人才。

本文以自动化专业中重要的单片机与嵌入式系统课程为例，结合我校自动化专业教育教学的实际情况，探讨在工程教育专业认证背景下如何设计和改进课程体系，以支撑可以适应新型人才培养模式的培养目标。

1 自动化专业培养目标的制订1.1 工程教育专业认证的要求工程教育专业认证要求人才培养必须有明确公开的毕业要求，这些毕业要求应能支撑培养目标的达成。